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

Какво представлява калкулаторът за енергия на активация?

Калкулаторът за енергия на активация намира енергията на активация (Ea) на химична реакция по две константи на скоростта, измерени при две различни температури. Той прилага уравнението на Арениус — стандартната формула, свързваща скоростта на реакцията с температурата. Потребителят въвежда константа на скоростта и температура за всяко от двете измервания, а инструментът изчислява Ea в килоджаули на мол (kJ/mol), като паралелно показва стойностите в kcal/mol и електронволти на молекула.

Какво представлява енергията на активация?

Всяка химична реакция трябва да преодолее енергийна бариера, преди реагентите да се превърнат в продукти. На молекулите им е необходимо минимално количество енергия, за да разкъсат старите връзки и да започнат да образуват нови. Това минимално количество е енергията на активация, означавана с Ea.

Полезна нагледна представа е топка, намираща се в долина. За да достигне следващата долина, топката първо трябва да се изкачи до билото и да го премине. Височината на билото е енергията на активация. Дори реакция, която като цяло отделя енергия, се нуждае от този първоначален тласък, за да започне.

Нагряването на реакцията не прави билото по-ниско. То дава на повече молекули достатъчно енергия, за да го преминат, така че всяка секунда реагират повече от тях. Затова дори повишаване на температурата с 10 °C често ускорява реакцията осезаемо.

Уравнението на Арениус

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

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

Тук A е предекспоненциалният фактор (приблизително — колко често молекулите се сблъскват с подходяща ориентация), R е моларната газова константа, а T е абсолютната температура в келвини. От ревизията на SI от 2019 г. 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. Тази линия се нарича графика на Арениус, а калкулаторът построява такава графика от въведените точки.

Решен пример

Нека разгледаме реакция с константа на скоростта k₁ = 0,0025 s⁻¹ при T₁ = 300 K и k₂ = 0,035 s⁻¹ при T₂ = 350 K.

1. Съотношение на скоростите: k₂/k₁ = 0,035 / 0,0025 = 14
2. ln(14) ≈ 2,63906
3. 1/T₁ − 1/T₂ = 1/300 − 1/350 = 1/2100 ≈ 4,76190 × 10⁻⁴ K⁻¹
4. Ea = 8,31446 × 2,63906 / (4,76190 × 10⁻⁴) ≈ 46 079 J/mol ≈ 46,08 кДж/mol

Същият резултат е приблизително 11,01 kcal/mol или около 0,478 eV за всяка реагираща молекула.

Как се тълкува резултатът

Катализаторите и ензимите ускоряват реакциите, като осигуряват път с по-ниска бариера. Например каталазата понижава бариерата за разлагането на водородния пероксид от около 75 кДж/mol до приблизително 8 кДж/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% от истинската стойност при умерен температурен диапазон, ако реакцията следва поведението на Арениус. Измерването при повече от две температури и прокарването на права подобрява допълнително точността.

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

Основното изчисление представлява едноредово аритметично изчисление. По-долу същата формула е дадена на няколко езика, като всеки от тях връща Ea в kJ/mol.

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

Източници

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 (10-о изд.). 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