Calculator voor activeringsenergie (Arrhenius-vergelijking)

Wat is een calculator voor de activeringsenergie?

Een calculator voor de activeringsenergie bepaalt de activeringsenergie (Ea) van een chemische reactie uit twee snelheidsconstanten die bij twee verschillende temperaturen zijn gemeten. De calculator gebruikt de vergelijking van Arrhenius, de standaardformule die het verband legt tussen de reactiesnelheid en de temperatuur. De gebruiker voert voor elk van de twee metingen een snelheidsconstante en een temperatuur in. De tool berekent Ea in kilojoule per mol (kJ/mol), waarbij ook kcal/mol en elektronvolt per molecuul worden weergegeven.

Wat is activeringsenergie?

Elke chemische reactie moet een energiedrempel overschrijden voordat reagentia in producten kunnen veranderen. Moleculen hebben een minimale hoeveelheid energie nodig om oude bindingen te verbreken en nieuwe te beginnen vormen. Die minimale energie is de activeringsenergie, aangeduid met Ea.

Een bruikbaar beeld is een bal die in een dal ligt. Om het volgende dal te bereiken, moet de bal eerst een heuvel op en over een kam rollen. De hoogte van de kam is de activeringsenergie. Zelfs een reactie waarbij in totaal energie vrijkomt, heeft die eerste impuls nodig om op gang te komen.

Door een reactie te verwarmen wordt de kam niet lager. Meer moleculen krijgen daardoor genoeg energie om eroverheen te komen, zodat er elke seconde meer reageren. Daarom versnelt een verhoging van de temperatuur met zelfs 10 °C een reactie vaak merkbaar.

De vergelijking van Arrhenius

De vergelijking van Arrhenius beschrijft hoe de snelheidsconstante van een reactie, k, afhangt van de temperatuur:

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

Hier is A de pre-exponentiële factor (ongeveer hoe vaak moleculen botsen met een geschikte oriëntatie), R de molaire gasconstante en T de absolute temperatuur in kelvin. Sinds de 2019-herziening van het SI heeft R een exacte waarde: 8,31446261815324 J/(mol·K). Vaak wordt de afgeronde waarde 8,314 gebruikt. De calculator gebruikt de exacte waarde.

Met één meting van k kunnen A en Ea niet van elkaar worden onderscheiden, omdat één vergelijking twee onbekenden bevat. Metingen van k bij twee temperaturen lossen dit probleem op, omdat A wegvalt wanneer de twee vergelijkingen door elkaar worden gedeeld.

Formule voor de activeringsenergie

Als de vergelijking van Arrhenius bij T₂ wordt gedeeld door dezelfde vergelijking bij T₁ en vervolgens van beide kanten de natuurlijke logaritme wordt genomen, ontstaat:

$\ln\!\left(\frac{k_2}{k_1}\right) = \frac{E_a}{R}\left(\frac{1}{T_1} - \frac{1}{T_2}\right)$

Oplossen naar Ea geeft:

$E_a = \frac{R\,\ln(k_2/k_1)}{\dfrac{1}{T_1} - \dfrac{1}{T_2}}$

Dit is de vergelijking die de calculator evalueert. Temperaturen kunnen worden ingevoerd in kelvin, Celsius of Fahrenheit; de tool zet ze eerst om naar kelvin, omdat de formule alleen geldt voor absolute temperatuur. Het antwoord wordt gedeeld door 1000 om kJ/mol te rapporteren, en hetzelfde getal wordt ook weergegeven in kcal/mol en elektronvolt per molecuul.

Wanneer ln k wordt uitgezet tegen 1/T, liggen de twee metingen op een rechte lijn waarvan de helling gelijk is aan −Ea/R. Die lijn heet een Arrhenius-plot; de calculator tekent er een op basis van de ingevoerde punten.

Uitgewerkt voorbeeld

Neem een reactie met snelheidsconstante k₁ = 0,0025 s⁻¹ bij T₁ = 300 K, en k₂ = 0,035 s⁻¹ bij T₂ = 350 K.

1. Verhouding van de snelheidsconstanten: 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 kJ/mol

Hetzelfde resultaat is ongeveer 11,01 kcal/mol, of ongeveer 0,478 eV per reagerend molecuul.

Het resultaat lezen

Katalysatoren en enzymen versnellen reacties door een reactiepad met een lagere energiebarrière aan te bieden. Catalase verlaagt bijvoorbeeld de barrière voor de ontleding van waterstofperoxide van ongeveer 75 kJ/mol naar ongeveer 8 kJ/mol. Activeringsenergie is een eigenschap van het reactiepad, niet van de temperatuur. Door verwarming verandert hoeveel moleculen de barrière kunnen overwinnen, niet de hoogte van de barrière zelf.

Beperkingen van de tweepuntsmethode

De tweepuntsformule veronderstelt dat de reactie tussen de twee gemeten punten exact de vergelijking van Arrhenius volgt. Die veronderstelling kan tekortschieten:

• Gekromde grafieken. Denaturatie van enzymen, kwantumtunneling bij lage temperatuur en diffusiebeperkte reacties kunnen de lijn van ln k tegen 1/T allemaal krommen. Een gekromde grafiek betekent dat de tweepuntswaarde van Ea slechts een gemiddelde over dat temperatuurbereik is en geen vaste energiebarrière.
• Negatieve Ea. Eén afzonderlijke reactiestap kan geen negatieve activeringsenergie hebben. Een negatieve uitkomst betekent meestal dat de reactie uit meerdere stappen bestaat, waarbij een temperatuurgevoelig evenwicht het mechanisme verhult.
• Betere precisie. Door k bij vier of vijf temperaturen te meten en een lijn door alle punten te fitten, in plaats van slechts twee punten te gebruiken, krijgt men een betrouwbaardere helling en een nauwkeurigere Ea.

Voor een grondigere behandeling drukt de overgangstoestandstheorie de snelheidsconstante uit met de vergelijking van Eyring, $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, die de energiedrempel opsplitst in afzonderlijke bijdragen van enthalpie en entropie.

Een nauwkeurig resultaat verkrijgen

• Voer de temperaturen correct in. Celsius en Fahrenheit worden geaccepteerd, maar de onderliggende berekening werkt alleen met kelvin, dus de tool voert de omzetting automatisch uit. Deze stap met de hand vergeten is de meest voorkomende fout bij handmatige berekeningen.
• Kies temperaturen die voldoende ver uit elkaar liggen. Een verschil van 20–50 °C geeft een stabiel antwoord. Temperaturen die te dicht bij elkaar liggen, brengen de verhouding van de snelheidsconstanten k₂/k₁ in de buurt van 1, waardoor kleine fouten door de logaritme worden versterkt.
• Gebruik consistente eenheden voor de snelheidsconstanten. Alleen de verhouding k₂/k₁ komt in de formule voor, dus de eenheden van k vallen tegen elkaar weg. De invoervelden zijn gelabeld met s⁻¹, de eenheid van een eerste-ordesnelheidsconstante, maar M⁻¹s⁻¹ of een andere eenheid werkt ook, zolang beide constanten dezelfde eenheid hebben en afkomstig zijn van hetzelfde type meting.

Veelgestelde vragen

Verandert de activeringsenergie als de temperatuur wordt verhoogd? Nee. Ea blijft vast voor een bepaald reactiepad. Door de temperatuur te verhogen neemt de fractie van de moleculen met genoeg energie om de barrière te overwinnen toe, waardoor de reactie sneller verloopt.

Waarin verschilt activeringsenergie van de enthalpieverandering (ΔH)? Ea is de hoogte van de energiebarrière die een reactie moet overwinnen. ΔH is de nettoverandering in energie tussen reagentia en producten. Een reactie kan een hoge drempel hebben (grote Ea) en toch eindigen op een lager energieniveau dan waarop zij begon (negatieve ΔH, een exotherme reactie).

Waarom zijn voor de berekening twee temperaturen nodig in plaats van één? Eén snelheidsmeting levert één vergelijking met twee onbekenden, A en Ea. Met een tweede meting bij een andere temperatuur valt A weg uit de verhouding, waardoor een vergelijking overblijft die alleen naar Ea kan worden opgelost.

Welke temperatuureenheden kunnen worden ingevoerd? Kelvin, Celsius of Fahrenheit. De calculator zet de gekozen eenheid vóór het berekenen van Ea om naar kelvin, omdat de vergelijking van Arrhenius absolute temperatuur vereist.

Hoe nauwkeurig is de tweepuntsmethode? Over een gematigd temperatuurbereik komt de methode doorgaans binnen 5–10% van de werkelijke waarde, mits de reactie Arrheniusgedrag vertoont. Metingen bij meer dan twee temperaturen en het fitten van een lijn verbeteren de nauwkeurigheid verder.

Implementaties in code

De kernberekening bestaat uit één rekenregel. Hieronder staat dezelfde formule in verschillende programmeertalen; elk voorbeeld geeft Ea in 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

Referenties

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 (10e 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