第8章 路径积分:第三条道路

📑 目录

第8章 路径积分:第三条道路

故事场景:费曼的迷宫

2089年,普林斯顿理论物理研究所的地下档案室,年轻的研究生卢卡发现了一堆泛黄的笔记——理查德·费曼在1942年博士论文期间的手稿。其中一页画满了一个粒子从A点到B点的所有可能路径:直线、弧线、螺旋、甚至看似荒诞的折返轨迹。旁边潦草地写着:"粒子不是选择了一条路,而是走过了所有的路。"

卢卡花了三天三夜才理解这句话的深意。在量子力学中,一个电子从双缝的一侧到达屏幕,并不只是穿过左缝或右缝——它同时"探索"了所有可能的路径,每条路径贡献一个相位因子,最终的干涉图样是所有这些贡献的叠加。

这就是路径积分的核心思想,也是费曼对量子力学最具革命性的贡献。在Shankar的书中,这一章不是放在研究生级别的高级课程里,而是作为入门教材的一部分——因为他相信,路径积分是理解量子力学最直觉化的方式之一。


前置知识:变分法与路径积分的关系

从Euler-Lagrange方程到作用量原理

在经典力学中,描述运动有两种等价的方式:Newton的力方程和Lagrange的作用量原理。路径积分正是建立在后者之上。

对于一维系统,Lagrange量定义为动能减去势能:

L(x,x˙,t)=12mx˙2V(x)L(x, \dot{x}, t) = \frac{1}{2}m\dot{x}^2 - V(x)

作用量(Action)是Lagrange量沿路径的时间积分:

S[x(t)]=t1t2L(x,x˙,t)dtS[x(t)] = \int_{t_1}^{t_2} L(x, \dot{x}, t) \, dt

Hamilton原理(最小作用量原理):真实路径使作用量取极值(通常是极小值),即δS=0\delta S = 0

通过变分法,δS=0\delta S = 0给出Euler-Lagrange方程

ddtLx˙Lx=0\frac{d}{dt}\frac{\partial L}{\partial \dot{x}} - \frac{\partial L}{\partial x} = 0

对于L=12mx˙2V(x)L = \frac{1}{2}m\dot{x}^2 - V(x)

mx¨=dVdxm\ddot{x} = -\frac{dV}{dx}

这正是Newton方程!

变分法的核心思想

变分法是微积分在"函数空间"上的推广。普通微积分寻找使函数f(x)f(x)取极值的点x0x_0,变分法寻找使泛函S[x(t)]S[x(t)]取极值的函数xcl(t)x_{cl}(t)

关键步骤:

  1. 考虑路径的微小变形x(t)x(t)+δx(t)x(t) \to x(t) + \delta x(t),其中δx(t1)=δx(t2)=0\delta x(t_1) = \delta x(t_2) = 0(固定端点)。
  2. 计算作用量的变化δS=S[x+δx]S[x]\delta S = S[x + \delta x] - S[x]到一阶。
  3. 要求δS=0\delta S = 0对所有δx(t)\delta x(t)成立。
  4. 通过分部积分得到Euler-Lagrange方程。
graph TD
    subgraph 变分法与路径积分
    A[Newton方程] --> B["Euler-Lagrange方程"]
    C[Hamilton原理] --> B
    B --> D["δS = 0"]
    D --> E["经典路径 x_cl(t)"]
    E --> F[作用量 S_cl]
    F --> G[路径积分中的相位因子]
    G --> H["exp("\\"iS_cl/ℏ\\"")"]
    
    style C fill:#c8e6c9
    style E fill:#e1f5fe
    style H fill:#fff3e0
    end

作用量的几何意义

作用量不仅仅是数学工具,它具有深刻的物理意义:

  • 相空间面积:对于周期运动,SS与相空间轨迹包围的面积直接相关。
  • 光学类比:Fermat原理说光线走时间最短的路径,而经典力学的作用量原理是其力学版本。
  • 量子相位:在路径积分中,eiS/e^{iS/\hbar}是每条路径的"权重",作用量决定了量子干涉的相位。

从经典到量子的桥梁:经典力学只关心使δS=0\delta S = 0的极值路径,而量子力学让所有路径参与,每条路径贡献eiS/e^{iS/\hbar}。当0\hbar \to 0,快速振荡抹除非极值路径的贡献,只剩下经典路径。


从传播子到路径积分

传播子的严格定义

在薛定谔绘景中,态的时间演化由哈密顿量生成:

ψ(t)=eiHt/ψ(0)|\psi(t)\rangle = e^{-iHt/\hbar}|\psi(0)\rangle

在位置表象中,这转化为一个积分:

ψ(x,t)=K(x,t;x,0)ψ(x,0)dx\psi(x,t) = \int K(x,t;x',0)\psi(x',0)dx'

其中K(x,t;x,0)K(x,t;x',0)传播子(propagator),表示粒子在t=0t=0时刻从xx'出发,在tt时刻到达xx的概率幅。

传播子可以写成算符形式:

K(x,t;x,0)=xeiHt/xK(x,t;x',0) = \langle x|e^{-iHt/\hbar}|x'\rangle

graph TD
    subgraph 传播子
    A["初始态 ψ("\\"x',0\\"")"] --> B["传播子 K("\\"x,t;x',0\\"")"]
    B --> C["演化 ψ("\\"x,t\\"")"]
    D[物理意义] --> E["从x'到x的所有量子路径"]
    D --> F[每条路径贡献一个相位]
    F --> G["相位 = exp("\\"iS/ℏ\\"")"]
    
    style B fill:#c8e6c9
    style G fill:#fff3e0
    end

时间切片与路径分解:费曼的洞察

费曼的关键洞察是将时间分成无穷小段ϵ=t/N\epsilon = t/N,并在每段插入完备性关系xxdx=1\int |x\rangle\langle x|dx = 1

K(x,t;x,0)=xeiHϵ/eiHϵ/eiHϵ/xK(x,t;x',0) = \langle x|e^{-iH\epsilon/\hbar}e^{-iH\epsilon/\hbar}\cdots e^{-iH\epsilon/\hbar}|x'\rangle

在每两个相邻演化算符之间插入完备性:

K(x,t;x,0)=dx1dx2dxN1j=0N1xj+1eiHϵ/xjK(x,t;x',0) = \int dx_1 dx_2 \cdots dx_{N-1} \prod_{j=0}^{N-1} \langle x_{j+1}|e^{-iH\epsilon/\hbar}|x_j\rangle

其中x0=xx_0 = x'xN=xx_N = x

NN \to \inftyϵ0\epsilon \to 0时,这变成对所有连续路径x(t)x(t')的"积分"——粒子从(x,0)(x',0)(x,t)(x,t)的每条可能路径都贡献一个因子。

graph TD
    subgraph 时间切片
    A[时间 t 分成 N段] --> B[每段插入完备性]
    B --> C["∫dx₁...dx_{"\\"N-1\\""}"]
    C --> D["N→∞ 极限"]
    D --> E[所有连续路径的积分]
    E --> F[路径积分]
    
    style D fill:#e1f5fe
    style F fill:#c8e6c9
    end

每条路径的相位推导

对于无穷小时间步ϵ\epsilon,利用Baker-Campbell-Hausdorff公式(在ϵ0\epsilon \to 0时,eAeBeA+Be^{A}e^{B} \approx e^{A+B}[A,B][A,B]高阶小):

eiHϵ/=eiϵ(P22m+V(X))eiϵP22meiϵV(X)e^{-iH\epsilon/\hbar} = e^{-i\frac{\epsilon}{\hbar}(\frac{P^2}{2m} + V(X))} \approx e^{-i\frac{\epsilon P^2}{2m\hbar}}e^{-i\frac{\epsilon V(X)}{\hbar}}

计算矩阵元xj+1eiϵP22meiϵV(X)xj\langle x_{j+1}|e^{-i\frac{\epsilon P^2}{2m\hbar}}e^{-i\frac{\epsilon V(X)}{\hbar}}|x_j\rangle

插入动量完备性ppdp=1\int |p\rangle\langle p|dp = 1

xj+1eiϵP22mppeiϵV(X)xjdp=eiϵV(xj)eiϵp22meip(xj+1xj)/dp2π\langle x_{j+1}|e^{-i\frac{\epsilon P^2}{2m\hbar}}|p\rangle\langle p|e^{-i\frac{\epsilon V(X)}{\hbar}}|x_j\rangle dp = e^{-i\frac{\epsilon V(x_j)}{\hbar}}\int e^{-i\frac{\epsilon p^2}{2m\hbar}}e^{ip(x_{j+1}-x_j)/\hbar}\frac{dp}{2\pi\hbar}

高斯积分给出:

xj+1eiHϵ/xjm2πiϵexp[iϵ(m2(xj+1xjϵ)2V(xj))]\langle x_{j+1}|e^{-iH\epsilon/\hbar}|x_j\rangle \approx \sqrt{\frac{m}{2\pi i\hbar\epsilon}}\exp\left[\frac{i\epsilon}{\hbar}\left(\frac{m}{2}\left(\frac{x_{j+1}-x_j}{\epsilon}\right)^2 - V(x_j)\right)\right]

对所有时间步相乘,指数中的求和在ϵ0\epsilon \to 0时变为积分:

K(x,t;x,0)=D[x(t)]exp(i0tL(x,x˙)dt)K(x,t;x',0) = \int \mathcal{D}[x(t')] \exp\left(\frac{i}{\hbar}\int_0^t L(x,\dot{x})dt'\right)

其中L=12mx˙2V(x)L = \frac{1}{2}m\dot{x}^2 - V(x)拉格朗日量D[x(t)]\mathcal{D}[x(t')]表示对所有路径的"测度"。

graph LR
    subgraph 路径积分公式
    A[传播子 K] --> B["∫ D[\"x(t)\"] exp("\\"iS/ℏ\\"")"]
    B --> C["S = ∫ L dt"]
    C --> D["L = T - V"]
    D --> E[作用量 S]
    
    style B fill:#c8e6c9
    style E fill:#fff3e0
    end

与拉格朗日力学的深刻联系

最小作用量原理的量子版本

经典力学中,真实路径由最小作用量原理确定:

δS=δt1t2L(x,x˙)dt=0\delta S = \delta \int_{t_1}^{t_2} L(x,\dot{x})dt = 0

在路径积分中,所有路径都参与,但相位eiS/e^{iS/\hbar}使得作用量接近极值的路径贡献最相干——周围路径的相位快速振荡,相互抵消。只有极值点附近的路径"存活"下来。

graph TD
    subgraph 经典极限
    A["ℏ → 0"] --> B[相位快速振荡]
    B --> C[非极值路径抵消]
    C --> D["只有 δS=0 路径存活"]
    D --> E[最小作用量原理]
    E --> F["经典力学!"]
    
    style A fill:#fff3e0
    style F fill:#c8e6c9
    end

经典对应与驻相近似

0\hbar \to 0,或更准确地说,当作用量SS远大于\hbar时:

  • 只有满足δS=0\delta S = 0的经典路径有显著贡献。
  • 量子涨落被"压制"在经典轨道附近。
  • 路径积分退化为经典轨迹。

这解释了为什么宏观世界看起来是经典的:一个抛出的棒球的作用量约102710^{27}\hbar,只有极值路径周围的极小区域有相干贡献。

驻相近似(Stationary Phase Approximation):对于积分eif(x)/dx\int e^{if(x)/\hbar}dx,当0\hbar \to 0时,主要贡献来自f(x)=0f'(x) = 0的点(驻相点)。在路径积分中,f=Sf = Sf=0f' = 0就是δS=0\delta S = 0,即经典路径。


自由粒子的路径积分:从求和到解析

直接计算

自由粒子的拉格朗日量L=12mx˙2L = \frac{1}{2}m\dot{x}^2是最简单的案例。让我们一步步推导其路径积分。

经典路径是直线匀速运动:xcl(t)=x+xxttx_{cl}(t') = x' + \frac{x-x'}{t}t'

经典作用量:

Scl=0t12m(xxt)2dt=m(xx)22tS_{cl} = \int_0^t \frac{1}{2}m\left(\frac{x-x'}{t}\right)^2 dt' = \frac{m(x-x')^2}{2t}

将任意路径分解为经典路径加涨落x(t)=xcl(t)+y(t)x(t') = x_{cl}(t') + y(t'),其中y(0)=y(t)=0y(0) = y(t) = 0

由于自由粒子的Lagrange量对xx是二次的,作用量可以精确分解:

S[x(t)]=Scl+0t12my˙2dtS[x(t')] = S_{cl} + \int_0^t \frac{1}{2}m\dot{y}^2 dt'

(没有交叉项,因为经典路径使作用量取极值,一阶变分为零)

涨落积分D[y(t)]ei2my˙2dt\int \mathcal{D}[y(t')]e^{\frac{i}{2\hbar}\int m\dot{y}^2 dt'}与端点xxxx'无关,只依赖于时间tt。记为F(t)F(t)

因此:

K(x,t;x,0)=F(t)eiScl/=F(t)exp(im(xx)22t)K(x,t;x',0) = F(t)e^{iS_{cl}/\hbar} = F(t)\exp\left(\frac{im(x-x')^2}{2\hbar t}\right)

通过与t0t \to 0时的δ\delta函数极限匹配(Kδ(xx)K \to \delta(x-x')),确定:

F(t)=m2πitF(t) = \sqrt{\frac{m}{2\pi i\hbar t}}

最终传播子

K(x,t;x,0)=m2πitexp(im(xx)22t)K(x,t;x',0) = \sqrt{\frac{m}{2\pi i\hbar t}}\exp\left(\frac{im(x-x')^2}{2\hbar t}\right)

验证:这正是自由粒子薛定谔方程的格林函数!

graph TD
    subgraph 自由粒子传播子推导
    A["L = ½mẋ²"] --> B["经典路径: 直线"]
    B --> C["S_cl = m("\\"x-x'\\"")²/2t"]
    C --> D["路径分解: x = x_cl + y"]
    D --> E["涨落积分 F(t)"]
    E --> F["δ极限定 F(t)"]
    F --> G["K = √(m/2πiℏt) exp("\\"im(x-x'\\")²/2ℏt")"]
    G --> H[与薛定谔方程一致]
    
    style G fill:#c8e6c9
    style H fill:#e1f5fe
    end

数值例子:电子的自由传播

考虑一个电子(m=9.11×1031m = 9.11 \times 10^{-31} kg),初始时刻局域在x=0x' = 0处(δ\delta函数),求t=1t = 1 ns后传播到x=1x = 1 μm的概率密度。

K=9.11×10312πi×1.055×1034×109exp(i×9.11×1031×(106)22×1.055×1034×109)K = \sqrt{\frac{9.11 \times 10^{-31}}{2\pi i \times 1.055 \times 10^{-34} \times 10^{-9}}}\exp\left(\frac{i \times 9.11 \times 10^{-31} \times (10^{-6})^2}{2 \times 1.055 \times 10^{-34} \times 10^{-9}}\right)

计算振幅:

K=m2πt=9.11×10312π×1.055×1034×1091.37×10121.17×106 m1/2|K| = \sqrt{\frac{m}{2\pi\hbar t}} = \sqrt{\frac{9.11 \times 10^{-31}}{2\pi \times 1.055 \times 10^{-34} \times 10^{-9}}} \approx \sqrt{1.37 \times 10^{12}} \approx 1.17 \times 10^{6} \text{ m}^{-1/2}

指数中的相位:

m(xx)22t=9.11×1031×10122×1.055×1034×1094.32 rad\frac{m(x-x')^2}{2\hbar t} = \frac{9.11 \times 10^{-31} \times 10^{-12}}{2 \times 1.055 \times 10^{-34} \times 10^{-9}} \approx 4.32 \text{ rad}

概率密度:K2=m2πt1.37×1012|K|^2 = \frac{m}{2\pi\hbar t} \approx 1.37 \times 10^{12} m1^{-1}。注意这不是归一化的概率(对xx积分发散),因为初始δ\delta函数不是物理态。

对于物理的Gauss波包初始态,概率密度会随时间扩散,中心以群速度移动。


谐振子的路径积分:经典与涨落的分离

路径分解:经典 + 量子涨落

谐振子的拉格朗日量为L=12mx˙212mω2x2L = \frac{1}{2}m\dot{x}^2 - \frac{1}{2}m\omega^2x^2

关键技巧是将任意路径分解为经典路径加涨落

x(t)=xcl(t)+y(t)x(t') = x_{cl}(t') + y(t')

其中xcl(t)x_{cl}(t')满足经典运动方程x¨cl+ω2xcl=0\ddot{x}_{cl} + \omega^2 x_{cl} = 0,边界条件为xcl(0)=xx_{cl}(0) = x'xcl(t)=xx_{cl}(t) = x。涨落y(t)y(t')满足y(0)=y(t)=0y(0) = y(t) = 0

graph TD
    subgraph 路径分解
    A["任意路径 x(t)"] --> B["经典路径 x_cl(t)"]
    A --> C["涨落 y(t)"]
    B --> D[满足运动方程]
    C --> E["y(0)=y(t)=0"]
    D --> F[边界匹配]
    
    style B fill:#c8e6c9
    style C fill:#e1f5fe
    end

经典路径的显式解

经典运动方程x¨+ω2x=0\ddot{x} + \omega^2 x = 0的通解为:

xcl(t)=Acos(ωt)+Bsin(ωt)x_{cl}(t') = A\cos(\omega t') + B\sin(\omega t')

应用边界条件:

  • x_{cl}(0) = x' \implies A = x'
  • x_{cl}(t) = x \implies x'\cos(\omega t) + B\sin(\omega t) = x

解得:

B=xxcos(ωt)sin(ωt)B = \frac{x - x'\cos(\omega t)}{\sin(\omega t)}

因此:

xcl(t)=xcos(ωt)+xxcos(ωt)sin(ωt)sin(ωt)x_{cl}(t') = x'\cos(\omega t') + \frac{x - x'\cos(\omega t)}{\sin(\omega t)}\sin(\omega t')

作用量的分解

谐振子的Lagrange量是二次型,因此作用量的展开到二阶为:

S[x(t)]=Scl+120t(my˙2mω2y2)dtS[x(t')] = S_{cl} + \frac{1}{2}\int_0^t \left(m\dot{y}^2 - m\omega^2y^2\right)dt'

(同样没有一阶交叉项,因为经典路径使作用量取极值)

路径积分变为:

K=eiScl/D[y(t)]exp(i20t(my˙2mω2y2)dt)K = e^{iS_{cl}/\hbar} \int \mathcal{D}[y(t')]\exp\left(\frac{i}{2\hbar}\int_0^t(m\dot{y}^2 - m\omega^2y^2)dt'\right)

涨落积分只依赖于时间间隔tt,与端点xxxx'无关!

涨落行列式的计算

涨落积分是一个高斯路径积分。将涨落y(t)y(t')展开为Fourier级数:

y(t)=n=1cnsin(nπtt)y(t') = \sum_{n=1}^{\infty} c_n \sin\left(\frac{n\pi t'}{t}\right)

(满足Dirichlet边界条件y(0)=y(t)=0y(0) = y(t) = 0

作用量的涨落部分变为:

Sfl=12n=1(m2(nπt)2mω22)cn2t2S_{fl} = \frac{1}{2}\sum_{n=1}^{\infty}\left(\frac{m}{2}\left(\frac{n\pi}{t}\right)^2 - \frac{m\omega^2}{2}\right)c_n^2 \cdot \frac{t}{2}

对每个cnc_n做高斯积分,得到贡献1λn\propto \frac{1}{\sqrt{\lambda_n}},其中λn=m2(nπt)2mω22\lambda_n = \frac{m}{2}\left(\frac{n\pi}{t}\right)^2 - \frac{m\omega^2}{2}

无限乘积n=1λn\prod_{n=1}^{\infty}\lambda_n需要正规化。利用ζ\zeta函数正规化或比较已知结果(ω=0\omega = 0时应退化为自由粒子),得到涨落因子:

Fω(t)=mω2πisin(ωt)F_{\omega}(t) = \sqrt{\frac{m\omega}{2\pi i\hbar\sin(\omega t)}}

完整传播子与验证

经典作用量经过计算(需要一些三角恒等式技巧)为:

S_{cl} = \frac{m\omega}{2\sin(\omega t)}\left[(x^2 + x'^2)\cos(\omega t) - 2xx'\right]

因此完整传播子为:

K(x,t;x',0) = \sqrt{\frac{m\omega}{2\pi i\hbar\sin(\omega t)}}\exp\left(\frac{im\omega}{2\hbar\sin(\omega t)}\left[(x^2+x'^2)\cos(\omega t) - 2xx'\right]\right)

验证——ω0\omega \to 0极限

利用sin(ωt)ωt\sin(\omega t) \approx \omega tcos(ωt)1(ωt)22\cos(\omega t) \approx 1 - \frac{(\omega t)^2}{2}

S_{cl} \approx \frac{m\omega}{2\omega t}\left[(x^2+x'^2)\left(1-\frac{\omega^2t^2}{2}\right) - 2xx'\right] = \frac{m(x-x')^2}{2t} - \frac{m\omega^2 t(x^2+x'^2)}{4}

ω0\omega \to 0,第二项消失,Sclm(xx)22tS_{cl} \to \frac{m(x-x')^2}{2t}。同时mω2πisin(ωt)m2πit\frac{m\omega}{2\pi i\hbar\sin(\omega t)} \to \frac{m}{2\pi i\hbar t}

确实退化为自由粒子传播子!

graph LR
    subgraph 谐振子传播子
    A[经典作用量 S_cl] --> B["exp("\\"iS_cl/ℏ\\"")"]
    C[涨落积分] --> D["√(mω/2πiℏsin(ωt))"]
    B --> E[完整传播子]
    D --> E
    E --> F["ω→0: 退化为自由粒子"]
    
    style E fill:#c8e6c9
    style F fill:#e1f5fe
    end

数值例子:谐振子传播子的振荡

考虑一个电子在谐振子势中,ω=1015\omega = 10^{15} Hz(光频),m=9.11×1031m = 9.11 \times 10^{-31} kg。计算t=π/(2ω)t = \pi/(2\omega)时刻的传播子。

此时sin(ωt)=sin(π/2)=1\sin(\omega t) = \sin(\pi/2) = 1cos(ωt)=0\cos(\omega t) = 0

传播子:

K=mω2πiexp(imω2(2xx))=mω2πiexp(imωxx)K = \sqrt{\frac{m\omega}{2\pi i\hbar}}\exp\left(\frac{im\omega}{2\hbar}(-2xx')\right) = \sqrt{\frac{m\omega}{2\pi i\hbar}}\exp\left(-\frac{im\omega xx'}{\hbar}\right)

计算振幅:

K=mω2π=9.11×1031×10152π×1.055×10341.37×10181.17×109 m1|K| = \sqrt{\frac{m\omega}{2\pi\hbar}} = \sqrt{\frac{9.11 \times 10^{-31} \times 10^{15}}{2\pi \times 1.055 \times 10^{-34}}} \approx \sqrt{1.37 \times 10^{18}} \approx 1.17 \times 10^{9} \text{ m}^{-1}

物理解释:在这个特殊时刻,传播子类似于Fourier变换核——exp(imωxx/)\exp(-im\omega xx'/\hbar)δ\delta函数的"振荡版本"。实际上,t=π/(2ω)t = \pi/(2\omega)时,谐振子完成了1/4周期,位置空间波函数与动量空间波函数相互转换。


配分函数与统计力学:虚时间的桥梁

虚时间与路径积分

通过Wick转动tiβt \to -i\beta\hbar(即τ=it\tau = it),传播子与统计力学的配分函数联系起来:

Z=Tr(eβH)=dxK(x,iβ;x,0)Z = \mathrm{Tr}(e^{-\beta H}) = \int dx \, K(x,-i\beta\hbar;x,0)

这对应于在虚时间τ[0,β]\tau \in [0,\beta\hbar]上,起点和终点相同的路径积分——即周期性边界条件的路径。

graph TD
    subgraph 虚时间路径积分
    A["实时间: t"] --> B["虚时间: τ = it"]
    B --> C["周期性路径 x(0)=x(βℏ)"]
    C --> D[配分函数 Z]
    D --> E[统计力学]
    
    style C fill:#e1f5fe
    style E fill:#c8e6c9
    end

谐振子的配分函数

对于谐振子,将t=iβt = -i\beta\hbar代入传播子:

sin(ωt)=sin(iβω)=isinh(βω)\sin(\omega t) = \sin(-i\beta\hbar\omega) = -i\sinh(\beta\hbar\omega)

cos(ωt)=cos(iβω)=cosh(βω)\cos(\omega t) = \cos(-i\beta\hbar\omega) = \cosh(\beta\hbar\omega)

对角项K(x,iβ;x,0)K(x, -i\beta\hbar; x, 0)

K=mω2πi(isinh(βω))exp(imω2(isinh(βω))[2x2cosh(βω)2x2])K = \sqrt{\frac{m\omega}{2\pi i\hbar(-i\sinh(\beta\hbar\omega))}}\exp\left(\frac{im\omega}{2\hbar(-i\sinh(\beta\hbar\omega))}\left[2x^2\cosh(\beta\hbar\omega) - 2x^2\right]\right)

=mω2πsinh(βω)exp(mωx2tanh(βω2))= \sqrt{\frac{m\omega}{2\pi\hbar\sinh(\beta\hbar\omega)}}\exp\left(-\frac{m\omega x^2}{\hbar}\tanh(\frac{\beta\hbar\omega}{2})\right)

xx积分(高斯积分):

Z=Kdx=12sinh(βω/2)Z = \int_{-\infty}^{\infty} K \, dx = \frac{1}{2\sinh(\beta\hbar\omega/2)}

一致性验证:直接求和也能得到相同结果:

Z=n=0eβ(n+1/2)ω=eβω/211eβω=12sinh(βω/2)Z = \sum_{n=0}^{\infty}e^{-\beta(n+1/2)\hbar\omega} = e^{-\beta\hbar\omega/2}\frac{1}{1-e^{-\beta\hbar\omega}} = \frac{1}{2\sinh(\beta\hbar\omega/2)}

两种方法给出相同结果,说明虚时间路径积分与统计力学的等价性。

graph TD
    subgraph 谐振子配分函数
    A[路径积分 Z] --> B["1/(2sinh("\\"βℏω/2\\""))"]
    C[直接求和 Z] --> D["Σ exp("\\"-βE_n\\"")"]
    D --> E[相同结果]
    B --> E
    
    style E fill:#c8e6c9
    end

数值例子:光频谐振子的热激发

考虑一个光频谐振子,ω=2π×5×1014\omega = 2\pi \times 5 \times 10^{14} Hz(绿光),在室温T=300T = 300 K下:

βω=ωkBT=6.626×1034×5×10141.38×1023×30080\beta\hbar\omega = \frac{\hbar\omega}{k_B T} = \frac{6.626 \times 10^{-34} \times 5 \times 10^{14}}{1.38 \times 10^{-23} \times 300} \approx 80

这是一个很大的数!因此:

Z=12sinh(40)e404.2×1018Z = \frac{1}{2\sinh(40)} \approx e^{-40} \approx 4.2 \times 10^{-18}

平均量子数:

n=1eβω1e800\langle n \rangle = \frac{1}{e^{\beta\hbar\omega} - 1} \approx e^{-80} \approx 0

物理结论:室温下光频谐振子几乎完全处于基态(n=0n=0)。这是因为光子的能量ω2\hbar\omega \approx 2 eV远大于热能kBT0.025k_B T \approx 0.025 eV。

对比一个微波谐振子,ω=2π×1010\omega = 2\pi \times 10^{10} Hz:

βω=6.626×1034×10101.38×1023×3001.6×103\beta\hbar\omega = \frac{6.626 \times 10^{-34} \times 10^{10}}{1.38 \times 10^{-23} \times 300} \approx 1.6 \times 10^{-3}

n11.6×103625\langle n \rangle \approx \frac{1}{1.6 \times 10^{-3}} \approx 625

微波谐振子在室温下有大量热激发!


费曼图的雏形:微扰展开的几何化

微扰展开的路径积分形式

考虑一个带有微扰V(x)V(x)的系统,路径积分可以展开为微扰级数:

D[x(t)]ei(S0V)dt=D[x(t)]eiS0dt(1iVdt+12(i)2(Vdt)2+)\int \mathcal{D}[x(t')]e^{\frac{i}{\hbar}\int(S_0 - V)dt'} = \int \mathcal{D}[x(t')]e^{\frac{i}{\hbar}\int S_0 dt'}\left(1 - \frac{i}{\hbar}\int Vdt' + \frac{1}{2}\left(\frac{i}{\hbar}\right)^2\left(\int Vdt'\right)^2 + \cdots\right)

其中S0S_0是自由(或可解)部分的作用量。

每一项对应于在不同时间点"感受"微扰的路径。这些项可以用费曼图表示,其中顶点代表VV的作用,线代表自由传播子。

graph TD
    subgraph 微扰展开
    A["exp("\\"-i/ℏ ∫Vdt\\"")"] --> B[1]
    A --> C["-i/ℏ ∫Vdt"]
    A --> D["(-i/ℏ)²/2 (∫Vdt)²"]
    B --> E["零阶: 自由传播"]
    C --> F["一阶: 单顶点"]
    D --> G["二阶: 双顶点"]
    
    style E fill:#c8e6c9
    style F fill:#e1f5fe
    style G fill:#fff3e0
    end

一阶微扰的显式计算

以自由粒子为例,S0S_0对应自由传播子K0K_0,微扰V(x)V(x)的一阶修正为:

K(1)(x,t;x,0)=i0tdt1dx1K0(x,t;x1,t1)V(x1)K0(x1,t1;x,0)K^{(1)}(x,t;x',0) = -\frac{i}{\hbar}\int_0^t dt_1 \int_{-\infty}^{\infty} dx_1 K_0(x,t;x_1,t_1)V(x_1)K_0(x_1,t_1;x',0)

物理解释:粒子从xx'自由传播到x1x_1,在t1t_1时刻与势V(x1)V(x_1)"相互作用",再自由传播到xx。对所有可能的"相互作用点"(x1,t1)(x_1, t_1)积分。

这可以画成费曼图:一条线从(x,0)(x',0)(x1,t1)(x_1,t_1)(自由传播),一个"叉"表示VV的作用,再一条线从(x1,t1)(x_1,t_1)(x,t)(x,t)

二阶项涉及两个相互作用点,可以画成两条线连接三个点。

路径积分的优势

路径积分方法特别适合处理:

  1. 非微扰问题:如瞬子(instanton)隧穿,不能用微扰论处理。
  2. 规范场论:路径积分自然保持规范对称性。
  3. 统计力学:虚时间路径积分直接连接量子与统计。
  4. 路径依赖的量:如Aharonov-Bohm效应,自然包含在路径积分中。
graph TD
    subgraph 路径积分优势
    A[路径积分] --> B[非微扰隧穿]
    A --> C[规范场论]
    A --> D[统计力学联系]
    A --> E[量子引力]
    A --> F[路径依赖效应]
    
    B --> G[瞬子解]
    C --> H[规范对称性保持]
    D --> I[虚时间周期性]
    
    style A fill:#fff3e0
    style G fill:#c8e6c9
    end

路径积分的数学挑战与严格化

测度问题

严格定义D[x(t)]\int \mathcal{D}[x(t')]是一个深刻的数学问题。Wiener测度(用于虚时间/统计力学)有严格的数学基础,但实时间路径积分的振荡积分eiS/e^{iS/\hbar}需要更精细的处理。

方法包括:

  • 格点化:将时间离散化,取极限。
  • 解析延拓:先做虚时间计算,再解析延拓回实时间。
  • 渐进方法:用驻相近似(stationary phase approximation)。
graph TD
    subgraph 数学严格化
    A["路径积分 ∫D[x]"] --> B[格点化极限]
    A --> C[Wiener测度虚时间]
    A --> D[解析延拓]
    A --> E[驻相近似]
    
    B --> F["时间切片 N→∞"]
    C --> G[概率测度]
    D --> H["统计→量子"]
    E --> I["ℏ→0 展开"]
    
    style B fill:#e1f5fe
    style C fill:#c8e6c9
    end

与算符方法的等价性

路径积分与标准的薛定谔/海森堡绘景是等价的,只是不同数学包装:

薛定谔绘景路径积分
态矢$\psi\rangle$
算符X,PX,P经典变量x,px,p
对易子[X,P]=i[X,P]=i\hbar路径积分的测度
本征值问题极值作用量
微扰论(Dyson级数)费曼图展开
graph TD
    subgraph 三种表述等价性
    A[薛定谔绘景] --> B[微分方程]
    C[海森堡绘景] --> D[算符运动]
    E[路径积分] --> F[积分形式]
    A --> G[等价]
    C --> G
    E --> G
    
    style G fill:#c8e6c9
    end

本章总结

graph TD
    subgraph 第8章知识地图
    A[路径积分] --> B[传播子 K]
    A --> C[时间切片]
    A --> D[经典极限]
    A --> E[自由粒子]
    A --> F[谐振子]
    A --> G[配分函数]
    A --> H[费曼图]
    A --> I[数学严格化]
    
    B --> B1["K = ∫D[x]exp("\\"iS/ℏ\\"")"]
    B --> B2[所有路径叠加]
    
    C --> C1["N→∞ 极限"]
    C --> C2[完备性插入]
    
    D --> D1["ℏ→0 时"]
    D --> D2[极值路径主导]
    D --> D3[最小作用量原理]
    
    E --> E1[高斯积分]
    E --> E2["经典路径 = 直线"]
    E --> E3[与薛定谔方程一致]
    
    F --> F1[路径分解]
    F --> F2["经典 + 涨落"]
    F --> F3[涨落行列式]
    F --> F4["ω→0 退化验证"]
    
    G --> G1[虚时间 Wick转动]
    G --> G2[周期性路径]
    G --> G3[统计力学联系]
    G --> G4[热激发数值]
    
    H --> H1[微扰展开]
    H --> H2["顶点 = V"]
    H --> H3["线 = 自由传播子"]
    
    I --> I1[格点化]
    I --> I2[Wiener测度]
    I --> I3[解析延拓]
    
    style A fill:#fff3e0
    style D fill:#c8e6c9
    style G fill:#e1f5fe
    end

这一章我们走了量子力学的"第三条道路"——从传播子到路径积分,从所有路径的叠加到经典极限的最小作用量。

费曼的路径积分不是薛定谔方程的替代品,而是它的另一副面孔。在这副面孔中,粒子不是神秘的波函数,而是走过所有可能路径的探险者;相位不是抽象的复数,而是作用量的影子。

Shankar将这一章放在入门级别,因为他相信这是理解量子力学最物理的方式。当你看到谐振子的传播子如何从经典路径和量子涨落中优雅地浮现,当你看到虚时间如何将量子与统计统一,你会明白:路径积分不是高深的数学技巧,而是自然本身的语言。

路径积分也是量子场论的天然语言——规范对称性、费曼图、甚至量子引力,都在路径积分的框架中找到了最自然的表述。


练习与思考

问题1:从自由粒子的路径积分传播子出发,推导高斯波包在时间演化下的扩散行为。证明波包宽度满足(ΔX)2(t)=(ΔX)2(0)+(ΔP)2(0)t2/m2(\Delta X)^2(t) = (\Delta X)^2(0) + (\Delta P)^2(0)t^2/m^2

问题2:对于谐振子路径积分,详细推导涨落行列式的计算。证明涨落算符(d2dτ2ω2)\left(-\frac{d^2}{d\tau^2} - \omega^2\right)在Dirichlet边界条件下的特征值为λn=n2π2t2ω2\lambda_n = \frac{n^2\pi^2}{t^2} - \omega^2,并由此推导涨落因子F(t)F(t)

问题3:考虑一个电子在双势阱V(x)=12mω2(xa)2V(x) = \frac{1}{2}m\omega^2(|x|-a)^2(两个谐振子势阱对称放置)。用虚时间路径积分和瞬子(instanton)方法估计基态能级的分裂。讨论这种隧穿导致的能级分裂如何随势垒高度和宽度变化。

问题4:一个质量m=106m = 10^{-6} kg的宏观粒子,初始制备在x=0x = 0、宽度σ=1010\sigma = 10^{-10} m的Gauss波包。计算t=1t = 1 s后的波包宽度。如果用路径积分方法计算其传播到x=1x = 1 mm处的概率密度,经典路径的作用量是多少?相位因子eiS/e^{iS/\hbar}会振荡多少次?


"粒子不是选择了一条路,而是走过了所有的路——只是有些路,走得更响亮。"