离子束衰减模型与数学物理推导

在离子注入机（如 Axcelis GSD）中，该公式用于描述带正电的离子穿过气体区域时，因捕获电子发生电荷交换中和而导致的束流衰减。

1. 物理理论基础

• 碰撞截面 ($\sigma$)： 将单个气体分子抽象为一个有效碰撞靶面积 $\sigma$（单位：$\text{m}^2$）。当离子穿过该区域时，会与气体原子发生库仑引力相互作用并捕获电子。
• 气体靶密度 ($n$)： 单位体积内的气体分子数量为 $n$（单位：$\text{m}^{-3}$）。
• 微元中和概率： 当离子束穿过厚度为 $dx$ 的微元气体层时，单个离子发生电荷交换中和的物理概率为： $$P(\text{neutralization}) = \sigma \cdot n \cdot dx$$

2. 数学推导过程（微分方程法）

1. 建立微分变化量方程：
   设在路径位置 $x$ 处，尚未被中和的离子数量为 $N(x)$。当离子束穿过无限小距离 $dx$ 时，损失的离子数量 $dN$ 正比于当前离子数及中和概率： $$dN = – N(x) \cdot \sigma \cdot n \cdot dx$$ （负号表示未中和离子数量随传输距离增加而减少）
2. 分离变量：
   $$\frac{dN}{N(x)} = – \sigma \cdot n \cdot dx$$
3. 定积分与边界条件：
   初始条件 $x = 0$ 时未中和离子数为 $N(0) = N_0$；末态条件为 $N(x)$。假设 $\sigma$ 与 $n$ 沿路径均匀恒定，对两侧进行定积分： $$\int_{N_0}^{N(x)} \frac{1}{N’} \, dN’ = \int_{0}^{x} – \sigma \cdot n \, dx’$$ $$\ln(N(x)) – \ln(N_0) = – \sigma \cdot n \cdot x$$ $$\ln\left(\frac{N(x)}{N_0}\right) = – \sigma \cdot n \cdot x$$
4. 指数求解：
   两边取自然指数 $\exp$，导出衰减公式： $$N(x) = N_0 \cdot \exp(-\sigma \cdot n \cdot x)$$

3. 等价的物理与统计模型

• 朗伯-比尔定律 (Beer-Lambert Law)： 光强/射线穿过介质时的衰减公式 $I = I_0 e^{-\alpha x}$ 与此完全同构（$\alpha = \sigma n$ 为衰减系数）。
• 平均自由程 ($\lambda$)： 定义平均中和自由程为 $\lambda = \frac{1}{\sigma \cdot n}$，代入后公式可写为： $$N(x) = N_0 \cdot \exp\left(-\frac{x}{\lambda}\right)$$
• 泊松过程 (Poisson Process)： 单次碰撞为独立随机事件，未发生中和碰撞的概率服从泊松分布的零事件概率 $P(K=0) = e^{-\mu}$，其中期望碰撞次数 $\mu = \sigma n x$。

4. 在注入机压力补偿 ($pcomp$) 中的工程应用

利用理想气体状态方程 $n = \frac{p}{k_B T}$ 将微观密度转化为宏观压力 $p$：

$$N(x) = N_0 \cdot \exp\left(-\frac{\sigma \cdot x}{k_B T} \cdot p\right)$$

Axcelis GSD 系统在此模型基础上定义了压力补偿系数 $pcomp$，实时计算真实束流 $I_0$：

$$I_0 = I_m \cdot \exp\left( \frac{p}{p_0} \cdot \ln\left(1 + \frac{pcomp}{100}\right) \right)$$