# CF-LIBS基础介绍LIBS的基本概念

> **本文定位**：解释免定标定量分析的基本思想、核心假设、计算顺序与局限。
>
> **关联阅读**：[等离子体物理与假设依据]({{< relref "blog/LIBS-Physics/index.md" >}}) · [计算步骤与实现细节]({{< relref "blog/CF-LIBS-details/index.md" >}})。

---

## 详细介绍一下 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% 的测量不确定度。