电子阻止本领：决定离子射程与深度分布的关键参数

1. 一个决定离子射程却总被一笔带过的参数

1.1 我为什么开始认真研究它

前几年给一个客户做 200 keV 质子辐照石墨烯的损伤评估，按他的预期，这个能量在材料里的平均射程应该在微米量级。我拿到模拟结果后却怎么都对不上实验里的层错分布，最后发现问题出在电子阻止本领（electronic stopping power）的低能区处理上。那段时间我把 SRIM 的帮助文档翻了个遍，又借了半天的 CASP 和 CasP 的论文，才意识到一个事实：和核阻止本领相比，电子阻止本领在绝大多数文献里都是“被默认存在”的东西，但一旦你的实验牵涉到精确深度、多层膜剖面或辐照剂量匹配，它反而是最先跳出来捣乱的一个参数。

这也是我今天想认真聊它的原因。电子阻止本领不是什么冷门概念，做离子注入、离子束分析（RBS）、加速器质谱、空间器件单粒子效应模拟的人都绕不开它。它描述的是入射离子在靶材料中穿行时，通过激发、电离靶原子电子而损失能量的速率。射线与物质相互作用里，能量损失机制只有两大分支：一个贡献给靶原子核的弹性碰撞，也就是核阻止本领；另一个贡献给电子系统，就是电子阻止本领。后者在多数能区里占大头，尤其在 MeV 量级离子分析中，几乎决定了你能探测到的深度和能量分辨。

1.2 精确定义和工程单位

用专业一点的说法，电子阻止本领定义为入射离子在单位路径长度上传递给靶电子系统的能量，通常写作 S_e = -(dE/dx)_e，或者除以靶原子体密度 N，得到截面形式的阻止截面 ε_e = S_e / N。这里加负号是因为离子在减速，能量在减少。你会看到不同的文献使用不同的单位，这是新人最容易懵的地方。

工程上常见单位有三个。第一是 SRIM 常用的 eV/(10^15 atoms/cm²)，它把能量损失除以面密度而不是线密度；第二是 MeV/(mg/cm²)，常用于离子束分析里的阻止因子表；第三是 keV/nm、eV/Å 这类直观的线性能量损失单位。它们之间不是差一个密度就完事，还要注意原子密度的量纲。举个例子，硅的原子密度约为 5×10^22 atoms/cm³，厚度 1 Å 对应的面堆积就是 5×10^14 atoms/cm²。所以如果 SRIM 给出某个化合物的阻止本领是 1 eV/(10^15 atoms/cm²)，换算成体密度下的能量损失就是约 0.5 eV/Å。许多人在对照文献数据时直接拿数字比大小，忽略了靶材料密度不同，结果差了接近一个数量级。

提示：在写论文或做模拟参数表时，务必把单位写完整。比如“S_e (Si) = 0.64 eV/Å at 1 MeV H”和“S_e (Si) = 0.64 keV/µm at 1 MeV H”其实是同一个量，但如果你只写数字，别人根本没法复现。

2. 电子阻止的物理图像：从高速库仑激发到电子气体阻力

2.1 高能区：Bethe-Bloch 框架与速度标尺

理解电子阻止本领，关键是抓住一个物理标尺：入射离子的速度与靶原子中电子轨道速度（用玻尔速度 v0 ≈ 2.18×10^6 m/s 作参照）的相对大小。

当入射离子速度远大于靶电子轨道速度时，入射离子可以看作一个高速飞行的带电粒子，靶电子在它产生的电磁场中被“踢”了一下，获得能量。这个过程可以用 Bethe-Bloch 公式描述。非相对论近似下，阻止本领大致正比于 Z_eff² / v²，再乘一个包含靶平均激发能的 ln 项。入射离子的等效电荷 Z_eff 不是恒定的，它需要考虑离子的电荷态和在靶内被电子俘获、剥离的平衡过程。高能区里，速度越高，单位路径上能量损失反而越小，这种“高速时减速反而慢”的反直觉行为，正是 1/ v² 项导致的。

在具体实验里，这个区域对应兆电子伏量级的氢、氦离子。比如 1 MeV 质子入射到硅中，速度大约是玻尔速度的 2 倍左右，电子阻止在各处的贡献已经明显压过核阻止。RBS 分析中常见的 2 MeV He⁺ 也主要落在 Bethe 框架的适用范围。此时电子阻止本领的变化比较平滑，实验数据与理论计算都能对到几个百分点内，所以平时工作中反而不太纠结这一区间的精度。

2.2 低能区：Lindhard-Scharff 与“拖拽”图像

一旦入射离子速度降到与电子轨道速度可比甚至更低，事情就变了。这时离子的库仑势会持续扰动周围的电子气体，形成了一个类似“粘滞阻力”的电子拖拽过程。低能区的图像通常用自由电子气体模型来描述，最有名的就是 Lindhard-Scharff 理论。这个理论有一个重要结论：电子阻止本领与入射离子速度近似成正比，S_e ≈ b v。换算成离子能量，因为 v 与 E^(1/2) 成正比，所以 S_e ≈ k E^(1/2)。这就是为什么许多软件在低能区给出随能量单调上升的电子阻止曲线，而不是高能区那样先下降再上升。

需要注意，这里的“低能”并不一定指绝对能量低。对质子来说，几 keV 到几十 keV 就已经进入低能区；但对金、银这类重离子，几十 keV/u 甚至几百 keV/u 时仍然可能处于复杂的电荷交换区。因为重离子的有效电荷态在靶内反复变化，电子阻止曲线会出现所谓“Z₁ 振荡”——阻止能力随原子序数周期性起伏，而不是单调变化。这背后是电子壳层结构的周期性影响：闭壳层离子更难被电离，等效电荷低，电子阻止也低。

2.3 中间区域：经验数据库存在的理由

最麻烦的是中间区域，也就是速度接近玻尔速度的前后那一段。这里量子力学微扰论不好用，自由电子气体模型也不完全适用，理论计算从第一性原理出发会很吃力。于是大家开始用实验测量来构筑数据库，最常用的就是 Ziegler 和 Anderson 的经验拟合曲线，后来被整合进 SRIM 的电子阻止数据库里。

这一段的电子阻止本领曲线往往有一个峰，峰值位置大致在离子速度 v ≈ Z_1^(2/3) v0 附近。为什么用这个标度？它本质上是入射离子与靶电子壳层速度匹配的判据。低于峰值时曲线上升，高于峰值时曲线下降。许多工程人员只知道“看曲线峰值”，但不知道这个峰值的位置可以从离子原子序数和玻尔速度推算。当你需要预测一个新材料的阻止本领时，先按这个标度估算峰值位置，再去查阅数据库，会更有方向感。

3. 模拟计算时如何正确使用电子阻止本领

3.1 SRIM里电子阻止的参与路径

SRIM 是很多人接触离子-物质相互作用时用的第一个工具，它的完整名称是 Stopping and Range of Ions in Matter。在 SRIM 中，电子阻止本领并不是一个孤立的输入参数，它通过 TRIM 的蒙特卡罗步长计算深度参与射程、损伤分布和背散射粒子能量谱的模拟。

TRIM 每一步的核心逻辑是这样的：先根据当前离子剩余能量查表得到核阻止本领和电子阻止本领，二者相加就是总阻止本领；然后用随机方式决定一次碰撞的能量转移和偏转方向，其中核碰撞用 ZBL 势描述弹性散射，而电子能量损失则在两次核碰撞之间连续消耗。这种“连续慢化+离散碰撞”的处理方式，在离子射程远大于原子间距时是合理的。但要注意，SRIM 的电子阻止数据库对低能区重离子的处理并不完美，很多情况下它是对实验数据的经验拟合，外推到没有实验点的组合时误差可能超过 20%。

我在实际使用中的做法是：先跑一次“SR Module”看电子阻止、核阻止和总阻止随能量的曲线，确认自己关心的能量点在哪个区域，然后再跑 TRIM。很多人上来直接跑 TRIM，最后从输出文件里看到 stopping 曲线，才发现参数设错了能量范围，属于本末倒置。

3.2 实验校正的几个关键点

电子阻止本领的模拟值能否直接用，你需要问自己三个问题：靶是不是单质？是不是非晶？能量区域有没有实验数据？

对单质非晶靶（比如无定形碳、多晶硅），SRIM 的精度通常还不错。但对化合物靶，情况会复杂一些。SRIM 默认按 Bragg 加和规则处理化合物，即把化合物看作各组成元素的线性加权平均。这个规则在能量较高时可行，在低能区却经常失效，因为化合物的化学键会改变电子态分布，进而影响低能电子阻止。比如 SiC、GaN 这类宽带隙半导体，直接按 Bragg 规则算出来的低能电子阻止和实验值可能有明显偏差。文献里有人提出用“核价带修正”来改进，但 SRIM 本身并不自动完成这种修正。

另一个容易被忽略的点是靶的密度。TRIM 输入里有一个靶密度，许多人直接填固体材料的宏观密度，但在纳米薄膜或多孔样品中，真实密度可能低得多。电子阻止本领与单位体积内电子数成正比，密度差 20%，射程就差 20% 量级。更隐蔽的是，同一靶材在不同工艺下的密度也会变化，比如不同方法生长的 SiO₂ 薄膜，密度可以从 2.1 变到 2.3 g/cm³。

3.3 单位换算与容易犯的错

我见过不止一次，有人在实验报告里写“电子阻止本领为 0.64”，既没有单位，也没有说明是针对什么入射离子和靶材。这种数据完全无法使用。下面这张表是常用换算关系，建议收藏：

这些换算看起来简单，但是只要有一个单位看错，结果就全乱。尤其是从 SRIM 输出文件复制数据时，它默认的单位经常是 eV/(10¹⁵ atoms/cm²)，而你在论文里引用时通常需要 keV/µm。我习惯在模拟之前先把所需的单位体系定下来，直接在 SRIM 的 stopping 输出界面选好，避免后续手工转换时引入低级错误。

4. 实际应用中的电子阻止效应：不得不算的几笔账

4.1 离子注入深度分布与射程尾部

离子注入的深度分布主要由射程决定，而射程绝大多数情况下由电子阻止贡献的总能损主导。以 MeV 离子注入为例，几百微米的注入深度几乎完全来自电子阻止对高能段能量的消耗，核阻止只在离子接近停止时才起作用。这带来一个直接后果：如果你电子阻止查表值偏大，常规射程预计会偏浅；偏小，则偏深。对浅层注入（几 keV 到几十 keV），核阻止与电子阻止的贡献相当，二者误差还会相互耦合。

一个记忆中的案例是低能硼注入硅。20 keV 的硼离子，在硅里的电子阻止本领处于低能区，S_e ≈ k E^(1/2)，很多人只关心核阻止，因为损伤峰主要由核碰撞产生。可一旦要精确模拟 B 的深度分布尾部，电子阻止的误差会改变注入离子的最终停止位置，边际浓度对器件漏电流又非常敏感。这种情况下，与其用 SRIM 默认值，不如先找同能量范围内硼注入硅的实验射程数据反推电子阻止，再回填到模拟中。

4.2 辐射损伤区域与电子激发耦合

电子阻止本身虽然不直接产生永久性晶格位移，但它和损伤的演化是耦合的。高能离子在靶中留下一条电子激发、电离的轨迹，局部电子温度骤升，然后通过电子-声子耦合加热晶格。这条“热峰”路径在实验上已被反复验证：同样的入射离子，在绝缘材料中的损伤形貌往往与纯核碰撞预测不同，就是因为电子阻止的能量在纳米尺度内被转换成了晶格热振动，造成局部无序甚至相变。

我的一个深刻体会来自重离子辐照氧化物陶瓷。用几百 keV 的 Au 离子辐照，核阻止峰值深度和最终的空位分布峰值并不完全重合，电子激发在近表面区域带来了额外的缺陷回复。你如果只盯着核阻止曲线做分析，很难解释近表面空位浓度为什么比模拟低。这种情况下把电子阻止的“非热”贡献单独列出来建模，比简单加和有用得多。

4.3 离子束分析中能量损失的另一面：缺陷深度剖面

在 RBS 和 ERDA 分析中，入射离子和出射离子在路径上的电子阻止决定了能量-深度刻度。我们通常用表面能近似计算深度标尺，假设能量损失率在分析窗口内近似恒定，这个近似在小深度范围内是合理的。可一旦测量几微米深的剖面，电子阻止随能量的变化就不能忽略，需要用数值积分逐层迭代。

反向出射的离子也会损失能量，而且由于出射路径能量更低，电子阻止可能更大。所以深度定标的精度，其实取决于你对电子阻止曲线在整个能量区间的掌握程度。我在处理 RBS 深度-浓度剖面时，都会同时跑两组计算：一组用 SRIM 的电子阻止表格，另一组用实验校正后的表格，看两者造成的深度差异是否在可接受范围内。如果差异太大，我会选择直接在目标样品上做多角度 RBS 自洽标定，而不是继续依赖默认参数。

5. 关于电子阻止本领的几条实践经验与边界

5.1 低速极限下不要默认忽略

很多人看到“电子阻止本领”这个词，第一反应是“高能才需要管它”。但实际上在非常低的能量下，电子阻止仍然存在，并且可能与核阻止同量级。比如几个 keV/u 的轻离子在固体中，电子气体阻力不再是线性正比，而是进入所谓“阈值效应”区域：当离子速度低于某个能量阈值时，电子阻止快速趋近于零，因为靶中电子的平均激发能已经高于入射离子可提供的最大能量转移。但这个阈值不是突然消失，中间有一段非常陡的斜坡，如果你使用的软件没有正确描述这一段，射程计算会显著偏差。

在模拟极低能离子（如本底注入、表面改性）时，我会专门检查电子阻止在 1 keV/u 以下是不是呈合理的下降趋势。有些经验拟合公式在这一区域会继续外推成线性上升，这是不物理的，必须人工截断或者换用更多实验点验证的数据库。

5.2 化合物靶材的“简单加和”陷阱

前文提到 Bragg 加和对化合物不一定准，这里再展开两句。对高能区，Bragg 规则误差往往在百分之几以内；对低能区，误差可能上升到百分之几十。原因不复杂：阻止本领依赖靶中电子的激发能谱，而激发能谱受化学键直接影响。氧化物、氮化物、碳化物中电子局域化程度差异巨大，价电子对大阻止本领的贡献在低能区尤其显著。

如果实验条件允许，最好用自由薄膜或已知厚度标样测量透射能损，直接反推电子阻止。比如用自支撑的 Si₃N₄ 或碳膜做透射实验，简单、直观，又能为整个模拟提供锚点。没有实验条件时，至少要在论文中注明“电子阻止按 SRIM 默认处理，可能带来约 X% 误差”，让人心里有数。

5.3 当电子阻止“看起来不对”的时候

如果你发现模拟的射程和实验差得离谱，先别急着怀疑数据库。按优先级排查：第一，单位对不对；第二，靶密度对不对；第三，入射离子能量有没有写错（比如把 keV 和 MeV 搞混）；第四，靶原子序数和质量数有没有填对；第五，才轮到电子阻止数据库的精度问题。

还有一个我自己的习惯：不要只依赖 SRIM 一个来源。国际原子能机构（IAEA）的带电粒子相互作用数据库、NIST 的 PSTAR/ASTAR，以及学术文献里针对特定离子-靶组合的测量数据，都值得对照。两个独立数据源如果对上了，你的模拟结论才算真正站得住。电子阻止本领这个参数，难点不在于定义一个公式，而在于你知道什么时候该相信模型、什么时候该相信实验、什么时候该老老实实做实验。只要把这一层搞清楚，离子束相关的深度问题就至少解决了一半。