# LIBS与等离子物理介绍LIBS中的等离子体基础知识

> **本文定位**：介绍激光烧蚀、等离子体演化、光谱辐射与诊断，并讨论 CF-LIBS 的物理前提。
>
> **关联阅读**：[数据分析与算法路线]({{< relref "blog/LIBS/index.md" >}}) · [CF-LIBS 入门]({{< relref "blog/CF-LIBS/index.md" >}}) · [CF-LIBS 的工程实现]({{< relref "blog/CF-LIBS-details/index.md" >}})。

---

### 激光烧蚀与等离子体演化

#### 一、概述

**激光诱导击穿光谱（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) 迭代自洽法

1. 先假设光学薄，得到初始浓度估计
2. 用初始浓度计算理论自吸收因子
3. 用校正后的强度重新计算浓度
4. 迭代至收敛

###### (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 主要挑战

1. **LTE假设的有效性**：LIBS等离子体的LTE程度仍有争议，尤其在早期和低密度区域
2. **自吸收校正的不确定性**：现有校正方法各有局限，尚无通用解决方案
3. **轻元素（H, C, N, O）分析困难**：谱线少、自吸收严重、空气中干扰大
4. **痕量元素精度不足**：信噪比低，统计误差大
5. **自动化与鲁棒性**：谱线识别、参数选择的自动化程度仍需提高

##### 9.2 前沿研究方向

- **第一性原理等离子体模拟**：用流体/粒子模拟精确预测等离子体状态，替代经验假设
- **深度学习端到端定量**：直接从原始光谱预测浓度，隐式学习物理模型偏差
- **多模态融合**：CF-LIBS + Raman + XRF 联合分析，互补验证
- **微型化/便携式CF-LIBS**：面向现场应用的紧凑系统设计
- **标准不确定度评估框架**：建立CF-LIBS结果的完整不确定度传播链（GUM框架）

---

#### 十、总结

CF-LIBS 的物理本质可以概括为一句话：

> **在LTE、光学薄、均匀性和化学计量烧蚀四个假设下，等离子体发射光谱是样品元素组成的完备编码；通过Boltzmann/Saha方程解码这一编码，并利用归一化约束封闭方程组，即可在无外部标样的条件下获得绝对浓度。**

CF-LIBS 不是要取代传统定标LIBS，而是在**标样不可得或不适用**的场景下提供一条基于第一性原理的定量分析路径。其精度虽不及定标法，但其**普适性和灵活性**使其在行星探测、文化遗产、应急检测等独特场景中不可替代。随着等离子体物理理解的深化、自吸收校正方法的进步以及人工智能技术的融入，CF-LIBS 正朝着更高精度、更强鲁棒性和更广适用性的方向持续发展。