Kalkulačka aktivačnej energie (Arrheniova rovnica)

Čo kalkulačka robí

Reaktácie sa zrýchľujú, keď ich zohrejete. Tento nástroj zmení to obvyklé pozorovanie na číslo: zmerajte, ako rýchlo sa reaktácia odohrá pri dvoch teplotách, a povie vám aktivačnú energiu reaktácie (Ea) — veľkosť energetickej "vyvýšeniny", ktorú musia prekonať molekuly, aby reagovali.

Potrebujete štyri hodnoty: dve rýchlostné konštanty (k₁, k₂) a dve teploty, pri ktorých boli meraané (T₁, T₂). Zadajte teploty v Kelvinoch, stupňoch Celzia alebo Fahrenheita — kalkulačka ich pre vás prevedie. Výsledok je v kJ/mol, štandardnej jednotke v chémii.

Aktivačná energia jednoduchšie povedané

Predstavte si guľu sediacu v údolí. Aby ste ju zvalili do ďalšieho údolia, musíte ju najprv pretlačiť cez chrbát. Aktivačná energia je výška tohto chrbáta. Aj reaktácie, ktoré celkovo uvoľňujú energiu, stále potrebujú tento počiatočný impulz na prerušenie počiatočných väzieb.

Zohriatie reaktácie nezníži chrbát — dá viac molekulám dostatok energie na jeho prekonanie, takže viac z nich reaguje každú sekundu. To je dôvod, prečo nárast o 10 °C môže výrazne urýchliť veci.

Vzorec

Arrheniova rovnica prepája rýchlostnú konštantu k s teplotou:

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

kde A je preexponenciálny faktor, R = 8.314 J/(mol·K) je plynová konštanta, a T je absolútna teplota v Kelvinoch.

Zmerajte k pri dvoch teplotách a neznáme A sa zrušia, keď vezmete pomer:

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

Toto je presne to, čo kalkulačka vyhodnocuje (potom delí 1000, aby hlásila kJ/mol). Geometricky je to sklon priamky na Arrheniover grafe ln k v závislosti od 1/T: sklon sa rovná −Ea/R. Nástroj vykreslí tento graf z vašich dvoch bodov.

Pracovný príklad

S k₁ = 0,0025 s⁻¹ pri 300 K a 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}$

Získanie spoľahlivej odpovede

• Interne používajte Kelvin. Celsius alebo Fahrenheit sú v poriadku na zadanie — nástroj konvertuje — ale matematika funguje len na absolútnej teplote. Zabudnutie na to je najčastejšou chybou.
• Rozprostrite vaše teploty. Medzera 20–50 °C dáva stabilné výsledky. Body príliš blízko pri sebe (pomer k₂/k₁ blízko 1) spôsobia chybu, pretože berieš logaritmus čísla blízko 1.
• Udržujte jednotky konzistentné. Pretože používate pomer k₂/k₁, jednotky k sa zrušia — s⁻¹, M⁻¹s⁻¹, čokoľvek — pokiaľ obe konštanty používajú tú istú jednotku a pochádzajú z rovnakej kinetickej analýzy.

Čítanie výsledku

Katalyzátory a enzýmy pracujú tým, že ponúkajú cestu s nižšou prekážkou. Kataláza napríklad znižuje prekážku rozkladu peroxidu vodíka z ~75 kJ/mol na približne 8 kJ/mol. Aktivačná energia je vlastnosť cesty, nie teploty — zohriatie zmení, koľko molekúl prekážku prekonajú, nie výšku prekážky.

Za metódou dvoch bodov

Vzorec dvoch bodov predpokladá dokonalé Arrheniovo správanie medzi vašimi bodmi. Keď to tak nie je:

• Non-Arrheniova krivka sa objavuje s denaturáciou enzýmu, kvantovým tunelovacím efektom pri nízkej teplote alebo reaktáciami limitovanými difúziou. Zakrivená ln k versus 1/T graf je signálom — potom je hodnota dvoch bodov efektívny priemer, nie pevná prekážka.
• Záporné Ea je nemožné pre elementárny krok. Ak to vypočítate, mechanizmus má viaceré kroky s teplotne citlivou pred-rovnováhou.
• Pre presnosť, teória prechodného stavu používa Eyringovu rovnicu $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, rozdelením prekážky na entalpiu a entropiu. Pre najlepšiu experimentálnu presnosť zmerajte k pri 4–5 teplotách a prispôsobte sklon Arrhenioveho grafu lineárnou regresiou namiesto použitia len dvoch bodov.

Code implementations

The core is one line of arithmetic. Ordered 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

FAQ

Does temperature change the activation energy? No. Ea is fixed for a given pathway. Raising the temperature increases the fraction of molecules that clear the barrier, which is why the rate goes up.

How is Ea different from ΔH? Ea is the height of the hill you must climb to react; ΔH is the net change in altitude from reactants to products. A reaction can have a tall barrier (large Ea) yet finish lower than it started (negative ΔH, exothermic).

Why do I need two temperatures? One measurement gives one equation with two unknowns (Ea and A). A second measurement lets A cancel so you can solve for Ea directly.

How accurate is the two-point method? Typically within 5–10% over a moderate temperature range, assuming linear Arrhenius behaviour. For more precision, fit ln k vs 1/T across 4–5 temperatures.

References

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