Aktivaatioenergia-laskin (Arrhenius-yhtälö)

Mikä on aktivaatioenergia-laskin?

Aktivaatioenergia-laskin määrittää kemiallisen reaktion aktivaatioenergian (Ea) kahdessa eri lämpötilassa mitattujen kahden nopeusvakion perusteella. Se käyttää Arrheniuksen yhtälöä, joka on reaktion nopeuden ja lämpötilan välinen standardiyhteys. Käyttäjä syöttää kummallekin mittaukselle nopeusvakion ja lämpötilan, ja työkalu ratkaisee Ea:n yksikössä kilojoulea moolia kohti (kJ/mol); sen rinnalla näytetään myös yksiköt kcal/mol ja elektronivolttia molekyyliä kohti.

Mitä aktivaatioenergia on?

Jokaisen kemiallisen reaktion on ylitettävä energiavalli, ennen kuin lähtöaineet voivat muuttua tuotteiksi. Molekyylit tarvitsevat vähimmäismäärän energiaa vanhojen sidosten katkaisemiseen ja uusien muodostamisen aloittamiseen. Tätä vähimmäisenergiaa kutsutaan aktivaatioenergiaksi, ja se merkitään Ea.

Hyödyllinen mielikuva on laaksossa lepäävä pallo. Päästäkseen seuraavaan laaksoon pallon on ensin noustava harjanteen yli. Harjanteen korkeus on aktivaatioenergia. Vaikka reaktio kokonaisuudessaan vapauttaisi energiaa, se tarvitsee silti alkuun tämän sysäyksen.

Reaktion kuumentaminen ei madalla harjannetta. Se antaa useammille molekyyleille riittävästi energiaa sen ylittämiseen, joten niitä reagoi enemmän joka sekunti. Siksi jopa lämpötilan nostaminen 10 °C:lla nopeuttaa usein reaktiota huomattavasti.

Arrheniuksen yhtälö

Arrheniuksen yhtälö kuvaa, miten reaktion nopeusvakio k riippuu lämpötilasta:

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

Tässä A on esieksponentiaalinen tekijä (likimäärin se, kuinka usein molekyylit törmäävät sopivassa orientaatiossa), R on moolinen kaasuvakio ja T on absoluuttinen lämpötila kelvineinä. SI-järjestelmän 2019 tehdyn tarkistuksen jälkeen R:llä on ollut täsmällinen arvo: 8,31446261815324 J/(mol·K). Se ilmoitetaan usein pyöristettynä arvoon 8,314. Laskin käyttää täsmällistä arvoa.

Yksi k:n mittaus ei voi erottaa A:ta Ea:sta, koska yhdessä yhtälössä on kaksi tuntematonta. Kun k mitataan kahdessa lämpötilassa, ongelma ratkeaa, koska A supistuu pois yhtälöiden jakamisen yhteydessä.

Aktivaatioenergian kaava

Kun Arrheniuksen yhtälö jaetaan lämpötilassa T₂ samalla yhtälöllä lämpötilassa T₁ ja tämän jälkeen otetaan yhtälön molemmilta puolilta luonnollinen logaritmi, saadaan:

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

Ratkaisemalla Ea:n suhteen:

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

Tätä yhtälöä laskin käyttää. Lämpötilat voidaan syöttää kelvineinä, celsiusasteina tai fahrenheitasteina; työkalu muuntaa ne ensin kelvineiksi, koska kaava pätee vain absoluuttiseen lämpötilaan. Tulos jaetaan luvulla 1000, jotta se voidaan ilmoittaa yksikössä kJ/mol, ja sama luku näytetään myös yksiköissä kcal/mol ja elektronivolttia molekyyliä kohti.

Kun ln k esitetään suureen 1/T funktiona, nämä kaksi mittausta sijoittuvat suoralle, jonka kulmakerroin on −Ea/R. Tätä suoraa kutsutaan Arrheniuksen kuvaajaksi, ja laskin piirtää sen syötettyjen pisteiden perusteella.

Laskuesimerkki

Tarkastellaan reaktiota, jonka nopeusvakio on k₁ = 0,0025 s⁻¹ lämpötilassa T₁ = 300 K ja k₂ = 0,035 s⁻¹ lämpötilassa T₂ = 350 K.

1. Nopeuksien suhde: 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

Sama tulos on noin 11,01 kcal/mol eli suunnilleen 0,478 eV kutakin reagoivaa molekyyliä kohti.

Tuloksen tulkitseminen

Katalyytit ja entsyymit nopeuttavat reaktioita tarjoamalla reaktioreitin, jonka energiavalli on matalampi. Esimerkiksi katalaasi madaltaa vetyperoksidin hajoamisen vallin noin arvosta 75 kJ/mol suunnilleen arvoon 8 kJ/mol. Aktivaatioenergia kuuluu reaktioreittiin, ei lämpötilaan. Kuumentaminen muuttaa niiden molekyylien määrää, jotka ylittävät vallin, mutta ei itse vallin korkeutta.

Kahden pisteen menetelmän rajoitukset

Kahden pisteen kaava olettaa, että reaktio noudattaa Arrheniuksen yhtälöä täsmällisesti näiden kahden mittauspisteen välillä. Oletus voi pettää seuraavissa tilanteissa:

• Kaarevat kuvaajat. Entsyymien denaturoituminen, kvanttitunnelointi matalassa lämpötilassa ja diffuusion rajoittamat reaktiot voivat kaikki taivuttaa ln k:n ja suureen 1/T välistä suoraa. Kaareva kuvaaja tarkoittaa, että kahden pisteen Ea on vain keskiarvo kyseisellä lämpötila-alueella, ei kiinteä valli.
• Negatiivinen Ea. Yksittäisellä reaktiovaiheella ei voi olla negatiivista aktivaatioenergiaa. Negatiivinen tulos tarkoittaa yleensä, että reaktio koostuu useista vaiheista ja mekanismiin kätkeytyy lämpötilaherkkä tasapaino.
• Parempi tarkkuus. Kun k mitataan neljässä tai viidessä lämpötilassa ja kaikkien pisteiden kautta sovitetaan suora kahden pisteen käyttämisen sijaan, saadaan luotettavampi kulmakerroin ja tarkempi Ea.

Systemaattisemmassa käsittelyssä siirtymätilateoria esittää nopeusvakion Eyringin yhtälön avulla, $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, joka jakaa energiavallin erillisiksi entalpia- ja entropiatermeiksi.

Tarkan tuloksen saaminen

• Syötä lämpötilat oikein. Celsius- ja Fahrenheit-asteet hyväksytään, mutta taustalla oleva matematiikka toimii vain kelvineillä, joten työkalu muuntaa ne automaattisesti. Tämän vaiheen unohtaminen käsin laskettaessa on yleisin virhelähde.
• Käytä riittävää lämpötilaeroa. Väli 20–50 °C tuottaa vakaan tuloksen. Liian lähellä toisiaan olevat lämpötilat vievät nopeuksien suhdetta k₂/k₁ kohti arvoa 1, jolloin logaritmi voimistaa pieniä virheitä.
• Pidä nopeusvakioiden yksiköt samoina. Kaavassa esiintyy vain suhde k₂/k₁, joten k:n yksiköt supistuvat pois. Syöttökenttien yksikkönä on s⁻¹, joka on ensimmäisen kertaluvun nopeusvakion yksikkö, mutta M⁻¹s⁻¹ tai mikä tahansa muu yksikkö käy, kunhan molemmilla vakioilla on sama yksikkö ja ne on saatu samantyyppisestä mittauksesta.

Usein kysyttyä

Muuttuuko aktivaatioenergia, kun lämpötilaa nostetaan? Ei. Ea pysyy vakiona tietyllä reaktioreitillä. Lämpötilan nostaminen kasvattaa niiden molekyylien osuutta, joilla on riittävästi energiaa vallin ylittämiseen, minkä vuoksi reaktio nopeutuu.

Miten aktivaatioenergia eroaa entalpian muutoksesta (ΔH)? Ea on korkeus, joka reaktion on noustava energiavallin yli. ΔH on lähtöaineiden ja tuotteiden välinen energian nettomuutos. Reaktiolla voi olla korkea valli (suuri Ea) ja sen energia voi silti lopulta olla lähtötilaa pienempi (negatiivinen ΔH, eksoterminen reaktio).

Miksi laskennassa tarvitaan kaksi lämpötilaa yhden sijaan? Yksi nopeusmittaus antaa yhden yhtälön, jossa on kaksi tuntematonta, A ja Ea. Toinen mittaus eri lämpötilassa mahdollistaa A:n supistumisen pois suhteesta, jolloin jäljelle jää yhtälö, josta Ea voidaan ratkaista yksin.

Mitä lämpötilayksiköitä voidaan syöttää? Kelvinit, celsiusasteet tai fahrenheitasteet. Laskin muuntaa valitun yksikön kelvineiksi ennen Ea:n laskemista, koska Arrheniuksen yhtälö edellyttää absoluuttista lämpötilaa.

Kuinka tarkka kahden pisteen menetelmä on? Se antaa tyypillisesti tuloksen 5–10%:n sisällä todellisesta arvosta kohtuullisella lämpötila-alueella, jos reaktio käyttäytyy Arrheniuksen yhtälön mukaisesti. Useamman kuin kahden lämpötilan mittaaminen ja suoran sovittaminen parantavat tarkkuutta entisestään.

Kooditoteutukset

Peruslasku on yksi aritmeettinen rivi. Sama kaava esitetään alla useilla ohjelmointikielillä, ja jokainen palauttaa Ea:n yksikössä 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

Lähteet

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 (10. painos). 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