📚 课程目录第4章非线性波动方程

4、非线性波动方程：从麦克斯韦方程组推导非线性波动方程、慢变包络近似（SVEA）、耦合波方程的基本形式

非线性光学与线性光学的根本区别在于：介质的极化强度 P 不再与电场强度 E 保持简单的线性关系。当入射光强足够高（例如来自锁模激光器的超短脉冲或聚焦后的连续光），介质中会感应出非线性极化项 P⁽²⁾、P⁽³⁾ 等，这些项成为新频率光波（和频、差频、倍频、参量振荡等）的“源”。要定量描述这些波的产生与相互作用，必须从麦克斯韦方程组出发，推导出适用于非线性介质的波动方程，并在实际计算中引入慢变包络近似（SVEA），最终得到描述不同频率分量之间能量交换的耦合波方程。本节将沿着这一逻辑链条，逐步展开推导过程，并给出工程应用中的关键判据与常见误区。

4.1 核心原理：从麦克斯韦方程组到非线性波动方程

非线性光学的一切理论计算都建立在经典电磁理论基础之上。我们首先写出无自由电荷、无自由电流、非磁性介质中的麦克斯韦方程组：

∇ · D = 0
∇ · B = 0
∇ × E = -∂B/∂t
∇ × H = ∂D/∂t

其中电位移矢量 D = ε₀E + P，磁感应强度 B = μ₀H（非磁性介质）。将 D 的表达式代入 ∇ × H 方程，再对两边取旋度，并利用 ∇ × (∇ × E) = ∇(∇·E) - ∇²E，在均匀介质中 ∇·E = 0，即可得到：

∇²E - (1/c²) ∂²E/∂t² = μ₀ ∂²P/∂t²

这里 c = 1/√(ε₀μ₀) 是真空光速。关键一步在于将极化强度 P 分解为线性部分 P^(L) 和非线性部分 P^(NL)：

P = P^(L) + P^(NL)
P^(L) = ε₀χ^(1) · E
P^(NL) = ε₀χ^(2) : EE + ε₀χ^(3) ⋮ EEE + …

将线性部分并入左边，定义介质中的介电常数 ε(ω) = ε₀[1 + χ^(1)(ω)]，折射率 n(ω) = √(ε(ω)/ε₀)，则波动方程改写为：

∇²E - (n²/c²) ∂²E/∂t² = μ₀ ∂²P^(NL)/∂t²

这个方程就是非线性波动方程的基本形式。它的物理意义非常明确：左边描述光波在线性介质中的传播（色散、衍射），右边非线性极化项 P^(NL) 作为“源”驱动新频率分量的产生。如果没有右边，方程退化为线性波动方程，各频率分量独立传播；有了右边，不同频率的波之间就发生了能量耦合。

工程要点：在实际计算中，我们通常将电场写成多个频率分量的叠加，每个分量具有独立的振幅、相位和传播常数。非线性极化项则会产生频率组合（如 ω₁ + ω₂、2ω₁、ω₁ - ω₂ 等），这些组合频率项正是新光波产生的根源。例如在倍频过程中，基频光 ω 通过 χ⁽²⁾ 产生 2ω 的极化分量，该分量作为源项驱动倍频光的增长。

4.2 慢变包络近似（SVEA）的引入与判据

非线性波动方程是一个二阶偏微分方程，直接求解非常困难。但在绝大多数实际场景中，光波的振幅和相位在空间上变化缓慢——在一个波长量级的距离内，振幅的相对变化远小于 1。这就是慢变包络近似的物理基础。

设沿 z 方向传播的平面波电场为：

E(z,t) = A(z) e^(i(kz - ωt)) + c.c.

其中 A(z) 是慢变的复振幅包络，k = nω/c 是波数。将上式代入非线性波动方程，并利用 SVEA 条件：

|∂²A/∂z²| ≪ |2ik ∂A/∂z|

即振幅的二阶空间导数远小于一阶导数乘以波数。这个条件等价于振幅在传播方向上的变化尺度远大于波长。在此近似下，二阶导数项可以忽略，波动方程降阶为一阶常微分方程：

∂A/∂z = (iμ₀ω²)/(2k) · P^(NL)(z) e^(-ikz)

SVEA 的适用性判据通常用“增益长度”或“相互作用长度”来衡量。对于典型非线性晶体（如 BBO、LBO），相互作用长度 L 通常为毫米到厘米量级，而波长在可见光波段为数百纳米，因此 L/λ ~ 10⁴，SVEA 条件严格满足。但在超短脉冲（飞秒量级）与强聚焦条件下，需要考虑更高阶的修正项。

物理量	典型值	SVEA 适用性
连续波倍频（CW SHG）	L = 1 cm, λ = 1 μm	完全适用
飞秒脉冲参量放大	L = 1 mm, λ = 800 nm	需谨慎，考虑高阶色散
波导中的非线性相互作用	L = 1 cm, 有效模面积小	适用，但需考虑模式重叠
强聚焦球面波	焦深 ~ 10 μm	不适用，需采用更精确模型

4.3 耦合波方程的基本形式与推导

耦合波方程是描述多个频率分量之间能量交换的定量工具。以最常见的三波混频（和频产生 SFG、差频产生 DFG、光学参量放大 OPA）为例，考虑三个频率分量 ω₁、ω₂、ω₃，满足能量守恒 ω₃ = ω₁ + ω₂。将总电场写为三个分量的叠加：

E(z,t) = Σ A_j(z) e^(i(k_j z - ω_j t)) + c.c.,  j = 1,2,3

将非线性极化 P^(NL) = ε₀χ⁽²⁾E² 展开，保留与各频率分量相位匹配的项。例如，对于 ω₃ 分量，其源项来自 ω₁ 和 ω₂ 的和频：

P^(NL)(ω₃) = ε₀χ⁽²⁾ A₁ A₂ e^(i(k₁+k₂)z)

代入 SVEA 后的波动方程，得到：

dA₃/dz = (iω₃²)/(2k₃c²) · χ⁽²⁾ A₁ A₂ e^(iΔkz)

其中 Δk = k₁ + k₂ - k₃ 是相位失配量。类似地，可以得到 A₁ 和 A₂ 的方程。写成对称形式：

dA₁/dz = iκ A₃ A₂* e^(-iΔkz)
dA₂/dz = iκ A₃ A₁* e^(-iΔkz)
dA₃/dz = iκ A₁ A₂ e^(iΔkz)

这里 κ 是耦合系数，与非线性极化率 χ⁽²⁾、频率和折射率有关：κ = (ω_j χ⁽²⁾)/(2 n_j c)。这组一阶常微分方程组就是耦合波方程的基本形式。它清晰地揭示了三个波之间的能量流动关系：A₃ 的增长依赖于 A₁ 和 A₂ 的乘积，而 A₁ 和 A₂ 的变化又反过来受 A₃ 的驱动。

记忆技巧：耦合波方程的结构可以类比“三体互动”——每个波的振幅变化率正比于另外两个波的振幅乘积，相位因子 e^(±iΔkz) 决定了能量流动的方向和效率。当 Δk = 0（完美相位匹配）时，能量持续从泵浦光流向信号光和闲频光；当 Δk ≠ 0 时，能量会在传播过程中来回振荡，转换效率受到限制。

4.4 案例拆解：倍频过程中的耦合波分析

以最典型的倍频（SHG）为例，ω₁ = ω₂ = ω，ω₃ = 2ω。此时耦合波方程简化为两个方程：

dA_ω/dz = iκ A_2ω A_ω* e^(-iΔkz)
dA_2ω/dz = iκ A_ω² e^(iΔkz)

其中 Δk = 2k_ω - k_2ω。假设泵浦光强度远大于倍频光（小信号近似），A_ω 近似不变，则第二个方程可直接积分：

A_2ω(L) = iκ A_ω² L · sinc(ΔkL/2) e^(iΔkL/2)

倍频效率 η = |A_2ω(L)|² / |A_ω(0)|² = (κ²L²) · sinc²(ΔkL/2)。这个结果揭示了两个重要结论：

相位匹配至关重要：当 Δk = 0 时，sinc²(0) = 1，效率随 L² 增长；当 Δk ≠ 0 时，效率随 L 振荡，最大值仅为 1/(Δk/2)²，远小于相位匹配情形。
相干长度概念：定义相干长度 L_c = π/Δk，超过这个长度后，能量开始从倍频光回流到基频光，因此单晶长度不宜超过 L_c。

实际工程中，常用角度调谐或温度调谐来实现相位匹配。例如 KDP 晶体在 1064 nm 倍频时，通过旋转晶体角度使 n_ω 和 n_2ω 满足 n_2ω = n_ω，从而 Δk = 0。若角度偏差 0.1°，Δk 就会显著增大，倍频效率可能下降一个数量级。

常见误区提醒：初学者常误以为“晶体越长，倍频效率越高”。实际上，只有在完美相位匹配条件下效率才随 L² 增长。一旦存在相位失配，晶体长度超过相干长度后效率反而下降。此外，SVEA 忽略了二阶导数项，但在强聚焦或超短脉冲场景下，衍射项和色散项可能不可忽略，需要引入更完整的非线性薛定谔方程（NLSE）来描述。

4.5 常见误区与总结复盘

在本节的学习中，以下几个误区需要特别警惕：

混淆线性与非线性极化：线性极化只改变折射率（色散），不产生新频率；非线性极化才是新频率的源。计算时必须明确区分 P^(L) 和 P^(NL) 的贡献。
误用 SVEA 条件：SVEA 要求振幅变化尺度远大于波长。在波导或微腔中，由于模式限制，振幅变化可能很快，此时需要检查 SVEA 是否仍然成立。
忽略相位匹配：耦合波方程中的 e^(iΔkz) 因子是决定能量流向的关键。没有相位匹配，非线性过程效率极低，甚至无法观测。
混淆耦合系数 κ 的单位：κ 的单位是 m⁻¹·V⁻¹（国际单位制），不同教材可能采用高斯单位制，数值差异很大，计算时务必统一单位体系。