第7章 谐振子:算符方法的胜利
故事场景:量子世界的音叉
在2162年的维也纳量子声学实验室,研究员索菲亚正在调试一台"量子音叉"——一个由超导电路构成的纳米谐振子。当她将系统冷却到接近绝对零度时,示波器上出现了奇怪的脉冲序列:单个的能量量子,像一颗颗珍珠,从基态开始逐个注入系统。
"每一颗珍珠都是一个ℏω的能量包,"她的合作者——一位来自印度的理论物理学家解释道,"这就是数算符N=a†a的本征值。谐振子不是连续变大的,它是’计数’的。"
这个场景捕捉了谐振子在量子力学中的核心地位——它不仅是最优雅的可解模型,更是量子场论的基础。从弹簧上的质点到电磁场的光子,从晶格振动到分子振动,谐振子的身影无处不在。Shankar在这里展示了他最拿手的代数解法——升降算符a和a†的舞蹈,将微分方程转化为代数游戏,将复杂的波函数转化为简单的数态。
前置知识:升降算符的思想来源
从因式分解到量子化
升降算符(ladder operators)的思想并非凭空出现——它有着深厚的数学历史。早在量子力学诞生之前,数学家们就发现,某些微分方程可以通过因式分解(factorization)方法简化。
考虑一维谐振子的薛定谔方程:
−2mℏ2dx2d2ψ+21mω2x2ψ=Eψ
这个方程可以"因式分解"为两个一阶算符的乘积。定义:
a=√2ℏmω(x+mωℏdxd),a†=√2ℏmω(x−mωℏdxd)
则薛定谔方程变为:
ℏω(a†a+21)ψ=Eψ
这就像一个二次多项式E=ℏω(N+21)被分解为"根"的形式。
历史上的因式分解方法
因式分解方法最早由Schrödinger本人在1940年系统化,用于求解各类可解势的本征值问题。后来Infeld和Hull在1951年将这种方法发展为系统的"因式分解方法"(Factorization Method),证明几乎所有可解的一维势都可以写成升降算符的形式。
谐振子的特殊性在于:它的升降算符之间的对易关系最简单——[a,a†]=1。这个关系与一维谐振子的势能形式直接相关:V(x)∝x2是唯一使得升降算符产生等间距能级的势。
从谐振子到量子场论
谐振子的代数结构是量子场论的基石:
- 电磁场:每个模式(k,λ)对应一个谐振子,akλ和akλ†分别是光子的湮灭和产生算符。
- 声子:晶格振动的量子化直接产生声子产生/湮灭算符。
- 量子电动力学:所有粒子都可以视为某种场的激发,而场激发的产生和湮灭遵循谐振子的代数。
graph TD
subgraph 升降算符的思想谱系
A[经典因式分解] --> B["Schrödinger 1940"]
B --> C["Infeld-Hull方法"]
C --> D["谐振子 a, a†"]
D --> E["[a, a†] = 1"]
E --> F["量子场论产生/湮灭算符"]
F --> G[光子]
F --> H[声子]
F --> I["夸克/胶子激发"]
style D fill:#c8e6c9
style E fill:#fff3e0
style F fill:#e1f5fe
end
谐振子的代数解法:一场算符的舞蹈
为什么谐振子如此重要?
谐振子是量子力学中最重要的模型,原因有三:
普适性:任何势能在极小值附近都可以泰勒展开为谐振子形式:V(x)≈V(x0)+21V′′(x0)(x−x0)2。
解析可解:谐振子的能级和本征函数有闭合形式,是测试新方法的基准。
场论基础:量子场论中的产生和湮灭算符直接继承自谐振子的升降算符。
graph TD
subgraph 谐振子的重要性
A[谐振子模型] --> B[分子振动]
A --> C[晶格振动声子]
A --> D[电磁场光子]
A --> E[量子场论]
A --> F[任何势阱的近似]
style A fill:#fff3e0
style E fill:#c8e6c9
end Hamilton量的重写与升降算符的完整推导
谐振子的Hamilton量为:
H=2mP2+21mω2X2
Shankar的招牌技巧是引入升降算符(或称产生/湮灭算符)。让我们完整推导这个过程。
首先,为了简化,引入特征长度x0=√mωℏ和特征动量p0=√ℏmω,定义无量纲算符:
X^=x0X,P^=p0P
则[X^,P^]=x0p01[X,P]=ℏiℏ=i,且Hamilton量变为:
H=2ℏω(P^2+X^2)
现在,类比经典复数z=x+ip的模平方∣z∣2=x2+p2,我们尝试"因式分解"X^2+P^2:
(X^−iP^)(X^+iP^)=X^2+P^2+i(X^P^−P^X^)=X^2+P^2−1
(注意对易子[X^,P^]=i,所以i(X^P^−P^X^)=i⋅i=−1)
因此:
X^2+P^2=(X^−iP^)(X^+iP^)+1
定义升降算符:
a=√21(X^+iP^)=√2ℏmω(X+mωiP)
a†=√21(X^−iP^)=√2ℏmω(X−mωiP)
这些算符不是厄米的(a≠a†),但它们的对易关系极其简洁:
[a,a†]=21[X^+iP^,X^−iP^]=21(−i[X^,P^]+i[P^,X^])=21(−i⋅i+i⋅(−i))=1
用升降算符表示Hamilton量
从定义反解出X^和P^:
X^=√21(a+a†),P^=√2i(a†−a)
计算X^2+P^2:
X^2=21(a2+aa†+a†a+(a†)2)
P^2=−21(a2−aa†−a†a+(a†)2)
相加:
X^2+P^2=aa†+a†a=2a†a+1
(利用aa†=a†a+1)
因此:
H=2ℏω(2a†a+1)=ℏω(a†a+21)=ℏω(N+21)
其中N=a†a是数算符。
关键洞察:能量本征值问题转化为寻找N的本征值。
graph LR
subgraph 升降算符
A[位置 X] --> C["a 和 a†"]
B[动量 P] --> C
C --> D["a = √(mω/2ℏ)(X + iP/mω)"]
C --> E["a† = √(mω/2ℏ)(X - iP/mω)"]
D --> F["[a, a†] = 1"]
E --> F
style F fill:#fff3e0
end 数算符的代数结构与本征态的显式构造
数算符N=a†a的本征态记为∣n⟩,对应本征值n:
N∣n⟩=n∣n⟩
从对易关系可以证明升降算符的作用:
证明a是降算符:
计算Na∣n⟩:
Na=a†aa=a(a†a)−a=a(N−1)
因此:
Na∣n⟩=a(N−1)∣n⟩=a(n−1)∣n⟩=(n−1)a∣n⟩
这说明a∣n⟩是N的本征态,本征值为n−1。因此a∣n⟩正比于∣n−1⟩:
a∣n⟩=Cn∣n−1⟩
归一化条件给出∣Cn∣2=⟨n∣a†a∣n⟩=n,因此:
a∣n⟩=√n∣n−1⟩
证明a†是升算符:
类似地,Na†=a†(N+1),因此:
a†∣n⟩=√n+1∣n+1⟩
graph TD
subgraph 升降算符作用
A["|n⟩"] --> B["a|n⟩ = √n |n-1⟩"]
A --> C["a†|n⟩ = √(n+1) |n+1⟩"]
B --> D["能量降低 ℏω"]
C --> E["能量升高 ℏω"]
style B fill:#ffccbc
style C fill:#c8e6c9
end 能级与本征态的完整推导
由于能量必须非负,降算符不能无限作用——必须存在一个基态∣0⟩使得a∣0⟩=0。
基态能量为:
E0=ℏω⟨0∣(a†a+21)∣0⟩=21ℏω
这就是著名的零点能——即使谐振子处于最低能态,它仍有非零能量。
激发态通过重复应用升算符得到:
∣n⟩=√n!(a†)n∣0⟩,En=(n+21)ℏω
能级间距:ΔE=En+1−En=ℏω,与n无关。这是谐振子最独特的性质——等间距能级。
graph LR
subgraph 能级结构
A["|0⟩"] --> B["E₀ = ½ℏω"]
A --> C["a†"]
C --> D["|1⟩"]
D --> E["E₁ = ³⁄₂ℏω"]
D --> F["a†"]
F --> G["|2⟩"]
G --> H["E₂ = ⁵⁄₂ℏω"]
style B fill:#c8e6c9
style E fill:#e1f5fe
style H fill:#fff3e0
end 基态波函数与激发态波函数的显式构造
从a∣0⟩=0在位置表象中的形式:
√2ℏmω(x+mωℏdxd)ψ0(x)=0
这是一个一阶微分方程:
dxdψ0=−ℏmωxψ0
解为高斯型基态波函数:
ψ0(x)=(πℏmω)1/4exp(−2ℏmωx2)
激发态波函数通过升算符作用得到。由于a†∝(x−mωℏdxd),反复作用产生厄米多项式:
ψn(x)=√2nn!1(πℏmω)1/4Hn(√ℏmωx)exp(−2ℏmωx2)
其中Hn是厄米多项式,满足H0(y)=1,H1(y)=2y,H2(y)=4y2−2,等等。
数值例子:光镊中的纳米粒子
考虑一个硅纳米球(m=10−15 kg,直径约100 nm)被光镊捕获在谐振子势中,陷阱频率ω=2π×104 Hz(10 kHz):
特征长度:
x0=√mωℏ=√10−15×2π×1041.055×10−34≈1.3×10−12 m=1.3 pm
零点能:
E0=21ℏω=21×1.055×10−34×2π×104≈3.3×10−30 J≈2×10−11 eV
基态位置展宽:
ΔX=√2mωℏ=√2x0≈0.92 pm
物理意义:在室温(kBT≈25 meV)下,这个纳米粒子会处于极高的激发态(n∼1012),量子效应完全被热噪声淹没。只有冷却到mK温度,才能看到量子化的迹象。
物理洞察:零点能与不确定性
基态能量E0=21ℏω不是数学 artifact,而是不确定性原理的必然结果。
对于谐振子:
ΔX=√2mωℏ,ΔP=√2ℏmω
乘积满足最小不确定性:ΔXΔP=2ℏ。
基态能量可以写成:
E0=2m(ΔP)2+21mω2(ΔX)2=4ℏω+4ℏω=2ℏω
物理图像:谐振子即使在"静止"时也在振动——这是一种量子涨落,是真空本身的不确定性的体现。
graph TD
subgraph 零点能来源
A[不确定性原理] --> B["ΔXΔP ≥ ℏ/2"]
B --> C[不能同时精确为零]
C --> D[位置涨落]
C --> E[动量涨落]
D --> F["势能 ½mω²(ΔX)²"]
E --> G["动能 (ΔP)²/2m"]
F --> H["零点能 E₀ = ℏω/2"]
G --> H
style H fill:#c8e6c9
end
相干态:最接近经典的量子态
定义与性质
相干态∣α⟩是湮灭算符的本征态:
a∣α⟩=α∣α⟩
其中α是任意复数。用数态展开:
∣α⟩=e−∣α∣2/2n=0∑∞√n!αn∣n⟩
关键特性:相干态中粒子数服从泊松分布:
P(n)=∣⟨n∣α⟩∣2=n!∣α∣2ne−∣α∣2
因此⟨n⟩=∣α∣2,Δn=∣α∣,相对涨落⟨n⟩Δn=∣α∣1=√⟨n⟩1。
graph TD
subgraph 相干态特性
A["|α⟩"] --> B["a|α⟩ = α|α⟩"]
A --> C[数态展开]
A --> D[泊松分布]
A --> E[最小不确定波包]
A --> F[时间演化保持形状]
B --> G["经典场振幅 α"]
D --> H["Δn = |α|"]
E --> I["ΔXΔP = ℏ/2"]
style E fill:#c8e6c9
style F fill:#e1f5fe
end 相干态的时间演化
相干态在时间演化下保持为相干态,只是相位旋转:
∣α(t)⟩=∣αe−iωt⟩
这意味着:
- 波包不扩散:形状保持最小不确定性。
- 质心作经典运动:⟨X(t)⟩=√mω2ℏ∣α∣cos(ωt−ϕ),其中α=∣α∣eiϕ。
- 最接近经典:相干态是量子态中最"经典"的态。
数值例子:相干态的物理尺度
考虑一个质量m=1 g的宏观谐振子(ω=1 Hz),制备在相干态∣α=106⟩:
平均量子数:⟨n⟩=1012
位置期望值:
⟨X⟩=√mω2ℏ∣α∣cos(ωt)=√10−3×12×1.055×10−34×106cos(t)≈1.45×10−11 mcos(t)
等等,这个振幅只有10−11 m,太小了。这是因为α的"自然单位"太小了。
实际上,如果我们想让宏观粒子有宏观振幅(如A=1 cm),需要:
A = \sqrt{\frac{2\hbar}{m\omega}}|\alpha| \implies |\alpha| = A\sqrt{\frac{m\omega}{2\hbar}} = 10^{-2}\sqrt{\frac{10^{-3}}{2 \times 1.055 \times 10^{-34}}} \approx 6.9 \times 10^{14}
这是一个巨大的量子数!⟨n⟩≈4.8×1029。
相对涨落:Δn/⟨n⟩=1/∣α∣≈1.4×10−15,完全不可测量。
物理结论:宏观相干态的能量不确定度相对于其平均能量小到无法察觉,这就是为什么宏观世界看起来是经典的。
相干态在位置表象中的波函数
相干态在位置表象中的波函数仍然是高斯型,但其中心随时间振荡:
ψα(x,t)=(πℏmω)1/4exp(−2ℏmω(x−⟨X(t)⟩)2+ℏi⟨P(t)⟩x+iθ(t))
其中⟨X(t)⟩和⟨P(t)⟩遵循经典谐振子运动方程:
⟨X(t)⟩=X0cos(ωt)+mωP0sin(ωt)
⟨P(t)⟩=P0cos(ωt)−mωX0sin(ωt)
这就是相干态被称为"最接近经典的量子态"的原因——它的期望位置和动量严格遵循经典方程。
物理实现:激光与微波
相干态在实验中无处不在:
- 激光:理想激光器的输出就是相干态,α与光场振幅成正比。
- 微波:超导量子电路中的微波谐振器通常制备在相干态。
- 力学振子:光镊悬浮的纳米粒子可以被冷却到量子基态,然后激发到相干态。
graph TD
subgraph 相干态实现
A["相干态 |α⟩"] --> B[激光]
A --> C[微波谐振器]
A --> D[纳米机械振子]
A --> E["玻色-爱因斯坦凝聚"]
B --> F[光场振幅]
C --> G[超导量子比特耦合]
D --> H[光镊悬浮粒子]
style A fill:#fff3e0
style B fill:#e1f5fe
style C fill:#c8e6c9
end
数表象与矩阵表示
算符的矩阵形式
在数态{∣n⟩}构成的表象中,算符有简单的矩阵形式:
数算符:
N=⎝⎜⎜⎛000⋮010⋮002⋮⋯⋯⋯⋱⎠⎟⎟⎞
降算符(非厄米,下对角线):
a=⎝⎜⎜⎛000⋮√100⋮0√20⋮⋯⋯⋯⋱⎠⎟⎟⎞
升算符(上对角线):
a†=⎝⎜⎜⎛0√10⋮00√2⋮000⋮⋯⋯⋯⋱⎠⎟⎟⎞
graph TD
subgraph 矩阵表示
A["数表象 {|n⟩}"] --> B[N 对角矩阵]
A --> C[a 下对角矩阵]
A --> D["a† 上对角矩阵"]
B --> E[本征值 n]
C --> F["元素 √n"]
D --> G["元素 √(n+1)"]
style B fill:#c8e6c9
style C fill:#ffccbc
style D fill:#e1f5fe
end 位置与动量的矩阵表示
从升降算符的定义反解:
X=√2mωℏ(a+a†),P=i√2ℏmω(a†−a)
这些算符在数表象中也是三对角矩阵,反映了X和P只能连接相邻的数态。
本章总结
graph TD
subgraph 第7章知识地图
A[谐振子] --> B["升降算符 a, a†"]
A --> C["数算符 N = a†a"]
A --> D["能级 E_n = (n+½)ℏω"]
A --> E[基态与激发态波函数]
A --> F[相干态]
A --> G[数表象矩阵]
B --> B1["[a, a†] = 1"]
B --> B2[代数推导]
C --> C1["本征值 n = 0,1,2,..."]
C --> C2[能级等间距]
D --> D1["零点能 ℏω/2"]
D --> D2[不确定性原理来源]
E --> E1[高斯基态]
E --> E2[厄米多项式激发态]
F --> F1["a|α⟩ = α|α⟩"]
F --> F2[最小不确定波包]
F --> F3[经典运动]
F --> F4[泊松分布]
G --> G1[对角化 N]
G --> G2[a 下对角]
G --> G3["a† 上对角"]
style A fill:#fff3e0
style B fill:#c8e6c9
style F fill:#e1f5fe
end这一章展示了谐振子的代数解法——算符方法的优雅与力量。升降算符a和a†的代数结构不仅给出了能级的优雅推导,还为量子场论奠定了基础。
相干态∣α⟩告诉我们,量子世界可以如此接近经典——激光的相干光、超导电路的微波,都是相干态的化身。数表象中的矩阵形式则让我们看到,整个无穷维Hilbert空间可以用简洁的递推关系操纵。
Shankar的"走廊谈话"风格在这里发挥到极致:谐振子不再是一堆微分方程的解,而是一场算符的舞蹈,一场数与概率的对话。
练习与思考
问题1:证明谐振子的相干态∣α⟩在位置表象中的波函数仍然是高斯型,但其中心随时间振荡。计算ΔX(t)和ΔP(t),证明它们保持最小不确定性。
问题2:对于谐振子,计算a2和(a†)2在数态∣n⟩中的期望值。利用这些结果推导X2和P2的期望值,并验证能量本征值。
问题3:考虑一个非谐振子V(x)=21mω2x2+λx4(λ>0)。用升降算符表示微扰λX4,计算基态能量的一阶修正。讨论当λ增大时,能级间距如何变化?这与纯谐振子的等间距能级有何本质区别?
问题4:一个光镊捕获的纳米粒子(m=10−15 kg,ω=2π×105 Hz)被冷却到基态。随后施加一个脉冲,制备出相干态∣α=3⟩。计算:(a) 位置期望值的最大振幅;(b) 该相干态中测得n=0的概率;© 如果测量粒子数,最可能测得哪个值?
"谐振子是量子力学的C大调音阶——所有更复杂的旋律,都是它的变奏。"