线性光学基础回顾：从麦克斯韦方程组到波传播

非线性光学研究的是强光与物质相互作用时产生的频率变换、折射率改变、参量放大等现象。要理解这些非线性效应，首先必须牢固掌握线性光学的基本框架。本章不是简单复习，而是从“为非线性引入做铺垫”的视角，重新梳理麦克斯韦方程组、线性极化率、折射率色散以及线性介质中的波传播。我们将特别关注哪些线性概念在非线性情形下会被修正、哪些假设会被打破，从而为后续章节的χ⁽²⁾、χ⁽³⁾效应建立清晰的参照系。

2.1 麦克斯韦方程组：光与物质相互作用的总框架

麦克斯韦方程组是经典电磁理论的基石，也是描述光在介质中传播的出发点。在非线性光学中，我们仍然从这组方程出发，但会修改本构关系（即D与E、B与H之间的关系）。线性光学中，电极化强度P与电场E成正比，比例系数为ε₀χ⁽¹⁾；而在非线性光学中，P与E之间呈现幂级数关系，这正是所有非线性频率变换的数学根源。

在均匀、各向同性、非磁性的线性介质中，麦克斯韦方程组可以简化为波动方程：

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

其中P = ε₀χ⁽¹⁾E。将P代入后，得到标准波动方程：

∇²E − (n²/c²)∂²E/∂t² = 0

这里n² = 1 + χ⁽¹⁾。这个方程的解是平面波E = E₀e^{i(kz−ωt)}，其中波矢k = nω/c。值得注意的是，在非线性光学中，波动方程右侧的极化项不再与E成简单正比，而是包含E²、E³等高阶项，这些高阶项相当于“源项”，会辐射出新的频率分量。

2.2 线性极化率与介电常数：频率响应的微观起源

线性极化率χ⁽¹⁾(ω)描述了介质对频率为ω的电场的响应。从微观角度看，它源于电子云在外场作用下的畸变、离子位移以及偶极子取向等机制。χ⁽¹⁾(ω)是复数，实部对应折射率，虚部对应吸收（或增益）。在非线性光学中，我们经常需要知道χ⁽¹⁾在基频和倍频处的取值，因为它们决定了相位匹配条件。

介电常数ε(ω) = ε₀[1 + χ⁽¹⁾(ω)]，折射率n(ω) = √[ε(ω)/ε₀] = √[1 + χ⁽¹⁾(ω)]。对于大多数光学介质，χ⁽¹⁾远小于1（例如石英玻璃χ⁽¹⁾≈0.5），因此n≈1+χ⁽¹⁾/2。这个近似在非线性光学中经常使用，但要注意在强吸收区域（如半导体带边）不成立。

2.3 折射率与色散：从相速度到群速度

折射率n(ω)随频率变化的现象称为色散。在非线性光学中，色散扮演着双重角色：一方面，它导致不同频率的波以不同速度传播，破坏相位匹配；另一方面，通过设计色散曲线（如光子晶体光纤中的反常色散），可以实现超连续谱产生和孤子传输。理解色散的关键是区分相速度v_p = c/n和群速度v_g = c/(n + ω·dn/dω)。

在非线性光学实验中，我们经常遇到以下色散相关参数：

• 群速度色散（GVD）：β₂ = d²k/dω²，描述脉冲展宽或压缩。β₂ > 0为正常色散，β₂ < 0为反常色散。
• 色散长度：L_D = T₀²/|β₂|，其中T₀为脉冲宽度。当介质长度L ≫ L_D时，色散效应显著。
• 非线性长度：L_NL = 1/(γP₀)，其中γ为非线性系数，P₀为峰值功率。当L_NL < L_D时，非线性效应主导。

2.4 线性介质中的波传播：相位积累与干涉

在均匀线性介质中，平面波传播距离L后，相位积累为Δφ = nk₀L，其中k₀ = ω/c为真空波矢。这个简单的相位关系是理解非线性光学中相位匹配的基础。当两束频率不同的光在非线性晶体中相互作用时，产生的极化波以和频或差频频率振荡，其相位取决于参与波的相位之和或之差。

以倍频为例：基频光E(ω)在晶体中传播，产生非线性极化P(2ω) = ε₀χ⁽²⁾E²(ω)。这个极化波以相速度v_p(2ω) = c/n(2ω)传播，同时辐射出倍频光。如果倍频光的相速度与极化波的相速度一致（即n(2ω) = n(ω)），则各点辐射的倍频光相干叠加，输出功率随晶体长度平方增长。否则，倍频光在相干长度L_c = π/Δk内增长，之后因相位失配而衰减。

倍频效率（小信号近似）：
η = (2ω²d_eff²L²)/(ε₀c³n²(ω)n(2ω)) · sinc²(ΔkL/2) · I(ω)

其中 d_eff 为有效非线性系数，I(ω) 为基频光强。
当 Δk = 0 时，sinc²(0) = 1，效率最高。

在真实实验中，我们还需要考虑高斯光束的聚焦、走离效应（walk-off）以及晶体吸收。走离效应发生在双折射晶体中，当基频光和倍频光的偏振方向不同时，它们的能量传播方向不一致，导致有效相互作用长度缩短。解决方法是采用非临界相位匹配（θ = 90°）或使用周期极化晶体。

2.5 案例拆解：从线性色散预测到倍频相位匹配设计

我们以一块5 mm长的BBO（β-BaB₂O₄）晶体为例，设计1064 nm → 532 nm的倍频实验。首先从线性光学参数出发：

1. 查表获取折射率：BBO在1064 nm处n_o = 1.655，n_e = 1.543；在532 nm处n_o = 1.675，n_e = 1.555（Sellmeier方程计算值）。
2. 计算相位匹配角：对于I类（o+o→e）匹配，需要n_e(2ω, θ) = n_o(ω)。利用n_e(θ) = n_o·n_e/√(n_o²sin²θ + n_e²cos²θ)，解得θ ≈ 22.8°。
3. 计算有效非线性系数：d_eff = d_22·cosθ·cos3φ（对于BBO，d_22 ≈ 2.2 pm/V）。取φ = 0°，d_eff ≈ 2.0 pm/V。
4. 估算转换效率：设基频光峰值功率10 kW，聚焦光斑半径50 μm，晶体长度5 mm。代入小信号公式，η ≈ 0.8%·sinc²(ΔkL/2)。若角度调谐精确，ΔkL/2 ≈ 0，η ≈ 0.8%。
5. 检查走离角：BBO在532 nm处的走离角约3.2°，对于5 mm晶体，走离位移约280 μm，大于光斑直径100 μm，因此实际效率会下降。改进方案：采用非临界相位匹配（θ=90°）或缩短晶体。

这个案例说明，线性光学参数（折射率、色散、走离角）直接决定了非线性过程的效率。在实验前，务必使用Sellmeier方程精确计算各波长处的折射率，并评估走离和群速度失配的影响。

2.6 从线性到非线性：极化率展开的桥梁

线性光学中，我们假设P = ε₀χ⁽¹⁾E，这意味着介质对光的响应与光强无关。当光强足够大（通常>10⁸ W/cm²），高阶项开始显现。将极化强度展开为幂级数：

P = ε₀(χ⁽¹⁾E + χ⁽²⁾E² + χ⁽³⁾E³ + …)

其中χ⁽²⁾的量级约为10⁻¹² m/V，χ⁽³⁾约为10⁻²² m²/V²。对于连续波激光，即使聚焦到10 μm光斑，电场强度也仅约10⁶ V/m，χ⁽²⁾E²项约为χ⁽¹⁾E的10⁻⁶倍，可以忽略。但对于调Q脉冲激光（峰值功率MW级），电场可达10⁸ V/m，非线性项变得显著。

从线性到非线性的过渡还体现在折射率上：n = n₀ + n₂I，其中n₂ = 3χ⁽³⁾/(4ε₀cn₀²)。对于石英光纤，n₂ ≈ 2.6×10⁻²⁰ m²/W，这意味着在100 μm²纤芯中传输1 W功率时，折射率变化仅约10⁻¹⁰，但经过数十公里积累后，会产生显著的自相位调制和孤子效应。