第VIII章:碰撞问题——量子散射的艺术

📑 目录

第VIII章:碰撞问题——量子散射的艺术

"碰撞不是终结,而是另一种形式的相遇。在经典物理中,碰撞是轨迹的改变;在量子力学中,散射是概率振幅的干涉艺术。"


前置知识:散射问题的经典图像

在深入量子散射理论之前,我们需要建立一个清晰的经典参照系。经典碰撞理论不仅是量子散射的历史前身,更是理解量子结果的直觉基础。

1. 经典散射的基本框架

考虑一个质量为 mm 的粒子以初速度 v0v_0 从远处入射,受到一个中心势场 V(r)V(r) 的作用后偏转。经典力学中,散射由以下要素刻画:

碰撞参数(impact parameter)bb:入射粒子如果没有受力,将与力心的最短距离。

散射角 θ\theta:入射方向与出射方向之间的夹角。

守恒量

  • 能量守恒:E=12mv02=12mr˙2+L22mr2+V(r)E = \frac{1}{2}mv_0^2 = \frac{1}{2}m\dot{r}^2 + \frac{L^2}{2mr^2} + V(r)
  • 角动量守恒:L=mv0b=mr2ϕ˙L = mv_0 b = mr^2\dot{\phi}

其中 LL 是角动量,ϕ\phi 是极角。

轨道方程:利用 u=1/ru = 1/r,轨道满足Binet方程:

d2udϕ2+u=mL2u2F(1/u)\frac{d^2u}{d\phi^2} + u = -\frac{m}{L^2u^2}F(1/u)

其中 F(r)=dV/drF(r) = -dV/dr 是力。

2. 散射角与碰撞参数的关系

对于中心势场,散射角由以下积分给出:

θ(b)=π2rminbdrr21V(r)/Eb2/r2\theta(b) = \pi - 2\int_{r_{min}}^{\infty} \frac{b dr}{r^2\sqrt{1 - V(r)/E - b^2/r^2}}

其中 rminr_{min} 是最近距离,由根号内表达式为零确定。

物理意义

  • 碰撞参数 bb 越小,粒子越接近力心,散射角越大
  • b=0b = 0(正碰)时,θ=π\theta = \pi(粒子被反弹回来)
  • bb \to \infty 时,θ0\theta \to 0(粒子不受偏转)

3. 经典截面

经典力学中,所有具有相同碰撞参数的粒子具有相同的散射角。因此,我们可以定义微分截面:

考虑一束均匀入射的粒子流,单位面积单位时间的粒子数为 jinj_{in}。碰撞参数在 [b,b+db][b, b+db] 环内的粒子数为 dN=jin2πbdbdN = j_{in} \cdot 2\pi b\, db。这些粒子散射到角度 [θ,θ+dθ][\theta, \theta+d\theta] 内。

经典微分截面

(dσdΩ)cl=bsinθdbdθ\left(\frac{d\sigma}{d\Omega}\right)_{cl} = \frac{b}{\sin\theta} \left|\frac{db}{d\theta}\right|

关键特征:经典截面只依赖于 b(θ)b(\theta) 的函数形式。如果存在多个 bb 给出相同的 θ\theta,需要对所有贡献求和。

彩虹效应与 glory:当 db/dθ=0db/d\theta = 0(即 θ(b)\theta(b) 取极值)时,经典截面发散——这是经典力学向量子力学过渡的信号。在量子力学中,这些发散被干涉效应平滑。

graph TD
    A["经典散射框架"] --> B["碰撞参数 b"]
    A --> C["散射角 θ(b)"]
    A --> D["守恒量: E, L"]
    
    B --> E["b越小 → θ越大"]
    C --> F["db/dθ=0 → 截面发散
彩虹/ glory"]
    
    D --> G["轨道积分:
θ = π - 2∫..."]
    G --> H["Binet方程:
d²u/dφ² + u = -mF/(L²u²)"]
    
    A --> I["经典微分截面:
dσ/dΩ = b|db/dθ|/sinθ"]
    
    style A fill:#e3f2fd
    style E fill:#c8e6c9
    style F fill:#ffcdd2
    style I fill:#90caf9

4. 经典到量子的过渡

量子力学修正了经典图像的几个核心方面:

  1. 无确定轨迹:粒子不具有同时确定的 rrpp,碰撞参数 bb 不再是一个经典确定的量
  2. 波包扩散:入射波包在传播中扩散,"碰撞参数"变成一个概率分布
  3. 干涉效应:不同"路径"(不同经典碰撞参数)的量子振幅相互干涉,导致彩虹发散被平滑
  4. 隧道效应:量子粒子可以穿透经典禁戒的势垒,进入 E<V(r)E < V(r) 的区域
  5. 共振:量子粒子可以暂时"困在"势场中形成准束缚态——这在经典力学中不存在

Dirac在《量子力学原理》第VIII章中,正是从这种经典图像出发,逐步构建量子散射的全新语言。


引言:星际尘埃中的相遇

2203年,"深空五号"探测器在穿越奥尔特云时,搭载的首席物理学家苏明遇到了一个前所未见的挑战。飞船的传感器捕获到一个奇异的现象:一束高能质子流在穿越一片稀薄的氢气云后,并未按照经典轨道力学所预言的那样发生偏转——它们的出射角度呈现出一种概率性的分布,仿佛每一颗粒子都在进行一场量子层面的"掷骰子"。更诡异的是,在特定能量下,散射强度出现了尖锐的峰值,随后又迅速跌落,如同一场精心编排的量子共振交响乐。

苏明想起了Dirac《量子力学原理》中的第VIII章。在经典物理中,碰撞只是一个轨迹的改变;但在量子力学中,散射是一种概率振幅的干涉艺术。Dirac写道,碰撞问题的核心不在于追踪粒子的轨迹,而在于计算一个粒子从入射态转变为出射态的振幅。这要求一种全新的数学语言——动量表象中的变换理论

飞船上的量子计算机开始运行Dirac的散射算法。苏明屏住呼吸:这些计算将决定探测器是否能安全穿越前方那片神秘的、密度不均的尘埃带。


一、从经典碰撞到量子散射:范式的转换

1.1 经典碰撞的图像

在牛顿力学中,两个粒子的碰撞由守恒定律完全决定:

  • 动量守恒m1v1+m2v2=m1v1+m2v2m_1\mathbf{v}_1 + m_2\mathbf{v}_2 = m_1\mathbf{v}_1' + m_2\mathbf{v}_2'
  • 能量守恒12m1v12+12m2v22=12m1(v1)2+12m2(v2)2\frac{1}{2}m_1v_1^2 + \frac{1}{2}m_2v_2^2 = \frac{1}{2}m_1(v_1')^2 + \frac{1}{2}m_2(v_2')^2

给定入射条件,出射速度和方向是唯一确定的。轨迹可以被精确计算,碰撞参数(impact parameter)与散射角之间存在一一对应的关系。

graph LR
    subgraph 经典碰撞
        A["入射粒子
动量 pᵢ"] --> B["碰撞参数 b"]
        B --> C["势场 V(r)"]
        C --> D["出射粒子
动量 p_f"]
        D --> E{"散射角 θ"}
        E -->|确定| F[唯一轨迹]
    end
    style A fill:#1a1a2e,color:#fff
    style B fill:#16213e,color:#fff
    style C fill:#0f3460,color:#fff
    style D fill:#16213e,color:#fff
    style E fill:#533483,color:#fff
    style F fill:#e94560,color:#fff

1.2 量子碰撞的本质

在量子力学中,粒子不具有同时确定的动量和位置。"碰撞"不再是一个发生在特定时空点的局域事件,而是一个入射波包与势场相互作用后演化为出射波包的过程。

Dirac在第VIII章中强调:量子碰撞问题的核心是散射系数(scattering coefficient)——一个复数,其模方给出散射概率。这不是对经典碰撞的"修正",而是一种根本不同的描述方式。

量子碰撞由以下要素刻画:

  1. 入射态:具有确定动量 p\mathbf{p} 的平面波 p|\mathbf{p}\rangle
  2. 相互作用区域:势场 V(r)V(\mathbf{r}) 存在的有限空间范围
  3. 出射态:沿不同方向传播的球面波叠加
  4. 散射振幅 f(θ,ϕ)f(\theta, \phi):决定出射方向上的概率分布
graph LR
    subgraph 量子散射
        A1["入射平面波
|p⟩"] --> B1["势场 V(r)"]
        B1 --> C1["出射球面波
|p'⟩"]
        C1 --> D1["散射振幅
f("θ,φ")"])]
        D1 --> E1["概率分布
|f|²dΩ"]
    end
    style A1 fill:#1a1a2e,color:#fff
    style B1 fill:#16213e,color:#fff
    style C1 fill:#0f3460,color:#fff
    style D1 fill:#533483,color:#fff
    style E1 fill:#e94560,color:#fff

二、散射系数的定义与物理意义

2.1 渐近态与散射算符

Dirac引入了一个精妙的概念:散射算符(Scattering Operator)S矩阵。在相互作用绘景中(回顾第VII章),如果相互作用在遥远的过去和未来都消失,那么系统的演化可以从 t=t = -\infty 的"入射自由态"映射到 t=+t = +\infty 的"出射自由态":

ψout=Sψin|\psi_{\text{out}}\rangle = S |\psi_{\text{in}}\rangle

S矩阵的矩阵元:

Sfi=fSiS_{fi} = \langle f | S | i \rangle

就是从初态 i|i\rangle 到末态 f|f\rangle总散射振幅(包含透射和反射的所有贡献)。

S矩阵具有两个关键性质:

  1. 幺正性SS=SS=IS^\dagger S = SS^\dagger = I——概率守恒
  2. 能量守恒SS 只连接相同能量的态——[S,H0]=0[S, H_0] = 0

2.2 微分截面与总截面

实验上可测量的量是散射截面(scattering cross section)。考虑一束入射粒子流(单位面积单位时间的粒子数)jinj_{\text{in}},在出射方向 (θ,ϕ)(\theta, \phi) 上,单位立体角 dΩd\Omega 内检测到的散射粒子数 dndn 满足:

dndΩ=jindσdΩ\frac{dn}{d\Omega} = j_{\text{in}} \cdot \frac{d\sigma}{d\Omega}

微分截面定义为:

dσdΩ=f(θ,ϕ)2\frac{d\sigma}{d\Omega} = |f(\theta, \phi)|^2

其中 f(θ,ϕ)f(\theta, \phi)散射振幅,与S矩阵的关系为:

Sfi=δfi2πiδ(EfEi)TfiS_{fi} = \delta_{fi} - 2\pi i \delta(E_f - E_i) T_{fi}

这里 TfiT_{fi}跃迁矩阵元(T-matrix element),而:

f(θ,ϕ)=m2π2Tfif(\theta, \phi) = -\frac{m}{2\pi\hbar^2} T_{fi}

总截面是所有方向的积分:

σtot=f(θ,ϕ)2dΩ\sigma_{\text{tot}} = \int |f(\theta, \phi)|^2 d\Omega

总截面具有面积的量纲,其物理意义是:入射粒子"看到"的有效靶面积。

graph TB
    A["入射流密度 j_in"] --> B[靶粒子]
    B --> C{"散射角 θ"}
    C --> D["出射粒子 dn/dΩ"]
    D --> E["dσ/dΩ = |f(θ,φ)|²"]²]
    E --> F["σ_tot = ∫|f|²dΩ"]
    style A fill:#1a1a2e,color:#fff
    style B fill:#16213e,color:#fff
    style C fill:#0f3460,color:#fff
    style D fill:#533483,color:#fff
    style E fill:#e94560,color:#fff
    style F fill:#e94560,color:#fff

三、动量表象中的散射:Dirac的变换方法

3.1 为什么选择动量表象?

Dirac在第VIII章中旗帜鲜明地选择了动量表象(momentum representation)作为处理碰撞问题的主战场。这不是偶然:

  1. 自由粒子的本征态在动量表象中是对角的——p\mathbf{p} 就是量子数
  2. 碰撞前后,粒子在远离相互作用区时表现为自由粒子——动量是好量子数
  3. 能量守恒天然地在动量表象中表达——E=p2/2mE = p^2/2m
  4. 散射振幅直接与S矩阵的矩阵元关联

在位置表象中,自由粒子的薛定谔方程为:

22m2ψ(r)=Eψ(r)-\frac{\hbar^2}{2m}\nabla^2 \psi(\mathbf{r}) = E \psi(\mathbf{r})

解为平面波 ψp(r)=1(2π)3/2eipr/\psi_{\mathbf{p}}(\mathbf{r}) = \frac{1}{(2\pi\hbar)^{3/2}} e^{i\mathbf{p}\cdot\mathbf{r}/\hbar}

Dirac的特色方法:Dirac不直接在位置表象中求解微分方程,而是在动量表象中将散射问题转化为一个积分方程。这是他的变换理论的自然延伸:既然动量表象中自由粒子是对角的,那么整个问题就简化为研究势场如何在动量本征态之间"耦合"。

3.2 Lippmann-Schwinger方程

Dirac的变换理论自然导出了散射的积分方程形式。在动量表象中,引入完整的格林函数(Green’s function)或传播子 G0(E)G_0(E)

G0(E)=1EH0+iϵG_0(E) = \frac{1}{E - H_0 + i\epsilon}

其中 ϵ0+\epsilon \to 0^+ 是一个无穷小正数,它保证了因果性——只向未来传播。

散射态 ψ(+)|\psi^{(+)}\rangle 满足著名的Lippmann-Schwinger方程

ψ(+)=ϕ+1EH0+iϵVψ(+)|\psi^{(+)}\rangle = |\phi\rangle + \frac{1}{E - H_0 + i\epsilon} V |\psi^{(+)}\rangle

其中:

  • ϕ|\phi\rangle 是入射自由态(具有相同能量 EE
  • ψ(+)|\psi^{(+)}\rangle 是包含出射球面波边界条件的完整散射态
  • 第二项是势场引起的"修正",以积分方程的递归形式出现

这个方程的深刻之处在于:它将微分方程转化为积分方程,自然包含了边界条件+iϵ+i\epsilon 的选择确保了出射波(而非入射波)的物理要求。

graph LR
    subgraph "Lippmann-Schwinger方程"
        A["|ψ⁽⁺⁾⟩"] ==> B["|φ⟩"]
        A ==> C["+ G₀(E+iε)V|ψ⁽⁺⁾⟩"]
        C ==> D[递归结构]
        D ==> E[Born级数展开]
    end
    subgraph 物理图像
        F[入射平面波]
        G[势场散射]
        H[出射球面波]
        F --> G --> H
    end
    style A fill:#1a1a2e,color:#fff
    style B fill:#16213e,color:#fff
    style C fill:#0f3460,color:#fff
    style D fill:#533483,color:#fff
    style E fill:#e94560,color:#fff
    style F fill:#1a1a2e,color:#fff
    style G fill:#16213e,color:#fff
    style H fill:#e94560,color:#fff

3.3 Born近似

如果势场 VV "足够弱",可以对Lippmann-Schwinger方程进行迭代展开:

ψ(+)=ϕ+G0Vϕ+G0VG0Vϕ+|\psi^{(+)}\rangle = |\phi\rangle + G_0 V |\phi\rangle + G_0 V G_0 V |\phi\rangle + \cdots

这就是Born级数。在一级Born近似下,散射振幅为:

f(1)(q)=m2π2pfVpi=m2π2V~(q)f^{(1)}(\mathbf{q}) = -\frac{m}{2\pi\hbar^2} \langle \mathbf{p}_f | V | \mathbf{p}_i \rangle = -\frac{m}{2\pi\hbar^2} \tilde{V}(\mathbf{q})

其中:

  • q=pfpi\mathbf{q} = \mathbf{p}_f - \mathbf{p}_i动量转移
  • V~(q)=d3reiqr/V(r)\tilde{V}(\mathbf{q}) = \int d^3r \, e^{-i\mathbf{q}\cdot\mathbf{r}/\hbar} V(\mathbf{r}) 是势场的傅里叶变换

Born近似的物理直觉:一级散射振幅正比于势场的傅里叶分量,其动量转移由入射和出射动量的差决定。大动量转移对应于势场的精细空间结构,小动量转移对应于大尺度结构。

graph LR
    A["势场 V(r)"] -->|傅里叶变换| B["Ṽ(q)"]
    B --> C["散射振幅
f ∝ Ṽ(q)"]
    C --> D["微分截面
dσ/dΩ = |f|²"]
    A2["短程势
如硬核"] --> B2["宽傅里叶谱
大角度散射"]
    A3["长程势
如库仑势"] --> B3["窄傅里叶谱
小角度散射"]
    style A fill:#1a1a2e,color:#fff
    style B fill:#16213e,color:#fff
    style C fill:#0f3460,color:#fff
    style D fill:#e94560,color:#fff
    style A2 fill:#1a1a2e,color:#fff
    style B2 fill:#533483,color:#fff
    style A3 fill:#1a1a2e,color:#fff
    style B3 fill:#533483,color:#fff

四、色散散射:势场的指纹

4.1 什么是色散散射?

"色散"(Dispersion)这个词在物理学中有悠久历史。在光学中,色散指不同波长的光在介质中以不同速度传播;在量子散射中,色散散射指的是不同能量的入射粒子具有不同的散射行为

从Born近似的结果可以看出:

f(θ)V~(q)=V~(2psinθ2)f(\theta) \propto \tilde{V}(q) = \tilde{V}\left(2p\sin\frac{\theta}{2}\right)

由于 p=2mEp = \sqrt{2mE},对于给定的散射角 θ\theta,动量转移 qq 随能量 EE 变化。因此,即使势场固定,不同能量的粒子"感受"到的也是势场的不同傅里叶分量。

4.2 Ramsauer-Townsend效应

在低能电子-稀有气体原子散射中,实验观测到一个惊人的现象:在特定能量(约0.7 eV对于氩气)附近,散射截面突然跌落到接近零——电子仿佛"透明"地穿过了原子!

这个效应无法用经典力学解释,但在量子力学中是自然的结果。考虑S波(l=0l=0)散射的相移 δ0(E)\delta_0(E)

σtot=4πk2sin2δ0\sigma_{\text{tot}} = \frac{4\pi}{k^2} \sin^2 \delta_0

δ0=nπ\delta_0 = n\pinn 为整数)时,sin2δ0=0\sin^2\delta_0 = 0,总截面为零。这在物理上意味着:入射波和散射波发生了破坏性干涉,使得出射波完全恢复为平面波

Ramsauer-Townsend效应是量子干涉在散射中的最美展示之一。

graph TD
    A["入射波 ψ_in"] --> B[原子势场]
    B --> C["散射波 ψ_sc"]
    C --> D["出射波 ψ_out = ψ_in + ψ_sc"]
    D --> E{"ψ_sc的相位?"}
    E -->|"δ₀ = nπ"| F["相消干涉
σ_tot → 0"]
    E -->|"δ₀ = (n+½)π"| G["相长干涉
σ_tot最大"]
    style A fill:#1a1a2e,color:#fff
    style B fill:#16213e,color:#fff
    style C fill:#0f3460,color:#fff
    style D fill:#533483,color:#fff
    style E fill:#e94560,color:#fff
    style F fill:#e94560,color:#fff
    style G fill:#e94560,color:#fff

五、共振散射:量子世界的共鸣

5.1 Breit-Wigner公式

在第VIII章中,Dirac讨论了散射截面在特定能量下出现尖锐峰的现象——共振散射。这在核物理和原子物理中极为常见:一个入射粒子与靶核形成复合态(compound state),短暂存在后衰变。

共振附近的截面可以用Breit-Wigner公式描述:

σl(E)=4πk2(2l+1)Γ2/4(EER)2+Γ2/4\sigma_l(E) = \frac{4\pi}{k^2} (2l+1) \frac{\Gamma^2/4}{(E - E_R)^2 + \Gamma^2/4}

其中:

  • ERE_R共振能量
  • Γ\Gamma共振宽度(半高全宽,FWHM)
  • τ=/Γ\tau = \hbar/\Gamma 是复合态的寿命

这个公式揭示了量子力学中最深刻的关系之一:

能量-时间不确定性:共振越尖锐(Γ\Gamma 越小),寿命越长(τ\tau 越大)。ΔEΔt\Delta E \cdot \Delta t \sim \hbar

5.2 复平面上的极点

Dirac的变换理论为共振提供了更深刻的数学视角。将散射振幅解析延拓到复能量平面上,共振对应于S矩阵在下半平面的极点:

E=ERiΓ2E = E_R - i\frac{\Gamma}{2}

这个复数极点意味着:

  • 实部 ERE_R:共振的能量位置
  • 虚部 Γ/2-\Gamma/2:指数衰减的速率

当入射能量接近 ERE_R 时,极点对散射振幅的贡献最大,产生尖锐的共振峰。这类似于一个被驱动的阻尼谐振子:当驱动频率接近固有频率时,振幅最大。

graph TB
    subgraph 复能量平面
        A["实轴: 物理能量"]
        B["E = E_R - iΓ/2
S矩阵极点"]
        C["上半平面: 因果性禁止"]
        D["下半平面: 共振态"]
    end
    B -->|接近实轴| E["尖锐共振
长寿命"]
    B -->|远离实轴| F["宽共振
短寿命"]
    style A fill:#1a1a2e,color:#fff
    style B fill:#e94560,color:#fff
    style C fill:#16213e,color:#fff
    style D fill:#0f3460,color:#fff
    style E fill:#e94560,color:#fff
    style F fill:#533483,color:#fff

5.3 散射的相移分析

对于球对称势场(V=V(r)V = V(r)),角动量守恒,不同 ll 的分波可以独立处理。散射振幅展开为:

f(θ)=12ikl=0(2l+1)(e2iδl1)Pl(cosθ)f(\theta) = \frac{1}{2ik} \sum_{l=0}^{\infty} (2l+1) (e^{2i\delta_l} - 1) P_l(\cos\theta)

其中:

  • δl\delta_l 是第 ll 分波的相移
  • PlP_l 是勒让德多项式
  • e2iδle^{2i\delta_l} 是S矩阵的 ll 分波本征值

相移的物理意义:散射只改变了各分波的相位(而非振幅),因为势场是弹性的(能量守恒)。相移 δl\delta_l 捕捉了势场如何"延迟"或"提前"了各分波的出射波。

在低能极限(ka1ka \ll 1aa 是势场范围),只有S波(l=0l=0)贡献:

σtot4πk2sin2δ0\sigma_{\text{tot}} \approx \frac{4\pi}{k^2} \sin^2\delta_0

这正是Ramsauer-Townsend效应和共振散射的低能版本。


六、发射与吸收:碰撞的时间反演

6.1 吸收作为逆散射

散射的逆过程是什么?如果散射是入射粒子与靶相互作用后改变方向,那么逆过程就是靶系统吸收一个粒子,从某个初态跃迁到激发态。

Dirac在第VIII章末尾指出:从变换理论的观点,发射和吸收是同一枚硬币的两面。S矩阵的幺正性要求:

fSfi2=1\sum_f |S_{fi}|^2 = 1

这意味着:如果存在从 i|i\ranglef|f\rangle 的散射过程,必然存在从 f|f\ranglei|i\rangle 的逆过程(吸收)。

6.2 细致平衡原理

更精确地说,在热平衡状态下,细致平衡原理要求:

P(if)Ni=P(fi)NfP(i \to f) \cdot N_i = P(f \to i) \cdot N_f

其中 NiN_iNfN_f 是初态和末态的占据数。这个原理由S矩阵的幺正性和时间反演对称性共同保证。

在光的发射与吸收中,细致平衡导出了爱因斯坦著名的A、B系数关系:一个能自发发射光子的系统,也必然能够吸收相同频率的光子。

6.3 发射与吸收的量子描述

Dirac在第VIII章中将发射和吸收统一在散射框架中:

  • 吸收:光子 k|k\rangle + 原子 g|g\rangle \to 原子 e|e\rangle(光子被"吞没")
  • 发射:原子 e|e\rangle \to 光子 k|k\rangle + 原子 g|g\rangle(光子被"产生")

从S矩阵的角度,这两个过程由同一组T矩阵元描述:

\langle e; 0 | T | g; k \rangle \quad \text{(吸收)} \langle g; k | T | e; 0 \rangle \quad \text{(发射)}

时间反演对称性要求这两个矩阵元在形式上密切相关。

graph LR
    subgraph 散射与吸收的统一
        A["散射
|p⟩ → |p'⟩"] -->|时间反演| B["逆散射
|p'⟩ → |p⟩"]
        B -->|等价| C["吸收
原子+光子 → 激发态"]
        C -->|逆过程| D["发射
激发态 → 原子+光子"]
    end
    style A fill:#1a1a2e,color:#fff
    style B fill:#16213e,color:#fff
    style C fill:#0f3460,color:#fff
    style D fill:#e94560,color:#fff

数值例子

例1:Born近似下Yukawa势的散射截面

问题:高能电子被中性原子散射,原子产生的势可以用Yukawa势近似:

V(r)=V0eμrrV(r) = -V_0 \frac{e^{-\mu r}}{r}

其中 V0=Ze24πϵ0V_0 = \frac{Ze^2}{4\pi\epsilon_0}μ=1/a0\mu = 1/a_0a0a_0 是玻尔半径)。用一级Born近似计算微分散射截面,并对高能极限(q/a0q \gg \hbar/a_0)和低能极限(q/a0q \ll \hbar/a_0)进行讨论。

解答

参数

  • V0=Ze24πϵ0=Z×1.44V_0 = \frac{Ze^2}{4\pi\epsilon_0} = Z \times 1.44 eV·nm =Z×2.307×1028= Z \times 2.307 \times 10^{-28} J·m
  • μ=1/a0=1/(0.529\mu = 1/a_0 = 1/(0.529 Å)=1.89×1010) = 1.89 \times 10^{10} m1^{-1}
  • Z=1Z = 1(氢原子),入射电子能量 E=100E = 100 eV

Born近似散射振幅

f(1)(q)=m2π2V~(q)f^{(1)}(\mathbf{q}) = -\frac{m}{2\pi\hbar^2} \tilde{V}(\mathbf{q})

计算Yukawa势的傅里叶变换:

V~(q)=d3reiqr/(V0eμrr)\tilde{V}(\mathbf{q}) = \int d^3r \, e^{-i\mathbf{q}\cdot\mathbf{r}/\hbar} \left(-V_0\frac{e^{-\mu r}}{r}\right)

利用标准积分公式(球坐标,选择 zz 轴沿 q\mathbf{q} 方向):

d3reμrreiqr/=4π2q2+(μ)2\int d^3r \frac{e^{-\mu r}}{r} e^{-i\mathbf{q}\cdot\mathbf{r}/\hbar} = \frac{4\pi\hbar^2}{q^2 + (\hbar\mu)^2}

因此:

V~(q)=V04π2q2+(μ)2\tilde{V}(q) = -V_0 \cdot \frac{4\pi\hbar^2}{q^2 + (\hbar\mu)^2}

散射振幅:

f(1)(q)=m2π2(4π2V0q2+(μ)2)=2mV0q2+(μ)2f^{(1)}(q) = -\frac{m}{2\pi\hbar^2} \cdot \left(-\frac{4\pi\hbar^2 V_0}{q^2 + (\hbar\mu)^2}\right) = \frac{2mV_0}{q^2 + (\hbar\mu)^2}

微分截面:

dσdΩ=f(1)2=4m2V02[q2+(μ)2]2\frac{d\sigma}{d\Omega} = |f^{(1)}|^2 = \frac{4m^2V_0^2}{[q^2 + (\hbar\mu)^2]^2}

动量转移

q=pfpi=2psinθ2q = |\mathbf{p}_f - \mathbf{p}_i| = 2p\sin\frac{\theta}{2}

其中 p=2mEp = \sqrt{2mE}。对 E=100E = 100 eV 的电子:

p=2×9.11×1031×100×1.602×1019p = \sqrt{2 \times 9.11 \times 10^{-31} \times 100 \times 1.602 \times 10^{-19}}

= \sqrt{2.92 \times 10^{-47}} = 5.40 \times 10^{-24}\text{ kg·m/s} \hbar\mu = 1.055 \times 10^{-34} \times 1.89 \times 10^{10} = 1.99 \times 10^{-24}\text{ kg·m/s}

注意 ppμ\hbar\mu 是同量级的。

小角度散射θ1\theta \ll 1qμq \ll \hbar\mu):

dσdΩ4m2V02(μ)4=4m2V024μ4\frac{d\sigma}{d\Omega} \approx \frac{4m^2V_0^2}{(\hbar\mu)^4} = \frac{4m^2V_0^2}{\hbar^4\mu^4}

这是常数(与角度无关)——小角度散射是各向同性的。

代入数值:

  • m=9.11×1031m = 9.11 \times 10^{-31} kg
  • V0=2.307×1028V_0 = 2.307 \times 10^{-28} J·m
  • =1.055×1034\hbar = 1.055 \times 10^{-34} J·s
  • μ=1.89×1010\mu = 1.89 \times 10^{10} m1^{-1}

dσdΩθ=0=4×(9.11×1031)2×(2.307×1028)2(1.055×1034)4×(1.89×1010)4\frac{d\sigma}{d\Omega}\bigg|_{\theta=0} = \frac{4 \times (9.11 \times 10^{-31})^2 \times (2.307 \times 10^{-28})^2}{(1.055 \times 10^{-34})^4 \times (1.89 \times 10^{10})^4}

=4×8.30×1061×5.32×10561.24×10136×1.28×1041= \frac{4 \times 8.30 \times 10^{-61} \times 5.32 \times 10^{-56}}{1.24 \times 10^{-136} \times 1.28 \times 10^{41}}

= \frac{1.77 \times 10^{-115}}{1.59 \times 10^{-95}} = 1.11 \times 10^{-20}\text{ m}^2 = 1.11 \times 10^{4}\text{ 靶(barn)}

(1 barn =1028= 10^{-28} m2^2

大角度散射θπ\theta \approx \piq2pq \approx 2p):

q = 2p = 1.08 \times 10^{-23}\text{ kg·m/s} q^2 = 1.17 \times 10^{-46}\text{ (kg·m/s)}^2 (\hbar\mu)^2 = (1.99 \times 10^{-24})^2 = 3.96 \times 10^{-48}\text{ (kg·m/s)}^2

注意 q2(μ)2q^2 \gg (\hbar\mu)^2,所以:

dσdΩθ=π4m2V02q4=4m2V02(2p)4=m2V024p4\frac{d\sigma}{d\Omega}\bigg|_{\theta=\pi} \approx \frac{4m^2V_0^2}{q^4} = \frac{4m^2V_0^2}{(2p)^4} = \frac{m^2V_0^2}{4p^4}

=8.30×1061×5.32×10564×(2.92×1047)2=4.42×101163.41×1093=1.30×1023 m2=13 barn= \frac{8.30 \times 10^{-61} \times 5.32 \times 10^{-56}}{4 \times (2.92 \times 10^{-47})^2} = \frac{4.42 \times 10^{-116}}{3.41 \times 10^{-93}} = 1.30 \times 10^{-23}\text{ m}^2 = 13\text{ barn}

物理讨论

  • 小角度截面(~10⁴ barn)远大于大角度截面(~10 barn)——这是屏蔽库仑势的典型特征
  • Yukawa势的有限力程(1/μ=a0\sim 1/\mu = a_0)使得大角度散射被强烈抑制,因为高能粒子只有非常靠近核才能感受到未屏蔽的库仑势
  • EE \to \inftypp \to \infty),截面整体趋于零——高能粒子"感受不到"短程势的存在
  • 如果令 μ0\mu \to 0(纯库仑势),Born近似给出Rutherford公式:dσ/dΩ1/q4d\sigma/d\Omega \propto 1/q^4,但在小角度发散——这是长程势的特征

例2:共振散射的Breit-Wigner参数

问题:低能中子被原子核散射,在 ER=1.2E_R = 1.2 eV 处观测到一个尖锐的共振峰。实验测得共振处的总截面为 σres=2.5×1020\sigma_{res} = 2.5 \times 10^{-20} m²(25000 barn)。估算共振宽度 Γ\Gamma 和复合核的寿命 τ\tau。并计算在 E=1.25E = 1.25 eV(略高于共振能量)处的总截面。

解答

已知参数

  • ER=1.2E_R = 1.2 eV =1.2×1.602×1019=1.922×1019= 1.2 \times 1.602 \times 10^{-19} = 1.922 \times 10^{-19} J
  • σres=2.5×1020\sigma_{res} = 2.5 \times 10^{-20}
  • 中子质量 mn=1.675×1027m_n = 1.675 \times 10^{-27} kg

波数计算

kR=2mnER2=2×1.675×1027×1.922×1019(1.055×1034)2k_R = \sqrt{\frac{2m_n E_R}{\hbar^2}} = \sqrt{\frac{2 \times 1.675 \times 10^{-27} \times 1.922 \times 10^{-19}}{(1.055 \times 10^{-34})^2}}

=6.44×10461.11×1067=5.80×1021=7.62×1010 m1= \sqrt{\frac{6.44 \times 10^{-46}}{1.11 \times 10^{-67}}} = \sqrt{5.80 \times 10^{21}} = 7.62 \times 10^{10}\text{ m}^{-1}

共振宽度估算

在共振峰处(E=ERE = E_R),Breit-Wigner公式给出最大截面:

σmax=4πkR2(2l+1)\sigma_{max} = \frac{4\pi}{k_R^2} (2l+1)

假设S波共振(l=0l=02l+1=12l+1=1):

\sigma_{max}^{(理论)} = \frac{4\pi}{(7.62 \times 10^{10})^2} = \frac{12.57}{5.81 \times 10^{21}} = 2.16 \times 10^{-21}\text{ m}^2

但实验测得 σres=2.5×1020\sigma_{res} = 2.5 \times 10^{-20} m²,比理论最大值大——这意味着:

  1. 要么不是纯弹性散射(有吸收/非弹性通道)
  2. 要么 l>0l > 0

实际上,核物理中的共振通常涉及多个通道,Breit-Wigner公式需要修正。但为简化,我们假设这是弹性散射且用等效公式:

σres=4πkR2Γ2/4(EER)2+Γ2/4E=ER=4πkR2\sigma_{res} = \frac{4\pi}{k_R^2} \frac{\Gamma^2/4}{(E-E_R)^2 + \Gamma^2/4}\bigg|_{E=E_R} = \frac{4\pi}{k_R^2}

这与理论最大值一致。如果实验值更大,可能是多个共振叠加或非弹性通道贡献。我们改用另一种方式:假设实验截面峰值就是 4π/k24\pi/k^2,然后 kk 略有不同。

反过来,如果接受 σres=2.5×1020\sigma_{res} = 2.5 \times 10^{-20} m² 为峰值,则对应的等效 keffk_{eff}

4πkeff2=2.5×1020\frac{4\pi}{k_{eff}^2} = 2.5 \times 10^{-20}

keff=4π2.5×1020=5.03×1020=2.24×1010 m1k_{eff} = \sqrt{\frac{4\pi}{2.5 \times 10^{-20}}} = \sqrt{5.03 \times 10^{20}} = 2.24 \times 10^{10}\text{ m}^{-1}

这与计算的 kR=7.62×1010k_R = 7.62 \times 10^{10} m⁻¹ 不同,说明共振可能涉及 l>0l > 0 或非弹性散射。

更实际的估算:在核物理中,中子共振宽度 Γ\Gamma 通常在 meV 到 eV 量级。让我们从截面公式反推:

假设在 E=ERE = E_R 处,截面达到峰值。在 E=ER±Γ/2E = E_R \pm \Gamma/2 处,截面降到一半:

σ(ER±Γ/2)=12σres\sigma(E_R \pm \Gamma/2) = \frac{1}{2}\sigma_{res}

如果实验能分辨共振峰,测得半高宽 Γexp0.05\Gamma_{exp} \approx 0.05 eV(典型值),则:

寿命

τ=Γ=1.055×10340.05×1.602×1019=1.055×10348.01×1021\tau = \frac{\hbar}{\Gamma} = \frac{1.055 \times 10^{-34}}{0.05 \times 1.602 \times 10^{-19}} = \frac{1.055 \times 10^{-34}}{8.01 \times 10^{-21}}

τ=1.32×1014 s=13.2 fs\tau = 1.32 \times 10^{-14}\text{ s} = 13.2\text{ fs}

E=1.25E = 1.25 eV 处的截面

ΔE=EER=0.05 eV\Delta E = E - E_R = 0.05\text{ eV}

如果 Γ=0.05\Gamma = 0.05 eV,则 ΔE=Γ\Delta E = \Gamma

σ(E)=σresΓ2/4(ΔE)2+Γ2/4=2.5×1020(0.05)2/4(0.05)2+(0.05)2/4\sigma(E) = \sigma_{res} \cdot \frac{\Gamma^2/4}{(\Delta E)^2 + \Gamma^2/4} = 2.5 \times 10^{-20} \cdot \frac{(0.05)^2/4}{(0.05)^2 + (0.05)^2/4}

=2.5×10206.25×1042.5×103+6.25×104=2.5×10200.0006250.003125= 2.5 \times 10^{-20} \cdot \frac{6.25 \times 10^{-4}}{2.5 \times 10^{-3} + 6.25 \times 10^{-4}} = 2.5 \times 10^{-20} \cdot \frac{0.000625}{0.003125}

=2.5×1020×0.2=5.0×1021 m2=5000 barn= 2.5 \times 10^{-20} \times 0.2 = 5.0 \times 10^{-21}\text{ m}^2 = 5000\text{ barn}

物理讨论

  • 复合核寿命约13飞秒——这是核反应的特征时间尺度
  • 在共振峰附近偏离 1个宽度(ΔE=Γ\Delta E = \Gamma),截面降到峰值的20%
  • 核共振宽度与原子光谱线宽(如氢原子跃迁的 natural linewidth 107\sim 10^{-7} eV)相比极其巨大——核相互作用远强于电磁相互作用
  • 中子共振在核反应堆设计中有重要应用:在特定能量窗口,中子截面极大,影响链式反应效率

本章总结

Dirac的第VIII章将碰撞问题从经典轨道的语言彻底转换为量子振幅的语言。我们学习了:

  1. 经典散射图像:碰撞参数 bb、散射角 θ(b)\theta(b)、经典微分截面 dσ/dΩ=bdb/dθ/sinθd\sigma/d\Omega = b|db/d\theta|/\sin\theta。当 db/dθ=0db/d\theta = 0 时,经典截面发散——这是量子干涉必须介入的信号。

  2. 散射系数与S矩阵:碰撞问题被重新定义为从入射自由态到出射自由态的幺正变换,概率守恒和能量守恒内嵌于S矩阵的结构之中。

  3. 动量表象的优势:Dirac选择动量表象作为散射的主战场——自由粒子在动量表象中是对角的,散射振幅直接与势场的傅里叶变换关联,Born近似提供了计算弱势场散射的系统方法。

  4. Lippmann-Schwinger方程:将散射问题转化为积分方程,+iϵ+i\epsilon 的因果性选择确保了出射波边界条件的正确实现。

  5. 色散与共振:Ramsauer-Townsend效应展示了量子干涉在散射中的威力,Breit-Wigner公式揭示了共振的寿命-能量宽度关系 τ=/Γ\tau = \hbar/\Gamma,复能量平面上的极点为共振提供了数学本质。

  6. 分波分析与相移:球对称势场中的散射分解为角动量本征态的独立贡献,相移 δl\delta_l 捕捉了势场对各分波的相位影响。低能极限下只有S波贡献。

  7. 发射与吸收的对偶:S矩阵的幺正性将散射、吸收和发射统一到同一个变换理论框架中,细致平衡原理是时间反演对称性的宏观体现。

mindmap
  root("(#quot;碰撞问题
第VIII章#quot;)")
    经典基础
      碰撞参数 b
      散射角 θ(b)
      经典截面
      彩虹与 glory
    散射系数
      S矩阵
      幺正性
      能量守恒
      T矩阵
    动量表象
      Born近似
      傅里叶变换
      Lippmann-Schwinger方程
      Dirac的变换方法
    色散散射
      Ramsauer-Townsend效应
      相消干涉
      低能极限
    共振散射
      Breit-Wigner公式
      复能量极点
      寿命与宽度
      核共振
    分波分析
      相移δ_l
      角动量分解
      S波主导
      勒让德多项式
    发射与吸收
      细致平衡
      时间反演
      爱因斯坦系数
      散射-吸收统一

练习与思考

1. Born近似下的Yukawa势散射

考虑屏蔽库仑势(Yukawa势):V(r)=V0eμrrV(r) = V_0 \frac{e^{-\mu r}}{r}。用一级Born近似计算散射振幅 f(q)f(\mathbf{q}),并讨论当 μ0\mu \to 0 时(纯库仑势)的极限行为。为什么纯库仑势的散射截面会发散?这与卢瑟福实验有什么关系?

2. 相移与束缚态的关系

Levinson定理指出:S波相移在零能量处的值 δ0(0)\delta_0(0) 与势场支持的束缚态数目 nbn_b 满足 nb=δ0(0)/πn_b = \delta_0(0)/\pi(在一定条件下)。请从物理上解释为什么束缚态的存在会"消耗"相移?提示:考虑一个从吸引势逐渐加深的过程,当一个虚能级穿过零点变成束缚态时,相移会发生什么变化?

3. 共振散射与不稳定粒子

在现代粒子物理中,许多"粒子"(如 ρ\rho 介子、Δ\Delta 重子)实际上是共振态。它们具有确定的质量 MM 和衰变宽度 Γ\Gamma,因此寿命 τ=/Γ\tau = \hbar/\Gamma。请从S矩阵的极点结构出发,讨论为什么这些共振态可以被视为"粒子"——它们与稳定粒子的本质区别是什么?在量子场论中,这种描述如何自然涌现?

4. 经典与量子截面的对比

对于硬球势:

V(r)={r<a0r>aV(r) = \begin{cases} \infty & r < a \\ 0 & r > a \end{cases}

(a) 计算经典微分截面。提示:碰撞参数 bab \leq a 的粒子被散射,散射角与 bb 的关系为 θ=2arccos(b/a)\theta = 2\arccos(b/a)

(b) 用S波(l=0l=0)近似计算低能量子总截面,并与经典截面 σcl=πa2\sigma_{cl} = \pi a^2 比较。

© 讨论为什么低能量子截面是经典截面的4倍(这是量子力学的著名结果)。


"散射不是偏转,而是波函数的重新编码。"

✍️ "深空五号"的数据还在回传。苏明的散射模型成功预测了那片尘埃带的穿透概率——在特定能量窗口,飞船几乎可以自由穿行;在另一些能量,则必须绕行。量子力学从不给出确定的轨迹,但它给了人类在不确定性中导航的精确工具。🖤