色散与波包

Learning objectives / 学习目标

完成本课后，你应当能够：

1. 用 Fourier 积分把局域初态分解成波数分量，并解释 \(\Delta x\Delta k\) 的互易性。
2. 从具体介质方程或离散模型求 dispersion relation \(\omega(k)\)，识别 nondispersive、normal/anomalous dispersion 与截止区。
3. 推导 \(v_p=\omega/k\) 和 \(v_g=d\omega/dk\)，说明它们各自追踪的几何对象和适用条件。
4. 对窄带波包将 \(\omega(k)\) 展开到二阶，定量预测包络中心和宽度随时间变化。
5. 区分 phase、group、energy、signal 和 front velocity，判断超光速或负群速度为何不自动违反因果性。
6. 用量纲、极限、对称性和数值尺度检查离散链、水波及一般波包结果。

Prerequisites / 先修知识

• 行波函数、\(k\)、\(\omega\) 与相位。
• 波动方程及初值问题。
• 线性叠加、拍与调制。
• Gaussian 积分的结论、Taylor 展开和复指数的实部表示。

Scope and conventions / 范围与约定

• 先讨论一维、均匀、时间不变的线性系统；实物理场取复表示的实部。
• Fourier 约定为

\[\psi(x,t)=\frac{1}{2\pi}\int_{-\infty}^{\infty} \tilde\psi(k)e^{i[kx-\omega(k)t]}\,dk.\]

• 取正频率分支；在周期介质中波数只需取一个 Brillouin zone，\(k\) 与 crystal momentum 的细节留到凝聚态课程。
• \(v_g\) 代表窄带、弱损耗、弱失真的包络速度；只有在合适的无损系统中才常与能量速度相同。
• 电磁信号的相对论因果界为真空光速 \(c\)；非相对论有效介质方程的前沿需回到其微观模型判断。

Notation / 符号与单位

Knowledge 1 — Fourier synthesis of localized waves / 局域波的傅里叶合成

Knowledge ID: P02-U04-L03-K01

1.1 单色波不能局域

\(e^{ikx}\) 在整个空间都有相同模长，不能表示有限脉冲。局域波必须叠加一段波数范围：

\[\psi(x,0)=\frac{1}{2\pi}\int\tilde\psi(k)e^{ikx}\,dk.\]

若 \(\tilde\psi(k)\) 只在 \(k_0\) 附近非零，\(\psi\) 是快速载波 \(e^{ik_0x}\) 乘慢变包络。空间越窄，谱通常越宽；这不是量子专属，而是 Fourier 变换的数学结构。

图 1：局域 Gaussian 包络由一簇相近波数的无限平面波合成。

1.2 Gaussian 的宽度互易

取归一化强度具有标准差 \(\sigma_{x0}\) 的 Gaussian 包络：

\[\psi(x,0)\propto \exp\left(-\frac{x^2}{4\sigma_{x0}^2}\right)e^{ik_0x}.\]

其波数谱为

\[\tilde\psi(k)\propto \exp[-\sigma_{x0}^2(k-k_0)^2],\]

从而强度谱标准差满足

\[\boxed{\sigma_x\sigma_k=\frac12}.\]

其他形状一般给不等式 \(\sigma_x\sigma_k\ge1/2\)；常数依 Fourier 约定与“振幅宽度/强度宽度”的定义，引用数值前必须声明约定。

图 2：窄空间包络对应宽波数谱；两图横轴量纲不同，宽度不能直接以图上厘米比较。

Worked Example 1 — Gaussian 脉冲的最小谱宽

已知： 某一维 Gaussian 波包的强度空间标准差 \(\sigma_x=0.250\ \mathrm m\)，采用本课 Fourier 约定。
研究对象与系统： 均匀线性介质中的初始包络；尚不讨论时间演化。
参考系与坐标： 实验室系，\(x\) 沿传播方向。
模型理由： 精确 Gaussian 达到 Fourier 宽度下界。

\[\sigma_k=\frac{1}{2\sigma_x} =\frac{1}{0.500\ \mathrm m} =2.00\ \mathrm{rad,m^{-1}}.\]

三位有效数字。\(\sigma_x\sigma_k=0.500\)，量纲为 1。若要求更窄的 \(\sigma_x\)，谱宽必增大；非 Gaussian 脉冲在相同 \(\sigma_x\) 下通常需要不小于此值的 \(\sigma_k\)。这里的 rad 是相位记账，量纲仍为 1。

Knowledge 2 — Dispersion relations / 色散关系

Knowledge ID: P02-U04-L03-K02

2.1 色散关系来自方程，不是额外公式

把平面波试探解 \(\psi=\psi_0e^{i(kx-\omega t)}\) 代入控制方程，得到 \(\omega\) 与 \(k\) 的关系，称 dispersion relation。例如理想 nondispersive wave equation 给

\[\omega^2=v^2k^2,\]

正频率右行分支为 \(\omega=vk\)，所有波数以同一速度传播，任意波形不展宽。若 \(\omega(k)\) 弯曲，不同 \(k\) 分量相位积累速率不同，介质具有色散。

图 3：同一坐标中示意三类 \(\omega(k)\)；曲线形状必须由具体介质方程决定。

2.2 离散单原子链

质量 \(m\)、间距 \(a\)、最近邻弹簧劲度 \(K\) 的无限链满足

\[m\ddot u_n=K(u_{n+1}-2u_n+u_{n-1}).\]

令 \(u_n=Ue^{i(kna-\omega t)}\)，得到

\[\boxed{\omega(k)=2\sqrt{\frac Km}\left\lvert\sin\frac{ka}{2}\right\rvert}.\]

长波极限 \(ka\ll1\) 时 \(\omega\simeq a\sqrt{K/m}\lvert k\rvert\)，恢复连续 nondispersive wave equation；接近 Brillouin-zone 边界 \(ka=\pi\) 时曲线变平，群速度趋零。

图 4：离散链的 acoustic branch；周期结构给出有限独立波数区间。

Worked Example 2 — 离散链的相速度与群速度

已知： \(a=0.200\ \mathrm m\)，\(\sqrt{K/m}=50.0\ \mathrm{s^{-1}}\)，考察 \(ka=\pi/2\) 的正 \(k\) 模。
研究对象与系统： 无限、无损、最近邻单原子链。
参考系与坐标： 链静止，\(+x\) 沿增大 \(n\)。
模型理由： 平移对称且线性，Bloch 平面波是正常模。

\[\omega=2(50.0)\sin\frac{\pi}{4}=70.7\ \mathrm{rad,s^{-1}}, \qquad k=\frac{\pi}{2a}=7.85\ \mathrm{rad,m^{-1}}.\]

因此

\[v_p=\frac{\omega}{k}=9.00\ \mathrm{m,s^{-1}},\]

\[v_g=\frac{d\omega}{dk} =a\sqrt{\frac Km}\cos\frac{ka}{2} =7.07\ \mathrm{m,s^{-1}}.\]

三位有效数字。长波速度为 \(a\sqrt{K/m}=10.0\ \mathrm{m,s^{-1}}\)，两结果均不超过它且 \(v_g<v_p\)，与曲线变平一致。量纲由 \(a\times\mathrm{s^{-1}}\) 得速度。

Knowledge 3 — Phase and group velocity / 相速度与群速度

Knowledge ID: P02-U04-L03-K03

3.1 两个速度来自曲线的割线与切线

单色相位 \(kx-\omega t=\text{constant}\) 的移动速度为

\[\boxed{v_p=\frac{\omega}{k}}.\]

对中心波数 \(k_0\) 附近的窄带波包，令 \(q=k-k_0\)，一阶展开

\[\omega(k)\simeq\omega_0+\omega'_0q.\]

代回 Fourier 积分后，包络变为 \(F(x-\omega'_0t)\)，故

\[\boxed{v_g=\left.\frac{d\omega}{dk}\right|_{k_0}}.\]

\(v_p\) 是从原点到曲线上一点的割线斜率；\(v_g\) 是该点切线斜率。

图 5：非线性色散曲线上相速度与群速度通常不同。

3.2 两近频分量的包络

两列相近的波相加：

\[\cos(k_1x-\omega_1t)+\cos(k_2x-\omega_2t) =2\cos\left(\frac{\Delta kx-\Delta\omega t}{2}\right) \cos(\bar kx-\bar\omega t).\]

快载波以 \(\bar\omega/\bar k\) 移动，慢包络以 \(\Delta\omega/\Delta k\to d\omega/dk\) 移动。若谱不窄、介质有强吸收或 \(\omega(k)\) 在谱宽内不光滑，一个单一 \(v_g\) 不能忠实描述整个脉冲。

图 6：载波波峰可穿过包络；这不代表介质质点或能量以相速度穿越整个波包。

Worked Example 3 — 深水重力波

已知： 深水表面重力波近似满足 \(\omega^2=gk\)，取波长 \(\lambda=10.0\ \mathrm m\)，\(g=9.81\ \mathrm{m,s^{-2}}\)。
研究对象与系统： 深度远大于 \(\lambda\) 的小振幅、无黏水面波。
参考系与坐标： 水体平均静止，\(+x\) 沿波列传播。
模型理由： 忽略表面张力与有限水深修正。

\[k=\frac{2\pi}{\lambda}=0.628\ \mathrm{rad,m^{-1}}, \qquad v_p=\frac{\omega}{k}=\sqrt{\frac gk}=3.95\ \mathrm{m,s^{-1}}.\]

由 \(\omega=\sqrt{gk}\)，

\[v_g=\frac{d\omega}{dk}=\frac12\sqrt{\frac gk}=1.98\ \mathrm{m,s^{-1}}.\]

三位有效数字。\(v_g=v_p/2\) 与 \(\omega(k)\) 向下弯曲一致。极限上波长越长，两速度越大，但当水深不再远大于波长时必须改用有限水深色散关系，不能把该趋势外推到 \(k\to0\)。

Knowledge 4 — Wave-packet spreading / 波包展宽

Knowledge ID: P02-U04-L03-K04

4.1 曲率控制一阶展宽

在 \(k_0\) 附近展开到二阶：

\[\omega(k)=\omega_0+v_gq+\frac12\omega''_0q^2+O(q^3), \qquad q=k-k_0.\]

常数项推进载波相位，一阶项平移包络，二阶项使不同 \(q\) 获得非线性相位而改变形状。若 \(\omega''_0=0\) 且高阶项可忽略，波包不因色散展宽。

图 7：色散波包的中心按 \(v_g\) 前进，同时宽度增加；图中为随中心移动的相对坐标示意。

4.2 Gaussian 宽度公式

对本课定义的初始 Gaussian，若二阶展开足够，强度标准差演化为

\[\boxed{\sigma_x(t)=\sigma_{x0} \sqrt{1+\left(\frac{\omega''_0t}{2\sigma_{x0}^2}\right)^2}}.\]

特征色散时间

\[t_D=\frac{2\sigma_{x0}^2}{\lvert\omega''_0\rvert}.\]

窄波包的 \(\sigma_k\) 大，往往更容易采样到显著曲率并展宽。若有吸收，谱幅也随 \(k\) 改变，展宽不再只由实 \(\omega''\) 描述。

图 8：短时宽度近似不变，\(t\gg t_D\) 时近似线性增长。

Worked Example 4 — Gaussian 波包传播与展宽

已知： \(\sigma_{x0}=0.500\ \mathrm m\)，\(v_g=3.00\ \mathrm{m,s^{-1}}\)，\(\omega''_0=0.120\ \mathrm{m^2,s^{-1}}\)，传播 \(t=20.0\ \mathrm s\)。
研究对象与系统： 无损一维均匀介质中的窄带 Gaussian 波包。
参考系与坐标： 初始包络中心在 \(x=0\)，\(+x\) 沿 \(v_g\)。
模型理由： 谱内二阶 Taylor 展开足够，高阶色散忽略。

中心位置为

\[x_c=v_gt=(3.00)(20.0)=60.0\ \mathrm m.\]

宽度因子为

\[\frac{\omega''_0t}{2\sigma_{x0}^2} =\frac{(0.120)(20.0)}{2(0.500)^2}=4.80,\]

故

\[\sigma_x=0.500\sqrt{1+4.80^2}=2.45\ \mathrm m.\]

三位有效数字。\(\omega''t\) 与 \(\sigma_x^2\) 都是 \(\mathrm{m^2}\)，根号内无量纲。若 \(\omega''_0\to0\)，宽度回到 \(0.500\ \mathrm m\)；正负曲率在纯二阶 Gaussian 宽度中给相同展宽，但啁啾相位符号不同。

Knowledge 5 — Causality and signal velocity / 因果性与信号速度

Knowledge ID: P02-U04-L03-K05

5.1 五种“速度”不能自动等同

在透明、弱色散的窄带区，\(v_g\) 常近似能量与信号速度；在强吸收或 anomalous dispersion 区，脉冲被重塑，峰值可呈负或超光速群速度，不能把峰的位移当作新信息超越因果前沿。

图 9：重塑后的峰可以领先参考峰，但仍位于 causal front 之后；只有前沿与可操作信息到达时间直接检验因果性。

5.2 因果性来自响应函数的时间顺序

线性介质的响应可写为

\[P(t)=\int_{-\infty}^{\infty}\chi(\tau)E(t-\tau)\,d\tau.\]

因果介质要求 \(\chi(\tau)=0\) 对所有 \(\tau<0\)，即输出不能依赖未来输入。这一条件使频域响应的实部与虚部受到 Kramers–Kronig 型联系，并保证“异常群速度”必须伴随色散、吸收/增益或波形重塑。对相对论电磁场，新信息前沿不能超过 \(c\)。

图 10：因果 impulse response 在激励前严格为零；后续振铃可改变峰值位置，却不能提前产生响应。

一个重要的判断边界

“\(v_g\) 是信息速度”只在窄带、弱失真、可明确追踪同一包络时是良好近似。若群速度大于 \(c\)、为负或谱内变化剧烈，第一反应应是检查：

1. 波包是否窄带；
2. 介质是否吸收或增益；
3. 输入与输出峰是否真是同一特征；
4. 新信息前沿何时首次超过噪声；
5. 微观理论的因果传播上限是什么。

Lesson synthesis / 全课整合

从初态到可观测信号

1. 对初始局域场作 Fourier 分解，得到 \(\tilde\psi(k)\) 与实际谱宽。
2. 从介质方程求每个允许分支的 \(\omega(k)\)，不要先猜速度。
3. 用 \(\omega/k\) 追踪相位，用 \(d\omega/dk\) 追踪窄带包络中心。
4. 用 \(d^2\omega/dk^2\) 和更高阶项估计展宽、啁啾与失真。
5. 若有复 \(k\) 或复 \(\omega\)，同时处理吸收/增益造成的谱重权。
6. 用 causal response 或初值问题确定前沿，再讨论能量和信息到达。

Formula-and-decision table / 公式—决策表

Model boundaries / 模型边界

• 窄带近似： 需要 \(\sigma_k\) 足够小，使 \(\omega(k)\) 在谱内可由低阶 Taylor 展开描述。
• 均匀介质： 空间变化介质需 ray/WKB、mode conversion 或散射分析，简单全局 \(k\) 不一定存在。
• 无损假设： 吸收使 \(k\) 或 \(\omega\) 复数，幅度谱与相位谱同时改变。
• Gaussian 公式： 初始啁啾、非 Gaussian 尾和高阶色散会改变宽度公式与不对称性。
• 连续模型截止： 离散链在 \(ka\sim1\) 偏离 continuum；不能把长波波速外推到 zone edge。
• 因果层级： 有效非相对论方程可能产生数学上立即出现的极小尾部；可观测因果界应由更完整微观理论与探测阈值共同解释。

Misconception clinic / 易错概念诊所

1. “一个波包只有一个波长。” 局域化要求一簇 \(k\)；只能定义中心波长和谱宽。
2. “有色散就是有耗散。” 色散描述相位与频率关系；耗散描述振幅衰减。因果介质中二者有关联，但概念不相同。
3. “\(v_p\) 是介质粒子速度。” \(v_p\) 是相位图样速度，粒子振动速度由场振幅和时间导数决定。
4. “\(v_g\) 永远是信息速度。” 仅在窄带弱重塑情形近似成立；异常区必须看前沿与编码。
5. “\(\omega''<0\) 会使波包收窄。” 对无初始啁啾 Gaussian，宽度依 \(\omega''^2\)；符号决定啁啾方向，不单独决定宽度增减。
6. “群速度为零意味着所有质点静止。” 它表示包络/能量传输的一阶速度为零；正常模内仍可强烈振动。

Practice and retrieval / 练习与检索

1. 局域化（K01）

为什么严格单色平面波不能表示有限长度脉冲？

2. 宽度互易（K01）

本课约定下 Gaussian 包络 \(\sigma_x=2.00\ \mathrm m\)，求 \(\sigma_k\)。

3. 求色散（K02）

把 \(\psi=e^{i(kx-\omega t)}\) 代入 \(\psi_{tt}+\Omega^2\psi-c^2\psi_{xx}=0\)，求正频率分支。

4. 长波极限（K02）

说明单原子链在正波数分支 \(0<ka\ll1\) 时的长波极限。

5. 相速与群速（K03）

若 \(\omega=\alpha k^2\) 且 \(k>0\)，求 \(v_p,v_g\)。

6. 读图（K03）

在 \(\omega(k)\) 的局部极大点，群速度是多少？相速度是否必为零？

7. 展宽时间（K04）

\(\sigma_{x0}=1.00\ \mathrm m\)，\(\lvert\omega''\rvert=0.400\ \mathrm{m^2,s^{-1}}\)。求 \(t_D\)。

8. 无色散检查（K04）

若 \(\omega=vk\)，一个任意初始包络如何演化？

9. 异常群速度（K05）

观察到输出脉冲峰比真空传播的参考峰更早到达，能否据此断言超光速传信？

10. 综合判据（K02–K05）

某谱很宽的脉冲跨越 \(\omega(k)\) 的拐点。为何用中心点单一 \(v_g\) 预测输出可能失败？

Exit ticket / 课末检验

1. 给一条弯曲的 \(\omega(k)\) 图，在任一点用图形标出 \(v_p\)、\(v_g\) 与决定 Gaussian 一阶展宽的量。
2. 用不超过三句话解释：为什么 \(v_g>c\) 的实验结果与“可用新信息超光速发送消息”不是同一命题？

References and source notes / 参考资料与来源说明

1. MIT OCW 8.03SC, Lecture 13 — Dispersive Medium, Phase Velocity, Group Velocity — 色散曲线、相速度和群速度主线。
2. MIT OCW 8.03SC, Lecture 14 — Fourier Transform, AM Radio — Fourier 合成、包络与频谱。
3. MIT OCW 8.03SC, Textbook Chapter 10 — Signals and Fourier Analysis — 局域化、带宽与 wave packet formalism。
4. Caltech Physics 12a — Waves, Dispersion and Group Velocity — 非线性色散、波包与调制的本科推导。
5. MIT 18.336 — Numerical Methods for PDEs, Lecture 8 summary (official course page) — \(\omega(k)\) 展开、群速度、色散展宽及耗散情形警告。
6. MIT OCW 8.07, Lecture Note 17 — Group Velocity and Anomalous Dispersion — 异常群速度与信号速度因果边界。
7. Duke Physics — The Fast-Light Debate — Sommerfeld–Brillouin 前沿、脉冲重塑与信息速度的物理说明。

所有图依据本课方程原创绘制。来源与链接于 2026-08-23 实际打开核对。MIT Open Learning 的入门材料把群速度概括为信息速度；本文根据 MIT 8.07 与 Duke 的强色散/吸收讨论，将该说法严格限制在窄带弱失真区域。

Source-to-content map

复习与拓展

• Unit 2.5 把一般波包速度应用到声脉冲、超声与 Doppler 测量。
• Unit 2.6 将系统研究离散链、Bloch 波、cutoff、band gap 与非线性包络。
• Part V 的光学色散、群延迟和脉冲传播会复用本课因果边界。
• Part VII 的量子自由波包将使用相同 Fourier 数学，但概率解释不同。