Title here
Summary here
LIBS与等离子物理
August 12, 2026 in LIBS 与光谱分析6 minutes
介绍LIBS中的等离子体基础知识
本文定位：介绍激光烧蚀、等离子体演化、光谱辐射与诊断，并讨论 CF-LIBS 的物理前提。
关联阅读： 数据分析与算法路线 · CF-LIBS 入门 · CF-LIBS 的工程实现。
激光烧蚀与等离子体演化
一、概述
激光诱导击穿光谱（Laser-Induced Breakdown Spectroscopy, LIBS） 是一种基于原子发射光谱的元素分析技术。其核心物理过程是：高能脉冲激光聚焦于样品表面，使微量材料瞬间经历加热→熔化→汽化→电离→等离子体形成→辐射冷却的完整物理过程，最终通过分析等离子体发射的特征光谱实现元素的定性与定量分析。
LIBS被誉为"未来化学分析之星"，具有无需样品前处理、多元素同时检测、适用于固/液/气全形态样品、可远程非接触测量等独特优势。
二、激光烧蚀的物理过程
2.1 激光-物质相互作用的基本阶段
激光烧蚀是一个多时间尺度、多物理场耦合的复杂过程，可分为以下几个阶段：
| 时间尺度 | 物理过程 | 典型温度 |
|---|---|---|
| < 1 ps | 多光子电离、电子激发 | — |
| 1 ps ~ 1 ns | 电子-晶格能量传递、加热 | 10³ ~ 10⁴ K |
| 1 ~ 10 ns | 熔化、汽化、初始等离子体形成 | 10⁴ K |
| 10 ns ~ 1 μs | 等离子体膨胀、逆轫致辐射吸收 | 10⁴ ~ 10⁵ K |
| 1 μs ~ 1 ms | 等离子体冷却、复合、特征辐射 | 10³ ~ 10⁴ K |
2.2 激光能量吸收机制
激光能量被材料吸收主要通过以下机制：
- 带间/带内吸收（金属中自由电子吸收）
- 多光子电离（MPI）：在强激光场下，原子同时吸收多个光子跨越电离能
- 逆轫致辐射吸收（Inverse Bremsstrahlung, IB）：自由电子在离子库仑场中吸收光子能量，这是等离子体形成后最主要的激光能量吸收机制
逆轫致辐射吸收系数为：
$$\alpha_{IB} \propto \frac{n_e \cdot n_i \cdot Z^2}{\sqrt{T_e} \cdot \lambda^2}$$其中 $n_e$ 为电子密度，$n_i$ 为离子密度，$Z$ 为离子电荷态，$T_e$ 为电子温度，$\lambda$ 为激光波长。这表明短波长激光（如266 nm）的等离子体耦合效率更高。
2.3 烧蚀阈值
不同材料存在不同的激光烧蚀阈值（Fluence Threshold），即开始产生可观测烧蚀坑的最低激光能量密度。典型值：
- 金属：约 0.1 ~ 1 J/cm²
- 半导体：约 0.5 ~ 2 J/cm²
- 绝缘体/陶瓷：约 1 ~ 10 J/cm²
当激光辐照度超过 ~1 GW/cm² 时，通常可产生足够密度的等离子体用于LIBS分析。
三、等离子体的形成机制（击穿过程）
3.1 初始种子电子
等离子体形成需要初始自由电子（“种子电子”），来源包括：
- 多光子电离（MPI）：在激光焦点处，电场强度极高（>10⁸ V/m），原子通过同时吸收多个光子实现电离。这是纳秒激光烧蚀中初始电子的主要来源。
- 热发射：材料表面受热后发射电子。
- 缺陷/杂质处的场致发射。
3.2 雪崩电离（Avalanche Ionization）
初始种子电子在激光电场中加速获得动能，当其动能超过原子电离能时，通过碰撞电离产生新的自由电子：
$$e^- + A \rightarrow 2e^- + A^+$$新产生的电子继续被加速并碰撞电离其他原子，电子数目呈指数增长——这就是"雪崩电离"或"级联电离"。
电子密度增长可近似为：
$$n_e(t) = n_{e,0} \cdot \exp(\alpha \cdot t)$$其中 $\alpha$ 为电离速率系数，与激光强度、波长、气体压强等因素相关。
3.3 光学击穿（Optical Breakdown）
当电子密度达到临界值时（约 $10^{16} \sim 10^{18}$ cm⁻³），等离子体变得光学不透明，强烈吸收激光能量，产生明亮的闪光和冲击波——这就是光学击穿。此时等离子体温度可瞬间达到 10⁴ ~ 10⁵ K，电子密度达 $10^{17} \sim 10^{19}$ cm⁻³。
3.4 等离子体屏蔽效应
一旦等离子体形成，它会吸收和散射后续入射激光，减少到达样品表面的激光能量。这一"等离子体屏蔽"效应是LIBS中影响烧蚀效率和定量分析的重要因素。
四、等离子体的时空演化
4.1 早期阶段（0 ~ 100 ns）：等离子体形成与膨胀
- 等离子体在激光脉冲持续期间快速形成。
- 电子温度极高（$T_e \sim 10^4 \sim 10^5$ K），但离子/原子温度较低（$T_i < T_e$），处于非平衡态。
- 等离子体以超声速向外膨胀，形成激光支持吸收波（LSA波）或激光支持爆轰波（LSD波）。
- 此阶段辐射以**连续谱（轫致辐射 + 复合辐射）**为主，光谱信噪比极差。
4.2 中期阶段（0.1 ~ 10 μs）：等离子体冷却与特征辐射
- 等离子体绝热膨胀冷却，温度降至 $10^3 \sim 10^4$ K。
- 电子-离子复合过程产生离子谱线（如 Fe II、Ca II）。
- 随着温度进一步降低，中性原子谱线（Fe I、Ca I 等）逐渐增强。
- 这是LIBS分析的最佳采集窗口（通常延迟 0.5 ~ 5 μs）。
4.3 晚期阶段（> 10 μs）：复合与消散
- 等离子体温度降至数千K以下。
- 原子/离子大量复合，辐射强度急剧衰减。
- 等离子体最终消散为中性气体和微粒。
4.4 等离子体形态
- 在大气环境中，等离子体呈扁平半球形，受环境气体约束，膨胀受抑制。
- 在真空中，等离子体呈细长椭球形，自由膨胀，寿命更长。
- 在液体中，等离子体被液体包围形成空化气泡，经历多次膨胀-坍缩振荡。
等离子体性质、辐射与诊断
五、等离子体的热力学与动力学特征
5.1 关键等离子体参数
| 参数 | 典型值（LIBS等离子体） | 说明 |
|---|---|---|
| 电子温度 $T_e$ | 0.5 ~ 3 eV（~6000 ~ 35000 K） | 早期高，晚期低 |
| 电子密度 $n_e$ | 10¹⁶ ~ 10¹⁹ cm⁻³ | 随时间快速衰减 |
| 离子温度 $T_i$ | 略低于 $T_e$ | 逐渐趋于平衡 |
| 等离子体尺寸 | 0.1 ~ 2 mm | 取决于激光能量和聚焦条件 |
| 寿命 | ~1 μs（大气）/ ~10 μs（真空） | 环境气体显著影响 |
5.2 局部热力学平衡（LTE）
LIBS定量分析的核心前提是等离子体处于局部热力学平衡（LTE）。LTE要求：
- 粒子速度分布服从 Maxwell-Boltzmann 分布
- 激发态布居服从 Boltzmann 分布
- 电离平衡服从 Saha-Boltzmann 方程
McWhirter 判据——LTE成立的最低电子密度条件：
$$n_e \geq 1.6 \times 10^{12} \cdot T_e^{1/2} \cdot (\Delta E)^3 \quad \text{(cm}^{-3}\text{)}$$其中 $\Delta E$ 为跃迁能级差（eV），$T_e$ 为电子温度（K）。
在LIBS等离子体的早期（< 1 μs），电子密度通常满足LTE条件；但在等离子体边缘和晚期，可能偏离LTE。
5.3 Boltzmann 分布与 Saha 方程
Boltzmann 分布（激发态布居）：
$$\frac{N_k}{N} = \frac{g_k}{U(T)} \exp\left(-\frac{E_k}{k_B T}\right)$$Saha 方程（电离平衡）：
$$\frac{n_{z+1} \cdot n_e}{n_z} = \frac{2U_{z+1}(T)}{U_z(T)} \left(\frac{2\pi m_e k_B T}{h^2}\right)^{3/2} \exp\left(-\frac{E_{ion}}{k_B T}\right)$$这些关系是LIBS定量分析（如Boltzmann图法、Saha-Boltzmann图法）的理论基础。
六、等离子体辐射机制
6.1 连续辐射（早期主导）
- 轫致辐射（Bremsstrahlung）：自由电子在离子库仑场中减速时辐射连续谱光子。
- 复合辐射（Recombination Radiation）：自由电子被离子俘获至束缚态时辐射光子。
连续辐射强度：
$$I_{cont} \propto n_e \cdot n_i \cdot Z^2 \cdot T_e^{-1/2} \cdot \exp\left(-\frac{h\nu}{k_B T_e}\right)$$6.2 线辐射（中后期主导）
- 原子谱线：中性原子（如 Fe I 371.99 nm）的电子从高能级跃迁至低能级。
- 离子谱线：一次电离离子（如 Fe II 259.94 nm）的电子跃迁。
- 谱线波长由能级差决定：$\lambda = hc / \Delta E$
6.3 谱线展宽机制
LIBS谱线并非理想的δ函数，而是存在多种展宽：
| 展宽机制 | 物理原因 | 线型 | 典型贡献 |
|---|---|---|---|
| 自然展宽 | 能级有限寿命（测不准原理） | Lorentz | 极小（~10⁻⁵ nm） |
| 多普勒展宽 | 粒子热运动 | Gaussian | 0.001 ~ 0.01 nm |
| Stark 展宽 | 电子/离子微扰 | Lorentz | 0.01 ~ 0.1 nm（LIBS中主要） |
| 压力展宽 | 粒子碰撞 | Lorentz | 大气压下显著 |
| 自吸收展宽 | 等离子体中心辐射被外层同种原子吸收 | 中心凹陷 | 高浓度时严重 |
Stark展宽是LIBS中最重要的展宽机制，其半高全宽（FWHM）与电子密度直接相关：
$$\Delta \lambda_{Stark} = 2\omega \left(\frac{n_e}{10^{16}}\right) + 3.5A \left(\frac{n_e}{10^{16}}\right)^{1/4}$$其中 $\omega$ 为电子碰撞参数，$A$ 为离子碰撞参数。这为电子密度诊断提供了依据。
七、LIBS光谱分析原理
7.1 定性分析
每种元素具有独特的能级结构，因此发射的谱线波长是元素的"指纹"。通过将观测到的谱线波长与标准数据库（如 NIST Atomic Spectra Database）比对，即可确定样品中含有哪些元素。
7.2 定量分析
谱线强度与元素浓度之间的关系由 Boltzmann 方程给出：
$$I_{ki} = \frac{hc}{4\pi \lambda_{ki}} \cdot A_{ki} \cdot \frac{g_k}{U(T)} \cdot N_s \cdot C \cdot \exp\left(-\frac{E_k}{k_B T}\right)$$其中 $A_{ki}$ 为Einstein自发辐射系数，$g_k$ 为统计权重，$U(T)$ 为配分函数，$N_s$ 为总粒子数密度，$C$ 为元素浓度。
7.3 定标曲线与基体效应
实际定量分析中，通常采用定标曲线法（Calibration Curve）：配制一系列已知浓度的标准样品，测量谱线强度，建立强度-浓度关系。
但LIBS面临严重的基体效应（Matrix Effect）：样品的物理化学性质（硬度、热导率、反射率、成分等）会影响烧蚀量和等离子体特性，导致同一元素在不同基体中产生不同的谱线强度。克服基体效应的方法包括：
- 内标法
- 免定标法（CF-LIBS）
- 化学计量学/机器学习方法
- 双脉冲LIBS
7.4 自吸收效应
当等离子体中某元素浓度较高时，中心区域发射的特征光子会被外层较冷的同种基态原子重新吸收，导致：
- 谱线中心强度下降（甚至出现自蚀，即中心凹陷）
- 谱线展宽
- 定量分析出现非线性偏差
自吸收是LIBS定量分析中最重要的系统误差来源之一。
八、等离子体诊断方法
8.1 电子温度诊断
- Boltzmann 图法（Boltzmann Plot）：选取同一元素的多条谱线，以 $\ln(I_{ki}\lambda_{ki}/g_kA_{ki})$ 对 $E_k$ 作图，斜率即为 $-1/(k_BT_e)$。
- Saha-Boltzmann 图法：同时使用原子线和离子线，可同时获得温度和电子密度。
- 双谱线比值法：利用两条谱线的强度比推算温度。
8.2 电子密度诊断
- Stark 展宽法：测量已知Stark展宽参数的谱线（如 H_α 656.3 nm）的FWHM，反推 $n_e$。
- Saha 方程法：结合温度和电离平衡计算。
8.3 时空分辨诊断
- 使用ICCD（增强型电荷耦合器件） 或iCCD相机，配合可调延迟门控，获得等离子体的时间分辨光谱和空间分辨图像。
- 典型设置：延迟时间 0.1 ~ 10 μs，门宽 0.5 ~ 5 μs。
九、影响等离子体特性的关键因素
9.1 激光参数
| 参数 | 影响 |
|---|---|
| 波长 | 短波长（266 nm）吸收效率更高，等离子体更小更热；长波长（1064 nm）烧蚀量更大 |
| 脉冲能量 | 能量越高，烧蚀量越大，等离子体温度越高，但过高会导致等离子体屏蔽加剧 |
| 脉宽 | 纳秒（ns）脉冲：热效应为主；飞秒（fs）脉冲：非热烧蚀，热影响区极小 |
| 聚焦光斑 | 光斑越小，能量密度越高，但烧蚀量减少 |
| 重复频率 | 高重频下前一个脉冲的等离子体/颗粒可能影响下一个脉冲 |
9.2 环境条件
- 大气 vs 真空：大气中等离子体被约束、冷却更快、寿命更短；真空中自由膨胀、寿命更长。
- 环境气体种类：Ar、He、N₂ 等惰性/活性气体影响等离子体温度和谱线强度。Ar 常用于增强LIBS信号。
- 气压：低压下等离子体膨胀更大，自吸收减轻。
9.3 样品性质
- 材料的热导率、比热、反射率、硬度等直接影响烧蚀效率。
- 基体效应是LIBS定量分析的最大挑战。
LIBS 增强技术与方程速览
十、先进LIBS技术中的等离子体物理
10.1 双脉冲LIBS（DP-LIBS）
使用两个时间间隔可调的激光脉冲：
- 预烧蚀模式：第一个脉冲预烧蚀/激发样品，第二个脉冲进一步加热等离子体 → 信号增强 2~100 倍。
- 再加热模式：第一个脉冲产生等离子体，第二个脉冲在等离子体膨胀后重新加热 → 减少自吸收。
物理机制包括：环境气体稀薄化（降低等离子体约束）、等离子体再加热、烧蚀量增加等。
10.2 飞秒LIBS（fs-LIBS）
- 飞秒脉冲（< 1 ps）远短于电子-晶格热弛豫时间，实现非热烧蚀。
- 几乎没有热影响区和熔融飞溅。
- 等离子体初始温度更高，但连续辐射背景更低。
- 定量分析性能更优（基体效应减轻）。
10.3 磁场约束LIBS
在等离子体区域施加外部磁场（通常 0.1 ~ 1 T），利用 Lorentz 力约束等离子体膨胀：
- 提高等离子体温度和电子密度
- 延长等离子体寿命
- 增强谱线强度，改善信噪比
10.4 纳米颗粒增强LIBS（NELIBS）
在样品表面沉积金属纳米颗粒（如 Au、Ag 纳米颗粒），利用局域表面等离子体共振（LSPR） 增强局部电磁场，降低烧蚀阈值，增强LIBS信号。
10.5 激光诱导荧光（LIF）辅助LIBS
在LIBS等离子体冷却后，用可调谐激光选择性激发特定元素的原子，通过检测荧光信号大幅提高灵敏度和选择性。
十一、LIBS等离子体物理的核心方程总结
| 物理量 | 方程/关系 |
|---|---|
| 激发态布居 | Boltzmann 分布 |
| 电离平衡 | Saha 方程 |
| 谱线强度 | $I \propto A_{ki} g_k N \exp(-E_k/k_BT) / U(T)$ |
| Stark 展宽 → $n_e$ | $\Delta\lambda_{FWHM} = f(n_e, T_e)$ |
| Boltzmann 图 → $T_e$ | $\ln(I\lambda/gA)$ vs $E_k$ 的斜率 |
| LTE 判据 | McWhirter 条件 |
| 逆轫致辐射吸收 | $\alpha_{IB} \propto n_e n_i Z^2 / (\sqrt{T_e}\lambda^2)$ |
十二、总结
激光烧蚀与LIBS等离子体物理是一个涵盖激光-物质相互作用、非平衡等离子体动力学、原子分子物理、辐射传输、流体力学等多学科交叉的复杂领域。其核心物理图像可概括为：
高能激光 → 多光子/雪崩电离 → 光学击穿 → 高温高密度等离子体 → 绝热膨胀冷却 → 特征辐射 → 光谱分析
深入理解等离子体的形成、演化和辐射机制，是优化LIBS分析性能（提高灵敏度、降低检测限、克服基体效应和自吸收）的物理基础。当前研究前沿包括飞秒LIBS、双/多脉冲LIBS、磁场/电场约束、纳米增强、机器学习辅助定量分析等方向，不断推动LIBS技术向更高精度、更广应用发展。
CF-LIBS 原理、假设与基本方程
一、CF-LIBS 概述
1.1 什么是 CF-LIBS？
CF-LIBS（Calibration-Free LIBS），又称"免定标LIBS"或"无标样LIBS"，是一种不依赖标准样品定标曲线即可实现元素定量分析的LIBS方法。它直接从等离子体发射光谱出发，基于等离子体物理模型和原子光谱数据库，通过自洽计算反演出样品中各元素的绝对浓度。
1.2 为什么需要 CF-LIBS？
传统LIBS定量分析面临两大瓶颈：
| 问题 | 说明 |
|---|---|
| 基体效应 | 不同样品的物理化学性质差异导致烧蚀量和等离子体特性不同，同一元素在不同基体中产生不同的定标曲线 |
| 标准样品匮乏 | 许多实际样品（如考古文物、外星岩石、工业废料）没有匹配的标准参考物质 |
CF-LIBS 的核心思想是：既然等离子体本身携带了样品成分的全部信息，那么只要正确解读这些信息，就不需要外部定标。
1.3 CF-LIBS 的发展里程碑
- 1999年：Ciucci, Corsi, Legnaioli, Tognoni, Palleschi 等人首次系统提出CF-LIBS理论框架（Applied Spectroscopy, 1999）。
- 2000s：逐步完善自吸收校正、非LTE修正等算法。
- 2010s至今：结合机器学习、双脉冲、时空分辨等技术不断优化。
二、CF-LIBS 的物理基础与核心假设
CF-LIBS 的数学推导建立在以下四个基本假设之上：
假设 1：局部热力学平衡（LTE）
等离子体在观测时间窗口内处于LTE状态，即：
- 激发态布居服从 Boltzmann 分布
- 电离平衡服从 Saha-Boltzmann 方程
- 粒子速度服从 Maxwell-Boltzmann 分布
⚠️ 这是最关键的假设，也是CF-LIBS误差的主要来源。 实际LIBS等离子体仅在有限时空范围内近似满足LTE。
假设 2：等离子体光学薄（Optically Thin）
等离子体对自身辐射透明，不存在自吸收效应。即每条谱线的强度正比于该跃迁上能级的粒子数密度。
⚠️ 实际等离子体中，强共振线几乎总是存在不同程度的自吸收。自吸收校正是CF-LIBS实用化的关键步骤。
假设 3：等离子体空间均匀性
在观测体积内，等离子体的温度 $T$、电子密度 $n_e$ 和各元素浓度是均匀的。
⚠️ 实际等离子体存在显著的径向和轴向梯度。通常通过选择合适的时间延迟和空间积分来近似满足。
假设 4：化学计量烧蚀（Stoichiometric Ablation）
烧蚀产生的等离子体中各元素的相对比例与原始样品中的比例一致，即不存在分馏效应（Fractionation）。
⚠️ 对于多组分样品，不同元素的挥发性差异可能导致烧蚀产物偏离原始配比。飞秒激光和适当的光束整形有助于减轻此效应。
三、CF-LIBS 的核心物理方程
3.1 谱线强度表达式
在LTE和光学薄条件下，一条从能级 $k$ 到能级 $i$ 的原子/离子谱线强度为：
$$I_{ki} = \frac{hc}{4\pi \lambda_{ki}} A_{ki} \frac{g_k}{U_s(T)} N_s C_s \exp\left(-\frac{E_k}{k_B T}\right)$$其中：
- $h$：普朗克常数
- $c$：光速
- $\lambda_{ki}$：谱线波长
- $A_{ki}$：Einstein自发辐射系数
- $g_k$：上能级统计权重
- $U_s(T)$：物种 $s$ 的配分函数
- $N_s$：物种 $s$ 的总粒子数密度
- $C_s$：元素 $s$ 的质量分数（待求量）
- $E_k$：上能级能量
- $k_B$：Boltzmann常数
- $T$：等离子体温度
3.2 归一化条件
所有元素的质量分数之和为1：
$$\sum_s C_s = 1$$这一约束条件是CF-LIBS能够求解绝对浓度的关键——它将相对测量转化为绝对定量。
3.3 Saha-Boltzmann 方程（连接原子线与离子线）
对于同一元素的不同电离态，LTE下满足：
$$\frac{n_{z+1}}{n_z} = \frac{2 U_{z+1}(T)}{n_e U_z(T)} \left(\frac{2\pi m_e k_B T}{h^2}\right)^{3/2} \exp\left(-\frac{E_{ion,z}}{k_B T}\right)$$这使得我们可以用同一套温度参数同时处理原子线和离子线，大大增加可用谱线数量和分析可靠性。
CF-LIBS 计算流程与关键问题
四、CF-LIBS 的标准算法流程
4.1 经典 Boltzmann 图法（Ciucci et al., 1999）
Step 1：获取等离子体温度 $T$
对每种元素 $s$，选取多条谱线，构建 Boltzmann 图：
$$y = \ln\left(\frac{I_{ki} \lambda_{ki}}{g_k A_{ki}}\right) = -\frac{E_k}{k_B T} + \ln\left(\frac{hc N_s C_s}{4\pi U_s(T)}\right)$$以 $y$ 对 $E_k$ 作线性拟合，斜率 $m = -1/(k_B T)$ 给出温度 $T$。
理想情况下，所有元素的Boltzmann图应给出相同的温度。若不一致，说明LTE假设不成立或存在自吸收。
Step 2：计算各元素的相对浓度
从Boltzmann图的截距 $q_s$ 中提取：
$$q_s = \ln\left(\frac{hc N_s C_s}{4\pi U_s(T)}\right)$$由于 $N_s \propto C_s / M_s$（$M_s$ 为摩尔质量），可得：
$$F_s = C_s \cdot K = \frac{U_s(T)}{M_s} \exp(q_s)$$其中 $K$ 是与实验几何、采集效率等相关的仪器常数（对所有元素相同）。
Step 3：归一化获得绝对浓度
利用 $\sum C_s = 1$：
$$C_s = \frac{F_s / M_s'}{\sum_j (F_j / M_j')}$$注：具体归一化形式取决于采用质量分数还是摩尔分数，以及是否考虑电离态修正。上述为简化表述。
Step 4：电离态修正
若只使用了原子线（或只使用了离子线），需通过Saha方程将部分浓度修正到总元素浓度：
$$C_s^{total} = C_s^{(z)} \left[1 + \frac{n_{z+1}}{n_z} + \frac{n_{z-1}}{n_z}\right]$$这需要已知电子密度 $n_e$（通常由Stark展宽获得）。
4.2 算法流程图
原始光谱
│
▼
谱线识别 & 峰值提取
│
▼
自吸收校正 ──────────────┐
│ │
▼ │
Boltzmann 图 → T │
│ │
▼ │
Stark 展宽 → nₑ │
│ │
▼ │
Saha 方程 → 电离态修正 │
│ │
▼ │
计算 Fₛ (各元素) │
│ │
▼ │
归一化 ΣCₛ = 1 │
│ │
▼ │
输出各元素绝对浓度 ◄───────┘五、关键物理问题的深入讨论
5.1 LTE 验证与偏离处理
McWhirter 判据
LTE成立的必要条件（但非充分条件）：
$$n_e \geq 1.6 \times 10^{12} \cdot T_e^{1/2} \cdot (\Delta E)^3 \quad [\text{cm}^{-3}]$$例如，对于 $T_e = 10000$ K、$\Delta E = 3$ eV 的跃迁，要求 $n_e \geq 4.3 \times 10^{15}$ cm⁻³。
更严格的判据
McWhirter判据仅保证碰撞过程主导辐射过程，还需验证：
- Boltzmann 图线性度：非线性表明非LTE或自吸收
- 多元素温度一致性：不同元素给出的 $T$ 应一致
- Saha-Boltzmann 图一致性：原子线和离子线应落在同一直线上
- 时间演化检验：在不同延迟时间下结果应稳定
非LTE修正
当LTE不完全成立时，可采用：
- 碰撞辐射模型（CR Model）：显式求解各能级的速率方程
- 两温度模型：分别描述电子温度和重粒子温度
- 部分LTE（pLTE）：高能级满足LTE，低能级偏离
5.2 自吸收校正
自吸收是CF-LIBS最大的系统误差源。校正方法包括：
(a) C-sigma 图法（Gornushkin et al.）
定义实验坐标 $(x, y)$：
$$x = \frac{g_k A_{ki} \lambda_{ki}}{U(T)} \exp\left(-\frac{E_k}{k_B T}\right), \quad y = I_{ki} \lambda_{ki}$$在光学薄极限下，$y \propto x$（直线过原点）。自吸收使数据点偏离直线向下弯曲。通过拟合 “C-sigma曲线”（考虑自吸收的理论曲线）可同时获得浓度和自吸收程度。
(b) 双线比值法
选取同一元素的两条谱线，其自吸收程度不同（如一条强共振线、一条弱非共振线），通过比值反推自吸收因子。
(c) 迭代自洽法
- 先假设光学薄，得到初始浓度估计
- 用初始浓度计算理论自吸收因子
- 用校正后的强度重新计算浓度
- 迭代至收敛
(d) 基于辐射传输模型的校正
求解等离子体中的辐射传输方程：
$$I_\nu = S_\nu \left[1 - \exp(-\tau_\nu)\right]$$其中源函数 $S_\nu$ 在LTE下等于Planck函数，光学深度 $\tau_\nu = \kappa_\nu L$。通过拟合实测谱线轮廓反推 $\tau_\nu$ 和真实发射强度。
5.3 等离子体不均匀性的处理
多层模型
将等离子体分为若干同心壳层，每层具有独立的 $T$、$n_e$、$C_s$，求解多层辐射传输。
Abel 逆变换
对于轴对称等离子体，通过侧向投影光谱的Abel逆变换重建径向分布。
时空门控优化
选择合适的ICCD延迟时间和门宽，使观测窗口内的等离子体尽可能均匀。通常在等离子体演化的"平台期"（约1~3 μs延迟）采集。
5.4 化学计量烧蚀的验证与修正
- 飞秒激光：非热烧蚀机制大幅减少分馏
- 旋转/扫描样品：平均化表面不均匀性
- 内标验证：用已知含量的元素验证烧蚀化学计量性
- 烧蚀动力学模型：考虑各元素蒸气压差异进行修正
CF-LIBS 改进、评估与应用
六、CF-LIBS 的改进算法与变体
6.1 Saha-Boltzmann 图法（SB-CF-LIBS）
同时使用原子线和离子线构建统一的Saha-Boltzmann图：
$$\ln\left(\frac{I_{ki}\lambda_{ki}}{g_k A_{ki}} \cdot \frac{1}{\Phi_z(T,n_e)}\right) = -\frac{E_k^*}{k_B T} + \text{const}$$其中 $\Phi_z$ 为电离态修正因子，$E_k^*$ 为等效能量。优点：
- 增加可用谱线数量
- 自动包含电离平衡信息
- 对温度误差更鲁棒
6.2 单谱线CF-LIBS
仅需一条谱线即可完成定量（配合已知的等离子体参数）。适用于在线/实时监测场景，但精度较低。
6.3 迭代CF-LIBS
将温度、电子密度、浓度作为耦合变量进行全局迭代优化，而非分步求解。提高了自洽性和精度。
6.4 机器学习辅助CF-LIBS
- 用神经网络学习"光谱→浓度"的非线性映射，隐式补偿非LTE、自吸收、基体效应
- 用遗传算法优化CF-LIBS中的自由参数
- 用PCA/PLS预处理光谱，提取与浓度相关的特征
6.5 双脉冲CF-LIBS
双脉冲配置可改善LTE条件和等离子体均匀性，使CF-LIBS的基本假设更好地成立。
七、CF-LIBS 的性能评估
7.1 典型精度
| 样品类型 | 主量元素误差 | 微量/痕量元素误差 | 备注 |
|---|---|---|---|
| 合金/金属 | 1 ~ 5% (相对) | 10 ~ 30% | 最佳应用场景 |
| 地质样品 | 3 ~ 10% | 20 ~ 50% | 基体复杂 |
| 生物/有机样品 | 5 ~ 15% | >50% | 轻元素困难 |
| 液体/气溶胶 | 5 ~ 20% | — | 信号波动大 |
7.2 与传统定标LIBS的比较
| 指标 | 定标LIBS | CF-LIBS |
|---|---|---|
| 准确度 | ★★★★★ | ★★★☆☆ |
| 精密度 | ★★★★☆ | ★★★☆☆ |
| 无需标样 | ✗ | ✓ |
| 适用新样品 | 需重新定标 | 直接分析 |
| 计算复杂度 | 低 | 高 |
| 对物理假设的依赖 | 弱 | 强 |
7.3 检测限
CF-LIBS的检测限通常比定标LIBS差1~2个数量级，因为缺乏定标曲线的统计平均效应。典型检测限：
- 主量元素：~0.01 wt%
- 微量元素：~10 ~ 100 ppm
- 痕量元素：>100 ppm（较难）
八、CF-LIBS 的典型应用领域
| 领域 | 应用实例 | CF-LIBS优势 |
|---|---|---|
| 文化遗产 | 青铜器、陶瓷釉料、壁画颜料成分分析 | 无损/微损，无标样可用 |
| 行星探测 | NASA Curiosity/Perseverance火星车ChemCam/SuperCam | 无法携带大量标样 |
| 冶金/回收 | 废钢分拣、合金牌号鉴定 | 快速在线，样品种类多变 |
| 环境 | 土壤重金属、大气颗粒物 | 现场原位分析 |
| 核工业 | 放射性废物表征 | 远程操作，标样制备危险 |
| 法医 | 玻璃碎片、油漆涂层比对 | 微量样品，无需破坏 |
九、当前挑战与前沿方向
9.1 主要挑战
- LTE假设的有效性：LIBS等离子体的LTE程度仍有争议，尤其在早期和低密度区域
- 自吸收校正的不确定性：现有校正方法各有局限，尚无通用解决方案
- 轻元素（H, C, N, O）分析困难：谱线少、自吸收严重、空气中干扰大
- 痕量元素精度不足：信噪比低，统计误差大
- 自动化与鲁棒性：谱线识别、参数选择的自动化程度仍需提高
9.2 前沿研究方向
- 第一性原理等离子体模拟：用流体/粒子模拟精确预测等离子体状态，替代经验假设
- 深度学习端到端定量：直接从原始光谱预测浓度，隐式学习物理模型偏差
- 多模态融合：CF-LIBS + Raman + XRF 联合分析，互补验证
- 微型化/便携式CF-LIBS：面向现场应用的紧凑系统设计
- 标准不确定度评估框架：建立CF-LIBS结果的完整不确定度传播链（GUM框架）
十、总结
CF-LIBS 的物理本质可以概括为一句话：
在LTE、光学薄、均匀性和化学计量烧蚀四个假设下，等离子体发射光谱是样品元素组成的完备编码；通过Boltzmann/Saha方程解码这一编码，并利用归一化约束封闭方程组，即可在无外部标样的条件下获得绝对浓度。
CF-LIBS 不是要取代传统定标LIBS，而是在标样不可得或不适用的场景下提供一条基于第一性原理的定量分析路径。其精度虽不及定标法，但其普适性和灵活性使其在行星探测、文化遗产、应急检测等独特场景中不可替代。随着等离子体物理理解的深化、自吸收校正方法的进步以及人工智能技术的融入，CF-LIBS 正朝着更高精度、更强鲁棒性和更广适用性的方向持续发展。