介绍LIBS的基本概念

本文定位：解释免定标定量分析的基本思想、核心假设、计算顺序与局限。

关联阅读： 等离子体物理与假设依据 · 计算步骤与实现细节。

---

详细介绍一下 CF-LIBS 的原理和步骤

无定标激光诱导击穿光谱技术（Calibration-Free Laser-Induced Breakdown Spectroscopy, 简称 CF-LIBS）是传统 LIBS 技术的一种重要演进。

传统的 LIBS 技术在进行定量分析时，往往需要依赖与待测样品基体高度相似的**标准参考物质（定标样品）**来建立定标曲线。而 CF-LIBS 的核心优势在于：它通过纯粹的等离子体物理模型和光谱学理论，直接从光谱数据中反推元素的绝对或相对含量，完全摆脱了对标准样品的依赖。

这使得 CF-LIBS 在探测未知基体（如火星探测、深海探测、极端环境原位分析）时具有不可替代的价值。

以下是 CF-LIBS 的物理原理与具体操作步骤的详细拆解：

一、CF-LIBS 的核心物理原理

CF-LIBS 的成立依赖于三个基本物理假设：

1. 局部热力学平衡（LTE, Local Thermodynamic Equilibrium）： 假设等离子体内部的电子、离子和中性原子之间碰撞足够频繁，整个等离子体可以用一个统一的**激发温度（$T$）**来描述。
2. 光学薄（Optically Thin）： 等离子体内部产生的光子可以完全穿透等离子体逃逸出来，不存在或仅存在极微弱的自吸收（Self-absorption）效应。
3. 化学计量烧蚀（Stoichiometric Ablation）： 激光打出的等离子体成分，与原本固态样品的元素组成比例完全一致。

基于以上假设，某一元素 $s$ 的原子（或离子）从激发态 $k$ 跃迁到低能态 $i$ 产生的光谱线强度 $I_{\lambda}^{ki}$，服从玻尔兹曼分布（Boltzmann Distribution）：

$$ I_{\lambda}^{ki} = F \cdot C_s \cdot \frac{A_{ki} g_k}{U_s(T)} \exp\left(-\frac{E_k}{k_B T}\right) $$

其中：

• $F$：实验系统常数（与激光能量、收集效率、等离子体体积等有关，通常未知）。
• $C_s$：元素 $s$ 在等离子体中的浓度（即我们所求的目标）。
• $A_{ki}$：跃迁概率（可查数据库）。
• $g_k$：上能级统计权重（可查数据库）。
• $E_k$：上能级能量（可查数据库）。
• $k_B$：玻尔兹曼常数。
• $T$：等离子体温度。
• $U_s(T)$：温度 $T$ 下该元素的配分函数（Partition Function）。

数学转化（玻尔兹曼图法）：

将上式两边取自然对数，可转化为一条直线方程：

$$ \ln\left(\frac{I_{\lambda}^{ki}}{A_{ki} g_k}\right) = -\frac{1}{k_B T} E_k + \ln\left(\frac{F \cdot C_s}{U_s(T)}\right) $$

这是一个标准的 $y = mx + q$ 形式。

• 斜率 $m = -\frac{1}{k_B T}$：我们可以通过拟合同一元素的多条特征谱线，计算出直线斜率，进而得出等离子体温度 $T$。
• 截距 $q = \ln\left(\frac{F \cdot C_s}{U_s(T)}\right)$：利用截距，我们可以解出带有未知系统常数 $F$ 的浓度项 $F \cdot C_s$。

二、CF-LIBS 的具体计算步骤

在实际编程或数据处理中，实施 CF-LIBS 通常需要遵循以下严格的流程：

1. 光谱数据采集与预处理（提取真实线强） 使用光谱仪采集激光诱导等离子体的发射光谱。接下来必须进行去基线（扣除连续背景辐射，Bremsstrahlung）、寻峰以及积分计算（计算谱线轮廓下的面积，代表真实的绝对光强 $I_{\lambda}^{ki}$，而非单纯的峰值）。

2. 查阅光谱学参数库 针对选出的各元素特征谱线，访问 NIST（美国国家标准与技术研究院）原子光谱数据库，获取每条谱线对应的 $A_{ki}$、$g_k$ 和 $E_k$ 等物理参数。

3. 计算等离子体温度 (T) 和电子密度 (Ne) 利用同一种元素（如 Fe 原子）的多条谱线，绘制玻尔兹曼图并进行线性拟合，通过斜率求出等离子体温度 $T$。同时，通常利用氢的 $H_\alpha$ 线（656.2 nm）或特定展宽明显的谱线，通过测量其斯塔克展宽（Stark Broadening）来计算电子密度 $N_e$。

4. 结合 Saha 方程计算全状态浓度 等离子体中同一种元素会以**中性原子（I）和一价离子（II）**两种状态存在。根据截距 $q$ 算出的仅是某种电离态的浓度。我们需要利用 Saha-Boltzmann 方程，结合刚才算出的温度 $T$ 和电子密度 $N_e$：

   $$ \frac{N_e \cdot C_{s, II}}{C_{s, I}} = 2 \frac{U_{s, II}(T)}{U_{s, I}(T)} \left(\frac{m_e k_B T}{2\pi \hbar^2}\right)^{3/2} \exp\left(-\frac{E_{ion}}{k_B T}\right) $$

   由此得出每种元素在等离子体中的总浓度量 $C_s^{total} = C_{s,I} + C_{s,II}$（此时仍包含未知因子 $F$）。

5. 归一化求解绝对浓度（消除实验常数 F） 由于 $F$ 难以直接测量，CF-LIBS 引入了归一化原则：假设等离子体中所有元素的质量百分比之和为 1（或 100%）。

   $$ \sum_{s} C_s^{total} = 1 $$

   将各元素包含 $F$ 的浓度代入该方程，由于所有元素共享同一个系统常数 $F$，$F$ 将被作为一个公共系数提取并解出。解出 $F$ 后，再代回各元素的公式中，即可得到各元素的绝对浓度（质量分数或摩尔分数）。

三、局限与挑战

尽管 CF-LIBS 理论完美，但在实际应用中仍面临几个核心挑战：

1. 轻元素的遗漏： 归一化要求探测到样品中的"所有元素"。如果样品中含有大量无法在特定光谱范围内被探测到的元素（例如 C、H、O 等轻元素或卤族元素），归一化总和将产生严重偏差。
2. 自吸收效应： 共振线（往往是强度最高的线）常常存在自吸收，违背了"光学薄"假设，导致高浓度元素计算值偏低。
3. 数据库精度： 计算极度依赖 NIST 数据库中的 $A_{ki}$ 跃迁概率，但许多元素的跃迁概率参数本身存在 10% - 20% 的测量不确定度。