高斯光束原理：腰斑、瑞利长度与光束质量 M²
从“激光为什么不是一束完美的平行光”说起
很多人第一次用激光笔时都会有个疑问：激光不是“平行光”吗？为什么光点会随着距离越来越远而越变越大、越来越暗？
答案是：激光确实有很强的方向性，但它不是数学意义上的理想平行光。任何有限尺寸的光束在传播过程中都会因为衍射而自然发散——这是波动光学的必然结果，与激光器做得好不好无关。物理学家发现，激光器输出中一种最常见、也最“标准”的光束形态，可以用非常简洁的数学精确描述：它的横向强度分布是高斯函数，传播时的发散行为可以逐点预测。这种光束就叫高斯光束（Gaussian beam）。
本文要回答的核心问题有三个：描述高斯光束需要哪几个参数，它们之间如何互相制约；真实激光束与理想高斯光束的差距用什么量衡量；以及如何用一个 2×2 的矩阵，快速算出透镜聚焦后的光斑大小。
什么是高斯光束：一条横向强度分布
把一束激光用光束质量分析仪（beam profiler）拍下来，会看到横截面的光斑亮度从中心向边缘平滑衰减，近似呈高斯函数形状。设距光轴的距离为 $r$，归一化光强 $I(r)$ 可以写成：
$$
I(r)=I_0e^{-2r^2/w^2}
$$
其中 $I_0$ 是中心峰值强度，$w$ 是一个有明确物理意义的特征半径：当 $r=w$ 时，光强下降到峰值的 $e^{-2}\approx 13.5%$。这个 $w$ 被称为光束的 $1/e^2$ 半径（也叫束宽），光斑直径记作 $2w$。之所以用 $1/e^2$ 而不是更常见的半高全宽（FWHM）作为“光束尺寸”，是因为高斯光束传播理论的很多公式用 $w$ 表达最简洁，且 $r<w$ 的圆内集中了约 86.5% 的总功率（由 $\int_0^w e^{-2r^2/w^2}2\pi rdr$ 积分可得 $1-e^{-2}\approx 0.865$）。
图 1：高斯光束横截面光强分布。(a) 归一化光强随径向距离 $r/w_0$ 的变化，$w_0$ 为 $1/e^2$ 半径，该处光强降至峰值约 13.5%；(b) 束腰处光斑的二维强度分布，虚线圆标出 $1/e^2$ 半径。
沿传播方向看，束宽 $w$ 不是常数：光束在某处最细，这个最小束宽位置称为束腰（beam waist），记作 $w_0$。束腰通常是激光器输出镜处或透镜焦点附近。整个高斯光束的形状，就是由束腰大小 $w_0$ 和波长 $\lambda$ 决定的——下面一节看它们如何决定传播行为。
传播中的三个关键量：腰斑、瑞利长度与发散角
高斯光束沿传播轴（记作 $z$）的束宽演化有精确的解析式：
$$
w(z)=w_0\sqrt{1+\left(\frac{z}{z_R}\right)^2}
$$
式中出现了一个新的长度量 $z_R$，称为瑞利长度（Rayleigh length）：
$$
z_R=\frac{\pi w_0^2}{\lambda}
$$
$z_R$ 的物理意义非常直观：在 $|z|\le z_R$ 的范围内，束宽最多只增大到 $\sqrt{2}w_0$，也就是光斑面积翻倍；工程上通常把 $2z_R$ 当作光束的“准直距离”。而一旦 $z\gg z_R$，上式近似为 $w(z)\approx w_0z/z_R$，束宽随距离线性增大，相当于一条以束腰为顶点、半角为 $\theta$ 的圆锥：
$$
\theta=\frac{\lambda}{\pi w_0}
$$
这就是高斯光束的远场发散半角。公式里最值得注意的关系是：发散角与腰斑成反比。想要把光束聚焦得很细（$w_0$ 小），必然换来更快的发散；想要光束传得很远（发散小），就必须先把腰斑做大。这是衍射极限给出的硬约束，不是工艺问题。
表 1：按公式计算的高斯光束典型参数示例（$\theta=\lambda/(\pi w_0)$，$z_R=\pi w_0^2/\lambda$）
| 波长 $\lambda$ | 腰斑 $w_0$ | 瑞利长度 $z_R$ | 远场发散半角 $\theta$ |
|---|---|---|---|
| 633 nm（氦氖激光） | 0.5 mm | 约 1.24 m | 约 0.40 mrad |
| 1064 nm（光纤激光） | 1 mm | 约 2.95 m | 约 0.34 mrad |
| 633 nm（聚焦光束） | 5 µm | 约 124 µm | 约 40 mrad（2.3°） |
表 1 的最后一行的意义值得多说一句：同样一束 633 nm 激光，腰斑从 0.5 mm 细到 5 µm，瑞利长度从 1.24 m 缩短到约 124 µm，发散角却放大了 100 倍。这说明“聚焦”和“准直”本质上是同一枚硬币的两面——透镜把光束压得越细，焦点前后能保持细束的范围就越短。
图 2：高斯光束沿传播方向的束宽演化（$\lambda=633$ nm，$w_0=0.5$ mm）。束腰处束宽最小；$\pm z_R$ 范围内束宽不超过 $\sqrt{2}w_0$，近似准直；远离束腰后束宽渐近于两条直线，其夹角即远场发散角 $\theta$。
为什么“恰好”是高斯形：傍轴波动方程的基模
高斯光束不是人为选定的拟合函数，而是波动方程在傍轴近似下的自然解。从亥姆霍兹方程出发，假设光束沿 $z$ 轴传播且横向尺度远大于波长、波矢与 $z$ 轴夹角很小（这就是傍轴近似（paraxial approximation）），可以得到傍轴波动方程；它存在一族正交解——厄米-高斯（Hermite–Gaussian）模，其中最低阶的 TEM00 模（TEM = transverse electromagnetic mode，横电磁模）横向光强分布正是上面写的高斯函数。
这解释了为什么这么多激光器都输出高斯光束：稳定谐振腔的本征模本来就接近厄米-高斯模，而基模 TEM00 的光束质量最好，是激光器设计中最常被选中的模式。也解释了为什么高斯光束的公式如此“好用”：q 参数、ABCD 矩阵这些工具都是围绕它建立的。
高斯光束还有一个容易被忽略的相位细节：它穿过束腰（焦点）时，波前会额外经历一个相位变化，称为 Gouy 相位，$\zeta(z)=\arctan(z/z_R)$。对大多数强度计算它不影响结果，但在干涉、谐振腔模式和锁模等场景中，这个 $\pi$ 量级的相位差会实实在在地改变物理现象（比如谐振腔的纵模间隔）。复振幅完整表达式为：
$$
U(r,z)=\frac{w_0}{w(z)}\exp\left(-\frac{r^2}{w^2(z)}\right)\exp\left[-ikz-\frac{ikr^2}{2R(z)}-i\zeta(z)\right]
$$
其中 $R(z)$ 是波前曲率半径，$k=2\pi n/\lambda$ 为波数。公式虽长，结构却很清晰：第一项是振幅与横向高斯包络，第二项里 $e^{-ikz}$ 是平面波相位，$e^{-ikr^2/(2R)}$ 描述波前弯曲，$e^{-i\zeta}$ 是 Gouy 相位。
真实激光束与光束质量 M²
实验室里真正从激光器出来的光束，横截面强度分布往往不是完美的单高斯函数，可能存在高阶模、像差或衍射导致的畸变，因此发散得比同腰斑的理想高斯光束更快。工程上用光束质量因子 M²（M-squared factor，亦写作 $M^2$）来量化这种差距：理想基模高斯光束的 $M^2=1$，真实光束 $M^2>1$，$M^2$ 越大，相同腰斑下发散越快。发散角公式相应修正为：
$$
\theta=\frac{M^2\lambda}{\pi w_0}
$$
M² 的测量有一套国际标准（ISO 11146）：沿传播方向测量一系列位置上的束宽，得到束宽随距离的演化曲线（caustic），再用双曲线拟合提取腰斑位置、腰斑大小和 M²。注意，M² 和腰斑 $w_0$ 并不是独立给出的——拟合时二者一起确定。
表 2：相同腰斑与波长下，不同 M² 的光束对比（按 $\theta=M^2\lambda/(\pi w_0)$ 计算）
| 光束质量 $M^2$ | 等效瑞利长度 | 远场发散半角 $\theta$（$w_0=1$ mm，$\lambda=1064$ nm） |
|---|---|---|
| 1.0（理想高斯） | $z_R\approx 2.95$ m | 约 0.34 mrad |
| 1.3（常见固体激光器） | $z_R/M^2\approx 2.27$ m | 约 0.44 mrad |
| 1.5（中等光束质量） | $z_R/M^2\approx 1.97$ m | 约 0.51 mrad |
表 2 注：非理想光束的等效瑞利长度按 $z_R/M^2$ 估计，发散角按 $M^2\lambda/(\pi w_0)$ 计算。
图 3：相同腰斑与波长下，理想高斯光束（$M^2=1$）与 $M^2=1.5$ 光束的束宽演化对比（$\lambda=1064$ nm，$w_0=1$ mm）。M² 越大，等效瑞利长度越短，远场发散越快。
M² 的实用价值在于“可测、可比”：激光器厂商会在规格书里标出典型 M²，用户据此判断这束光能聚焦到多小、能准直传多远，而不必逐台测量整个光场。需要说明的是，$M^2$ 描述的是光束的聚焦/发散能力，不直接等于输出功率或单色性，选购激光器时它们要分开看。
透镜与 ABCD 变换：如何计算聚焦光斑
把高斯光束送进透镜、想预知焦点处的光斑大小，最常用的工具是复光束参数 q（complex beam parameter）和 ABCD 矩阵。q 参数把束宽 $w(z)$ 与波前曲率半径 $R(z)$ 打包成一个复数：
$$
\frac{1}{q(z)}=\frac{1}{R(z)}-i\frac{2}{kw^2(z)}
$$
在束腰处波前是平面（$R\to\infty$），于是 $1/q$ 只剩虚部，$q$ 为纯虚数，其模等于 $z_R$。任意光学系统（透镜、反射镜、自由空间段）都可以用一个 2×2 矩阵 $\begin{bmatrix}A&B\cr C&D\end{bmatrix}$ 表示，高斯光束穿过该系统的 q 参数变换满足 ABCD 定律：
$$
q_2=\frac{Aq_1+B}{Cq_1+D}
$$
多个元件顺序排列时，系统矩阵等于各元件矩阵按顺序相乘，然后只做一次 q 变换即可。这让“光路里有五块透镜，焦点光斑多大”这类问题变成纯粹的矩阵乘法，不用逐段做衍射积分。例如一束在透镜处腰斑为 $w_{\mathrm{in}}$ 的准直光束，经焦距 $f$ 的薄透镜聚焦，出射腰斑约为：
$$
w_0’\approx\frac{\lambda f}{\piw_{\mathrm{in}}}
$$
也就是说，聚焦光斑与波长和焦距成正比、与入射束宽成反比——入射光束越粗，聚焦越细。这与“发散角-腰斑”的反比关系同源，都是衍射极限的不同侧面。
图 4：透镜对高斯光束的聚焦示例（$\lambda=1064$ nm，入射腰斑 1 mm，$f=100$ mm）。(a) 聚焦光路示意，出射束腰位于焦点附近；(b) 用 ABCD 定律算出的束宽演化，出射腰斑约 34 µm。
实际工程中，入射光束往往并非恰好“腰在透镜上”，这时用 ABCD 定律先让光束自由传播一段距离、再经过透镜矩阵，就能严格算出出射腰斑的位置和大小。q 参数与 ABCD 方法不仅适用于透镜，也适用于谐振腔稳定性分析、光纤耦合效率估算等场景，是高斯光束工程化的“万能钥匙”。
意义与展望
高斯光束理论把“激光束长什么样、怎么传播”这个看似复杂的问题，压缩成了腰斑、瑞利长度、发散角、M² 和 q 参数这几个量，让光学工程师可以像算几何光学一样快速估算光束行为。它在现实世界的作用远不止理论：激光加工需要把焦点光斑和瑞利长度匹配到材料厚度；长距离光通信要用透镜把出射光束耦合进单模光纤（单模光纤的基模与高斯光束非常接近，耦合效率通常可达 80% 以上）；光学镊子利用聚焦高斯光束的强度梯度捕获微小粒子，光力的大小直接由束腰附近的光强分布决定。即使是近年来热门的贝塞尔光束、涡旋光束，也常常以高斯光束为包络（如 Bessel–Gauss 光束），高斯理论是理解它们的出发点。
从原理验证走向实际应用，仍有不少待推进的方向：真实光束的 M² 测量受采样位置、孔径效应影响，需要标准化流程；聚焦到极细腰斑时，傍轴近似的适用边界开始显现，需要矢量光束理论修正；以及如何用自适应光学等手段实时补偿光束质量退化。这些问题的共同基础，仍然是我们在这篇文章里建立的高斯光束语言。
参考资料
- Rüdiger Paschotta, “Gaussian beams,” RP Photonics Encyclopedia. https://www.rp-photonics.com/gaussian_beams.html
- Rüdiger Paschotta, “Beam quality,” RP Photonics Encyclopedia. https://www.rp-photonics.com/beam_quality.html
- D. Martin-Sanchez, J. Li, E. Z. Zhang, P. C. Beard, J. A. Guggenheim, “ABCD transfer matrix model of Gaussian beam propagation in plano-concave optical microresonators,” Optics Express 31(10), 16523 (2023). DOI: https://doi.org/10.1364/OE.484212（开放全文：https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11146662/）
- J. Gong, S. Zhang, G. Duan, L. Qi, Y. Yang, “Optical force exerted on the two dimensional transition-metal dichalcogenide coated dielectric particle by Gaussian beam,” Heliyon 9(6), e14314 (2023). DOI: https://doi.org/10.1016/j.heliyon.2023.e14314（开放全文：https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10015242/）