Aktivaatioenergia-laskin (Arrhenius-yhtälö)

Mitä tämä laskin tekee

Reaktiot nopeutuvat, kun niitä kuumennetaan. Tämä työkalu muuntaa tämän päivittäisen havainnon numeroksi: mittaa, kuinka nopeasti reaktio etenee kahdessa eri lämpötilassa, ja se kertoo reaktion aktivaatioenergialle (Ea) — energian "kukkulan" koko, jonka molekyylit joutuvat ylittämään reagoidakseen.

Tarvitset neljä arvoa: kaksi nopeusvakiota (k₁, k₂) ja kaksi lämpötilaa, joissa ne mitattiin (T₁, T₂). Syötä lämpötilat Kelvineinä, Celsiuksina tai Fahrenheiteinä — laskin muuntaa puolestasi. Tulos saadaan kJ/mol, kemian standardi yksikkönä.

Aktivaatioenergia yksinkertaisesti

Kuvittele pallo laaksossa. Vierittääksesi sen seuraavaan laaksoon, sinun on ensin työnnettävä se ylös harjanteen yli. Aktivaatioenergia on tuon harjanteen korkeus. Jopa reaktiot, jotka vapauttavat energiaa kokonaisuudessaan, tarvitsevat tämän alkuperäisen työnnön lähtöisiltään katkaisemiseksi.

Reaktion kuumennus ei madalla harjannetta — se antaa useammille molekyyleille tarpeeksi energiaa sen ylittämiseksi, joten useammat niistä reagoivat joka sekunti. Siksi 10 °C nousu voi huomattavasti nopeuttaa asioita.

Kaava

Arrhenius-yhtälö yhdistää nopeusvakion k lämpötilaan:

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

missä A on pre-eksponentiaalinen tekijä, R = 8.314 J/(mol·K) on kaasuvakio, ja T on absoluuttinen lämpötila Kelvineinä.

Mittaa k kahdessa lämpötilassa ja tuntematon A poistuu, kun otat suhteen:

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

Tämä on juuri se, mitä laskin laskee (sitten jakaa 1000:lla raportoidakseen kJ/mol). Geometrisesti se on suoran viivan kaltevuus Arrhenius-kuvaajalla ln k versus 1/T: kaltevuus on yhtä kuin −Ea/R. Työkalu piirtää tuon kuvaajan kahdesta pisteestäsi.

Työesimerkki

Jos k₁ = 0.0025 s⁻¹ lämpötilassa 300 K ja k₂ = 0.035 s⁻¹ lämpötilassa 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}$

Luotettavan vastauksen saaminen

• Käytä Kelviniä sisäisesti. Celsius tai Fahrenheit sopivat syötettäväksi — laskin muuntaa — mutta matematiikka toimii vain absoluuttisella lämpötilalla. Tämän unohtaminen on yleisin virhe.
• Jaa lämpötilasi. 20–50 °C ero antaa vakaita tuloksia. Pisteet liian lähellä toisistaan (k₂/k₁-suhde lähellä 1:tä) paisuttavat virhettä, koska otat logaritmin luvusta lähellä 1:tä.
• Pidä yksiköt yhdenmukaisina. Koska käytät suhdetta k₂/k₁, k:n yksiköt poistuvat — s⁻¹, M⁻¹s⁻¹, mikä tahansa — kunhan molemmat vakiot käyttävät samaa yksikköä ja tulevat samasta kineettisesta analyysista.

Tuloksen lukeminen

Katalyytit ja entsyymit toimivat tarjoamalla matalamman esteen polun. Katalaaasi esimerkiksi pudottaa vedynperoksidi-hajoamisen esteen ~75 kJ/mol:sta noin 8 kJ/mol:iin. Aktivaatioenergia on polun ominaisuus, ei lämpötila — kuumennus muuttaa kuinka monta molekyyliä ylittää esteen, ei esteen korkeutta.

Kahden pisteen menetelmän yli

Kahden pisteen kaava olettaa täydellisen Arrhenius-käyttäytymisen pisteidesi välillä. Kun se ei ole lineaarinen:

• Ei-Arrhenius-kaarevuus ilmestyy entsyymin denaturoimisen, kvanttikaivannon tai diffuusio-rajoitettujen reaktioiden kanssa. Kaareva ln k vs 1/T kuvaaja on merkki — silloin kahden pisteen arvo on tehokas keskiarvo, ei kiinteä este.
• Negatiivinen Ea on mahdoton alkeelliselle askeleelle. Jos lasket sellaisen, mekanismissa on useita askelia lämpötila-herkkään esitasapainoon.
• Tarkuutta varten, siirtymätilateoria käyttää Eyring-yhtälöä $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, jakamalla esteen entalpian ja entropian osuuksiksi. Parasta kokeellista tarkkuutta varten mittaa k 4–5 lämpötilassa ja sovita Arrhenius-kuvaajan kaltevuus lineaarisella regressiolla kahden pisteen käyttämisen sijasta.

Kooditekniikat

Ydin on yksi aritmeettinen rivi. Järjestetty 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

UKK

Muuttaako lämpötila aktivaatioenergiaa? Ei. Ea on kiinteä tietylle polulle. Lämpötilan nostaminen lisää fraktioita molekyyleistä, jotka ylittävät esteen, minkä vuoksi nopeus nousee.

Miten Ea eroaa ΔH:sta? Ea on kukkulan korkeus, joka sinun on kiivettävä reagoidaksesi; ΔH on netto korkeus muutos reaktanteista tuotteisiin. Reaktiolla voi olla korkea este (suuri Ea) mutta lopeta alemmas kuin se alkoi (negatiivinen ΔH, eksoterminen).

Miksi tarvitsen kaksi lämpötilaa? Yksi mittaus antaa yhden yhtälön kahdella tuntemattomalla (Ea ja A). Toinen mittaus antaa A:lle peruuttaa niin voit ratkaista Ea:n suoraan.

Kuinka tarkka on kahden pisteen menetelmä? Yleensä 5–10 % sisällä kohtuullisen lämpötilaalueen, olettaen lineaarista Arrhenius-käyttäytymistä. Suuremmalle tarkkuudelle sovita ln k vs 1/T neljässä-viidessä lämpötilassa.

Viitteet

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