3. 扩散路径与过渡态理论

在上一节中，我们把扩散问题归结为“原子从一个稳定位置迁移到另一个稳定位置需要越过一个能量障碍”。本节要把这个直观图像严格化：为什么扩散系数会随温度呈指数变化？过渡态到底对应势能面上的什么几何对象？为什么NEB方法要沿着一条“最小能量路径”去搜索？理解这些，是后面正确设置NEB计算、判断结果是否可信的前提。

一、阿伦尼乌斯关系：扩散的宏观温度依赖

实验上，扩散系数 D 与温度 T 的关系通常写成阿伦尼乌斯形式：

D = D₀ · exp(−Eₐ / k_B T)

其中 D₀ 是指前因子，Eₐ 是宏观激活能，k_B 是玻尔兹曼常数。取对数后得到 ln D = ln D₀ − Eₐ/(k_B T)，因此若把 ln D 对 1/T 作图，应得到一条直线，斜率给出 −Eₐ/k_B。这正是实验上提取激活能的标准做法。

关键问题是：Eₐ 与我们在原子尺度上算出的势垒是什么关系？在简单情形下，Eₐ 就等于迁移势垒 E_m；如果材料中还存在缺陷形成能 E_f（例如空位机制中空位的生成），则 Eₐ = E_f + E_m。NEB计算直接给出的是 E_m，因此解释实验激活能时必须先明确扩散机制，不能把两者混为一谈。

二、过渡态理论：从势能面到跃迁速率

过渡态理论（TST）把原子迁移看成体系在势能面上从一个极小值（初态）越过鞍点到达另一个极小值（末态）的过程。在谐振近似下，跃迁速率可写成：

k = (k_B T / h) · (q‡ / q₀) · exp(−Eₐ / k_B T)

其中 q‡ 是过渡态（鞍点）的配分函数，q₀ 是初态的配分函数，h 是普朗克常数。指数项由势垒高度决定，指前因子由振动频率等熵效应决定。对固体中的原子跳跃，常把指前因子近似为 ν* ≈ ∏ν_i / ∏ν_i‡，即初态与鞍点振动频率之比。若只关心势垒，NEB给出 Eₐ；若还要预测绝对速率，则需进一步做频率计算。

三、鞍点与最小能量路径

在N维势能面上，初态和末态都是局部极小值。连接两者的路径中，能量最高点最低的那条路径称为最小能量路径（MEP）。MEP上的能量最高点就是鞍点，它沿反应坐标方向是极大值，而在所有垂直于反应坐标的方向上是极小值——这正是“一阶鞍点”的定义。

为什么必须找MEP而不是随便一条路径？因为扩散速率由最低势垒决定，原子实际迁移会优先走能量最低的通道。如果NEB收敛到一条非MEP路径，得到的势垒会偏高，进而低估扩散系数。因此，判断NEB结果是否可信，首先要看路径是否收敛到MEP，其次看鞍点是否为一阶鞍点。

四、NEB方法的基本假设与实现要点

NEB（Nudged Elastic Band）方法的核心思想是：在初态和末态之间插入若干中间构型（images），用弹簧力把它们串成一条“弹性带”，同时让每个image在垂直于路径方向受力，使整条带逐渐松弛到MEP。其基本假设包括：

• 路径连续且可参数化：初末态之间不存在其他更低的通道，MEP是单一路径。
• 弹簧力只沿路径切向：切向弹簧力维持image分布均匀，不干扰垂直方向的真实力。
• 垂直方向力来自真实势能面：每个image沿垂直于路径方向受力 F⊥ = −∇E + (∇E·τ̂)τ̂，其中 τ̂ 是路径切向单位矢量。
• 鞍点由最高image近似：收敛后能量最高的image接近鞍点，但通常还需climbing image（CI-NEB）精确定位。

实际操作中，一个典型NEB输入包含：初末态结构、image数量（一般5～9个）、弹簧常数（如5 eV/Ų）、收敛判据（如最大力 < 0.05 eV/Å）。以空位在Cu(111)表面迁移为例，若初末态间距约2.5 Å，取7个image，初始路径用线性插值生成。运行后观察各image能量，若最高点两侧能量对称下降，且最大力收敛，则路径可信。若最高image的力仍偏大，应增加image数量或改用CI-NEB。

五、从势垒到扩散系数：一个可复算的示例

假设NEB计算得到某体系迁移势垒 Eₐ = 0.65 eV，指前因子取 D₀ = 1×10⁻³ cm²/s（示例假设值）。在 T = 600 K 时：

k_B T = 8.617×10⁻⁵ eV/K × 600 K ≈ 0.0517 eV
Eₐ / k_B T = 0.65 / 0.0517 ≈ 12.57
D = 1×10⁻³ × exp(−12.57) ≈ 1×10⁻³ × 3.5×10⁻⁶ ≈ 3.5×10⁻⁹ cm²/s

若温度升到 900 K，k_B T ≈ 0.0776 eV，Eₐ/k_B T ≈ 8.38，D ≈ 1×10⁻³ × 2.3×10⁻⁴ ≈ 2.3×10⁻⁷ cm²/s。可以看到，温度提高300 K，扩散系数增大约66倍。这个数量级变化正是阿伦尼乌斯关系的直观体现，也说明势垒的微小误差会被指数放大——NEB计算中0.05 eV的偏差在600 K下就会导致扩散系数约2.6倍的差异。

因此，在工程计算中，不仅要关注NEB是否收敛，还要评估势垒对计算参数（交换关联泛函、k点网格、截断能、image数量）的敏感性。通常建议对关键体系做收敛性测试，确认势垒变化小于0.05 eV后再用于宏观扩散系数的预测。