第7章 谐振子:算符方法的胜利

📑 目录

第7章 谐振子:算符方法的胜利

故事场景:量子世界的音叉

在2162年的维也纳量子声学实验室,研究员索菲亚正在调试一台"量子音叉"——一个由超导电路构成的纳米谐振子。当她将系统冷却到接近绝对零度时,示波器上出现了奇怪的脉冲序列:单个的能量量子,像一颗颗珍珠,从基态开始逐个注入系统。

"每一颗珍珠都是一个ω\hbar\omega的能量包,"她的合作者——一位来自印度的理论物理学家解释道,"这就是数算符N=aaN = a^\dagger a的本征值。谐振子不是连续变大的,它是’计数’的。"

这个场景捕捉了谐振子在量子力学中的核心地位——它不仅是最优雅的可解模型,更是量子场论的基础。从弹簧上的质点到电磁场的光子,从晶格振动到分子振动,谐振子的身影无处不在。Shankar在这里展示了他最拿手的代数解法——升降算符aaaa^\dagger的舞蹈,将微分方程转化为代数游戏,将复杂的波函数转化为简单的数态。


前置知识:升降算符的思想来源

从因式分解到量子化

升降算符(ladder operators)的思想并非凭空出现——它有着深厚的数学历史。早在量子力学诞生之前,数学家们就发现,某些微分方程可以通过因式分解(factorization)方法简化。

考虑一维谐振子的薛定谔方程:

22md2ψdx2+12mω2x2ψ=Eψ-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} + \frac{1}{2}m\omega^2x^2\psi = E\psi

这个方程可以"因式分解"为两个一阶算符的乘积。定义:

a=mω2(x+mωddx),a=mω2(xmωddx)a = \sqrt{\frac{m\omega}{2\hbar}}\left(x + \frac{\hbar}{m\omega}\frac{d}{dx}\right), \quad a^\dagger = \sqrt{\frac{m\omega}{2\hbar}}\left(x - \frac{\hbar}{m\omega}\frac{d}{dx}\right)

则薛定谔方程变为:

ω(aa+12)ψ=Eψ\hbar\omega\left(a^\dagger a + \frac{1}{2}\right)\psi = E\psi

这就像一个二次多项式E=ω(N+12)E = \hbar\omega(N + \frac{1}{2})被分解为"根"的形式。

历史上的因式分解方法

因式分解方法最早由Schrödinger本人在1940年系统化,用于求解各类可解势的本征值问题。后来Infeld和Hull在1951年将这种方法发展为系统的"因式分解方法"(Factorization Method),证明几乎所有可解的一维势都可以写成升降算符的形式。

谐振子的特殊性在于:它的升降算符之间的对易关系最简单——[a,a]=1[a, a^\dagger] = 1。这个关系与一维谐振子的势能形式直接相关:V(x)x2V(x) \propto x^2是唯一使得升降算符产生等间距能级的势。

从谐振子到量子场论

谐振子的代数结构是量子场论的基石:

  • 电磁场:每个模式(k,λ)(k, \lambda)对应一个谐振子,akλa_{k\lambda}akλa^\dagger_{k\lambda}分别是光子的湮灭和产生算符。
  • 声子:晶格振动的量子化直接产生声子产生/湮灭算符。
  • 量子电动力学:所有粒子都可以视为某种场的激发,而场激发的产生和湮灭遵循谐振子的代数。
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

谐振子的代数解法:一场算符的舞蹈

为什么谐振子如此重要?

谐振子是量子力学中最重要的模型,原因有三:

  1. 普适性:任何势能在极小值附近都可以泰勒展开为谐振子形式:V(x)V(x0)+12V(x0)(xx0)2V(x) \approx V(x_0) + \frac{1}{2}V''(x_0)(x-x_0)^2

  2. 解析可解:谐振子的能级和本征函数有闭合形式,是测试新方法的基准。

  3. 场论基础:量子场论中的产生和湮灭算符直接继承自谐振子的升降算符。

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=P22m+12mω2X2H = \frac{P^2}{2m} + \frac{1}{2}m\omega^2X^2

Shankar的招牌技巧是引入升降算符(或称产生/湮灭算符)。让我们完整推导这个过程。

首先,为了简化,引入特征长度x0=mωx_0 = \sqrt{\frac{\hbar}{m\omega}}特征动量p0=mωp_0 = \sqrt{\hbar m\omega},定义无量纲算符:

X^=Xx0,P^=Pp0\hat{X} = \frac{X}{x_0}, \quad \hat{P} = \frac{P}{p_0}

[X^,P^]=1x0p0[X,P]=i=i[\hat{X}, \hat{P}] = \frac{1}{x_0 p_0}[X, P] = \frac{i\hbar}{\hbar} = i,且Hamilton量变为:

H=ω2(P^2+X^2)H = \frac{\hbar\omega}{2}(\hat{P}^2 + \hat{X}^2)

现在,类比经典复数z=x+ipz = x + ip的模平方z2=x2+p2|z|^2 = x^2 + p^2,我们尝试"因式分解"X^2+P^2\hat{X}^2 + \hat{P}^2

(X^iP^)(X^+iP^)=X^2+P^2+i(X^P^P^X^)=X^2+P^21(\hat{X} - i\hat{P})(\hat{X} + i\hat{P}) = \hat{X}^2 + \hat{P}^2 + i(\hat{X}\hat{P} - \hat{P}\hat{X}) = \hat{X}^2 + \hat{P}^2 - 1

(注意对易子[X^,P^]=i[\hat{X}, \hat{P}] = i,所以i(X^P^P^X^)=ii=1i(\hat{X}\hat{P} - \hat{P}\hat{X}) = i \cdot i = -1

因此:

X^2+P^2=(X^iP^)(X^+iP^)+1\hat{X}^2 + \hat{P}^2 = (\hat{X} - i\hat{P})(\hat{X} + i\hat{P}) + 1

定义升降算符:

a=12(X^+iP^)=mω2(X+iPmω)a = \frac{1}{\sqrt{2}}(\hat{X} + i\hat{P}) = \sqrt{\frac{m\omega}{2\hbar}}\left(X + \frac{iP}{m\omega}\right)

a=12(X^iP^)=mω2(XiPmω)a^\dagger = \frac{1}{\sqrt{2}}(\hat{X} - i\hat{P}) = \sqrt{\frac{m\omega}{2\hbar}}\left(X - \frac{iP}{m\omega}\right)

这些算符不是厄米的(aaa \neq a^\dagger),但它们的对易关系极其简洁:

[a,a]=12[X^+iP^,X^iP^]=12(i[X^,P^]+i[P^,X^])=12(ii+i(i))=1[a, a^\dagger] = \frac{1}{2}[\hat{X} + i\hat{P}, \hat{X} - i\hat{P}] = \frac{1}{2}(-i[\hat{X}, \hat{P}] + i[\hat{P}, \hat{X}]) = \frac{1}{2}(-i \cdot i + i \cdot (-i)) = 1

用升降算符表示Hamilton量

从定义反解出X^\hat{X}P^\hat{P}

X^=12(a+a),P^=i2(aa)\hat{X} = \frac{1}{\sqrt{2}}(a + a^\dagger), \quad \hat{P} = \frac{i}{\sqrt{2}}(a^\dagger - a)

计算X^2+P^2\hat{X}^2 + \hat{P}^2

X^2=12(a2+aa+aa+(a)2)\hat{X}^2 = \frac{1}{2}(a^2 + aa^\dagger + a^\dagger a + (a^\dagger)^2)

P^2=12(a2aaaa+(a)2)\hat{P}^2 = -\frac{1}{2}(a^2 - aa^\dagger - a^\dagger a + (a^\dagger)^2)

相加:

X^2+P^2=aa+aa=2aa+1\hat{X}^2 + \hat{P}^2 = aa^\dagger + a^\dagger a = 2a^\dagger a + 1

(利用aa=aa+1aa^\dagger = a^\dagger a + 1

因此:

H=ω2(2aa+1)=ω(aa+12)=ω(N+12)H = \frac{\hbar\omega}{2}(2a^\dagger a + 1) = \hbar\omega\left(a^\dagger a + \frac{1}{2}\right) = \hbar\omega\left(N + \frac{1}{2}\right)

其中N=aaN = a^\dagger a数算符

关键洞察:能量本征值问题转化为寻找NN的本征值。

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=aaN = a^\dagger a的本征态记为n|n\rangle,对应本征值nn

Nn=nnN|n\rangle = n|n\rangle

从对易关系可以证明升降算符的作用:

证明aa是降算符

计算NanNa|n\rangle

Na=aaa=a(aa)a=a(N1)Na = a^\dagger a a = a(a^\dagger a) - a = a(N - 1)

因此:

Nan=a(N1)n=a(n1)n=(n1)anNa|n\rangle = a(N-1)|n\rangle = a(n-1)|n\rangle = (n-1)a|n\rangle

这说明ana|n\rangleNN的本征态,本征值为n1n-1。因此ana|n\rangle正比于n1|n-1\rangle

an=Cnn1a|n\rangle = C_n|n-1\rangle

归一化条件给出Cn2=naan=n|C_n|^2 = \langle n|a^\dagger a|n\rangle = n,因此:

an=nn1a|n\rangle = \sqrt{n}|n-1\rangle

证明aa^\dagger是升算符

类似地,Na=a(N+1)Na^\dagger = a^\dagger(N+1),因此:

an=n+1n+1a^\dagger|n\rangle = \sqrt{n+1}|n+1\rangle

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|0\rangle使得a0=0a|0\rangle = 0

基态能量为:

E0=ω0(aa+12)0=12ωE_0 = \hbar\omega\langle 0|(a^\dagger a + \frac{1}{2})|0\rangle = \frac{1}{2}\hbar\omega

这就是著名的零点能——即使谐振子处于最低能态,它仍有非零能量。

激发态通过重复应用升算符得到:

n=(a)nn!0,En=(n+12)ω|n\rangle = \frac{(a^\dagger)^n}{\sqrt{n!}}|0\rangle, \quad E_n = \left(n + \frac{1}{2}\right)\hbar\omega

能级间距ΔE=En+1En=ω\Delta E = E_{n+1} - E_n = \hbar\omega,与nn无关。这是谐振子最独特的性质——等间距能级。

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

基态波函数与激发态波函数的显式构造

a0=0a|0\rangle = 0在位置表象中的形式:

mω2(x+mωddx)ψ0(x)=0\sqrt{\frac{m\omega}{2\hbar}}\left(x + \frac{\hbar}{m\omega}\frac{d}{dx}\right)\psi_0(x) = 0

这是一个一阶微分方程:

dψ0dx=mωxψ0\frac{d\psi_0}{dx} = -\frac{m\omega}{\hbar}x\psi_0

解为高斯型基态波函数:

ψ0(x)=(mωπ)1/4exp(mω2x2)\psi_0(x) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\exp\left(-\frac{m\omega}{2\hbar}x^2\right)

激发态波函数通过升算符作用得到。由于a(xmωddx)a^\dagger \propto (x - \frac{\hbar}{m\omega}\frac{d}{dx}),反复作用产生厄米多项式

ψn(x)=12nn!(mωπ)1/4Hn(mωx)exp(mω2x2)\psi_n(x) = \frac{1}{\sqrt{2^n n!}}\left(\frac{m\omega}{\pi\hbar}\right)^{1/4}H_n\left(\sqrt{\frac{m\omega}{\hbar}}x\right)\exp\left(-\frac{m\omega}{2\hbar}x^2\right)

其中HnH_n是厄米多项式,满足H0(y)=1H_0(y) = 1H1(y)=2yH_1(y) = 2yH2(y)=4y22H_2(y) = 4y^2 - 2,等等。

数值例子:光镊中的纳米粒子

考虑一个硅纳米球(m=1015m = 10^{-15} kg,直径约100 nm)被光镊捕获在谐振子势中,陷阱频率ω=2π×104\omega = 2\pi \times 10^4 Hz(10 kHz):

特征长度

x0=mω=1.055×10341015×2π×1041.3×1012 m=1.3 pmx_0 = \sqrt{\frac{\hbar}{m\omega}} = \sqrt{\frac{1.055 \times 10^{-34}}{10^{-15} \times 2\pi \times 10^4}} \approx 1.3 \times 10^{-12} \text{ m} = 1.3 \text{ pm}

零点能

E0=12ω=12×1.055×1034×2π×1043.3×1030 J2×1011 eVE_0 = \frac{1}{2}\hbar\omega = \frac{1}{2} \times 1.055 \times 10^{-34} \times 2\pi \times 10^4 \approx 3.3 \times 10^{-30} \text{ J} \approx 2 \times 10^{-11} \text{ eV}

基态位置展宽

ΔX=2mω=x020.92 pm\Delta X = \sqrt{\frac{\hbar}{2m\omega}} = \frac{x_0}{\sqrt{2}} \approx 0.92 \text{ pm}

物理意义:在室温(kBT25k_B T \approx 25 meV)下,这个纳米粒子会处于极高的激发态(n1012n \sim 10^{12}),量子效应完全被热噪声淹没。只有冷却到mK温度,才能看到量子化的迹象。

物理洞察:零点能与不确定性

基态能量E0=12ωE_0 = \frac{1}{2}\hbar\omega不是数学 artifact,而是不确定性原理的必然结果。

对于谐振子:

ΔX=2mω,ΔP=mω2\Delta X = \sqrt{\frac{\hbar}{2m\omega}}, \quad \Delta P = \sqrt{\frac{\hbar m\omega}{2}}

乘积满足最小不确定性:ΔXΔP=2\Delta X \Delta P = \frac{\hbar}{2}

基态能量可以写成:

E0=(ΔP)22m+12mω2(ΔX)2=ω4+ω4=ω2E_0 = \frac{(\Delta P)^2}{2m} + \frac{1}{2}m\omega^2(\Delta X)^2 = \frac{\hbar\omega}{4} + \frac{\hbar\omega}{4} = \frac{\hbar\omega}{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

相干态:最接近经典的量子态

定义与性质

相干态α|\alpha\rangle是湮灭算符的本征态:

aα=ααa|\alpha\rangle = \alpha|\alpha\rangle

其中α\alpha是任意复数。用数态展开:

α=eα2/2n=0αnn!n|\alpha\rangle = e^{-|\alpha|^2/2}\sum_{n=0}^{\infty}\frac{\alpha^n}{\sqrt{n!}}|n\rangle

关键特性:相干态中粒子数服从泊松分布

P(n)=nα2=α2neα2n!P(n) = |\langle n|\alpha\rangle|^2 = \frac{|\alpha|^{2n}e^{-|\alpha|^2}}{n!}

因此n=α2\langle n \rangle = |\alpha|^2Δn=α\Delta n = |\alpha|,相对涨落Δnn=1α=1n\frac{\Delta n}{\langle n \rangle} = \frac{1}{|\alpha|} = \frac{1}{\sqrt{\langle n \rangle}}

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)=αeiωt|\alpha(t)\rangle = |\alpha e^{-i\omega t}\rangle

这意味着:

  1. 波包不扩散:形状保持最小不确定性。
  2. 质心作经典运动X(t)=2mωαcos(ωtϕ)\langle X(t) \rangle = \sqrt{\frac{2\hbar}{m\omega}}|\alpha|\cos(\omega t - \phi),其中α=αeiϕ\alpha = |\alpha|e^{i\phi}
  3. 最接近经典:相干态是量子态中最"经典"的态。

数值例子:相干态的物理尺度

考虑一个质量m=1m = 1 g的宏观谐振子(ω=1\omega = 1 Hz),制备在相干态α=106|\alpha = 10^6\rangle

平均量子数n=1012\langle n \rangle = 10^{12}

位置期望值

X=2mωαcos(ωt)=2×1.055×1034103×1×106cos(t)1.45×1011 mcos(t)\langle X \rangle = \sqrt{\frac{2\hbar}{m\omega}}|\alpha|\cos(\omega t) = \sqrt{\frac{2 \times 1.055 \times 10^{-34}}{10^{-3} \times 1}} \times 10^6 \cos(t) \approx 1.45 \times 10^{-11} \text{ m} \cos(t)

等等,这个振幅只有101110^{-11} m,太小了。这是因为α\alpha的"自然单位"太小了。

实际上,如果我们想让宏观粒子有宏观振幅(如A=1A = 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}

这是一个巨大的量子数!n4.8×1029\langle n \rangle \approx 4.8 \times 10^{29}

相对涨落Δn/n=1/α1.4×1015\Delta n / \langle n \rangle = 1/|\alpha| \approx 1.4 \times 10^{-15},完全不可测量。

物理结论:宏观相干态的能量不确定度相对于其平均能量小到无法察觉,这就是为什么宏观世界看起来是经典的。

相干态在位置表象中的波函数

相干态在位置表象中的波函数仍然是高斯型,但其中心随时间振荡:

ψα(x,t)=(mωπ)1/4exp(mω2(xX(t))2+iP(t)x+iθ(t))\psi_\alpha(x,t) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\exp\left(-\frac{m\omega}{2\hbar}(x - \langle X(t)\rangle)^2 + \frac{i}{\hbar}\langle P(t)\rangle x + i\theta(t)\right)

其中X(t)\langle X(t)\rangleP(t)\langle P(t)\rangle遵循经典谐振子运动方程:

X(t)=X0cos(ωt)+P0mωsin(ωt)\langle X(t)\rangle = X_0\cos(\omega t) + \frac{P_0}{m\omega}\sin(\omega t)

P(t)=P0cos(ωt)mωX0sin(ωt)\langle P(t)\rangle = P_0\cos(\omega t) - m\omega X_0\sin(\omega t)

这就是相干态被称为"最接近经典的量子态"的原因——它的期望位置和动量严格遵循经典方程。

物理实现:激光与微波

相干态在实验中无处不在:

  • 激光:理想激光器的输出就是相干态,α\alpha与光场振幅成正比。
  • 微波:超导量子电路中的微波谐振器通常制备在相干态。
  • 力学振子:光镊悬浮的纳米粒子可以被冷却到量子基态,然后激发到相干态。
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\rangle\}构成的表象中,算符有简单的矩阵形式:

数算符

N=(000010002)N = \begin{pmatrix} 0 & 0 & 0 & \cdots \\ 0 & 1 & 0 & \cdots \\ 0 & 0 & 2 & \cdots \\ \vdots & \vdots & \vdots & \ddots \end{pmatrix}

降算符(非厄米,下对角线):

a=(010002000)a = \begin{pmatrix} 0 & \sqrt{1} & 0 & \cdots \\ 0 & 0 & \sqrt{2} & \cdots \\ 0 & 0 & 0 & \cdots \\ \vdots & \vdots & \vdots & \ddots \end{pmatrix}

升算符(上对角线):

a=(000100020)a^\dagger = \begin{pmatrix} 0 & 0 & 0 & \cdots \\ \sqrt{1} & 0 & 0 & \cdots \\ 0 & \sqrt{2} & 0 & \cdots \\ \vdots & \vdots & \vdots & \ddots \end{pmatrix}

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=imω2(aa)X = \sqrt{\frac{\hbar}{2m\omega}}(a + a^\dagger), \quad P = i\sqrt{\frac{\hbar m\omega}{2}}(a^\dagger - a)

这些算符在数表象中也是三对角矩阵,反映了XXPP只能连接相邻的数态。


本章总结

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

这一章展示了谐振子的代数解法——算符方法的优雅与力量。升降算符aaaa^\dagger的代数结构不仅给出了能级的优雅推导,还为量子场论奠定了基础。

相干态α|\alpha\rangle告诉我们,量子世界可以如此接近经典——激光的相干光、超导电路的微波,都是相干态的化身。数表象中的矩阵形式则让我们看到,整个无穷维Hilbert空间可以用简洁的递推关系操纵。

Shankar的"走廊谈话"风格在这里发挥到极致:谐振子不再是一堆微分方程的解,而是一场算符的舞蹈,一场数与概率的对话。


练习与思考

问题1:证明谐振子的相干态α|\alpha\rangle在位置表象中的波函数仍然是高斯型,但其中心随时间振荡。计算ΔX(t)\Delta X(t)ΔP(t)\Delta P(t),证明它们保持最小不确定性。

问题2:对于谐振子,计算a2a^2(a)2(a^\dagger)^2在数态n|n\rangle中的期望值。利用这些结果推导X2X^2P2P^2的期望值,并验证能量本征值。

问题3:考虑一个非谐振子V(x)=12mω2x2+λx4V(x) = \frac{1}{2}m\omega^2x^2 + \lambda x^4λ>0\lambda > 0)。用升降算符表示微扰λX4\lambda X^4,计算基态能量的一阶修正。讨论当λ\lambda增大时,能级间距如何变化?这与纯谐振子的等间距能级有何本质区别?

问题4:一个光镊捕获的纳米粒子(m=1015m = 10^{-15} kg,ω=2π×105\omega = 2\pi \times 10^5 Hz)被冷却到基态。随后施加一个脉冲,制备出相干态α=3|\alpha = 3\rangle。计算:(a) 位置期望值的最大振幅;(b) 该相干态中测得n=0n=0的概率;© 如果测量粒子数,最可能测得哪个值?


"谐振子是量子力学的C大调音阶——所有更复杂的旋律,都是它的变奏。"