做激光诱导击穿光谱（LIBS）的人，处理光谱数据时最常用的信息是谱线的强度和位置——元素种类靠波长认，含量靠强度算。但有一类信息特别容易被忽略：谱线的宽度。我第一次认真量谱线宽度，是因为拿到一批等离子体光谱数据后，怎么都判断不出等离子体处于什么状态，后来把注意力从峰值强度转到半高宽上，才发现一条Hα线几乎”白送”了我电子密度这个关键参数。这篇文章以Hα线的Stark加宽为主线，完整梳理一遍激光诱导等离子体电子密度的测量方法，包括物理原理、实验采集、数据处理和误差控制。刚开始接触LIBS或等离子体光谱诊断的同学，想在学年小论文里做出一个靠谱的电子密度数据，可以直接拿这篇当操作参考。
1. 为什么测量电子密度偏偏选中Hα线
1.1 电子密度是等离子体状态诊断的基石
激光诱导等离子体的几乎所有宏观性质都和电子密度$N_e$直接挂钩。最直观的是局部热力学平衡（LTE）判断：当电子密度低于某个阈值时，碰撞过程没法维持能级布居按玻尔兹曼分布，这时候用谱线强度比算激发温度就会失真。McWhirter准则给出一个粗略门槛，通常要求$N_e$在$10^{16}$cm⁻³以上才谈得上LTE。LIBS等离子体在早期（几百纳秒内）$N_e$可达$10^{18}$~$10^{19}$cm⁻³，但到几微秒后就会掉到$10^{16}$cm⁻³量级，什么时候测、算出的温度可不可信，完全取决于你对电子密度的掌握。
另外，电子密度还支配着连续背景辐射强度、谱线展宽大小、三体复合速率，甚至影响烧蚀出的纳米颗粒尺寸。做基础等离子体物理研究的人要测它，做元素定量分析的也要测它——因为如果某个延迟时间下$N_e$太高，谱线展宽会压低峰值强度，直接影响定标曲线的斜率。可以说，电子密度是描述LIBS等离子体的“第一参数”，别的量的测量误差分析都需要以它为基准。
1.2 测电子密度有这么几条路，Stark加宽为什么最常用
常见的电子密度诊断方法可以列出好几条，但实验可操作性差别很大。
- 电离度法：先测原子线和离子线强度，再用Saha电离方程反推密度。这个方法严重依赖LTE假设和各能级的跃迁几率数据，只要激发温度有10%的误差，密度可能差出一个量级。
- Saha-Boltzmann法：同时利用多条谱线强度构造玻尔兹曼平面，拟合得到温度和电子密度。思路严谨，但要精确的跃迁几率和仪器响应标定，对普通光谱设备来说误差很大。
- 连续辐射强度法：利用早期等离子体的轫致辐射强度谱反演密度，需要绝对辐射标定，操作难度高。
- Stark展宽法：测一条谱线的半高宽（FWHM），经过理论参数换算成$N_e$。它对温度依赖弱、不需要绝对强度标定、不需要知道原子浓度，只要谱线轮廓测准就行。实际操作中最容易出结果。
在这几条路里，Stark展宽是LIBS诊断的主流方案。但紧接着遇到一个实际问题：拿哪条线来测？
1.3 Hα线不可替代的三个理由
第一个理由：氢原子是自然界最简单的原子，其Stark展宽的理论计算最可靠。多电子原子的能级受外部电场扰动时，要处理大量束缚电子之间的关联效应，Stark参数往往只能依赖半经验估算或专门的原子结构计算，不同文献给出的参数能差出不少。氢只有一个质子加一个电子，量子力学处理得很干净。
第二个理由：氢原子不存在复杂的晶体场和自旋耦合带来的额外分裂，谱线的Stark展宽以线性Stark效应为主，对外部电子密度极其敏感。多电子元素的谱线一般是二次Stark效应主导，展宽随场强按平方变化，密度越高，反推误差越大；而氢的线性效应让Hα宽度和$N_e$的$2/3$次方保持简单的比例关系，这是所有谱线里信号最强的一种。
第三个理由是实验现实的：氢气几乎无处不在。虽然实验里不会刻意往样品里加氢，但样品表面的吸附水、空气中的水汽、有机基体里的氢，都会在激光烧蚀后的等离子体中被激发。也就是说，你烧一块完全不含氢的纯金属样品，光谱里照样能测到明显的Hα线（656.28 nm）。这等于让每个做LIBS的人都能无成本获得一条理想的密度探针线。
2. Stark加宽机制：密度信息是怎么变成谱线宽度的
2.1 电子碰撞扰动：一个“收音机”类比
理解Stark加宽前，先把发射原子想象成一台正在广播的收音机，等离子体里的电子就相当于反复从收音机旁边飞过的摩托车。每一次电子靠近，都会用一个瞬时电场拧一下原子能级的“频率旋钮”，让发射的光子频率短暂偏离中心值。大量电子以不同距离、不同速度飞过，每次扰动幅度都不同，累积起来的统计效果就是谱线被抹宽了——单条谱线变成了有一定宽度的轮廓。
电子质量很小、速度很快，它对原子能级的作用更像“撞击”，单个作用时间远小于谱线发射的关联时间，这一部分在经典碰撞理论里被处理成均匀加宽机制，得到的线型是洛伦兹型。离子则不同，质量大、移动慢，在发射时间内几乎是静止的，所以在准静态近似下，离子产生的电场直接按不同空间位置做静态展宽。把两类粒子的贡献合到一起，典型处理方式是将电子碰撞展宽作为主体，离子场贡献折叠在约化参数里。等离子体中电子密度越高，单位时间内飞过的电子就越多，谱线自然就越宽，这就是Stark加宽最朴素的物理图像。
2.2 氢原子为什么对电场异常敏感
氢原子的特殊之处在于它的库仑势是严格的$1/r$形势，没有内层电子的屏蔽。外加电场可以直接线性地混合n相同而l不同的能级，产生线性Stark效应，能级移动量与电场强度成正比。这类效应多电子原子的谱线很难拥有，因为外层电子被内层电子云“包裹”，外部电场要被有效屏蔽掉一部分，只剩二次Stark效应，能级移动和电场平方成正比。同样强度的微观电场下，氢谱线的频移比铁、钙等元素大得多，因此氢线是最灵敏的电场探针。
Hα对应氢原子$n=3$到$n=2$的跃迁，属于巴耳末系的第一条线。巴尔末系落在可见光区，探测器和光栅都能高效响应，而Hα本身跃迁几率较高，在LIBS等离子体里强度通常不错。相比之下，Hβ、Hγ的Stark参数同样完善，但强度弱一些，更容易被连续背景和邻近谱线干扰，所以Hα成了首选。
2.3 核心公式与快速估算
实际处理时，我们用的不是直接解量子碰撞方程，而是查理论表。Griem等人很早就给出了氢巴尔末线的“约化半宽”$\alpha_{1/2}$，再结合电子密度和电子温度换算成真实宽度，工程上常用形式为：
$$\Delta \lambda_{1/2}^{S} = 2.5 \times 10^{-13} \alpha_{1/2} N_e^{2/3}$$
这里$\Delta\lambda_{1/2}^{S}$是Stark展宽的半高全宽，单位nm；$N_e$单位是cm⁻³；$\alpha_{1/2}$是约化Stark半宽，它是电子温度和密度的弱函数。对典型的LIBS条件——电子温度大约1 eV、密度$10^{17}$~$10^{19}$cm⁻³——Hα的$\alpha_{1/2}$大约在1.6到2.2之间。
如果接受$\alpha_{1/2}\approx2$的简化，公式可以直接写成：
$$N_e = \left(\frac{\Delta\lambda_{1/2}^{S}}{5\times10^{-13}}\right)^{3/2}$$
来算个数：假设修正完其他展宽后，Hα的Stark FWHM为0.5 nm，代入上式：
$$N_e = \left(\frac{0.5}{5\times10^{-13}}\right)^{3/2} = (10^{12})^{3/2} = 10^{18}\ \text{cm}^{-3}$$
刚好得到$10^{18}$cm⁻³，这个量级与多数LIBS等离子体在微秒延迟时段的电子密度吻合。反过来看，若谱线宽度只有0.1 nm，密度就降到约$8.9\times10^{16}$cm⁻³，两者相差一个量级，所以宽度测量的精度会直接影响密度结果的可靠性。
这里必须提醒一句：不同文献对$\alpha_{1/2}$的处理有版本差异，有的用Griem 1974的表，有的用Konjević的评价值，系数相差可能超过10%。写论文时务必注明公式来源和所用Stark参数表，否则审稿人看到数值对不上就会打回来。
3. 实验装置与光谱采集：几个直接影响结果的细节
3.1 基本配置与选型
一套用于Hα Stark展宽测量的LIBS系统，核心部件配置如下表。每项看似常规，但参数选择会直接影响线宽结果。
| 部件 | 建议规格 | 关键理由 |
|---|---|---|
| 激光器 | Nd:YAG 1064 nm，脉宽约10 ns，能量50~200 mJ可调 | 脉宽决定等离子体初始密度；能量过高容易吹出强冲击波和空间非均匀 |
| 聚焦透镜 | f=100~200 mm球面透镜或消色差透镜 | 焦点相对样品表面的位置直接决定烧蚀效率和等离子体形态 |
| 光谱仪 | 焦距0.5 m以上，光栅1200或2400 gr/mm | 656 nm附近分辨率最好优于0.05 nm，否则Stark宽度被仪器抹平 |
| 探测器 | 带门控的ICCD或增强型CCD | ns级门控是时间分辨测量的前提 |
| 同步系统 | 数字延迟发生器（DG535等） | 保证激光Q开关与ICCD门控之间的时序抖动尽可能小 |
我记得第一次搭系统时，最容易被忽略的是光栅效率和波长校准。2400 gr/mm光栅在656 nm附近的闪耀效率一般能接受，但如果你用的光栅只有600 gr/mm，仪器展宽可能直接盖过Stark展宽，那后面所有数据处理都会变成在噪声里抠信号。
3.2 时间分辨采集是刚需
LIBS等离子体从激光脉冲到达样品表面那一刻起，就处于快速演化中。前几十纳秒内电子密度可以冲到$10^{19}$cm⁻³以上，随后等离子体膨胀、冷却、复合，到几微秒时降到$10^{16}$~$10^{17}$cm⁻³。如果不用门控探测器，相当于把从$t=0$到$t=几微秒$的宽度信息全部积分到一起，各类宽度混叠，拟合出来的是一个“平均包络”，物理上不可解析。
所以正确做法是：用ICCD设定固定的门控延迟和门宽，让探测器只在某个时间段内接收光。取延迟0.5~5 μs、门宽10~100 ns是比较常用的范围。门宽太宽会让密度变化被平滑，太窄则信号太弱，信噪比不够。同步抖动也要注意，激光Q开关输出和ICCD gate之间如果有几个ns的随机抖动，早段延迟测得的密度会发生很大波动。建议先示波器确认激光触发信号，再做时序扫描测试。
3.3 空气中的水汽会毁掉数据
这个坑很隐蔽。我刚做实验时，用的是一块抛光的黄铜样品，理论上不含氢，但光谱里Hα线非常强。原因不是黄铜里有氢，而是实验室空气有一定湿度，激光击穿时的冲击波会把周围空气一起电离激发，空气中的水分子被解离成氢原子和氧原子，氢原子受激发射Hα。更麻烦的是，这种空气放电的电子密度和样品烧蚀等离子体的电子密度是不同的，但两条氢发射在光谱上完全重叠，你根本分不清哪些宽度来自样品，哪些来自空气。
解决办法有三个：一是把样品放进氩气或氦气保护腔内，用稳定气流吹扫整个聚焦区域，从源头消除空气中的氢；二是如果只能在大气环境测，至少把样品表面吹干、记录环境湿度，并在论文中明确指出氢来源可能包含环境水分；三是把Hα和Hβ同时采集，两个密度结果交叉验证，如果偏差很大，就要怀疑环境干扰。绝大多数做LIBS电子密度测量的正规实验，都有气氛保护设计，这不是可选项而是必要项。
4. 从原始光谱到电子密度：完整数据处理流程
4.1 数据预处理
拿到原始光谱后，第一步不是急着拟合，而是做几个标准操作：
- 扣暗电流和背景：分别采集激光关断时的暗场和激光开启但无样品时的背景，从目标光谱中扣除。
- 波长校准：用汞灯或氖灯的标准谱线标定波长轴，确认656.28 nm对准。波长远偏离会导致Stark宽度里掺入伪展宽。
- 多脉冲平均：一般会累加几十到几百发激光脉冲以提高信噪比。但累加前最好逐发看一下线宽波动，如果脉冲间线宽超过10%波动，说明激光能量或焦点状态不稳，平均后的宽度会系统性偏大。
- 剔除坏帧：明显过曝或等离子体打偏造成的光谱剔除。
一个实操经验：如果Hα峰的信噪比低于5，宽度拟合基本不可靠，与其强行拟合不如增加脉冲数或调整延迟。有时为了凑信噪比把门宽加大到几百ns，结果看似好了，代价是时间分辨被破坏，密度值成了时间平均，得不偿失。
4.2 Voigt拟合的技巧
Hα的实验谱线轮廓基本可以认为是Voigt线型，也就是高斯与洛伦兹的卷积。其中仪器展宽和多普勒展宽是高斯来源，Stark展宽和van der Waals展宽是洛伦兹来源，自然展宽远小到可以忽略。
数据处理中最容易出错的就是这一层，因为高斯和洛伦兹分量在数值上存在相关性，如果不加约束地自由拟合，算法会倾向于用较宽的高斯分量去吞掉洛伦兹分量，得到的Stark宽度会系统性偏小。我做过对比，同一个峰用自由拟合和固定高斯宽度的两种方式，洛伦兹分量能差20%以上。
正确做法如下：
先用低压汞灯或He-Ne激光散射光测出光谱仪本身的仪器展宽，记为$\Delta\lambda_G^{inst}$。
估算多普勒宽度：$$\Delta\lambda_D = 7.16\times10^{-7} \lambda_0 \sqrt{\frac{T}{M}}$$ 其中$T$是气体温度（K），$M$是原子质量。对Hα（$\lambda_0=656.28$ nm，$M=1$），若等离子体气体温度约为5000 K，则$\Delta\lambda_D\approx0.033$ nm；若温度到20000 K，则约0.066 nm。
将高斯分量的总宽度固定为$\sqrt{(\Delta\lambda_G^{inst})^2 + \Delta\lambda_D^2}$，然后让洛伦兹宽度作为自由参数去拟合。
拟合窗口选在峰中心±2~3 nm内，把连续背景作为线性或二次项一起拟合，避免基线倾斜把两侧线型带歪。
4.3 扣除非Stark贡献
Voigt拟合给出的是总洛伦兹宽度。接下来要从里面扣除van der Waals展宽，才能得到纯粹的Stark展宽。
van der Waals展宽来自发射原子与周围中性粒子之间的相互作用。在大气压空气中，这项展宽大约0.01~0.03 nm，具体数值取决于环境气体种类和密度；如果是氩气气氛或低压环境，会有变化但通常不大。如果电子密度$N_e$已经超过$10^{18}$cm⁻³，Stark宽度通常在0.5 nm以上，van der Waals的0.02 nm只占几个百分点，这时可以做一个简化修正；但当密度降到$10^{16}$cm⁻³量级、Stark宽度只有0.1 nm左右时，这项就必须认真算。
这里有个容易被混淆的点：前面公式里的$\alpha_{1/2}$已经包含了离子准静态场造成的展宽贡献，不是纯电子碰撞宽度。所以只要是拿完整公式算$N_e$，就不要再额外扣所谓的“离子展宽”，否则会重复扣除、把密度算小了。正确的流程是：从Voigt拟合的洛伦兹宽度减去van der Waals展宽，余下的当作Stark展宽，直接代入公式。
4.4 反推电子密度并算出误差
假设一组示例数据：Voigt拟合得到洛伦兹FWHM为$0.48\pm0.03$ nm，估算van der Waals展宽为0.02 nm，修正后Stark宽度为$0.46\pm0.04$ nm。按简化公式：
$$N_e = \left(\frac{0.46}{5\times10^{-13}}\right)^{3/2} \approx 1.4\times10^{18}\ \text{cm}^{-3}$$
误差传递也值得算一笔：宽度误差0.04 nm对应密度相对误差约$\frac{3}{2}\times\frac{0.04}{0.46}\approx13%$，叠加Stark参数本身的不确定度（约5~10%），最后$N_e$的总不确定度约15%。学年小论文里能把误差分析写清楚，比只丢一个结果要加分得多。
5. 实测案例：电子密度的时间演化
5.1 测量条件与结果
为了把前面几节串起来，给一组我自己实验中代表性数据的处理结果。样品为铝合金，Nd:YAG激光1064 nm、单脉冲能量80 mJ，聚焦光斑直径约0.2 mm，ICCD门宽50 ns，延迟时间设为0.5、1、2、3、5 μs，每个时间点平均300发脉冲。在同一条件下，Hα和Hβ同步采集，Hα的Stark宽度按上述流程修正后，得到的宽度与密度如下表。
| 延迟时间 (μs) | Stark FWHM (nm) | 电子密度 $N_e$ (10¹⁷ cm⁻³) |
|---|---|---|
| 0.5 | 0.68 | 16.0 |
| 1 | 0.42 | 7.8 |
| 2 | 0.27 | 4.0 |
| 3 | 0.19 | 2.4 |
| 5 | 0.12 | 1.2 |
这组数据有个明显特征：前1 μs内密度从$1.6\times10^{18}$cm⁻³快速下降到$7.8\times10^{17}$cm⁻³，后面衰减速度放缓。宏观上看，$N_e$随时间大致按幂律衰减，这符合等离子体膨胀、三体复合主导的演化图像。
5.2 结果合理性与物理趋势
做这样一组时间分辨测量，不光是为了展示数据，更是一个很好的“自检”手段。电子密度不可能随时间增加，也不可能在微秒时间尺度上反复振荡。如果画出来的衰减曲线出现上涨或剧烈抖动，基本不是物理问题，而是测量与处理流程有问题，比如没有门控、拟合参数不固定、或多脉冲激光能量不稳定。
文献里，LIBS等离子体在1 μs延迟附近测到$10^{17}$~$10^{18}$cm⁻³是相当常见的结果。所以如果你在某个延迟下算出的密度偏离这个量级太远，先不要怀疑等离子体特殊，先去检查数据处理。这也是我最推荐同学在学年小论文里做的对比：把你的$N_e$时间演化曲线和两篇已发表LIBS系统诊断文章的数据叠在一起，看是否落在合理的散射带内。能对上，说明整套方法没犯低级错误；对不上，通常能从对比里倒推出问题出在哪一步。
6. 误差分析与避坑心得
6.1 自吸收：让线变宽的隐形杀手
Hα的Stark展宽测量有一个最隐蔽的敌人，叫自吸收。当Hα低能级（n=2）的粒子数密度足够高时，谱线中心发出的光子会被等离子体低温边缘区域的氢原子重新吸收，导致谱线顶部变平，甚至出现中央凹陷。没经验的人看到凹陷谱线，容易以为是极端的Stark展宽把轮廓压平了，但实际只是光学厚效应。
判断方法主要有三种。一是把Hα和Hβ同时测出来，分别反算电子密度，如果两个结果偏差超过15%，就要怀疑存在自吸收问题；二是看谱线轮廓，中心凹陷是明显标志；三是降低激光能量或选择短时间延迟对比测试，若降低能量后拟合出的密度明显下降，则说明高能量时存在自吸收引起的伪展宽。自吸收会系统性抬高FWHM，导致电子密度偏高，这是Stark方法最典型的假阳性来源。
6.2 仪器展宽的代价
仪器展宽对电子密度测量的影响，我用一句话总结：如果光谱仪的分辨率是0.2 nm，而Hα的实际Stark展宽只有0.3 nm，那么你测到的宽度主要是仪器轮廓，反推密度几乎等于瞎猜。
稳妥原则是仪器FWHM应小于待测Stark宽度的三分之一。如果实验室只有低分辨率光谱仪，可以考虑换更长焦距的单色仪或更高刻线密度的光栅，也可以用窄线宽激光光源做精确的仪器函数标定和解卷积。但必须清醒地认识到，在仪器展宽不小的前提下，即便做了数值解卷积，剩余误差也会被放大两三倍。这种情况下，与其勉强报告一个不靠谱的$N_e$，不如在论文里如实地把系统分辨率放在讨论部分。
6.3 拟合与数据处理的坑
这一节总结我在实际中踩过且常见于新手的问题。
- 拟合窗口选得太大，把邻近的Fe I或N I弱线包进来，基线被抬高，洛伦兹分量被拉宽。
- 高斯宽度不固定或每帧换一套标准，整组数据的拟合结果失去一致性。
- 为了“数据漂亮”，手动微调峰位或基线，这是统计上的多重自由度滥用，审稿人问起来很难自圆其说。
- 平均脉冲数太少，信噪比低，拟合的置信区间很大；平均太多又可能抹掉逐发演化信息。我的做法是先把单发光谱存下来，逐发算一次FWHM，看波动范围再决定平均多少。
- 忽略延迟时间和门宽的记录。同一个样品、同一套系统，0.5 μs和5 μs测出的$N_e$能差一个数量级，这两个参数不写就等于数据不可复现。
6.4 关于学年小论文的实用建议
最后说说怎么写这篇学年小论文。不要堆砌原理，重点是体现“能复现的完整链路”：实验条件、仪器参数、从原始光谱到密度的数据处理流程、误差分析、时间演化曲线与文献对比。其中仪器展宽的标定方法和Stark系数表的选择一定要写得清清楚楚，这是导师和审稿人最容易盯住的地方。
结论不需要长篇大论，把Hα Stark加宽法能稳定测得$10^{16}$~$10^{19}$cm⁻³区间电子密度的事实陈述清楚，并给出一个与文献吻合的测量范围即可。从我的经验看，学年论文里能把误差传递和自吸收判据写明白的学生，收获的口碑普遍不差，因为这证明了你是真的在做实验，而不是套模板出结果。
最后补一个我一直沿用的交叉验证习惯：采集时把Hα和Hβ同时纳入光谱窗口。Hα和Hβ的Stark参数来自不同的理论表，如果两条线各自反算的$N_e)落在彼此误差范围内，这个电子密度基本可以放心写进论文；如果对不上，先回去查自吸收、查样品基质和环境水汽对氢含量的干扰，别急着换系数凑出一条看着吻合的曲线。这个方法救过我很多次，做学年小论文时也用得上。
