活化能计算器 (阿伦尼乌斯方程)

什么是活化能计算器？

活化能计算器根据在两个不同温度下测得的两个速率常数，求出化学反应的活化能（Ea）。它运用阿伦尼乌斯方程，这是将反应速率与温度联系起来的标准公式。用户为两次测量分别输入一个速率常数和一个温度，工具便会求解出以千焦每摩尔（kJ/mol）为单位的Ea，并同时给出以千卡每摩尔（kcal/mol）和每分子电子伏特表示的数值。

什么是活化能？

每一个化学反应，在反应物转变为产物之前，都必须越过一道能量势垒。分子需要具备最低限度的能量，才能断开旧键并开始形成新键。这个最低能量就是活化能，记作Ea。

一个有用的比喻是一颗停在山谷中的球。要到达下一个山谷，这颗球必须先爬上并越过一道山脊。山脊的高度就相当于活化能。即使一个反应总体上会释放能量，它也仍然需要这样一次初始的推动才能启动。

给反应加热并不会使山脊变矮。它只是让更多的分子获得足以越过山脊的能量，因此每秒钟发生反应的分子数量增多。这就是为什么哪怕只把温度升高10°C，也常常能让反应速率明显加快。

阿伦尼乌斯方程

阿伦尼乌斯方程描述了反应速率常数k如何依赖于温度：

$k = A \cdot e^{-E_a/RT}$

这里的A是指前因子（大致表示分子以合适取向发生碰撞的频率），R是摩尔气体常数，T是以开尔文表示的绝对温度。自2019年国际单位制（SI）修订以来，R具有精确值：8.31446261815324 J/(mol·K)，通常被引用为四舍五入后的8.314。本计算器使用的是精确值。

单次测量k值无法将A与Ea分离开来，因为一个方程含有两个未知数。而在两个温度下分别测量k值就能解决这个问题，因为将两个方程相除时，A会被消去。

活化能公式

将T₂时的阿伦尼乌斯方程除以T₁时的同一方程，再对两边取自然对数，可得：

$\ln\!\left(\frac{k_2}{k_1}\right) = \frac{E_a}{R}\left(\frac{1}{T_1} - \frac{1}{T_2}\right)$

求解Ea：

$E_a = \frac{R\,\ln(k_2/k_1)}{\dfrac{1}{T_1} - \dfrac{1}{T_2}}$

这就是计算器所求解的方程。温度可以用开尔文、摄氏度或华氏度输入；工具会先将其换算为开尔文，因为该公式只对绝对温度成立。计算结果除以1000后以kJ/mol的形式给出，同一个数值也会以kcal/mol和每分子电子伏特的形式一并显示。

以ln k对1/T作图时，两次测量值会落在同一条直线上，该直线的斜率等于−Ea/R。这条直线称为阿伦尼乌斯图，计算器会根据输入的两点绘制出这样一条直线。

计算示例

以一个反应为例，其速率常数在T₁ = 300 K时为k₁ = 0.0025 s⁻¹，在T₂ = 350 K时为k₂ = 0.035 s⁻¹。

速率比：k₂/k₁ = 0.035 / 0.0025 = 14
ln(14) ≈ 2.63906
1/T₁ − 1/T₂ = 1/300 − 1/350 = 1/2100 ≈ 4.76190 × 10⁻⁴ K⁻¹
Ea = 8.31446 × 2.63906 / (4.76190 × 10⁻⁴) ≈ 46,079 J/mol ≈ 46.08 kJ/mol

同一结果约等于11.01 kcal/mol，或者说每个参与反应的分子约为0.478 eV。

如何解读结果

Ea（kJ/mol）	含义	典型示例
低于40	反应非常快，势垒很低	自由基反应、酸碱中和
40–100	室温下反应速率适中	大多数溶液相有机反应
100–200	反应缓慢，需要加热	狄尔斯-阿尔德反应、酯水解
高于200	若不强力加热则反应极慢	C–H键活化、固氮反应

催化剂和酶通过提供势垒更低的反应路径来加快反应速率。例如，过氧化氢酶能把过氧化氢分解反应的势垒从约75 kJ/mol降低到约8 kJ/mol。活化能取决于反应路径本身，而与温度无关。加热改变的是能够越过势垒的分子数量，而不是势垒本身的高度。

双点法的局限性

双点公式假设反应在两个测量点之间完全遵循阿伦尼乌斯方程。这一假设可能会失效：

曲线弯曲。 酶变性、低温下的量子隧穿效应，以及扩散限制反应，都可能使ln k对1/T的曲线发生弯曲。曲线弯曲意味着双点法算出的Ea只是该温度范围内的一个平均值，而不是一个固定的势垒。
负值Ea。 单一反应步骤不可能具有负的活化能。出现负值结果，通常意味着该反应包含多个步骤，机理中隐藏着一个对温度敏感的平衡过程。
更高的精度。 在四五个温度下测量k值并对所有数据点拟合一条直线，而不是只用两个点，能得到更可靠的斜率和更准确的Ea。

若需要更严谨的处理方式，过渡态理论用艾林方程 $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$ 来表示速率常数，该方程将能量势垒拆分为焓和熵各自的贡献。

如何得到准确的结果

正确输入温度。 系统接受摄氏度和华氏度，但底层运算只在开尔文下成立，因此工具会自动进行换算。手动计算时忘记这一步，是最常见的出错原因。
拉开两个温度的差距。 20–50°C的温差能给出较为稳定的结果。两个温度过于接近，会使速率比k₂/k₁趋近于1，微小误差经对数运算后会被放大。
保持速率常数单位一致。 公式中只用到比值k₂/k₁，因此k的单位会相互抵消。输入框标注的单位是s⁻¹，即一级反应速率常数的单位，但只要两个常数使用相同单位、来自同一类测量，使用M⁻¹s⁻¹或其他任何单位同样可行。

常见问题

升高温度会改变活化能吗？ 不会。对于给定的反应路径，Ea保持不变。升高温度会增加具备足够能量越过势垒的分子比例，这正是反应加快的原因。

活化能与焓变（ΔH）有什么区别？ Ea是反应必须攀越的势垒高度，ΔH是反应物与产物之间的净能量变化。一个反应可以拥有很高的势垒（Ea较大），却仍然以比起始状态更低的能量结束（ΔH为负，属于放热反应）。

为什么计算需要两个温度而不是一个？ 一次速率测量只能得到一个含有两个未知数（A和Ea）的方程。在另一个温度下再进行一次测量，就能让A在比值中被消去，从而得到一个只需求解Ea的方程。

可以输入哪些温度单位？ 开尔文、摄氏度或华氏度均可。计算器会先将所选单位换算为开尔文，再计算Ea，因为阿伦尼乌斯方程要求使用绝对温度。

双点法的准确度如何？ 在温度范围适中且反应符合阿伦尼乌斯行为的情况下，双点法通常能达到与真实值相差5%至10%以内的精度。测量两个以上的温度点并拟合直线，能进一步提高准确度。

代码实现

核心计算只是一行算式。以下用多种编程语言实现了同一公式，每种实现都返回以kJ/mol为单位的Ea。

1' Cells: A1=k1, A2=T1(K), A3=k2, A4=T2(K); result in kJ/mol
2=8.31446261815324*LN(A3/A1)/((1/A2)-(1/A4))/1000
3

1import math
2
3def activation_energy(k1, T1, k2, T2):
4    """Ea in kJ/mol from two rate constants at Kelvin temperatures T1, T2."""
5    return 8.31446261815324 * math.log(k2 / k1) / (1 / T1 - 1 / T2) / 1000
6
7print(f"{activation_energy(0.0025, 300, 0.035, 350):.2f} kJ/mol")
8

1function activationEnergy(k1, T1, k2, T2) {
2  // Ea in kJ/mol; T1, T2 in Kelvin
3  return (8.31446261815324 * Math.log(k2 / k1)) / (1 / T1 - 1 / T2) / 1000;
4}
5
6console.log(activationEnergy(0.0025, 300, 0.035, 350).toFixed(2), "kJ/mol");
7

1public class ActivationEnergy {
2    static double ea(double k1, double T1, double k2, double T2) {
3        return 8.31446261815324 * Math.log(k2 / k1) / (1 / T1 - 1 / T2) / 1000; // kJ/mol
4    }
5
6    public static void main(String[] args) {
7        System.out.printf("%.2f kJ/mol%n", ea(0.0025, 300, 0.035, 350));
8    }
9}
10

1using System;
2
3class ActivationEnergy {
4    static double Ea(double k1, double T1, double k2, double T2) =>
5        8.31446261815324 * Math.Log(k2 / k1) / (1 / T1 - 1 / T2) / 1000; // kJ/mol
6
7    static void Main() =>
8        Console.WriteLine($"{Ea(0.0025, 300, 0.035, 350):F2} kJ/mol");
9}
10

1#include <cmath>
2#include <cstdio>
3
4double ea(double k1, double T1, double k2, double T2) {
5    return 8.31446261815324 * std::log(k2 / k1) / (1 / T1 - 1 / T2) / 1000; // kJ/mol
6}
7
8int main() {
9    std::printf("%.2f kJ/mol\n", ea(0.0025, 300, 0.035, 350));
10}
11

1package main
2
3import (
4	"fmt"
5	"math"
6)
7
8func ea(k1, T1, k2, T2 float64) float64 {
9	return 8.31446261815324 * math.Log(k2/k1) / (1/T1 - 1/T2) / 1000 // kJ/mol
10}
11
12func main() {
13	fmt.Printf("%.2f kJ/mol\n", ea(0.0025, 300, 0.035, 350))
14}
15

1fn ea(k1: f64, t1: f64, k2: f64, t2: f64) -> f64 {
2    8.31446261815324 * (k2 / k1).ln() / (1.0 / t1 - 1.0 / t2) / 1000.0 // kJ/mol
3}
4
5fn main() {
6    println!("{:.2} kJ/mol", ea(0.0025, 300.0, 0.035, 350.0));
7}
8

1<?php
2function ea($k1, $T1, $k2, $T2) {
3    return 8.31446261815324 * log($k2 / $k1) / (1 / $T1 - 1 / $T2) / 1000; // kJ/mol
4}
5
6printf("%.2f kJ/mol\n", ea(0.0025, 300, 0.035, 350));
7

1def ea(k1, t1, k2, t2)
2  8.31446261815324 * Math.log(k2 / k1) / (1.0 / t1 - 1.0 / t2) / 1000 # kJ/mol
3end
4
5puts format("%.2f kJ/mol", ea(0.0025, 300, 0.035, 350))
6

1ea <- function(k1, T1, k2, T2) {
2  8.31446261815324 * log(k2 / k1) / (1 / T1 - 1 / T2) / 1000  # kJ/mol
3}
4
5cat(sprintf("%.2f kJ/mol\n", ea(0.0025, 300, 0.035, 350)))
6

参考资料

Laidler, K. J. (1984). "The development of the Arrhenius equation." J. Chem. Educ. 61(6), 494. https://doi.org/10.1021/ed061p494
Eyring, H. (1935). "The Activated Complex in Chemical Reactions." J. Chem. Phys. 3(2), 107. https://doi.org/10.1063/1.1749604
Atkins, P., & de Paula, J. (2014). Atkins' Physical Chemistry (10th ed.). Oxford University Press.
IUPAC. Compendium of Chemical Terminology (Gold Book), "activation energy." https://goldbook.iupac.org/terms/view/A00102
NIST. Fundamental Physical Constants — molar gas constant. https://physics.nist.gov/cgi-bin/cuu/Value?r

活化能计算器 (阿伦尼乌斯方程)