基模高斯光束的数学描述：核心参数全解析
搞光学仿真的人都知道，基模高斯光束是所有高阶模式的基础。说白了，你如果连 TEM₀₀ 模都算不明白，后面那些高阶模式根本玩不转。我个人习惯先把数学表达式吃透，再动手写代码——这样出错了也知道往哪个方向找。
复振幅表达式：一切从这里开始
基模高斯光束沿 z 轴传播时，它的复振幅可以写成这样：
E(r, z) = E₀ · [w₀ / w(z)] · exp[-r² / w²(z)] · exp[-i·k·z - i·k·r² / (2·R(z)) + i·ζ(z)]看着有点长对吧？别怕，我们拆开看：
- 振幅项：E₀ · [w₀ / w(z)] · exp[-r² / w²(z)] —— 描述光强在横截面上的高斯分布
- 相位项：exp[-i·k·z] —— 沿轴向的平面波相位
- 波前弯曲项：exp[-i·k·r² / (2·R(z))] —— 等相位面的曲率
- 古伊相位：exp[i·ζ(z)] —— 这个后面细说
关键点：r 是径向坐标，z 是传播距离。当 r = 0 时，振幅最大；当 r = w(z) 时，振幅降到中心的 1/e。
我在项目中遇到过一个问题：直接用这个公式算远场分布时，数值溢出。后来发现是 w(z) 在 z 很大时增长太快，需要做归一化处理。嗯，这个坑你们后面也会遇到。
束腰半径 w₀：光束最细的地方
束腰半径是高斯光束最核心的参数。它定义在 z = 0 处，是光束横截面上光强降到中心 1/e² 时的半径。
为什么会有一个最细的地方？因为高斯光束不是理想平面波，它在传播过程中会先汇聚再发散。束腰就是那个转折点。
| 参数 | 符号 | 物理意义 | 典型值（He-Ne激光器） |
|---|---|---|---|
| 束腰半径 | w₀ | z=0 处的光斑半径 | 0.5 mm |
| 波长 | λ | 激光波长 | 632.8 nm |
| 束腰位置 | z=0 | 光束最细截面 | — |
我的习惯：做仿真时，先把 w₀ 和 λ 确定下来。这两个参数定了，其他所有参数都能算出来。省心。
瑞利长度 z_R：近场和远场的分界线
瑞利长度的定义很简单：
z_R = π · w₀² / λ它把光束传播分成两个区域：
- |z| < z_R：近场区，光束近似平行，光斑变化不大
- |z| > z_R：远场区，光束开始明显发散
我记得第一次做激光加工仿真时，没注意瑞利长度，直接把近场公式用在远场，结果算出来的能量密度差了 3 倍。从那以后，我每次都会先算 z_R，再决定用哪个近似公式。
光斑半径 w(z)：光束怎么变粗
光斑半径随传播距离的变化关系：
w(z) = w₀ · √[1 + (z / z_R)²]这个公式很漂亮。当 z = 0 时，w(0) = w₀；当 z = z_R 时，w(z_R) = √2 · w₀；当 z → ∞ 时，w(z) 线性增长。
注意：w(z) 是振幅降到 1/e 的半径，不是光强降到 1/e² 的半径。这两个概念容易搞混。光强半径是振幅半径的 √2 倍。
发散角 θ：光束能打多远
远场发散角定义为：
θ = λ / (π · w₀)单位是弧度。从这个公式能看出什么？束腰越小，发散角越大。想聚焦得很细，代价就是光束很快散开。
我曾经设计过一个自由空间光通信系统，需要把光束打 10 公里远。算了一下，如果 w₀ 取 1 cm，发散角只有 0.02 mrad，到 10 km 处光斑直径才 40 cm。但如果 w₀ 取 1 mm，发散角就变成 0.2 mrad，到 10 km 处光斑直径 4 米——接收端根本收不到信号。
等相位面曲率半径 R(z)：波前长什么样
等相位面的曲率半径：
R(z) = z · [1 + (z_R / z)²]这个函数有几个特点：
- z = 0 时，R(0) → ∞，等相位面是平面（束腰处）
- z = z_R 时，R(z_R) = 2 · z_R，曲率半径最小
- z → ∞ 时，R(z) ≈ z，趋近于球面波
你想想看，为什么束腰处等相位面是平面？因为那里是光束最汇聚的地方，波前刚好被拉平了。过了束腰，波前开始向外弯曲。
古伊相位 ζ(z)：被忽略的相位偏移
这个参数很多人会忽略，但在干涉和共振腔计算中特别重要：
ζ(z) = arctan(z / z_R)从 -π/2 到 π/2 变化。它表示高斯光束相对于平面波的额外相位延迟。在 z = 0 处，ζ = 0；在 z = ±z_R 处，ζ = ±π/4；在 z → ±∞ 处，ζ → ±π/2。
实际影响：在激光谐振腔设计中，古伊相位决定了腔的纵模间隔。忽略它，你算出来的谐振频率可能差几十 MHz。
知识体系总览
下面这张图把基模高斯光束的核心参数和它们之间的关系串起来了：
参数之间的依赖关系
从上面的图和公式可以看出，所有参数都跟 w₀ 和 λ 有关。我建议你记住这个依赖链：
- w₀ 和 λ → 算出 z_R
- w₀ 和 z_R → 算出任意 z 处的 w(z)
- w₀ 和 λ → 算出 θ
- z 和 z_R → 算出 R(z) 和 ζ(z)
仿真小技巧：写代码时，我习惯把 w₀ 和 λ 作为输入参数，其他全部写成函数。这样改一个参数，所有结果自动更新。后面做参数扫描时特别方便。
常见误区提醒
最后，分享几个我踩过的坑：
- 单位混淆：w₀ 用 mm，λ 用 nm，算 z_R 时忘了换算。结果差了 6 个数量级。
- 近远场不分：在 z < z_R 的区域用远场近似公式，误差很大。
- 忽略古伊相位：做干涉仿真时，少了 ζ(z) 项，干涉条纹位置全错。
嗯，把这些参数吃透了，后面学高阶模式就轻松多了。记住，基模是所有高阶模的基石，花时间把这里搞明白，绝对值。