Aktiveringsenergiräknare (Arrhenius-ekvation)

Vad är en kalkylator för aktiveringsenergi?

En kalkylator för aktiveringsenergi beräknar aktiveringsenergin (Ea) för en kemisk reaktion utifrån två hastighetskonstanter som uppmätts vid två olika temperaturer. Den använder Arrhenius ekvation, den standardformel som kopplar en reaktions hastighet till temperaturen. Användaren anger en hastighetskonstant och en temperatur för var och en av de två mätningarna, och verktyget löser ut Ea i kilojoule per mol (kJ/mol), medan kcal/mol och elektronvolt per molekyl visas bredvid.

Vad är aktiveringsenergi?

Alla kemiska reaktioner måste ta sig över en energibarriär innan reaktanterna kan omvandlas till produkter. Molekyler behöver en minsta mängd energi för att bryta gamla bindningar och börja bilda nya. Denna minsta energi är aktiveringsenergin, som skrivs Ea.

En användbar bild är en boll som ligger i en dal. För att nå nästa dal måste bollen först rulla upp och över en ås. Åsens höjd är aktiveringsenergin. Även en reaktion som totalt sett avger energi behöver denna första knuff för att komma igång.

Uppvärmning av en reaktion gör inte åsen lägre. Den ger fler molekyler tillräckligt med energi för att ta sig över den, så att fler av dem reagerar varje sekund. Därför ökar en höjning av temperaturen med även 10 °C ofta reaktionshastigheten märkbart.

Arrhenius ekvation

Arrhenius ekvation beskriver hur en reaktions hastighetskonstant, k, beror på temperaturen:

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

Här är A den preexponentiella faktorn (ungefär hur ofta molekyler kolliderar i en gynnsam riktning), R är den molära gaskonstanten och T är den absoluta temperaturen i kelvin. Sedan 2019 års revidering av SI har R ett exakt värde: 8,31446261815324 J/(mol·K). Det anges ofta avrundat till 8,314. Kalkylatorn använder det exakta värdet.

En enda mätning av k kan inte skilja A från Ea, eftersom en ekvation har två obekanta. Genom att mäta k vid två temperaturer kan detta lösas, eftersom A försvinner när de två ekvationerna divideras med varandra.

Formel för aktiveringsenergi

Om Arrhenius ekvation vid T₂ divideras med samma ekvation vid T₁ och man sedan tar den naturliga logaritmen av båda leden får man:

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

Om man löser ut Ea får man:

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

Det är denna ekvation som kalkylatorn beräknar. Temperaturer kan anges i kelvin, Celsius eller Fahrenheit; verktyget omvandlar dem först till kelvin, eftersom formeln endast gäller för absolut temperatur. Svaret divideras med 1000 för att anges i kJ/mol, och samma tal visas även i kcal/mol och elektronvolt per molekyl.

Om de två mätningarna plottas som ln k mot 1/T ligger de på en rät linje vars lutning är −Ea/R. Denna graf kallas ett Arrheniusdiagram, och kalkylatorn ritar ett sådant utifrån de angivna punkterna.

Räkneexempel

Betrakta en reaktion med hastighetskonstanten k₁ = 0,0025 s⁻¹ vid T₁ = 300 K, och k₂ = 0,035 s⁻¹ vid T₂ = 350 K.

1. Hastighetskvot: 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

Samma resultat motsvarar ungefär 11,01 kcal/mol, eller cirka 0,478 eV för varje reagerande molekyl.

Så tolkas resultatet

Katalysatorer och enzymer påskyndar reaktioner genom att erbjuda en reaktionsväg med lägre barriär. Katalas sänker till exempel barriären för nedbrytning av väteperoxid från omkring 75 kJ/mol till ungefär 8 kJ/mol. Aktiveringsenergin hör till reaktionsvägen, inte till temperaturen. Uppvärmning ändrar hur många molekyler som tar sig över barriären, inte själva barriärens höjd.

Begränsningar hos tvåpunktsmetoden

Tvåpunktsformeln förutsätter att reaktionen följer Arrhenius ekvation exakt mellan de två uppmätta punkterna. Det antagandet kan brista:

• Böjda diagram. Enzymdenaturering, kvanttunnling vid låg temperatur och diffusionsbegränsade reaktioner kan alla ge upphov till en krökt linje för ln k som funktion av 1/T. En krökt kurva innebär att tvåpunktsvärdet för Ea bara är ett medelvärde över det aktuella temperaturintervallet, inte en fast barriär.
• Negativ Ea. Ett enskilt reaktionssteg kan inte ha negativ aktiveringsenergi. Ett negativt resultat betyder vanligtvis att reaktionen består av flera steg, där en temperaturkänslig jämvikt döljer sig i mekanismen.
• Bättre precision. Om k mäts vid fyra eller fem temperaturer och en linje anpassas till alla punkter, i stället för att bara använda två punkter, fås en tillförlitligare lutning och ett noggrannare Ea-värde.

För en mer rigorös behandling uttrycker övergångstillståndsteorin hastighetskonstanten med Eyrings ekvation, $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, som delar upp energibarriären i separata bidrag från entalpi och entropi.

Så får man ett noggrant resultat

• Ange temperaturerna korrekt. Celsius och Fahrenheit accepteras, men de underliggande beräkningarna utförs endast med temperaturer i kelvin, så verktyget omvandlar dem automatiskt. Att glömma detta steg vid manuell beräkning är den vanligaste felkällan.
• Ha tillräckligt stort avstånd mellan temperaturerna. Ett intervall på 20–50 °C ger ett stabilt svar. Temperaturer som ligger för nära varandra gör att hastighetskvoten k₂/k₁ närmar sig 1, och små fel förstärks av logaritmen.
• Använd konsekventa enheter för hastighetskonstanten. Endast kvoten k₂/k₁ ingår i formeln, så enheterna för k förkortas bort. Inmatningsfälten är märkta s⁻¹, enheten för en hastighetskonstant av första ordningen, men M⁻¹s⁻¹ eller någon annan enhet fungerar, så länge båda konstanterna använder samma enhet och kommer från samma typ av mätning.

Vanliga frågor

Ändrar en höjning av temperaturen aktiveringsenergin? Nej. Ea förblir konstant för en given reaktionsväg. En höjning av temperaturen ökar andelen molekyler som har tillräckligt med energi för att ta sig över barriären, vilket gör att reaktionen går snabbare.

Hur skiljer sig aktiveringsenergin från entalpiförändringen (ΔH)? Ea är höjden på den barriär som reaktionen måste ta sig över. ΔH är den totala energiförändringen mellan reaktanter och produkter. En reaktion kan ha en hög barriär (stort Ea) och ändå sluta på en lägre energinivå än den började på (negativ ΔH, en exoterm reaktion).

Varför kräver beräkningen två temperaturer i stället för en? En hastighetsmätning ger en ekvation med två obekanta, A och Ea. En andra mätning vid en annan temperatur gör att A försvinner ur kvoten, så att en ekvation som endast kan lösas för Ea återstår.

Vilka temperaturenheter kan anges? Kelvin, Celsius eller Fahrenheit. Kalkylatorn omvandlar den valda enheten till kelvin innan Ea beräknas, eftersom Arrhenius ekvation kräver absolut temperatur.

Hur noggrann är tvåpunktsmetoden? Den hamnar vanligtvis inom 5–10% av det sanna värdet över ett måttligt temperaturintervall, förutsatt att reaktionen följer Arrheniusbeteende. Noggrannheten förbättras ytterligare om fler än två temperaturer mäts och en linje anpassas.

Kodimplementationer

Grundberäkningen består av en enda aritmetisk rad. Samma formel visas nedan på flera språk, och varje implementation returnerar Ea i 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

Referenser

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