📚 课程目录第1章光的波动方程与亥姆霍兹方程

1. 光的波动方程与亥姆霍兹方程：从麦克斯韦到傍轴近似

好，我们正式开始。这一章是整门课的数学地基，说白了就是回答一个问题：光到底是怎么传播的？

我记得刚读研那会儿，导师丢给我一篇关于高阶贝塞尔光束的论文，上来就是亥姆霍兹方程。我当时心里直犯嘀咕——这玩意儿跟麦克斯韦方程组到底啥关系？后来自己推了一遍，才明白其中的逻辑链条。今天我就把这个链条给你拆开揉碎。

先别急着跳代码。咱们从最根本的麦克斯韦方程组说起。在无源、线性、各向同性的介质中，麦克斯韦方程组长这样：

∇ · E = 0
∇ · B = 0
∇ × E = -∂B/∂t
∇ × B = με ∂E/∂t

嗯，这里要注意：μ是磁导率，ε是介电常数。真空中μ₀ε₀ = 1/c²。

怎么推出波动方程？很简单——对第三个方程两边取旋度，再把第四个方程代进去。我习惯这么写：

∇ × (∇ × E) = -∂(∇ × B)/∂t = -με ∂²E/∂t²

左边用矢量恒等式 ∇×(∇×E) = ∇(∇·E) - ∇²E，而 ∇·E = 0，所以：

∇²E - με ∂²E/∂t² = 0

这就是波动方程。磁场B也满足同样的形式。光，本质上就是电磁波。

  核心结论：波动方程 ∇²E - (1/c²) ∂²E/∂t² = 0 是麦克斯韦方程组的直接推论，描述了光在空间和时间上的传播规律。

波动方程是时域的，处理起来有点麻烦。我们更常用的是频域形式——亥姆霍兹方程。

怎么做？假设电场是单色波：

E(r, t) = E(r) · exp(-iωt)

代入波动方程，时间导数项变成 -ω²，整理后得到：

∇²E(r) + k²E(r) = 0

其中 k = ω/c = 2π/λ 是波数。这就是标量亥姆霍兹方程。

我个人习惯把E(r)理解成复振幅——它包含了振幅和相位信息，但不包含时间振荡。这在模拟光束传播时非常方便。

小技巧：实际计算中，我们通常处理的是标量场（忽略偏振），用U(r)表示复振幅。亥姆霍兹方程简化为 ∇²U + k²U = 0。对于高斯光束，这个方程的解就是我们要用的。

亥姆霍兹方程虽然精确，但直接求解很困难。好在大多数激光光束满足一个条件：传播方向上的变化远慢于横向变化。这就是傍轴近似。

怎么理解？你想想看，一束激光沿着z轴传播，它的波前在z方向变化很慢，但在x-y平面变化很快。我们可以把复振幅写成：

U(r) = A(r) · exp(ikz)

其中A(r)是慢变包络。代入亥姆霍兹方程，忽略∂²A/∂z²项（因为变化慢），得到：

∇ₜ²A + 2ik ∂A/∂z = 0

这就是傍轴亥姆霍兹方程，也叫抛物型波动方程。∇ₜ² = ∂²/∂x² + ∂²/∂y² 是横向拉普拉斯算子。

避坑指南：我曾经在模拟紧聚焦光束时直接用了傍轴近似，结果发现焦点附近的场分布完全不对。后来才意识到——傍轴近似要求光束发散角很小（通常 < 30°）。对于高数值孔径的聚焦系统，必须用非傍轴方法。记住：傍轴近似不是万能的。

下面这张图把整个推导逻辑串起来了。我建议你把它存下来，后面每章都会用到。

为了让你看得更清楚，我把四个方程放在一起对比：

我的建议：刚开始做光束模拟时，直接从傍轴亥姆霍兹方程入手。等后面需要处理大角度或紧聚焦情况时，再回头用完整的亥姆霍兹方程。这样学习曲线最平滑。

这一章我们走完了从麦克斯韦到傍轴近似的完整推导链。你只需要记住三件事：