Kalkulátor aktivační energie (Arrheniusova rovnice)

Co tento kalkulátor dělá

Reakce se zrychlují při zahřívání. Tento nástroj změní každodenní pozorování v číslo: změřte, jak rychle probíhá reakce při dvou teplotách, a zjistíte aktivační energii (Ea) – velikost energetické „hory", kterou musí molekuly překonat, aby reagovaly.

Potřebujete čtyři hodnoty: dvě konstanty rychlosti (k₁, k₂) a dvě teploty, při kterých byly měřeny (T₁, T₂). Zadejte teploty v Kelvinech, Celsiech nebo Fahrenheitech – kalkulátor převede. Výsledek vychází v kJ/mol, standardní jednotka v chemii.

Aktivační energie v jednoduchých slovech

Představte si aktivační energii jako kopec. Aby jste valili míč do dalšího údolí, musíte jej nejprve vytlačit přes hřeben. Aktivační energie je výška hřebene. Ani reakce, které celkově uvolňují energii, nepotřebují tento počáteční impuls k přerušení startovních vazeb.

Zahřívání reakce neznižuje hřeben – dává více molekulám dost energie, aby jím překonaly, a tedy více z nich reaguje každou sekundu. Proto zvýšení o 10 °C může znatelně urychlit průběh.

Vzorec

Arrheniusova rovnice spojuje konstantu rychlosti k s teplotou:

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

kde A je pre-exponenciální faktor, R = 8,314 J/(mol·K) je plynová konstanta a T je absolutní teplota v Kelvinech.

Změřte k při dvou teplotách a neznámé A se ruší, když vezmete poměr:

$\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 přesně to, co kalkulátor vyhodnocuje (pak dělí 1000 pro hlášení v kJ/mol). Geometricky je to sklon přímky na Arrheniusův graf ln k versus 1/T: sklon se rovná −Ea/R. Nástroj kreslí tento graf z vašich dvou bodů.

Zpracovaný příklad

S k₁ = 0,0025 s⁻¹ při 300 K a k₂ = 0,035 s⁻¹ při 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}$

Získání spolehlivé odpovědi

• Používejte Kelvin interně. Celsius nebo Fahrenheit je v pořádku zadat – nástroj převádí – ale matematika funguje jen s absolutní teplotou. Zapomenutí je nejčastější chyba.
• Rozložte své teploty. Mezera 20–50 °C dává stabilní výsledky. Body příliš blízko sebe (poměr k₂/k₁ blízko 1) způsobují chyby, protože berete logaritmus čísla blízko 1.
• Udržujte jednotky konzistentní. Protože používáte poměr k₂/k₁, jednotky k se ruší – s⁻¹, M⁻¹s⁻¹, cokoliv – pokud obě konstanty používají stejnou jednotku a pocházejí z stejné kinetické analýzy.

Čtení výsledku

Katalyzátory a enzymy fungují nabídnutím cesty s nižší bariérou. Kataláza například snižuje bariéru pro rozklad peroxidu vodíku z ~75 kJ/mol na přibližně 8 kJ/mol. Aktivační energie je vlastnost dráhy, ne teploty – zahřívání změní, kolik molekul překoná bariéru, ne výšku bariéry.

Mimo dvouznamankovou metodu

Vzorec dvou bodů předpokládá dokonalé Arrheniusovo chování mezi vašimi body. Když tomu tak není:

• Nearrheniusovo zakřivení se objevuje s denaturací enzymů, kvantovým tunelováním při nízké teplotě nebo difúzí omezené reakce. Zakřivený graf ln k vs 1/T je známka – pak dvouznamankový hodnota je efektivní průměr, ne pevná bariéra.
• Negativní Ea je nemožný pro elementární krok. Pokud jej vypočítáte, mechanismus má více kroků s citlivostí na teplotu v pre-rovnováhu.
• Pro přesnost, teorie přechodového stavu používá Eyringovu rovnici $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, rozdělující bariéru na příspěvky entalpie a entropie. Pro nejlepší experimentální přesnost změřte k při 4–5 teplotách a přizpůsobte sklon Arrheniusova grafu lineární regresí místo pouhého použití dvou bodů.

Implementace kódu

Jádro je jeden řádek aritmetiky. Seřazeno 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

Často kladené otázky

Mění teplota aktivační energii? Ne. Ea je pevná pro danou dráhu. Zvýšení teploty zvyšuje podíl molekul, která překoná bariéru, a proto se rychlost zvýší.

Jak se liší Ea od ΔH? Ea je výška hory, kterou musíte zdolat, aby reagovaly; ΔH je čistá změna výšky od reaktantů k produktům. Reakce může mít vysokou bariéru (velký Ea), ale skončí nižší, než začala (negativní ΔH, exotermní).

Proč potřebuji dvě teploty? Jedno měření dává jednu rovnici se dvěma neznámými (Ea a A). Druhé měření umožňuje A zrušit, takže můžete vyřešit Ea přímo.

Jak přesná je dvouznamankový metoda? Typicky v rozmezí 5–10 % v mírném temperaturním rozsahu, za předpokladu lineárního Arrheniusova chování. Pro větší přesnost přizpůsobte ln k vs 1/T přes 4–5 teplot.

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