Kalkulator aktivacijske energije (Arrheniusova jednadžba)

Što ovaj kalkulator radi

Reakcije se ubrzavaju kada ih zagrijete. Ovaj alat pretvara tu svakodnevnu opažanja u broj: izmjerite kako brzo se reakcija odvija na dvjema temperaturama, i dobit ćete aktivacijsku energiju (Ea) te reakcije — veličinu energetske "brijege" koju molekule moraju prevaliti da bi došlo do reakcije.

Trebate četiri vrijednosti: dvije konstante brzine (k₁, k₂) i dvije temperature na kojima su izmjerene (T₁, T₂). Temperaturu unesite u Kelvinima, Celzijusima ili Fahrenheitima — kalkulator će je konvertirati za vas. Rezultat dolazi u kJ/mol, standardnoj jedinici u kemiji.

Aktivacijska energija jednostavnim jezikom

Zamislite loptu koja sjedna u dolini. Kako biste je progurlali u sljedeću dolinu, prvo je trebate gurnuti gore preko brijega. Aktivacijska energija je visina tog brijega. Čak i reakcije koje oslobađaju energiju trebaju ovaj početni udarac da bi prekinule početne veze.

Zagrijavanje reakcije ne snižava brijem — daje više molekula dovoljno energije da je pređu, pa više njih reagira svake sekunde. Zato porast od 10 °C može primjetno ubrzati stvari.

Formula

Arrheniusova jednadžba povezuje konstantu brzine k s temperaturom:

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

gdje je A faktor prije eksponencijala, R = 8.314 J/(mol·K) je plinska konstanta, a T je apsolutna temperatura u Kelvinima.

Izmjerite k na dvjema temperaturama i nepoznati A se poništi kada uzmete omjer:

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

To je točno što kalkulator izračunava (zatim dijeli s 1000 da bi dao rezultat u kJ/mol). Geometrijski, to je nagib ravne linije na Arrheniusovom grafikonu gdje je ln k na ordinati a 1/T na apscisi: nagib je jednak −Ea/R. Alat crta taj grafikon od vaših dvije točke.

Obrađeni primjer

S k₁ = 0.0025 s⁻¹ na 300 K i k₂ = 0.035 s⁻¹ na 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}$

Kako dobit pouzdane rezultate

• Koristite Kelvin interno. Celsius ili Fahrenheit su u redu za upis — alat će konvertirati — ali matematika funkcionira samo s apsolutnom temperaturom. Zaboraviti to je najčešća greška.
• Rasporedite temperature. Razlika od 20–50 °C daje stabilne rezultate. Točke koje su premalo razmaknute (omjer k₂/k₁ blizu 1) preuveličavaju grešku jer računate logaritam broja blizu 1.
• Održavajte jedinice konzistentnim. Budući da koristite omjer k₂/k₁, jedinice k se poništavaju — s⁻¹, M⁻¹s⁻¹, bilo što — dokle god obje konstante koriste istu jedinicu i dolaze iz iste kinematske analize.

Čitanje rezultata

Katalizatori i enzimi funkcioniraju tako što nude put s nižom barijerom. Katalaza, na primjer, snižava barijeru za razgradnju vodikovog peroksida s ~75 kJ/mol na otprilike 8 kJ/mol. Aktivacijska energija je svojstvo puta, ne temperature — zagrijavanje mijenja koliko molekula može preći barijeru, ne samu visinu brijega.

Izvan metode s dvjema točkama

Formula s dvjema točkama pretpostavlja savršeno Arrheniusovo ponašanje između vaših točaka. Kada to nije linearno:

• Arrheniusova zakrivljenost pojavljuje se s denaturacijom enzima, kvantnom tuneliranjem pri niskoj temperaturi ili difuzijskom ograničenom brzinom. Zakrivljeni ln k prema 1/T grafikon je pokazitelj — tada je vrijednost s dvjema točkama efektivni prosjek, ne fiksna barijera.
• Negativna Ea je nemoguća za elementarni korak. Ako je izračunate, mehanizam ima više koraka s temperaturno osjetljivom pre-ravnotežom.
• Za strogost, teorija prijelaznog stanja koristi Eyringovu jednadžbu $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, razdvajajući barijeru na entalpiju i entropiju doprinos. Za najbolju eksperimentalnu točnost, izmjerite k pri 4–5 temperatura i prilagodite nagib Arrheniusova grafikona linearnom regresijom umjesto korištenja samo dvjema točaka.

Implementacije u kodu

Jezgra je jedan red aritmetike. Poredane 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

FAQ

Mijenja li temperatura aktivacijsku energiju? Ne. Ea je fiksna za određeni put. Podizanje temperature povećava udio molekula koji mogu preći barijeru, što je razlog zašto se brzina povećava.

Kako se Ea razlikuje od ΔH? Ea je visina brijega koji trebate penjati da biste reagirali; ΔH je neto promjena nadmorske visine od reaktanata do produkata. Reakcija može imati visoku barijeru (veliki Ea) a završiti niže nego što je počela (negativna ΔH, egzotermna).

Zašto trebam dvije temperature? Jedno mjerenje daje jednu jednadžbu s dvije nepoznanice (Ea i A). Drugo mjerenje omogućava A da se poništi kako biste mogli direktno riješiti Ea.

Koliko je točna metoda s dvjema točkama? Obično unutar 5–10% preko umjerenog temperaturnog raspona, pretpostavljajući linearno Arrheniusovo ponašanje. Za više preciznosti, prilagodite ln k prema 1/T preko 4–5 temperatura.

Reference

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