Калькулятор энергии активации (уравнение Аррениуса)

Что делает этот калькулятор

Реакции ускоряются при нагревании. Этот инструмент превращает это повседневное наблюдение в число: измерьте скорость реакции при двух температурах, и он скажет вам энергию активации реакции (Ea) — размер энергетического "холма", который молекулы должны преодолеть, чтобы прореагировать.

Вам нужны четыре значения: две константы скорости (k₁, k₂) и две температуры, при которых они были измерены (T₁, T₂). Введите температуры в Кельвинах, градусах Цельсия или Фаренгейта — калькулятор преобразует их для вас. Результат выводится в кДж/моль, стандартную единицу в химии.

Энергия активации в простых словах

Представьте мяч, лежащий в долине. Чтобы катить его в следующую долину, сначала нужно толкнуть его вверх через гребень. Энергия активации — это высота этого гребня. Даже реакции, которые выделяют энергию в целом, все еще нуждаются в этом начальном толчке, чтобы разорвать исходные связи.

Нагревание реакции не снижает гребень — это дает большему количеству молекул достаточно энергии, чтобы преодолеть его, поэтому большее их количество реагирует каждую секунду. Вот почему подъем на 10 °C может заметно ускорить процесс.

Формула

Уравнение Аррениуса связывает константу скорости k с температурой:

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

где A — предэкспоненциальный множитель, R = 8.314 Дж/(моль·К) — газовая постоянная, и T — абсолютная температура в Кельвинах.

Измерьте k при двух температурах, и неизвестное A сокращается при взятии отношения:

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

Именно это калькулятор вычисляет (затем делит на 1000, чтобы получить кДж/моль). Геометрически это наклон прямой на диаграмме Аррениуса ln k относительно 1/T: наклон равен −Ea/R. Инструмент рисует эту диаграмму из двух ваших точек.

Решенный пример

С k₁ = 0.0025 s⁻¹ при 300 K и k₂ = 0.035 s⁻¹ при 350 K:

$E_a = \frac{8.314 \cdot \ln(0.035/0.0025)}{\tfrac{1}{300} - \tfrac{1}{350}} = \frac{8.314 \cdot 2.639}{4.76\times10^{-4}} \approx 46\,100 \text{ J/mol} = 46.1 \text{ kJ/mol}$

Получение надежного ответа

• Используйте Кельвин внутри. Цельсий или Фаренгейт подходят для ввода — инструмент преобразует — но математика работает только с абсолютной температурой. Забыть об этом — самая распространенная ошибка.
• Распределяйте ваши температуры. Разрыв 20–50 °C дает стабильные результаты. Точки, расположенные слишком близко (отношение k₂/k₁ близко к 1), вызывают ошибку, потому что вы берете логарифм числа близко к 1.
• Сохраняйте единицы согласованными. Потому что вы используете отношение k₂/k₁, единицы k сокращаются — s⁻¹, M⁻¹s⁻¹, что угодно — пока обе константы используют одну и ту же единицу и поступают из одного и того же кинетического анализа.

Чтение результата

Катализаторы и ферменты работают, предлагая путь с более низким барьером. Например, каталаза снижает барьер для разложения перекиси водорода с ~75 кДж/моль до примерно 8 кДж/моль. Энергия активации является свойством пути, а не температуры — нагревание изменяет, сколько молекул преодолевают барьер, но не высоту самого барьера.

За пределами двухточечного метода

Двухточечная формула предполагает совершенное поведение Аррениуса между вашими точками. Когда это не так:

• Кривизна, не подчиняющаяся Аррениусу, появляется при денатурации фермента, квантовом туннелировании при низкой температуре или реакциях, ограниченных диффузией. Кривая диаграмма ln k против 1/T — это признак — тогда значение двух точек является эффективным средним, а не фиксированным барьером.
• Отрицательная Ea невозможна для элементарного шага. Если вы вычислите одно, механизм имеет несколько шагов с температурно-чувствительным предравновесием.
• Для строгости, теория переходного состояния использует уравнение Айринга $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, разделяя барьер на энтальпию и энтропийные вклады. Для лучшей экспериментальной точности измерьте k при 4–5 температурах и подгоните наклон диаграммы Аррениуса с помощью линейной регрессии вместо использования только двух точек.

Реализации кода

Ядро — одна строка арифметики. Упорядочено Excel → Python → JavaScript → Java → C# → C++ → Go → Rust → PHP → Ruby → R.

1' Cells: A1=k1, A2=T1(K), A3=k2, A4=T2(K); result in kJ/mol
2=8.314*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.314 * 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.314 * 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.314 * 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.314 * 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.314 * 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.314 * 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.314 * (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.314 * 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.314 * 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.314 * 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

Часто задаваемые вопросы

Изменяет ли температура энергию активации? Нет. Ea фиксирована для данного пути. Повышение температуры увеличивает долю молекул, которые преодолевают барьер, что объясняет, почему скорость возрастает.

Чем Ea отличается от ΔH? Ea — это высота холма, который вы должны преодолеть, чтобы прореагировать; ΔH — это чистое изменение высоты от реагентов к продуктам. Реакция может иметь высокий барьер (большая Ea), но закончиться ниже, чем она началась (отрицательная ΔH, экзотермическая).

Почему мне нужны две температуры? Одно измерение дает одно уравнение с двумя неизвестными (Ea и A). Второе измерение позволяет A сократиться, чтобы вы могли решить для Ea напрямую.

Насколько точен двухточечный метод? Обычно в пределах 5–10% в умеренном диапазоне температур, предполагая линейное поведение Аррениуса. Для большей точности подгоните ln k против 1/T на 4–5 температурах.

Ссылки

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