3. 高阶厄米-高斯光束：厄米多项式与高斯函数的乘积

好，咱们进入第三讲。厄米-高斯光束，说白了就是给标准高斯光束加上了「空间结构」。我刚开始接触这个的时候，总觉得它很抽象，直到在实验室里亲眼看到那些奇形怪状的光斑——嗯，从那以后我就彻底理解了。

3.1 从数学形式说起

厄米-高斯光束的电场分布，可以写成这样：

E_{mn}(x, y, z) = E_0 · H_m(√2·x/w(z)) · H_n(√2·y/w(z)) · exp[-(x²+y²)/w²(z)] · exp[-ikz - iφ(z)]

这里面，H_m 和 H_n 就是厄米多项式。m 和 n 是模式指数，分别对应 x 方向和 y 方向。

我个人习惯把厄米多项式看作「振荡器」——它决定了光斑在横向上有多少个「瓣」。m=0 时，H_0=1，就是普通高斯；m=1 时，H_1=2x，出现一个节线；m=2 时，H_2=4x²-2，出现两个节线……你想想看，是不是越来越有意思了？

3.2 模式指数 (m,n) 的物理意义

m 和 n 到底代表什么？我用自己的话给你翻译一下：

• m：x 方向上的节线数量（暗纹数量）
• n：y 方向上的节线数量（暗纹数量）
• m+n：决定了光束的「阶数」，阶数越高，光束发散角越大

我在项目中遇到过一个问题：用 (2,2) 模式做激光加工，结果光斑太大，能量密度不够。后来换成 (0,0) 基模，效果立竿见影。所以啊，选对模式指数，比调功率还重要。

模式 (m,n)	节线数	光斑特征	典型应用
(0,0)	0	圆斑，中心最亮	光纤通信、激光打标
(1,0)	1 条竖线	左右两瓣	光学镊子、粒子操控
(0,1)	1 条横线	上下两瓣	光束整形
(1,1)	1 竖 + 1 横	四瓣，像四叶草	量子纠缠、信息编码
(2,0)	2 条竖线	三瓣横排	多通道成像

3.3 光斑图案特征——眼见为实

咱们用 Python 画几个典型模式，看看光斑长什么样。

import numpy as np
import matplotlib.pyplot as plt
from scipy.special import hermite

def hermite_gaussian(m, n, x, y, w=1.0):
    """计算厄米-高斯光束的振幅分布"""
    Hm = hermite(m)
    Hn = hermite(n)
    u = np.sqrt(2) * x / w
    v = np.sqrt(2) * y / w
    envelope = np.exp(-(x**2 + y**2) / w**2)
    return Hm(u) * Hn(v) * envelope

# 生成网格
x = np.linspace(-3, 3, 200)
y = np.linspace(-3, 3, 200)
X, Y = np.meshgrid(x, y)

# 计算 (1,1) 模式
E = hermite_gaussian(1, 1, X, Y)

# 画强度分布
plt.figure(figsize=(6, 5))
plt.imshow(np.abs(E)**2, extent=[-3,3,-3,3], cmap='hot')
plt.colorbar(label='Intensity')
plt.title('Hermite-Gaussian (1,1)')
plt.xlabel('x')
plt.ylabel('y')
plt.show()

运行这段代码，你会看到四个亮斑，中间有个十字暗纹。这就是 (1,1) 模式的典型特征。

3.4 知识结构图

下面这张图帮你理清本章的核心逻辑：

3.5 实战避坑指南

最后，分享几个我这些年积累的经验：

• 模式纯度问题： 实际激光器输出的高阶模式往往不纯，会有模式串扰。我曾经用空间光调制器生成 (2,2) 模式，结果拍出来的光斑总有基模成分。后来加了模式滤波器才解决。
• 数值计算精度： 厄米多项式在阶数高时数值会爆炸。比如 m=20 时，H_20(3) 的值能到 10^12 量级。这时候一定要用归一化版本，或者用递推关系计算。
• 光斑对称性： 如果实验中发现光斑不对称，先检查入射光是否对准。我遇到过好几次，以为是模式问题，结果只是镜子歪了。