Kalkulator aktivacijske energije (Arrheniusova enačba)

Kaj ta kalkulator počne

Reaktivnosti se pospešijo, ko jih ogrevate. To orodje pretvori opazovanje iz vsakdana v številko: merite, kako hitro poteka reaktivnost pri dveh temperaturah, in vam pove reaktivnostno aktivacijsko energijo (Ea) — velikost energijskega »hriba«, ki ga morajo molekule preplezati, da reagirajo.

Potrebujete štiri vrednosti: dve konstanti hitrosti (k₁, k₂) in dve temperaturi, pri katerih sta bili izmerjeni (T₁, T₂). Vnesite temperature v Kelvinih, Celziju ali Fahrenheitu — kalkulator jih pretvori za vas. Rezultat je v kJ/mol, standardni enoti v kemiji.

Aktivacijska energija na preprost način

Predstavljajte si žogo, ki sedi v dolini. Da bi jo valili v naslednje dolino, jo morate najprej potisniti čez greben. Aktivacijska energija je višina tega grebena. Tudi reaktivnosti, ki na splošno sprostijo energijo, potrebujejo to začetno potisnitev, da prekinejo začetne vezi.

Ogrevanje reaktivnosti ne znižuje grebena — več molekulam daje dovolj energije, da ga preseže, zato več njih reagira vsako sekundo. Zato 10 °C porast opazno pospeši stvari.

Formula

Arrheniusova enačba povezuje konstanto hitrosti k s temperaturo:

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

kjer je A predeksponentni faktor, R = 8.314 J/(mol·K) je gasna konstanta, in T je absolutna temperatura v Kelvinih.

Merite k pri dveh temperaturah in neznana A se izniči, ko vzamete razmerje:

$\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}}$

To je točno to, kar kalkulator ovrednoti (nato deli z 1000, da poroča v kJ/mol). Geometrično je to naklon premice na Arrheniusovem grafu ln k v primerjavi z 1/T: naklon je enak −Ea/R. Orodje narisuje ta graf iz vaših dveh točk.

Obdelani primer

Z k₁ = 0.0025 s⁻¹ pri 300 K in k₂ = 0.035 s⁻¹ pri 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}$

Zanesljiv odgovor

• Uporabite Kelvin interno. Celzij ali Fahrenheit sta dovolj za vnos — orodje pretvori — vendar matematika deluje samo s temperaturo. Pozaba tega je najpogostejša napaka.
• Razširite vaše temperature. Razlika 20–50 °C daje stabilne rezultate. Točke, ki so premalo narazen (razmerje k₂/k₁ blizu 1), povečajo napako, ker vzamete logaritem števila blizu 1.
• Ohranite dosledne enote. Ker uporabite razmerje k₂/k₁, se enote k izničijo — s⁻¹, M⁻¹s⁻¹, kaj koli — dokler obe konstanti uporabita isto enoto in sta iz iste kinetične analize.

Branje rezultata

Katalizatorji in encimi delujejo z zagotavljanjem poti z nižjo pregrado. Katalaza npr. znižuje pregrado za razgradnjo vodikove peroksida s ~75 kJ/mol na približno 8 kJ/mol. Aktivacijska energija je lastnost poti, ne temperature — ogrevanje spremeni, koliko molekul prečka pregrado, ne pa višine pregrade.

Izven metode dveh točk

Formula dveh točk predpostavlja popolno Arrheniusovo vedenje med vašimi točkami. Ko ni linearna:

• Nelinearna ukrivljenost Arrheniusa se pojavi pri denaturaciji encimov, kvantnem tuneliranju pri nizki temperaturi ali difuzijsko omejenih reaktivnostih. Ukrivljen graf ln k v primerjavi z 1/T to kaže — potem je vrednost dveh točk efektivno povprečje, ne fiksna pregrada.
• Negativna Ea je nemogoča za elementarni korak. Če jo izračunate, ima mehanizem več korakov s temperaturno občutljivo predravnovesjem.
• Za natančnost teorija prehodnega stanja uporablja Eyringovo enačbo $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, ki razdeli pregrado na prispevke entalpije in entropije. Za najboljšo eksperimentalno natančnost merite k pri 4–5 temperaturah in prilagodite naklon Arrheniusovega grafa z linearno regresijo namesto da bi uporabili le dve točki.

Implementacije kode

Jedro je ena vrstica aritmetike. Razporejeno 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

Pogosta vprašanja

Ali temperatura spremeni aktivacijsko energijo? Ne. Ea je fiksna za dano pot. Zvišanje temperature poveča delež molekul, ki presežejo pregrado, zato se hitrost poveča.

Kako se Ea razlikuje od ΔH? Ea je višina hriba, ki ga morate preplezati, da reagirate; ΔH je neto sprememba višine od reaktantov do produktov. Reaktivnost ima lahko visoko pregrado (velika Ea) in še vedno konča nižje, kot se je začela (negativna ΔH, eksotermna).

Zakaj potrebujem dve temperaturi? Ena meritev daje eno enačbo z dvema neznankama (Ea in A). Druga meritev pusti A, da se izniči, tako da lahko rešite Ea neposredno.

Kako natančna je metoda dveh točk? Običajno v 5–10% v zmjernem temperaturnem razponu, ob predpostavki linearnega Arrheniusovega vedenja. Za večjo natančnost prilagodite ln k v primerjavi z 1/T v 4–5 temperaturah.

Reference

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