2、光学基础回顾：几何光学与波动光学基础、高斯光束基本性质、ABCD矩阵方法简介

激光谐振腔的设计与模式分析，本质上是对光场在腔内往返传播行为的精确描述与控制。无论是设计一台稳定输出的He-Ne激光器，还是构建一台高功率固体激光谐振腔，工程师都必须首先掌握描述光场传播的三套语言：几何光学（光线追迹）、波动光学（衍射与干涉）以及高斯光束的复参数描述。本章作为整个课程的理论基石，将系统梳理这三套工具的核心概念、适用边界与相互联系，并重点引入ABCD矩阵这一贯穿谐振腔分析全过程的通用数学方法。理解了三者之间的层级关系——几何光学是波动光学在短波长极限下的近似，而高斯光束是波动方程在近轴条件下的本征解——你才能在实际设计中准确选择分析工具，避免因模型误用而导致的灾难性设计错误。

2.1 几何光学基础：光线传播与成像的矩阵化描述

2.1.1 核心原理：光线矢量与变换矩阵

几何光学将光能量传播抽象为携带位置和方向信息的光线。在近轴近似（光线与光轴夹角θ极小，满足sinθ≈θ）下，任意光线在某一横截面处的状态可用一个二维列矢量 [y, θ]ᵀ 完整描述，其中y为光线距光轴的高度，θ为光线与光轴的夹角（以弧度为单位，逆时针为正）。当光线通过一个光学元件或一段自由空间后，其输出矢量与输入矢量呈线性关系，该线性变换可用一个2×2矩阵表示，即ABCD矩阵：

[y₂]   [A  B] [y₁]
[θ₂] = [C  D] [θ₁]

矩阵元具有明确的物理意义：A为光线高度放大率（当θ₁=0时，A=y₂/y₁）；B为等效自由空间距离（当y₁=0时，B=y₂/θ₁）；C为光焦度（当y₁=0时，C=θ₂/θ₁）；D为角度放大率（当θ₁=0时，D=θ₂/θ₁）。对于无源介质，矩阵行列式满足AD-BC=1（折射率均匀时），这一性质在谐振腔稳定性分析中至关重要。

2.1.2 常用光学元件的ABCD矩阵

级联规则极为简单：若光线依次通过矩阵M₁、M₂、…、Mₙ，则总矩阵为 M_total = Mₙ × … × M₂ × M₁（注意矩阵乘法顺序与光线传播顺序相反）。这一规则使得复杂光学系统的分析被简化为矩阵连乘，为后续谐振腔往返矩阵的构建奠定了基础。

2.1.3 案例拆解：两镜腔的光线往返矩阵构建

考虑一个由曲率半径R₁和R₂的两个球面反射镜构成的线性谐振腔，腔长L。光线从镜M₁出发，向M₂传播，经M₂反射后返回M₁，再经M₁反射完成一个往返。单程矩阵为：

M_单程 = M_反射(M₂) × M_传播(L) = [1 0; -2/R₂ 1] × [1 L; 0 1] = [1 L; -2/R₂ 1-2L/R₂]

往返矩阵为单程矩阵的平方：M_往返 = M_单程 × M_单程。实际工程中，我们更关心光线在腔内往返多次后是否仍能保持在镜面附近，即稳定性问题。通过计算往返矩阵的迹（A+D），可以判断谐振腔的稳定性：当|A+D|<2时，光线在腔内往返不会逃逸，腔是稳定的；当|A+D|=2时处于临界状态；当|A+D|>2时腔不稳定。这一判据与后续将介绍的g参数判据完全等价。

2.2 波动光学基础：从亥姆霍兹方程到衍射积分

2.2.1 核心原理：光场复振幅与亥姆霍兹方程

波动光学将光视为电磁波，在均匀介质中，单色光场的复振幅U(x,y,z)满足亥姆霍兹方程：∇²U + k²U = 0，其中k=2π/λ为波数。对于沿z轴传播的光束，若引入慢变近似——即光场沿z方向的变化远慢于横向变化——可推导出近轴波动方程（抛物型方程）：

∂²U/∂x² + ∂²U/∂y² - 2ik·∂U/∂z = 0

这一方程是高斯光束理论的数学基础。其物理意义在于：光场在传播过程中，横向分布与纵向演化相互耦合，衍射效应导致光束横向展宽，而介质或透镜的聚焦效应则抑制展宽。当两者平衡时，就形成了稳定的横向模式——这正是谐振腔中本征模式存在的物理根源。

2.2.2 衍射积分与模式自再现

在谐振腔分析中，最常用的波动光学工具是菲涅耳-基尔霍夫衍射积分。设腔内某一参考面上的光场分布为U₁(x₁,y₁)，经过一个往返传播后到达同一参考面的光场为U₂(x₂,y₂)，则：

U₂(x₂,y₂) = ∬ K(x₂,y₂; x₁,y₁) · U₁(x₁,y₁) dx₁dy₁

其中K为传播核，包含了自由空间衍射和镜面反射的全部效应。谐振腔的本征模式要求光场在腔内往返一次后“自再现”——即U₂与U₁只差一个复常数因子γ：U₂ = γ·U₁。这构成了一个积分本征值问题，其本征函数即为腔的横模（如TEM₀₀、TEM₀₁等），本征值γ的模给出单程损耗，幅角给出单程相移。这一自再现条件将谐振腔模式问题转化为数学上的本征方程求解，是Fox-Li数值迭代法的理论基础。

2.2.3 常见误区：几何光学与波动光学的适用边界

• 误区一：用几何光学计算光斑尺寸。几何光学中光线没有横向尺寸概念，无法给出光斑半径w(z)的定量结果。当腔长与波长之比L/λ>10⁴时，衍射效应不可忽略，必须使用波动光学。
• 误区二：忽略衍射损耗。在短腔或大孔径腔中衍射损耗可能小于0.1%，但在长腔或小孔径腔（如波导激光器）中，衍射损耗可达百分之几甚至更高，直接影响激光阈值和效率。
• 误区三：将几何光学稳定性判据用于所有腔型。几何光学稳定性只回答“光线是否逃逸”，不回答“模式损耗多大”。对于高功率激光器，即使几何稳定，高阶模的衍射损耗也可能导致模式竞争问题。

2.3 高斯光束基本性质：谐振腔模式分析的基石

2.3.1 核心原理：基模高斯光束的完整参数体系

高斯光束是近轴波动方程的一组特解，其基模（TEM₀₀模）的复振幅分布为：

U(r,z) = U₀ · [w₀/w(z)] · exp[-r²/w²(z)] · exp[-ikz - ikr²/(2R(z)) + iφ(z)]

其中三个关键参数构成完整的描述体系：光斑半径 w(z) = w₀·√[1+(z/z_R)²]，描述光场横向能量集中程度；等相面曲率半径 R(z) = z·[1+(z_R/z)²]，描述波前弯曲程度；瑞利长度 z_R = πw₀²/λ，是区分近场（z<<z_R）与远场（z>>z_R）的特征尺度。此外，远场发散角 θ₀ = λ/(πw₀) 是衡量光束方向性的核心指标，共焦参数 b = 2z_R 在谐振腔理论中频繁出现。

2.3.2 高斯光束的q参数与传播规律

为了将高斯光束纳入ABCD矩阵框架，定义复光束参数q(z)：

1/q(z) = 1/R(z) - i·λ/(π·w²(z))

q参数的实部与等相面曲率半径相关，虚部与光斑尺寸相关。高斯光束通过任意ABCD光学系统后，q参数按双线性变换传播：

q₂ = (A·q₁ + B) / (C·q₁ + D)

这一公式是谐振腔模式分析的“万能钥匙”。它同时包含了光斑尺寸和波前曲率的信息，且与几何光学的光线矩阵完全兼容——只需将光线矢量替换为q参数即可。特别地，当光束在腔内往返一次后满足q = (Aq+B)/(Cq+D)时，即实现了模式自再现，由此可解出腔的本征q参数，进而得到光斑尺寸和等相面曲率。

2.3.3 案例拆解：平凹腔的基模参数计算

场景：设计一个氦氖激光器谐振腔，腔长L=300mm，输出镜为平面镜（R₁=∞），后腔镜为凹面镜（R₂=1000mm），波长λ=632.8nm。求基模光斑在平面镜和凹面镜上的尺寸。

分析过程：首先计算往返矩阵。从平面镜出发，经自由空间L到达凹面镜，反射后返回平面镜，再经平面镜反射。往返矩阵为：

M = [1 L; 0 1] × [1 0; -2/R₂ 1] × [1 L; 0 1] × [1 0; 0 1]（平面镜反射矩阵为单位阵）
  = [1-2L/R₂, 2L-2L²/R₂; -2/R₂, 1-2L/R₂]

代入L=0.3m，R₂=1.0m，得A=0.4，B=0.42，C=-2.0，D=0.4。稳定性判据|A+D|=0.8<2，腔稳定。利用自再现条件q = (Aq+B)/(Cq+D)，解得：

1/q = (D-A)/(2B) ± i·√[4-(A+D)²]/(2B)
     = 0 ± i·√[4-0.64]/(0.84) = ±i·2.18

取虚部为正的解，得w² = -λ/(π·Im(1/q)) = 632.8×10⁻⁹/(π×2.18) = 9.24×10⁻⁸ m²，故平面镜上光斑半径w₁≈0.304mm。利用q参数传播公式可求得凹面镜上的光斑半径w₂≈0.335mm。

结果验证：凹面镜上的等相面曲率半径应等于镜面曲率半径R₂=1000mm，计算R(z)=1/Re(1/q₂)可验证一致性。该设计下基模体积约为π·w₁·w₂·L/2≈0.096mm³，增益介质需与之匹配以获得高效模式重叠。

2.4 ABCD矩阵方法：统一框架与工程应用

2.4.1 核心原理：从光线到q参数的统一变换

ABCD矩阵方法的核心价值在于其统一性：同一套矩阵，既可用于几何光学中光线矢量[y, θ]ᵀ的变换，也可用于波动光学中高斯光束q参数的变换。这种统一并非巧合，而是源于近轴波动方程与几何光学方程在数学结构上的同源性。对于谐振腔分析，ABCD矩阵方法提供了从“腔结构参数”到“模式参数”的直接映射路径：

1. 构建往返矩阵：根据腔型结构（镜面曲率、腔长、腔内元件），按传播顺序级联各元件的ABCD矩阵，得到往返矩阵M。
2. 稳定性判断：计算矩阵迹的值|A+D|，与2比较判断腔的稳定性。
3. 求解本征q参数：利用自再现条件q=(Aq+B)/(Cq+D)，解二次方程得到腔内参考面上的q参数。
4. 提取模式参数：从q参数分离出光斑半径w(z)和等相面曲率半径R(z)。
5. 传播到任意位置：利用q参数传播公式，计算腔内任意位置的光斑尺寸。

2.4.2 知识结构表格：ABCD矩阵方法在谐振腔设计中的全流程应用

2.4.3 实战场景：含热透镜效应的固体激光器腔型设计

背景与约束：设计一台Nd:YAG固体激光器，泵浦功率10W时增益介质等效热透镜焦距f_th=500mm，腔长L=400mm，输出镜为平面镜，全反镜曲率半径R=1000mm。热透镜位于增益介质中心（距两镜各200mm）。

初始现象：若忽略热透镜效应，按冷腔设计（R=1000mm，L=400mm），计算得|A+D|=|2-2L/R|=|2-0.8|=1.2<2，腔稳定。但实际运行中激光输出功率不稳定，光斑模式劣化。

分析过程：将热透镜等效为焦距f=500mm的薄透镜插入腔内。往返矩阵需分段构建：从平面镜出发，传播200mm→热透镜→传播200mm→凹面镜反射→传播200mm→热透镜→传播200mm→平面镜反射。计算得：

M_往返 = [1 0; 0 1] × [1 200; 0 1] × [1 0; -1/500 1] × [1 200; 0 1] × [1 0; -2/1000 1] × [1 200; 0 1] × [1 0; -1/500 1] × [1 200; 0 1]

计算得A=0.28，D=0.28，|A+D|=0.56<2，腔仍稳定，但g参数从冷腔的g₁g₂=0.6变为0.28，工作点明显移动。进一步求解q参数，发现平面镜上光斑从冷腔的0.35mm增大到0.52mm，模式体积增大，但热透镜像散可能导致模式畸变。

实际操作与结果：在设计中预留热透镜补偿措施——将全反镜曲率调整为R=800mm，使冷腔g参数略微偏向非稳区（g₁g₂=0.75），当热透镜效应出现时，g参数向稳定区中心移动，实现“动态稳定”运行。实测输出功率稳定性从±8%改善至±2%。

风险与代价：动态稳定设计牺牲了冷腔时的模式体积，导致低泵浦功率时效率下降约5%。需在启动阶段采用渐变泵浦策略，避免冷腔非稳区造成的振荡建立困难。

2.4.4 常见误区与总结复盘

• 误区一：矩阵乘法顺序错误。级联矩阵时，先通过的元件矩阵在右侧，后通过的在左侧。建议在纸上画出光路传播顺序，再按“从右到左”书写矩阵乘积。
• 误区二：忽略矩阵行列式约束。对于无源介质，AD-BC=1是必须满足的。若计算结果不满足，说明某个元件矩阵写错了（如焦距符号、曲率半径正负号）。
• 误区三：将q参数传播公式与光线传播公式混用。光线矢量[y,θ]ᵀ适用于几何光学；q参数适用于高斯光束。两者虽然共用ABCD矩阵，但变换公式不同——光线是线性变换，q参数是双线性变换（含除法）。
• 误区四：认为稳定腔一定低损耗。稳定性只保证光线不逃逸，不保证衍射损耗低。对于小孔径腔或高阶模，衍射损耗可能很大，需用Fox-Li迭代或模式本征值计算精确评估。