4. 高斯光束基础:波动方程的高斯解、基模高斯光束的表达式与物理意义

好,我们进入第四章。前面三章我们把谐振腔的几何结构、稳定性判据和模式概念理了一遍。从这一章开始,要接触一个真正贯穿整个激光工程的核心工具——高斯光束。

说实话,我早年刚接触激光时,觉得谐振腔就是两个镜子来回反射,光在里面走直线。后来做第一个实际项目,要设计一个长距离自由空间光通信的发射端,才发现光根本不是直线传播的。它有个“腰”,会发散,会弯曲波前。嗯,高斯光束就是描述这个“真实光”的数学工具。

4.1 为什么是高斯?——从波动方程说起

激光的本质是电磁波,所以它必须满足麦克斯韦方程组导出的波动方程。在稳态、单色、标量近似下,我们通常处理的是亥姆霍兹方程:

∇²U + k²U = 0

其中 U 是复振幅,k = 2π/λ 是波数。

这个方程的解有很多种。平面波是一种解,球面波也是一种解。但激光谐振腔里出来的光,既不是严格的平面波,也不是严格的球面波。它有一个有限的光斑尺寸,能量集中在轴线附近,而且横向分布是高斯函数形状。

核心结论:在近轴近似条件下(即光线传播方向与光轴夹角很小),亥姆霍兹方程的一个特解就是高斯光束。这个解同时描述了光场的横向分布和纵向传播特性。

为什么会得到高斯分布?我个人理解是:谐振腔的衍射效应和镜面边界条件共同作用,使得稳定腔内的自再现模必然呈现高斯型横向场分布。你想想看,一个有限孔径的腔,光来回振荡无数次,只有那些能量最集中、衍射损耗最小的模式才能存活下来——这就是基模。

4.2 基模高斯光束的完整表达式

好,我们直接上公式。基模高斯光束的复振幅分布可以写成:

U(r, z) = A₀ · [ω₀/ω(z)] · exp[-r²/ω²(z)] · exp[-i(kz + kr²/2R(z) - ψ(z))]

这个式子看着长,其实拆开看就几个部分。我习惯把它分成振幅项和相位项来理解:

组成部分 表达式 物理含义
振幅因子 A₀ · ω₀/ω(z) 轴上振幅随传播距离变化,能量守恒
横向高斯包络 exp[-r²/ω²(z)] 光强在横截面上呈高斯分布
纵向相位 exp(-ikz) 沿传播方向的平面波相位
波前弯曲相位 exp[-ikr²/2R(z)] 等相位面弯曲,类似球面波
古依相移 exp[iψ(z)] 通过焦点时附加的π相位跃变

我的记忆技巧:高斯光束 = 高斯振幅包络 × 球面波前 × 古依相移。三个因子分别对应“光斑多大”、“波前多弯”、“相位多跳”。

4.3 关键参数及其物理意义

公式里出现了几个参数,它们不是独立的,而是通过一个核心量——瑞利长度 z_R——联系在一起。

4.3.1 光斑半径 ω(z)

光斑半径定义为振幅下降到轴上值的 1/e 处的径向距离。它的变化规律是:

ω(z) = ω₀ · √[1 + (z/z_R)²]

其中 ω₀ 是束腰半径,即 z=0 处的最小光斑。z_R = πω₀²/λ 是瑞利长度。

这个公式告诉我们两件事:

  • 在束腰附近(|z| < z_R),光斑几乎不变,光束近似平行
  • 远离束腰(|z| >> z_R),光斑线性增长,发散角 θ = λ/(πω₀)

我曾经在设计激光测距仪时踩过一个坑:选了一个小束腰的激光器想获得小光斑,结果发散角太大,到目标处光斑比脸盆还大。后来才明白,束腰和发散角是反比关系——你压得越小,它发散得越快。

4.3.2 波前曲率半径 R(z)

R(z) = z · [1 + (z_R/z)²]

这个参数描述等相位面的弯曲程度:

  • 在束腰处(z=0),R → ∞,波前是平面
  • 在远场(z >> z_R),R ≈ z,波前近似球面
  • 在 z = z_R 处,R = 2z_R,曲率半径最小

注意:高斯光束的波前曲率半径 R(z) 和球面波的曲率半径不是一回事。球面波的曲率半径就是到点源的距离,而高斯光束的 R(z) 在束腰处无穷大,在远场才趋近于 z。这个区别在做腔模匹配时特别重要。

4.3.3 古依相移 ψ(z)

ψ(z) = arctan(z/z_R)

这个相位项是高斯光束独有的。它表示高斯光束通过焦点时,相对于平面波多了一个 π 的相位跃变(从 -π/2 到 +π/2)。

我记得第一次做谐振腔设计时,死活算不对纵模间隔,后来发现就是忘了考虑古依相移。在短腔或近共焦腔中,这个相移对模式频率的影响不可忽略。

4.4 知识结构图:高斯光束参数关系

下面这张图把基模高斯光束的核心参数和它们之间的逻辑关系画出来了。我建议你把它当作一个思维导图来记。

基模高斯光束参数关系图 亥姆霍兹方程 近轴近似 基模高斯解 光斑半径 ω(z) 波前曲率 R(z) 古依相移 ψ(z) 束腰 ω₀ 瑞利长度 z_R 发散角 θ 远场 R(z)≈z 通过焦点 Δψ=π 三个参数通过瑞利长度 z_R = πω₀²/λ 相互关联,知道 ω₀ 和 λ 即可计算所有传播特性 束腰越小 → 瑞利长度越短 → 发散角越大(反比关系)

4.5 一个实用的例子:估算你的激光光斑

假设你手头有一个 He-Ne 激光器,波长 632.8 nm,束腰半径 ω₀ = 0.5 mm。我们来算几个关键值:

瑞利长度:z_R = π × (0.5×10⁻³)² / (632.8×10⁻⁹) ≈ 1.24 m
发散角:θ = 632.8×10⁻⁹ / (π × 0.5×10⁻³) ≈ 0.4 mrad
在 10 m 处的光斑半径:ω(10) = 0.5 × √[1 + (10/1.24)²] ≈ 4.04 mm

这意味着什么?你的激光在 10 米外,光斑直径已经接近 8 mm 了。如果你要做精密对准,这个发散量必须考虑进去。

工程经验:做自由空间光通信或激光雷达时,我一般先算瑞利长度。如果链路距离远大于 z_R,光束就进入远场区,光斑随距离线性增大。这时候要么用扩束镜增大束腰(减小发散角),要么接受光斑变大的事实去设计接收口径。

4.6 本章小结

基模高斯光束是激光工程中最基础也最常用的模型。记住三个核心参数就够了:

  • ω(z)——光斑半径,决定能量集中度
  • R(z)——波前曲率半径,决定如何与光学元件匹配
  • ψ(z)——古依相移,决定谐振腔的纵模频率

它们都通过瑞利长度 z_R 联系在一起。理解了这个关系,后面讲谐振腔模式匹配、ABCD 矩阵、高阶模等内容时,你会觉得顺理成章。