4.1 速率方程：连续运行的动力学骨架

4.1 速率方程：连续运行的动力学骨架

本节摘要：速率方程把激光器压缩成两个变量——反转密度 ΔN 与腔内光子数 φ——的一阶微分方程组。泵浦项推 ΔN 上升，受激辐射项把 ΔN 转成 φ，增益饱和自动刹车。解稳态得到阈值与输出功率（第 2 章公式的出处），解暂态得到弛豫振荡与尖峰（调Q的原型）。它是激光物理里性价比最高的一组方程。

前面三章的叙述都是"结果态"：反转建立了、增益钳制了、功率线性了。本节把镜头对准"过程"：这些状态怎么随时间建立、怎么被扰动、怎么失稳。工具只有两个耦合方程，但它们解释的现象横跨稳态效率到毫秒尖峰，值得逐项拆开。

方程组：两个变量讲完一台激光器

四能级系统的速率方程（面积分形式，φ 为腔内总光子数，V 为模体积）：

dΔN/dt = W_p(N_tot − ΔN) − ΔN/τ_f − (σc/V)·ΔN·φ
dφ/dt = σc·ΔN·φ − φ/τ_c + (自发辐射种子项)

三项读法：泵浦 W_p 以略小于 N_tot 的饱和幅度注入反转；τ_f 是上能级荧光寿命，反转自发漏光；最后一对 σc·ΔN·φ 是受激辐射——它让 ΔN 减少、φ 增长，是两个变量之间的"转账"。光子方程里 1/τ_c 是腔损耗（τ_c 是光子寿命），种子项说明激光起振不靠魔法、靠自发辐射漏进腔模的那一点光——呼应 1.2 节"火种"的比喻。

方程虽简，两个极限已经足够定大事。稳态极限（导数为零）给出 φ ∝ (P_pump − P_th)、ΔN = ΔN_th：增益钳制的出处。暂态（扰动 ΔN、φ 偏离稳态）线性化后得到阻尼振荡，即弛豫振荡，频率

f_r = (1/2π)·√[(1/τ_c)(P_pump/P_th − 1)/τ_f]

代入 Nd:YAG 典型值（τ_f = 230 μs，τ_c ≈ 20 ns，泵浦三倍阈值）得 f_r ≈ 360 kHz——与实测的功率噪声谱峰严丝合缝。公式的定性读法：荧光寿命越长、腔损耗越大、越贴近阈值，振荡越慢越狂野；反之越快越稳。光纤激光器 τ_c 短（长腔低损耗），振荡频率上移、阻尼加强，所以光纤器件的功率噪声谱比固体干净得多——同一条公式解释了两个家族的噪声性格差异。

数值演练：亲手复现尖峰

速率方程是非线性的，解析解只覆盖小扰动；想看大信号行为（开泵瞬间的尖峰）得数值积分。下面这段代码用最朴素的欧拉法（步长 0.1 ns，远小于两个时间常数），可完整复现"上电 → 首个巨尖峰 → 阻尼弛豫振荡 → 稳态"的全过程：

# rate_eq_sim.py —— 四能级激光器速率方程数值积分（欧拉法）
# 变量：dn = 反转密度(相对值), ph = 腔内光子数(相对值)
import math

tau_f = 230e-6   # 上能级荧光寿命 s (Nd:YAG)
tau_c = 20e-9    # 光子腔寿命 s
pump = 3.0       # 泵浦速率（3 倍阈值，单位：1/s 归一化后）
sigma_c = 1.0    # 受激项系数（归一化）

dt = 0.1e-9      # 步长 0.1 ns，远小于 tau_c，欧拉法稳定
t_end = 400e-6   # 模拟 400 微秒，覆盖寿命量级
n_steps = int(t_end / dt)

dn, ph = 0.0, 1e-12   # 初始：无反转，仅自发辐射种子
peak = 0.0
for i in range(n_steps):
    t = i * dt
    dn_new = dn + dt * (pump * (1 - dn) - dn / tau_f - sigma_c * dn * ph)
    ph_new = ph + dt * (sigma_c * dn * ph - ph / tau_c)
    if ph_new > peak:
        peak, t_peak = ph_new, t
    dn, ph = dn_new, ph_new
    if i % 200000 == 0:
        print(f"t = {t*1e6:7.1f} us | 反转 {dn:8.4f} | 光子 {ph:10.4e}")

print(f"首个尖峰出现在 t = {t_peak*1e6:.1f} us，峰值光子数 = {peak:.3e}")
# 预期输出形态：t 约 200 us 处出现陡峭尖峰（宽约百纳秒），
# 随后若干个间隔数十微秒的衰减振荡，最终稳定在稳态光子数。
​

跑一遍就能"看见"动力学：反转被泵浦慢推（微秒-毫秒尺度），光子近乎瞬间响应（纳秒尺度）——两个时间尺度差四个数量级，正是这种"慢积累、快释放"让阈值附近天生不安分。调Q的发明就是把这种不安分推到极端：先憋着不让光子长（高损耗腔），让反转攒到远超阈值，再瞬间放开。

阈值以上：增益钳制的动力学读法

第 2 章的"增益钳制"在方程语言里更精确：稳态时 dΔN/dt = 0 解出 ΔN_ss ≈ ΔN_th（当 φ 足够大，泵浦与受激项平衡），反转被"钉"在阈值附近，多出的泵浦功率全额转为光子输出。钳制是饱和刹车（2.1 节 I_sat）的稳态表达——两处讲的是同一件事的时域/稳态两副面孔。理解这点，输出镜透射率实验（换 T 后功率曲线平移）就能预判：T 增大 → τ_c 变短 → 阈值升高 → 斜率略增，工作点沿 P_out 直线族滑动。

尖峰：调Q的物理原型

开泵初期（低Q状态、未起振）反转自由积累，可能超过阈值数倍；一旦某个扰动（模式匹配突变、机械振动）突然降低损耗，光子在几十纳秒内爆发，把超额反转一泻而空——功率尖峰比稳态高几十倍。1960 年代早期所有红宝石激光器都饱受这种"不可控尖峰"之苦，工程师们后来干脆反过来利用它：主动制造"先憋后放"，这就是 4.3 节的调Q。历史上故障先于发明、故障启发发明，速率方程把这条工程史浓缩在两行微分方程里。

稳态解亲手解一遍：三行代数回收第 2 章

速率方程与第 2 章公式的关系，值得用一次完整求解钉死。稳态令两个导数为零：光子方程解出 φ_ss = (σc·ΔN_ss − 1/τ_c)·τ_c，要 φ > 0 即 ΔN_ss > 1/(στ_c) ≡ ΔN_th——阈值反转的直接出处（增益 σΔN 等于腔损耗 1/τc·(1/c)）。反转方程解出 ΔN_ss = W_p N_tot/(W_p + 1/τ_f + σcφ_ss)：当 φ_ss 较大，分母被受激项主导，ΔN_ss 被压回 ΔN_th 附近——增益钳制。两式联立消 ΔN，得 φ_ss ∝ (W_p/W_th − 1)：输出功率随泵浦线性增长、斜率由 σc 与输出耦合决定。第 2 章的三件套（阈值反转、增益钳制、线性输出）至此全部有了动力学出处——速率方程不是另一套理论，是同一套理论的"过程视角"。

FAQ：两个问题

问：三能级系统的速率方程要改哪些地方？ 两处：泵浦项的饱和上限从 N_tot 改成 N_tot/2（基态搬空一半后无法继续净泵浦），且吸收项 B₁₂ρN₁ 不再可忽略（基态始终有人）。结果是没有"干净的线性段"——三能级器件的输出曲线在阈值附近更缓、效率更低，正是 2.2 节三笔账的方程版。

问：弛豫振荡能彻底消除吗？ 阻尼可以加强（提高泵浦倍数、缩短腔、加饱和吸收体阻尼），但物理上无法归零——它是双变量耦合系统的固有模态。工程上真正"消除"的方式是跳出锅外：锁模激光器工作在深度饱和与快速恢复的动态里，弛豫振荡让位于锁模动力学；或用反馈主动抑制（AOM 内腔损耗调制），用于对噪声苛刻的计量光源。

本节要点回顾

• 两个变量（ΔN、φ）的耦合方程组是激光动力学的最小模型，稳态给阈值与斜率，暂态给弛豫振荡与尖峰；
• 弛豫振荡频率 f_r = (1/2π)√[(1/τ_c)(P_pump/P_th − 1)/τ_f]，Nd:YAG 典型数百 kHz，光纤激光阻尼更强噪声更干净；
• 数值积分揭示"慢积累（μs-ms）快释放（ns）"的双时间尺度结构，是理解一切脉冲技术的直觉基础；
• 增益钳制 = 饱和刹车的稳态表达；开泵尖峰是调Q的物理原型。