阿伦尼乌斯图为什么会弯？曲率来源、局部表观活化能与三类场景的判据

说明：本文华算科技主要介绍阿伦尼乌斯图曲率的物理来源、局部表观活化能的含义，以及有限温计算、多步反应和离子输运中的判断条件。

阿伦尼乌斯图的直线与曲率怎样定义？

在固定压力、组成和观测对象下，把不同温度的速率常数 k 或扩散系数 D 取对数，再以倒数温度为横轴，就得到阿伦尼乌斯图。倒数温度增大对应降温，常见热激活过程向右逐渐变慢。纵轴可以直接标 ln k，也可以把 k 放在对数刻度上；两种画法都把跨越多个数量级的变化展开。

经典形式 k=A exp[−E_a/(RT)] 在A 与 E_a 均为常数时给出直线：ln k=ln A−E_a/(RT)。R 是气体常数，T 用 K 表示，E_a 为摩尔能量。曲线弯曲表示不同位置的切线斜率发生变化；局部量 E_app(T)=−R d ln k/d(1/T) 随温度改变。单次原子跃迁用 eV 计能量时，R 换为玻尔兹曼常数 k_B。

钨自扩散的例子采用 T_m/T 作横轴，T_m 为熔点，右侧面板把左侧曲线的切线斜率单独画出。曲线局部斜率与全温区拟合斜率描述不同对象：前者对应某个温度附近的响应，后者压缩了整个拟合区间。计算与实验各自采用对应熔点归一化，斜率换成能量时还须计入 T_m。

以 x=1000/T 代替 1/T 时，ln k 的斜率为 −E_a/(1000R)；改用 log₁₀ k，斜率还多出 1/ln 10。坐标的常数缩放会改变斜率数值，但保留曲线的弯曲形状。用摄氏温度倒数代替绝对温度倒数，则改变了函数关系，即使原过程服从经典公式也可能弯曲。

电导数据常画 ln(σT)，因为独立载流子的 Nernst–Einstein 关系含有 1/T 因子。ln σ 与 ln(σT)相差 ln T，在载流子浓度与电荷固定、所用扩散量匹配时，前者得到的局部摩尔活化能比后者少 RT。不同论文若用了不同纵轴，斜率差异中就包含这项数学贡献。

同一迁移或反应机制为什么也能产生曲率？

过渡态理论写作 k=κ(k_BT/h) exp[−ΔG^‡/(RT)]，其中 h 为普朗克常数，κ 为传输系数。固定标准态、κ 近似恒定时，局部活化能为E_app=ΔH^‡+RT。即使活化焓保持恒定，显式的 T 前因子会使 ln k 对 1/T 出现小幅曲率；ln(k/T) 才消除了这项前因子。

对于平衡空位参与的晶体自扩散，D 包含空位出现的概率和已有空位的跳跃率。指数项涉及G_form(T)+G_mig(T)，还乘有跳距、关联因子和尝试频率。体心立方钨的有限温计算显示，单空位形成与迁移自由能均有非线性温度变化，原子始终按单空位机制迁移，扩散图仍然弯曲。

准谐近似让振动频率随晶格体积变化，但原子在给定体积下仍沿二次势能振动。非谐振动使高温采样离开这一区域，体相、空位和迁移过渡态的自由能修正随之不同。钨的形成自由能明显降低，迁移自由能则呈现更复杂的非单调变化；这里的曲率来自有限温自由能差。

钨在 3000 K 的轨迹快照中，体相、空位和过渡态的电子态密度在费米能级附近不同，热平均曲线也偏离静态晶格结果。电子占据由费米分布决定，费米能级附近可激发电子态的差异，会改变三个构型的电子熵与电子自由能。

DFT 静态 NEB 得到的是指定晶格和电子状态下的势能差。预测宽温区速率时，初态与过渡态的振动自由能、电子激发和热膨胀进入同一温度下的能量比较。机器学习势延长轨迹，再通过热力学积分和 DFT 校正计算自由能，可检验曲率是否在固定迁移机制内产生。

活化焓的温度导数为活化热容 ΔC_p^‡，即过渡态与初态的热容差。κ 恒定时，dE_app/dT=ΔC_p^‡+R；ΔC_p^‡ 的正负和大小规定活化焓怎样随温度变化。若涉及轻原子转移，低温隧穿还可改变 κ(T)，即使化学反应路径保持原样，速率会偏离经典热激活关系。

多步反应和材料状态变化怎样改变曲线？

设两个独立并行过程具有恒定前因子和能垒，且观测速率常数确为 k=k₁+k₂，则 E_app=w₁E₁+w₂E₂，w_i=k_i/k。各过程的速率权重随温度改变，能垒保持不变也能产生曲率。此模型中 ln k 对 1/T 的二阶导数为能垒加权方差除以 R²，数值非负；曲率方向仅适用于这些限定条件。

表面反应测得的 r 常含 k(T) 与覆盖度 θ(T) 的乘积。若 r=kθ，则覆盖度对温度的响应额外贡献 −R d ln θ/d(1/T)。Pt 上水煤气变换的模型对比给出速率、CO 覆盖度和速率控制度：忽略吸附物相互作用，会改变 CO 占位与空位比例，连同高低温速率斜率一起改变。

相邻位点共同参与水解离时，空位对出现概率也进入速率，单个位点覆盖度未必足够描述反应。串联步骤的等待时间则按另一种规则累加：在简化的不可逆循环中，周转时间可写成 τ=Σ(1/k_i)，总周转率为 1/τ。并行速率相加与串联时间相加，会对各步温度响应产生不同权重。

LGPS，即 Li₁₀Ge(PS₆)₂，沿 c 方向与 ab 平面的扩散分量在降温时变化不同。方向分辨的扩散系数能显示总 D 中各方向的贡献：高温与低温数据可能各自接近直线，而中间温区发生连续偏离，整条曲线呈现跨温区的弯曲。

这组 LGPS 模拟的高低温直线斜率相近，外推截距却存在明显差别，350—500 K 间的方向分量也改变。斜率相近而截距改变可由迁移熵和离子运动关联变化产生。子晶格转变的解释要具体对应 Li 位点分布和方向分量；不同组成、晶界与模型设置下，转变温区未必相同。

同一模型的 Li–P 径向分布在低温保留约 3.0 Å 与 3.75 Å 的两个峰，高温下逐渐展宽并合。Li 与骨架原子的距离分布反映可访问局域环境的变化。RDF 展宽本身也可能只是热振动增强；方向分辨扩散、位点占据和离子交叉关联的温度序列，进一步约束子晶格重排的解释。

在浓集离子体系中，某个离子的位移可能伴随邻近离子连续换位，离子交叉关联使电荷输运与自扩散具有不同温度响应。LGPS 的模拟显示低温离子运动关联减弱。该文以 H_R=D^*/D_σ 定义 Haven 比，自扩散与电荷扩散的比值随温度变化，室温电导的外推误差可能来自 D 和关联项。

怎样区分真实曲率与采样、传输和外推的影响？

重复温度序列可检验曲线偏离是否超出测量误差，直线拟合残差的连续偏向比几个点目测成弧更有信息：残差若随温度有规律地由正变负再变正，单一直线遗漏了温度响应；误差范围内随机散布则仍支持该温区的线性近似。两段直线的连接温度也应由额外结构或动力学信号约束。

Li₅PS₄Cl₂ 的例子直接计算不同温度的输运，再用高温点预测低温值；原图纵轴为 σT，红蓝线表示外推误差范围。直接低温点与外推线的偏离须相对于这一范围评价，远距离外推产生的偏差与真实曲率是两个不同来源。

低温 MD 中少量跳跃会使 D 的估计偏斜，未观察到跳跃的轨迹也无法直接取对数。独立轨迹的跳跃数、长时间 MSD 和分块误差限定低温点的可信范围。机器学习势还须覆盖低温位点与迁移构型；仅在高温训练集上复现能量和力，无法保证低温扩散统计。

催化实验中的测量温度若偏离活性位局部温度，放热引起的热点会改变作图横轴。升温后反应快于物质补给，传质控制会使表观速率对温度的响应减弱。流量、颗粒尺寸和反应物分压改变后曲率同步变化，支持传输或覆盖度参与；结构相变的解释则应对应同温区的衍射或局域结构信号。

假设两种模型在 800 K 给出相同速率、前因子均不随温度变化，但活化能相差 0.10 eV，则到 300 K 时速率预测相差约 11 倍，来自 exp[(0.10 eV/k_B)(1/300 K−1/800 K)]。这一估计使用固定能垒模型；表面反应的 300 K 速率还含有该温度下的位点占据概率。