Calculator voor activeringsenergie (Arrhenius-vergelijking)

Wat deze calculator doet

Reacties gaan sneller wanneer je ze verwarmt. Deze tool zet die alledaagse waarneming in een getal: meet hoe snel een reactie verloopt bij twee temperaturen, en je krijgt de activeringsenergie (Ea) van de reactie — de grootte van de energie "heuvel" die moleculen moeten beklimmen om te reageren.

Je hebt vier waarden nodig: twee reactiesnelheidsconstanten (k₁, k₂) en de twee temperaturen waarop ze zijn gemeten (T₁, T₂). Voer temperaturen in in Kelvin, Celsius of Fahrenheit — de calculator converteert voor je. Het resultaat is in kJ/mol, de standaardeenheid in scheikunde.

Activeringsenergie in eenvoudige termen

Stel je een bal voor die in een vallei ligt. Om hem in de volgende vallei te rollen, moet je hem eerst over een bergrug heen duwen. Activeringsenergie is de hoogte van die bergrug. Zelfs reacties die over het geheel genomen energie vrijgeven, hebben deze eerste push nodig om de beginbindingen te breken.

Een reactie verwarmen verlaagt de bergrug niet — het geeft meer moleculen genoeg energie om erover heen te gaan, dus meer van hen reageren per seconde. Daarom kan een stijging van 10 °C de dingen merkbaar versnellen.

De formule

De Arrhenius-vergelijking verbindt een reactiesnelheidsconstante k met temperatuur:

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

waarbij A de pre-exponentiële factor is, R = 8.314 J/(mol·K) de gasconstante, en T de absolute temperatuur in Kelvin is.

Meet k bij twee temperaturen en het onbekende A valt weg wanneer je de verhouding neemt:

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

Dit is precies wat de calculator berekent (vervolgens gedeeld door 1000 om kJ/mol te rapporteren). Geometrisch is het de helling van een rechte lijn op een Arrhenius-plot van ln k versus 1/T: de helling is gelijk aan −Ea/R. De tool tekent die plot op basis van je twee punten.

Uitgewerkt voorbeeld

Met k₁ = 0.0025 s⁻¹ bij 300 K en k₂ = 0.035 s⁻¹ bij 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}$

Een betrouwbaar antwoord krijgen

• Gebruik Kelvin intern. Celsius of Fahrenheit kun je typen — de tool converteert — maar de berekening werkt alleen met absolute temperatuur. Dit vergeten is de meest voorkomende fout.
• Zet je temperaturen uit elkaar. Een verschil van 20–50 °C geeft stabiele resultaten. Punten die te dicht bij elkaar liggen (een k₂/k₁-verhouding dicht bij 1) vergroten de fout, omdat je de logaritme van een getal dicht bij 1 neemt.
• Houd eenheden consistent. Omdat je de verhouding k₂/k₁ gebruikt, vallen de eenheden van k weg — s⁻¹, M⁻¹s⁻¹, wat dan ook — zolang beide constanten dezelfde eenheid gebruiken en uit dezelfde kinetische analyse komen.

Het resultaat interpreteren

Katalysatoren en enzymen werken door een pad met lagere barrière aan te bieden. Catalase verlaagt bijvoorbeeld de barrière voor waterstofperóxide-ontleding van ~75 kJ/mol naar ongeveer 8 kJ/mol. De activeringsenergie is een eigenschap van het pad, niet van de temperatuur — verwarming verandert hoeveel moleculen de barrière passeren, niet de hoogte van de barrière.

Voorbij de twee-punts-methode

De twee-punts-formule gaat ervan uit dat er perfect Arrhenius-gedrag tussen je punten is. Wanneer dat niet lineair is:

• Non-Arrhenius-kromming verschijnt bij enzymdenaturatie, quantum-tunneling bij lage temperatuur of diffusiebeperkte reacties. Een gekromde ln k vs 1/T-plot is het teken — dan is de twee-punts-waarde een effectief gemiddelde, geen vaste barrière.
• Een negatieve Ea is onmogelijk voor een elementaire stap. Als je er een berekent, heeft het mechanisme meerdere stappen met een temperatuursgevoelige voorraadevenwicht.
• Voor nauwkeurigheid gebruikt transitiestate-theorie de Eyring-vergelijking $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, waarbij de barrière wordt opgesplitst in enthalpie- en entropiebijdragen. Voor de beste experimentele nauwkeurigheid meet je k bij 4–5 temperaturen en pas je de helling van de Arrhenius-plot aan via lineaire regressie in plaats van alleen twee punten te gebruiken.

Code-implementaties

De kern is één regel rekenkunde. Gerangschikt 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

Veelgestelde vragen

Verandert temperatuur de activeringsenergie? Nee. Ea is vast voor een bepaald pad. Het verhogen van de temperatuur verhoogt het deel van moleculen dat de barrière passeert, daarom gaat het tarief omhoog.

Hoe verschilt Ea van ΔH? Ea is de hoogte van de heuvel die je moet beklimmen om te reageren; ΔH is de netto verandering in hoogte van reactanten naar producten. Een reactie kan een hoge barrière hebben (grote Ea) en toch lager eindigen dan ze begon (negatieve ΔH, exotherm).

Waarom heb ik twee temperaturen nodig? Eén meting geeft één vergelijking met twee onbekenden (Ea en A). Een tweede meting laat A wegvallen zodat je direct Ea kunt oplossen.

Hoe nauwkeurig is de twee-punts-methode? Typisch binnen 5–10% over een matig temperatuurbereik, ervan uitgaande dat het Arrhenius-gedrag lineair is. Voor meer precisie pas je ln k vs 1/T aan over 4–5 temperaturen.

Referenties

1. Laidler, K. J. (1984). "De ontwikkeling van de Arrhenius-vergelijking." J. Chem. Educ. 61(6), 494. https://doi.org/10.1021/ed061p494
2. Eyring, H. (1935). "Het geactiveerde complex in chemische reacties." 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