高斯型积分的微扰展开(一)
By 苏剑林 | 2015-02-14 | 32453位读者 |前段时间在研究费曼的路径积分理论,看到路径积分的微扰方法,也就是通过小参数展开的方式逐步逼近传播子。这样的技巧具有非常清晰的物理意义,有兴趣了解路径积分以及量子力学的读者,请去阅读费曼的《量子力学与路径积分》。然而从数学角度看来,这种逼近的技巧实际上非常粗糙,收敛范围和速度难以得到保证。事实上,数学上发展了各种各样的摄动技巧,来应对不同情况的微扰。下面我们研究积分
$$\int_{-\infty}^{+\infty} e^{-ax^2-\varepsilon x^4} dx\tag{1}$$
或者更一般地
$$\int_{-\infty}^{+\infty} e^{-ax^2-\varepsilon V(x)} dx\tag{2}$$
路径积分的级数展开比它稍微复杂一些,但是仍然是类似的形式。
简单直接的展开 #
为了计算上述积分,最简单直接的方法,就是把$e^{-\varepsilon V(x)}$作级数展开
$$e^{-\varepsilon V(x)}=\sum_{n=0}^{\infty} (-1)^n \frac{1}{n!}\varepsilon^n V^n (x)$$
从而
$$\int_{-\infty}^{+\infty} e^{-ax^2-\varepsilon V(x)} dx=\sum_{n=0}^{\infty} (-1)^n\frac{1}{n!}\varepsilon^n \int_{-\infty}^{+\infty} e^{-ax^2} V^n (x)dx$$
如果$e^{-ax^2}V^n (x)$的定积分容易计算,那么我们就可以得到原积分的一个逼近,比如对于$V(x)=x^4$的情况,我们有
$$\int_{-\infty}^{+\infty} e^{-ax^2-\varepsilon x^4} dx=\sqrt{\frac{\pi}{a}}\left(1-\frac{3\varepsilon}{4a^2}+\frac{105\varepsilon^2}{32a^4}-\frac{3465\varepsilon^3}{128a^6}+\dots\right)\tag{3}$$
该级数在$\varepsilon=a^2$的时候便完全失效了,而且就算$\varepsilon < a^2$,收敛速度也是很慢的。事实上,这是一个渐近级数,严格来讲,它的收敛区域就是$\varepsilon=0$的时候!当然,我们可以取它的前几项,来对比较大的(相对于0来说)$\varepsilon$进行近似计算。渐近级数的“渐近”就是说,取的项越少,有效计算的半径越大,能计算,但精度不一定高;取的项越多,自然越精确,但是有效半径越小。取无穷多项时,收敛半径就只有0了。
在指数中展开:一个封闭的解 #
在研究微分方程摄动展开的时候,为了得到周期解,我们不仅仅会把解关于小参数展开,而且还会把周期作关于小参数的展开。由此,笔者受到启发:能不能将原来的系数作参数展开呢?具体来讲,就是我们引入$A$
$$\int_{-\infty}^{+\infty} e^{-ax^2-\varepsilon x^4} dx=\int_{-\infty}^{+\infty} e^{-Ax^2-[(a-A)x^2+\varepsilon x^4]} dx$$
现在把方括号中的$(a-A)x^2+\varepsilon x^2$当作$V(x)$作一阶展开,得到
$$\int_{-\infty}^{+\infty} e^{-Ax^2} \left\{1-[(a-A)x^2+\varepsilon x^4]\right\}dx$$
计算
$$\int_{-\infty}^{+\infty} e^{-Ax^2}[(a-A)x^2+\varepsilon x^4]dx=\sqrt{\frac{\pi}{A}}\frac{2 A (a-A)+3 \varepsilon}{4 A^2}$$
接下来是关键的一步,选取适当的$A$,让这一项为$0$,也就是
$$2 A (a-A)+3 \varepsilon=0\,\to\, A=\frac{1}{2} \left(a+\sqrt{a^2+6 \varepsilon}\right)$$
所以
$$\int_{-\infty}^{+\infty} e^{-ax^2-\varepsilon x^4} dx\approx\int_{-\infty}^{+\infty} e^{-Ax^2} =\sqrt{\frac{\pi}{A}}=\sqrt{\frac{2\pi}{a+\sqrt{a^2+6 \varepsilon}}}\tag{4}$$
这是一个非常好的近似。有多好?好到难以置信!考虑最极端的情况:$a=0$,上式给出
$$\int_{-\infty}^{+\infty} e^{-\varepsilon x^4} dx\approx\sqrt{\frac{2\pi}{\sqrt{6 \varepsilon}}}$$
上式不仅给出了积分的一个近似值,而且就连尺度都是正确的,也就是说$\int_{-\infty}^{+\infty} e^{-\varepsilon x^4} dx$的精确值,就是一个常数乘上$\frac{1}{\varepsilon^{1/4}}$的形式!对于不同的$a,\varepsilon$组合,算出来的值比$(3)$式的一阶近似要好(有时甚至比它的高阶近似都要好)。另外一个极端情况是$\varepsilon\to+\infty$,此时积分的精确值为0,而$(4)$居然也可以给出这一结果,而$(3)$则再次完全失效了。
下表是精确值、取$(3)$式的一阶近似$\sqrt{\frac{\pi}{a}}\left(1-\frac{3\varepsilon}{4a^2}\right)$、$(4)$式的计算结果。从表中就可以看出,公式$(4)$的在所有场合都优于公式$(3)$的一阶近似。
$$\begin{array}{c|cccc}
\hline
\text{情形} & \varepsilon=1,a=10 & \varepsilon=0.1,a=1 & \varepsilon=1,a=1 & \varepsilon=1,a=0 \\
\hline
\text{精确值} & 0.556466 &1.67409 &1.36843 &1.8128 \\
\sqrt{\frac{\pi}{a}}\left(1-\frac{3\varepsilon}{4a^2}\right) & 0.556295 &1.63952 & 0.443113 & \text{负无穷} \\
\sqrt{\frac{2\pi}{a+\sqrt{a^2+6 \varepsilon}}} & 0.556402 &1.66558 & 1.31279 & 1.60159 \\
\hline
\end{array}$$
虽然$(4)$式的稳定性和精确度都较优,但它本质上仍然是一阶的近似式,后面我们会继续谈及它的高阶近似。
转载到请包括本文地址:https://kexue.fm/archives/3217
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (Feb. 14, 2015). 《高斯型积分的微扰展开(一) 》[Blog post]. Retrieved from https://kexue.fm/archives/3217
@online{kexuefm-3217,
title={高斯型积分的微扰展开(一)},
author={苏剑林},
year={2015},
month={Feb},
url={\url{https://kexue.fm/archives/3217}},
}
March 8th, 2018
对式(1)作变换$y=\sqrt{a}x$,可是使得$\varepsilon$吸收掉$a$的,这样只需讨论一个参数。
不吸收掉的原因是$a$可能为0,保留$a$是希望得到$a$和$\varepsilon$都可以任意的近似解。