企业官网建设流程全解析

简介：《专题一 激光技术》是一份面向光学、物理及光电子相关专业师生的课堂演示课件，适合激光原理入门学习与教学备课使用。课件围绕受激吸收、自发辐射与受激辐射三种跃迁展开，梳理了从爱因斯坦提出受激辐射概念到梅曼研制红宝石激光器的发展简史，并系统讲解粒子数反转、光学谐振腔与阈值条件等核心机制，同时对比三能级与四能级系统的差异，延伸至激光在精密加工、光纤通信、医疗手术与科研测量中的典型应用。压缩包内仅含1个PPT文件，大小约6.23MB，页面配有能级示意图与谐振腔结构图，便于直接用于课堂讲授或个人复习。目前已有85人学习下载，适合希望借助图示快速建立激光技术知识框架的初学者，也可作为课程汇报与考前梳理的参考资料。

1. 一份《激光技术》课件为什么值得逐页拆

前阵子帮一家做激光打标的客户排查工艺漂移，现场工程师对脉宽、重频、单脉冲能量这些参数如数家珍，但问到「为什么 Nd:YAG 要用 808 nm 的泵浦光」「腔长和纵模间隔到底什么关系」，回答就开始含糊了。回来我把手头这份《专题一 激光技术.ppt》重新翻了一遍，它属于「物理学与激光技术」的入门讲义，从 1916 年爱因斯坦提出受激辐射，讲到 1960 年梅曼做出第一台红宝石激光器、1962 年贝尔实验室研制成 He-Ne 气体激光器，再落到三能级/四能级系统、光学谐振腔与阈值条件。课件里没有一句工艺参数，却把工艺参数背后那套账算得清清楚楚。适合两类人：一类是要把激光原理讲给学生听的老师，另一类是天天调激光设备、却说不清增益从哪来的工程师。

2. 从受激吸收到粒子数反转：把三种跃迁写成速率方程

激光的全称是 Light Amplification by Stimulated Emission of Radiation，直译是「受激辐射的光放大」，1964 年 10 月钱学森建议译作「激光」。要把这句话讲透，光画能级图是不够的，得把三种跃迁写成一组装得进代码的方程，才知道泵浦功率为什么必须过某个坎。

2.1 三种跃迁过程与爱因斯坦系数的关系

课件把光子与原子的相互作用拆成三件事。处于基态的原子吸收一个能量恰好等于某两能级之差的光子，跃迁到高能态，这是受激吸收；处在激发态的原子不稳定，停留时间数量级只有 10⁻⁸ s，会自发跳回低能态并放出一个光子，这是自发辐射，霓虹灯就是典型例子，各原子独立发光，频率、相位、方向都不同，所以不相干；如果激发态原子在自发辐射之前，先被一个能量合适的外来光子「撞」了一下，它会被诱导跃迁，同时放出一个与入射光子频率、初相、传播方向完全相同的光子，这就是受激辐射——输入一个光子、输出两个全同光子，两个再诱发四个，链式反应就是光放大。

三者的竞争关系可以用爱因斯坦系数写成速率方程，简写为：

• 受激吸收速率：B12 · ρ(ν) · N1
• 受激辐射速率：B21 · ρ(ν) · N2
• 自发辐射速率：A21 · N2

其中ρ(ν)是辐射场能量密度，N1、N2分别是低、高能级粒子数密度。热平衡下g1B12 = g2B21，A21/B21 = 8πhν³/c³。这组关系说明一件反直觉的事：受激辐射的「本事」和受激吸收完全对称，谁能占上风，唯一取决于N2和N1谁大。

2.2 玻尔兹曼分布：常温下高能级几乎是空的

原子在各能级上的分布服从玻尔兹曼分布，N2/N1 = (g2/g1)·exp(-ΔE/kT)。常温T = 300 K时kT ≈ 0.0259 eV，只要能级差到电子伏特量级，指数项就是天文数字般的倒数。与其空口说「高能级没人」，不如直接把几种常见激光波长的比值算出来。

import numpy as np k_B = 8.617333262e-5 # 玻尔兹曼常数，eV/K HC = 1239.84193 # hc 的常用换算因子，eV·nm def delta_E_eV(wavelength_nm): """由波长反推能级差（eV）""" return HC / wavelength_nm def boltzmann_ratio(wavelength_nm, T=300.0, g2_over_g1=1.0): """给定波长与温度，返回上下能级粒子数之比 N2/N1""" dE = delta_E_eV(wavelength_nm) return g2_over_g1 * np.exp(-dE / (k_B * T)) # 依次为红宝石、He-Ne、Nd:YAG、CO2 激光器的特征波长 for lam in (694.3, 632.8, 1064.0, 10600.0): print(f"{lam:>8.1f} nm dE={delta_E_eV(lam):.3f} eV " f"N2/N1={boltzmann_ratio(lam):.3e}")

跑出来的结果很说明问题：694.3 nm 对应能级差 1.786 eV，N2/N1只有 10⁻³⁰ 量级；1064 nm 是 10⁻²⁰ 量级；只有到 10.6 μm 的中红外，比值才勉强爬到 10⁻²。这解释了为什么热平衡状态下受激吸收永远压着受激辐射，光穿过介质只会被吸收。要让受激辐射占优，必须从外界灌能量把粒子「搬」到高能级去，这个动作就是课件的「激励」或「泵浦」（Pump），搬完之后的非平衡状态就叫粒子数反转。

2.3 三能级与四能级：阈值差一倍还是差几个数量级

课件点出两个硬条件：一是能实现粒子数反转，二是工作物质必须存在亚稳态能级。亚稳态的寿命是毫秒量级，比普通激发态长五个数量级，粒子在这里「待得住」，才轮得到受激辐射来消耗它。不是所有物质都有亚稳态，所以不是所有物质都能当激光工作物质。

红宝石之所以要用闪光灯猛抽，是因为它的下能级就是基态，粒子数反转意味着要把一半以上的「住户」搬到楼上；而 Nd:YAG 是四能级，下能级靠无辐射跃迁几纳秒就排空，反转门槛低得多。这也是为什么实际设备里连续激光器几乎清一色是四能级体系，三能级体系更多出现在脉冲、高能量的场景中。

3. 光学谐振腔与阈值条件：增益怎么盖过损耗

有反转不等于出激光。课件里反复强调的第三个要素是光学谐振腔，以及控制它能否起振的阈值条件。这两块是整套原理里最容易讲成「玄学」的部分，实际上它完全可以写成几行公式和几行代码。

3.1 谐振腔做的三件事：方向性、放大、选频

课件描述的腔体结构很朴素：在激光器两端面对面装上两块反射率很高的平面镜，一块全反射，一块大部分反射、少部分透射，让激光从透射那一端射出来。这三块作用值得逐条拆：

第一是方向性。只有沿工作物质轴线传播的光才能在两面镜之间来回反射，偏离轴线的光走几次就逸出腔外，所以出射光束的方向性极好。

第二是延长工作物质。光在腔内往复，每经过一次工作物质就被放大一次，等价于把 20 cm 的晶体「拉」成几十米的有效增益长度。

第三是选频。受激辐射虽然起源于自发辐射，但自发辐射谱宽很大，直接放大出来的是多波长杂光。在反射镜上镀一层选频膜，只让某个波长的光满足高反射条件，其余波长反射一次就漏出腔外，单色性由此而来。

腔内的纵模间隔由腔长决定：Δν = c/(2nL)。取 Nd:YAG 腔长 0.2 m、折射率 1.82，Δν = 3×10⁸ / (2×1.82×0.2) ≈ 412 MHz；而 Nd:YAG 的增益线宽约 120 GHz，两者一除，腔内同时有近 300 个纵模能起振。要单频输出，就得再加标准具或让腔长短到 0.5 mm 量级。

3.2 阈值条件的推导：往返增益必须等于往返损耗

光在腔内走一圈，被增益介质放大、被介质本身吸收散射、又被两面镜子吸收和透射，能量有增有减。只有往返一周的增益大于总损耗，激光才起得来。往返一次的光强变化可以写成：

I_out = I_in · R1 · R2 · exp(2gL) · exp(-2αL)

其中R1、R2是两端反射率，g是增益系数（单位 cm⁻¹ 或 m⁻¹），α是工作物质内部损耗系数，L是腔长。令I_out = I_in解出阈值增益：

g_th = α + (1/2L)·ln(1/(R1R2))

再结合受激辐射截面σ和反转粒子数密度ΔN的关系g = σ·ΔN，就能反推出阈值反转粒子数ΔN_th = g_th / σ。这比「增益大于损耗」这句口头结论有用得多，因为它把器件参数直接和材料参数绑在了一起。

import numpy as np def threshold_gain(L_cm, R1, R2, alpha_cm): """阈值增益系数 g_th，单位 cm^-1""" return alpha_cm + np.log(1.0 / (R1 * R2)) / (2.0 * L_cm) def threshold_inversion(g_th, sigma_cm2): """由增益系数反推阈值反转粒子数密度，单位 cm^-3""" return g_th / sigma_cm2 L = 20.0 # 腔长，cm R1, R2 = 1.00, 0.95 # 全反镜、输出镜反射率 alpha = 0.005 # Nd:YAG 内部损耗系数，cm^-1 sigma = 2.8e-19 # Nd:YAG 受激辐射截面，cm^2 g_th = threshold_gain(L, R1, R2, alpha) print(f"g_th = {g_th:.5f} cm^-1") print(f"dN_th = {threshold_inversion(g_th, sigma):.3e} cm^-3") # 扫描输出镜反射率：耦合输出越强，阈值越高 for R in (0.99, 0.95, 0.90, 0.80, 0.70): print(f"R2={R:.2f} g_th={threshold_gain(L, R1, R, alpha):.5f} cm^-1")

这段代码里三个参数最需要留意。alpha用 cm⁻¹ 时，L也必须是 cm，单位混用是这类计算最常见的错误源。ln(1/(R1R2))这一项代表镜面耦合损耗，Nd:YAG 在 R2 从 0.95 降到 0.70 时，阈值增益从 0.00628 cm⁻¹ 涨到 0.01137 cm⁻¹，几乎翻倍，说明输出耦合不是越大越好——耦合太强会拉高阈值甚至压灭激光，太弱则激光出不来，存在一个最优反射率。

提示：实际工程里选输出镜反射率，一般先按上面的公式算一遍阈值增益，再让泵浦功率留出 2～3 倍余量，避免温度漂移或镜片污染把阈值顶穿。

3.3 稳定腔与临界腔的边界

腔型不只是「两块镜子」那么简单，两块镜的曲率半径R1、R2与腔长L共同决定腔的稳定性。引入g1 = 1 - L/R1、g2 = 1 - L/R2，满足0 ≤ g1·g2 ≤ 1的腔叫稳定腔，光能在腔内长期往返不发散；g1·g2 = 0或1是临界腔，对准误差极其敏感。平行平面腔（R1 = R2 = ∞，即g1·g2 = 1）属于临界腔，虽然模体积大、适合大功率，但装调时两面镜必须严格平行，稍有角度偏差损耗就飙升。共焦腔（R1 = R2 = L，g1 = g2 = 0）稳定性最好，对准容差大，中小功率固体激光器用得多。

4. 激光器分类与选型：波长、脉宽、光束质量三张账

课件给的分类法有两条线：按输出方式分脉冲与连续，按工作物质分气体、固体、液体、半导体、自由电子。这两条线在选型时是交叉起作用的——同一个波长，脉冲和连续对应的设备和工艺完全不同。

4.1 按工作物质分类的谱系

半导体激光器的地位需要单独说一句：它是唯一能直接在电注入下高效出光的类型，插头效率能到 50% 以上，而灯泵固体激光器通常不到 5%。所以现在绝大多数固体激光器已经不再用闪光灯，而是用半导体激光器阵列做泵浦源——课件里提到的「八十年代光电子学与半导体技术推动光纤激光器出现」，说的就是这条链路。

4.2 连续与脉冲：调Q 和锁模在时间轴上的位置

连续输出（CW）适合需要稳定热输入的工艺，比如薄板焊接、表面热处理。脉冲输出要的是峰值功率，同样的平均功率，把能量压缩到更短时间内再释放，峰值能高好几个数量级，对材料的烧蚀机制也从熔融变成气化，热影响区大幅缩小。

在单脉冲能量的量级上，可以把脉冲方式分成三层：

• 普通脉冲（自由振荡）：脉宽百微秒到毫秒，峰值千瓦级。
• 调 Q：腔内插入可饱和吸收体或电光开关，先让腔损耗很高把粒子攒在上能级，再瞬间把损耗降到最低，能量在极短时间内倾泻，脉宽压缩到纳秒，峰值功率到兆瓦级。
• 锁模：让腔内的各个纵模相位锁死，等效于把上百个频率分量合成一个短脉冲在腔内跑，脉宽做到皮秒甚至飞秒。

4.3 选型先算聚焦光斑，再算功率密度

波长决定材料吸收率，脉宽决定热影响区，而这两者最终都要通过聚焦光斑落到工件上。衍射极限下的焦斑直径可以按d = M²·4λf/(πD)估算，M²是光束质量因子（理想高斯光束为 1），f是聚焦镜焦距，D是入射到镜面的光斑直径。

import numpy as np def focus_diameter(wavelength_m, M2, f_m, D_m): """衍射极限聚焦光斑直径，单位 m""" return M2 * 4 * wavelength_m * f_m / (np.pi * D_m) def peak_power_density(energy_J, pulse_width_s, diameter_m): """高斯光束近似的峰值功率密度，单位 W/m^2""" area = np.pi * (diameter_m / 2) ** 2 return (energy_J / pulse_width_s) / area lam = 1064e-9 # 波长 d = focus_diameter(lam, M2=1.2, f_m=0.1, D_m=0.01) print(f"聚焦光斑直径 d = {d*1e6:.2f} um") p_peak = peak_power_density(1e-3, 20e-9, d) # 1 mJ / 20 ns print(f"峰值功率密度 = {p_peak/1e4:.2e} W/cm^2")

按这组参数算出来，焦斑约 1.63 μm，峰值功率密度落在 10¹² W/cm² 量级。这个数字要引起警觉：空气中的光学击穿阈值对纳秒脉冲通常在 10¹⁰～10¹¹ W/cm² 附近，功率密度再往上走，焦点附近会直接打火产生等离子体，不但材料没切好，镜片还可能被反噬。所以真正的工艺窗口不是「功率越高越好」，而是在击穿阈值之下找一个既能气化材料、又不产生等离子体屏蔽的区间。

注意：M²是个很容易被忽略的参数。同一台激光器，M²从 1.2 退化到 3，焦斑直径和功率密度都会跟着变，工艺参数表如果不带M²，换一台机器就可能完全不适用。

5. 速率方程数值实验：把静态参数跑成动态曲线

前面的阈值计算都是稳态的，只回答了「能不能起振」。实际调试更关心「多久起振、过冲多少、弛豫振荡几个周期」，这就需要把上能级粒子数和腔内光子数一起对时间积分。简化后的四能级速率方程组常见写法是：

dN2/dt = Rp - N2/τ_f - σ·c·φ·N2dφ/dt = σ·c·φ·N2 - φ/τ_c + β·N2/τ_f

Rp是泵浦速率，τ_f是上能级荧光寿命，τ_c是光子寿命，β是自发辐射耦合进腔模的比例——最后这项是起振的关键，没有它，φ从零出发永远涨不起来。

5.1 用 solve_ivp 求解并核对阈值

import numpy as np from scipy.integrate import solve_ivp c = 3.0e8 L, R1, R2 = 0.2, 1.0, 0.95 alpha = 0.5 # 内部损耗系数，m^-1（注意是 SI 单位） sigma = 2.8e-23 # 受激辐射截面，m^2 tau_f = 230e-6 # 上能级寿命，s beta = 1e-5 # 自发辐射耦合因子 Rp = 5.0e26 # 泵浦速率，m^-3 s^-1 delta = -np.log(R1 * R2) + 2 * alpha * L # 往返总损耗 tau_c = 2 * L / (c * delta) # 光子寿命 def laser_ode(t, y): N2, phi = y dN2 = Rp - N2 / tau_f - sigma * c * phi * N2 dphi = sigma * c * phi * N2 - phi / tau_c + beta * N2 / tau_f return [dN2, dphi] sol = solve_ivp(laser_ode, [0, 4e-6], [0.0, 1.0], method="LSODA", rtol=1e-8, atol=1e-6, dense_output=True) g_th = alpha + np.log(1 / (R1 * R2)) / (2 * L) print(f"tau_c={tau_c*1e9:.2f} ns 阈值反转 N2_th={g_th/sigma:.3e} m^-3") print(f"稳态 N2={sol.y[0, -1]:.3e} m^-3 稳态光子密度={sol.y[1, -1]:.3e} m^-3")

method="LSODA"是必须的：τ_f是微秒量级、τ_c是纳秒量级，两者差三四个数量级，方程属于刚性系统，用默认的 RK45 会为了稳定把步长压到极小，跑到天荒地老。跑完可以做一个自检——外加泵浦停止后，稳态反转粒子数应当被「钳」在阈值附近，这是因为起振后受激辐射会消耗掉多余的反转粒子，增益自动锁定在等于损耗的位置，这也是激光器能稳定输出的物理原因。

5.2 三个高频踩坑点

beta设成 0 是最常见的错误，结果就是曲线全程一条直线，光子数为零。物理上自发辐射是不可避免的，数值上给个 10⁻⁵～10⁻⁶ 的耦合因子就能在纳秒尺度内起振。

单位混用排在第二。alpha用 cm⁻¹ 而L用 m，tau_c会算错 100 倍，起振时间的量级也就跟着错。建议全部换到 SI 单位再算，最后输出时统一换算。

第三是时间窗口。想看弛豫振荡的完整阻尼过程，积分区间至少覆盖 10 倍上能级寿命；只想看起振前沿，几百纳秒就够。窗口设小了会误判为「振荡不衰减」，设大了又被刚性拖慢。验证手段也简单：把泵浦速率从阈值的一半扫到三倍，观察稳态光子密度是否在阈值处出现明显的拐点，拐点位置应当与g_th/σ算出的N2_th吻合。

本文还有配套的精品资源，点击获取