什么是阿伦尼乌斯公式？数学形式、线性拟合方法、活化能与前因子的物理来源及多尺度对应条件

说明：本文华算科技主要介绍阿伦尼乌斯公式的数学形式、线性拟合方法、活化能与前因子的物理来源，以及它与第一性原理反应能垒、分子动力学和实验表观速率之间的对应条件。

阿氏公式中k/A/Ea与T指何？

阿伦尼乌斯公式写作 k(T)=Aexp[−E_a/(RT)]。k 是某条动力学方程中定义的速率常数，T 是绝对温度，E_a 是由温度响应得到的活化能，A 是前因子。公式描述的是温度改变时 k 的变化幅度；反应物浓度、表面位点数或载流子浓度是否已经包含在 k 中，由采用的速率方程决定。

图1. NH₄Cl/NaNO₂ 反应的 ln K–1/T 直线；两种速率单位给出相同斜率和不同截距。DOI：10.1021/acsomega.2c07466

量纲必须随速率级数写明。一级反应的 k 常用 s⁻¹，二级反应可用 L·mol⁻¹·s⁻¹；A 与 k 具有相同量纲。采用摩尔活化能时，R=8.314 J·mol⁻¹·K⁻¹，E_a 与 RT 都以 J·mol⁻¹ 表示。单次原子跃迁若以 eV 计能量，指数项写成 E/(k_BT)，其中 k_B 常用 eV·K⁻¹。

在不同温度的比较中，指数项控制热激活过程对温度的敏感程度：A 固定时，E_a 增大使同一温度下的 k 降低；E_a 固定时，A 改变会平移整条速率曲线。A 可包含振动尝试频率、反应构型出现概率、活化熵或输运中的几何因子，并非一个脱离过程类型的通用常数。

取自然对数得到 ln k=ln A−(E_a/R)(1/T)。以 1/T 为横轴、ln k 为纵轴，直线斜率 m=−E_a/R，截距 b=ln A。若纵轴使用 log₁₀k，斜率改为 −E_a/(2.303R)。横轴向右对应温度降低，这一点决定了斜率的负号。

活化能的微分定义为 E_a(T)=RT²d(ln k)/dT=−Rd(ln k)/d(1/T)。直线对应温区内近似恒定的 E_a，曲线则对应随温度改变的局部斜率。电导、扩散或介电弛豫研究有时使用 ln(σT)、ln D 或 ln f 作纵轴；纵轴形式来自各自的输运方程，斜率换算应沿用论文给出的数学表达。

图1中的黑色和红色数据采用不同的速率单位，斜率几乎相同，截距相差一个常数。单位换算不会改动由斜率求得的 E_a，却会改变 A 的数值与单位。比较前因子时，速率定义、反应级数和单位必须一致。

温度序列如何给出阿氏斜率与截距？

每个温度点应从同一种动力学模型得到 k。浓度—时间曲线可按反应级数拟合，催化实验可由低转化区的微分速率获得表观速率常数，分子动力学可由事件计数或物种随时间变化反演速率。温度改变之外的状态量要保持可比，包括相组成、表面覆盖度、反应物分压、电极电位、载流子浓度定义与统计窗口。

固定转化程度对应的特征时间 t 常满足 t∝1/k，于是 ln t 对 1/T 的斜率为正，数值接近 E_a/R。这种写法依赖相同转化阈值和相同速率方程；转化阈值改变若引入诱导期、传热滞后或副反应，斜率会成为过程的表观量。

速率常数的误差会传递到 ln k，温度误差又会通过 1/T 改变横坐标。重复实验或独立轨迹可估计每个温度点的方差，随后采用带不确定度的线性回归。低温点间隔在 1/T 轴上被放大，恒温偏差和温度标定误差可能显著改变斜率，尤其是在温区狭窄或数据点较少时。

反应分子动力学给出的原始对象通常是一组轨迹。HOOH 分解模拟中，不同温度下 HOOH、H₂O、O₂、HOO 和 OH 的数量随时间变化；分析程序从这些曲线识别主要基元反应并拟合各自的速率常数。轨迹中的物种统计先于阿伦尼乌斯拟合，否则直线缺少明确的反应对象。

线性回归应报告温区、点数、斜率不确定度和残差。温区很窄时，斜率与截距高度相关，A 可能跨越多个数量级；活化能的不确定度由斜率标准误乘以 R 得到。较高的 R² 只表示该温区近似线性。残差若随 1/T 呈连续弯曲，单一直线会把温度依赖的动力学状态压缩成一组表观参数。

活化能前因子为何随模型与机制变？

过渡态理论把速率写成 k=κ(k_BT/h)exp(ΔS^‡/R)exp[−ΔH^‡/(RT)]。κ 是透射系数，k_BT/h 给出热频率尺度，ΔS^‡反映反应物到活化构型时可访问微观状态数的改变。在 ΔH^‡、ΔS^‡ 和 κ 随温度变化较小时，阿伦尼乌斯活化能近似满足 E_a≈ΔH^‡+RT。

静态 DFT 的 ΔE^‡、有限温度的 ΔG^‡ 与实验 E_a 分别属于电子能差、活化自由能和温度导数。零点能、振动熵、构型熵、溶剂重排、表面覆盖度及势垒回穿会使三者出现数值差异。把一个 NEB 电子能垒直接写成实验阿伦尼乌斯活化能，会遗漏这些项。

谐振近似下，反应物与过渡构型的振动频率进入配分函数比值。低频受限转动、表面平移和软晶格模式若按普通谐振子处理，ΔS^‡ 与 A 会对频率处理方式敏感。液相重排、固体中的非谐振动以及低温隧穿还会改变 κ 或 ΔG^‡(T)，由此产生超出简单常数前因子的温度响应。

前因子在不同过程中包含的物理量并不相同：气相双分子反应包含碰撞频率和取向，表面反应包含吸附态配分函数与活化熵，离子输运还会带入载流子数、跳距、关联因子和尝试频率。前因子承担指数项之外的速率尺度，其数值随速率常数定义改变。

Li⁺ 固态电解质中的 Meyer–Neldel 关系展示了 E_a 与前因子协同变化：某些组成沿近似直线分布，另一些组成因晶体对称或局域结构改变而偏离。较低的 E_a 若伴随较低的 A，指定温度下的电导提升幅度会被抵消；电导率仍由 Aexp[−E_a/(RT)] 的乘积确定。

材料计算中何状态使阿氏直线失效？

实验表观速率常把多个基元步骤、表面覆盖度和传质过程汇总为一个数值。Ni/Al₂O₃ 催化 CO₂ 甲烷化的热反应与介质阻挡放电反应具有不同 Arrhenius 斜率和 Eyring 截距；热反应在高温端出现传质影响，数据点开始偏离原有斜率。同一拟合区间内的控速状态应保持一致，拟合斜率才具有统一的动力学对象。

若实验采用总速率 r 代替速率常数，r=kΠa_i^n_i 中的活度 a_i、反应级数 n_i 和覆盖度都可能随温度改变。由 ln r 得到的是表观活化能，其中可含吸附平衡、载流子生成或传质的温度贡献。微分条件、转化率范围和归一化面积应随参数一并记录。

弯曲的 Arrhenius 图常对应温度升高时材料状态或速率表达式发生改变。固态离子导体可能出现载流子浓度变化、协同迁移、无序程度变化或相变；催化表面可能改变覆盖度、吸附构型和主导基元步骤；低温反应还可能受到量子隧穿影响。此时局部斜率定义的 E_a(T) 随温度变化，全温区拟合所得单个斜率会把转折两侧的温度响应合并。

折线或曲线可按材料状态划分温区。相变温度、玻璃化转变、表面重构、电导机制转换或扩散控制起点可作为分段位置，随后分别检验残差与斜率。分段拟合的依据是可识别的状态变化；仅为提高直线相关系数而任意缩短温区，会得到缺少物理归属的斜率。

DFT 或 NEB 给出指定原子路径的 ΔE^‡，振动频率与过渡态理论可计算 ΔG^‡(T) 和单步速率，有限温分子动力学提供事件统计，微观动力学模型把多步速率、覆盖度与反应物分压写入稳态方程。实验 Arrhenius 斜率对应这些过程共同形成的温度响应。相组成、覆盖度、反应路径和传质状态在拟合温区保持不变时，所得 E_a 与 ln A 才能标注为该状态下的阿伦尼乌斯参数。