Kalkulator energii aktywacji (równanie Arrheniusa)

Co robi ten kalkulator

Reakcje przyspieszają się podgrzewane. Narzędzie to przekształca codzienną obserwację w liczbę: zmierz szybkość reakcji w dwóch temperaturach, a uzyskasz energię aktywacji (Ea) — wielkość energetycznego "wzgórza", które muszą pokonać cząsteczki, aby zareaagować.

Potrzebujesz czterech wartości: dwie stałe szybkości (k₁, k₂) i dwie temperatury, w których je zmierzono (T₁, T₂). Wpisz temperatury w Kelwinach, Celsjuszu lub Fahrenheicie — kalkulator przelicza dla ciebie. Wynik podawany jest w kJ/mol, standardowej jednostce w chemii.

Energia aktywacji w prostych słowach

Wyobraź sobie piłkę leżącą w dolinie. Aby toczyć ją do następnej doliny, musisz ją najpierw pchnąć nad grzbietem. Energia aktywacji to wysokość tego grzbietu. Nawet reakcje, które ogólnie uwalniają energię, potrzebują tego początkowego pchnięcia, aby przerwać wiązania początkowe.

Podgrzewanie reakcji nie obniża grzbietu — daje więcej cząsteczkom wystarczającą energię, aby go pokonać, więc więcej z nich reaguje każdej sekundy. Dlatego wzrost o 10 °C może wyraźnie przyspieszyć rzeczy.

Wzór

Równanie Arrheniusa łączy stałą szybkości k z temperaturą:

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

gdzie A to współczynnik przedwykładniczy, R = 8.314 J/(mol·K) to stała gazowa, a T to temperatura bezwzględna w Kelwinach.

Zmierz k w dwóch temperaturach, a nieznany A anuluje się, gdy weźmiesz stosunek:

$\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 jest dokładnie to, co oblicza kalkulator (a następnie dzieli przez 1000, aby podać w kJ/mol). Geometrycznie jest to nachylenie linii prostej na wykresie Arrheniusa ln k względem 1/T: nachylenie równa się −Ea/R. Narzędzie rysuje ten wykres z twoich dwóch punktów.

Przykład rozwiązany

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

Uzyskanie wiarygodnej odpowiedzi

• Używaj Kelwina wewnętrznie. Celsjusz lub Fahrenheit są w porządku do wpisania — narzędzie konwertuje — ale matematyka działa tylko na temperaturze bezwzględnej. Zapomnienie tego to najczęstszy błąd.
• Rozsunęcie temperatur. Przedział 20–50 °C daje stabilne wyniki. Punkty zbyt blisko siebie (stosunek k₂/k₁ bliski 1) powodują błąd, ponieważ bierzesz log liczby bliskiej 1.
• Zachowaj spójność jednostek. Ponieważ używasz stosunku k₂/k₁, jednostki k się znoszą — s⁻¹, M⁻¹s⁻¹, cokolwiek — tak długo, jak oba stałe używają tej samej jednostki i pochodzą z tej samej analizy kinetycznej.

Czytanie wyniku

Katalizatory i enzymy działają, oferując ścieżkę o niższej bariery. Na przykład katalaza obniża barierę dla rozkładu nadtlenku wodoru z ~75 kJ/mol do około 8 kJ/mol. Energia aktywacji jest właściwością drogi, a nie temperatury — podgrzewanie zmienia, ile cząsteczek przechodzi przez barierę, a nie wysokość bariery.

Poza metodą dwupunktową

Wzór dwupunktowy zakłada doskonałe zachowanie Arrheniusa między twoimi punktami. Gdy nie jest liniowy:

• Krzywiznę non-Arrheniusa pojawia się przy denaturacji enzymu, tunelowaniu kwantowym w niskiej temperaturze lub reakcjach ograniczonych dyfuzją. Zakrzywiony wykres ln k vs 1/T to wskaźnik — wtedy wartość dwupunktowa to efektywna średnia, a nie stała bariera.
• Ujemny Ea jest niemożliwy dla kroku elementarnego. Jeśli go obliczysz, mechanizm ma wiele kroków z czułą temperaturowo pre-równowagą.
• Dla ścisłości, teoria stanu przejścia używa równania Eyringa $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, dzieląc barierę na wkład entalpii i entropii. Dla najlepszej dokładności eksperymentalnej, zmierz k w 4–5 temperaturach i dopasuj nachylenie wykresu Arrheniusa za pomocą regresji liniowej zamiast używać tylko dwóch punktów.

Implementacje kodu

Rdzeń to jedna linia arytmetyki. Uporządkowane 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

Często zadawane pytania

Czy temperatura zmienia energię aktywacji? Nie. Ea jest stała dla danej ścieżki. Podniesienie temperatury zwiększa ułamek cząsteczek, które przechodzą przez barierę, dlatego szybkość rośnie.

Jak Ea różni się od ΔH? Ea to wysokość wzgórza, które musisz wejść, aby zareaagować; ΔH to zmiana netto wysokości od reagentów do produktów. Reakcja może mieć wysoką barierę (duży Ea), ale skończyć się niżej niż się zaczęła (ujemna ΔH, egzotermiczna).

Dlaczego potrzebuję dwóch temperatur? Jeden pomiar daje jedno równanie z dwiema niewiadomymi (Ea i A). Drugi pomiar pozwala A anulować, aby móc rozwiązać Ea bezpośrednio.

Jak dokładna jest metoda dwupunktowa? Zazwyczaj w granicach 5–10% w umiarkowanym zakresie temperatur, zakładając liniowe zachowanie Arrheniusa. Aby uzyskać większą precyzję, dopasuj ln k vs 1/T w 4–5 temperaturach.

Referencje

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