Kalkulátor aktivační energie (Arrheniusova rovnice)

Co je kalkulačka aktivační energie?

Kalkulačka aktivační energie určuje aktivační energii (Ea) chemické reakce ze dvou rychlostních konstant naměřených při dvou různých teplotách. Používá Arrheniovu rovnici, standardní vztah mezi rychlostí reakce a teplotou. Uživatel zadá rychlostní konstantu a teplotu pro každé ze dvou měření a nástroj vypočítá Ea v kilojoulech na mol (kJ/mol); vedle této hodnoty zobrazí také kcal/mol a elektronvolty na molekulu.

Co je aktivační energie?

Každá chemická reakce musí překonat energetickou bariéru, než se reaktanty mohou změnit na produkty. Molekuly potřebují minimální množství energie k rozštěpení starých vazeb a zahájení tvorby nových. Tato minimální energie se nazývá aktivační energie a značí se Ea.

Užitečnou představu poskytuje kulička spočívající v údolí. Aby se dostala do dalšího údolí, musí nejprve vystoupat přes hřeben. Výška hřebene představuje aktivační energii. I reakce, při níž se celkově energie uvolňuje, potřebuje k zahájení počáteční impuls.

Zahřívání reakce hřeben nesnižuje. Dodává většímu počtu molekul dostatek energie k jeho překonání, takže každou sekundu zreaguje více molekul. Proto zvýšení teploty i o 10 °C často znatelně urychlí reakci.

Arrheniova rovnice

Arrheniova rovnice popisuje, jak závisí rychlostní konstanta reakce k na teplotě:

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

Zde A označuje předexponenciální faktor (přibližně četnost srážek molekul ve vhodné orientaci), R molární plynovou konstantu a T absolutní teplotu v kelvinech. Od 2019 revize SI má R přesnou hodnotu 8,31446261815324 J/(mol·K). Často se uvádí zaokrouhlená hodnota 8,314. Kalkulačka používá přesnou hodnotu.

Jediné měření k nedokáže oddělit A od Ea, protože jedna rovnice obsahuje dvě neznámé. Měření k při dvou teplotách tento problém řeší, protože při vydělení obou rovnic se A vykrátí.

Vzorec pro aktivační energii

Vydělením Arrheniovy rovnice pro T₂ stejnou rovnicí pro T₁ a následným zlogaritmováním obou stran přirozeným logaritmem získáme:

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

Po vyjádření Ea:

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

Toto je rovnice, kterou kalkulačka vyhodnocuje. Teploty lze zadat v kelvinech, stupních Celsia nebo stupních Fahrenheita; nástroj je nejprve převede na kelviny, protože vzorec platí pouze pro absolutní teplotu. Výsledek se dělí 1000, aby se uvedl v kJ/mol, a stejné číslo se zobrazí také v kcal/mol a elektronvoltech na molekulu.

Pokud se ln k vynese proti 1/T, leží obě měření na přímce, jejíž směrnice se rovná −Ea/R. Tento graf se nazývá Arrheniův graf a kalkulačka jej vykreslí ze zadaných bodů.

Příklad výpočtu

Uvažujme reakci s rychlostní konstantou k₁ = 0,0025 s⁻¹ při T₁ = 300 K a k₂ = 0,035 s⁻¹ při T₂ = 350 K.

1. Poměr rychlostí: 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

Stejný výsledek odpovídá přibližně 11,01 kcal/mol neboli zhruba 0,478 eV na každou reagující molekulu.

Jak výsledek interpretovat

Katalyzátory a enzymy urychlují reakce tím, že nabízejí cestu s nižší bariérou. Například kataláza snižuje bariéru rozkladu peroxidu vodíku přibližně z 75 kJ/mol na zhruba 8 kJ/mol. Aktivační energie náleží reakční dráze, nikoli teplotě. Zahřívání mění počet molekul, které bariéru překonají, nikoli samotnou výšku bariéry.

Omezení metody dvou bodů

Vzorec pro dva body předpokládá, že reakce se mezi dvěma naměřenými body přesně řídí Arrheniovou rovnicí. Tento předpoklad nemusí platit:

• Zakřivené grafy. Denaturace enzymů, kvantové tunelování při nízké teplotě a reakce omezené difuzí mohou zakřivit přímku ln k v závislosti na 1/T. Zakřivený graf znamená, že Ea vypočtená ze dvou bodů je pouze průměrem pro daný rozsah teplot, nikoli pevnou bariérou.
• Záporná Ea. Jediný krok reakce nemůže mít zápornou aktivační energii. Záporný výsledek obvykle znamená, že reakce probíhá v několika krocích a v mechanismu se skrývá rovnováha citlivá na teplotu.
• Vyšší přesnost. Měření k při čtyřech nebo pěti teplotách a proložení přímky všemi body namísto použití pouhých dvou bodů poskytne spolehlivější směrnici a přesnější Ea.

Pro rigoróznější popis vyjadřuje teorie přechodového stavu rychlostní konstantu Eyringovou rovnicí, $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, která rozděluje energetickou bariéru na samostatné entalpické a entropické příspěvky.

Jak získat přesný výsledek

• Zadávejte teploty správně. Stupně Celsia a Fahrenheita jsou přijímány, ale výpočty v pozadí fungují pouze s kelviny, takže je nástroj automaticky převádí. Opomenutí tohoto kroku při ručním výpočtu je nejčastějším zdrojem chyb.
• Zvolte dostatečný rozdíl mezi teplotami. Rozdíl 20–50 °C poskytuje stabilní výsledek. Příliš blízké teploty posouvají poměr rychlostí k₂/k₁ směrem k 1 a malé chyby se logaritmem zesilují.
• Používejte shodné jednotky rychlostní konstanty. Ve vzorci se objevuje pouze poměr k₂/k₁, takže jednotky k se vykrátí. Vstupní pole jsou označena jednotkou s⁻¹, což je jednotka rychlostní konstanty reakce prvního řádu, ale lze použít M⁻¹s⁻¹ nebo jakoukoli jinou jednotku, pokud obě konstanty používají stejnou jednotku a pocházejí ze stejného typu měření.

Často kladené otázky

Mění zvýšení teploty aktivační energii? Ne. Ea zůstává pro danou reakční dráhu konstantní. Zvýšení teploty zvětšuje podíl molekul s dostatkem energie k překonání bariéry, a proto se reakce urychlí.

Jak se aktivační energie liší od změny entalpie (ΔH)? Ea je výška bariéry, kterou musí reakce překonat. ΔH je celková změna energie mezi reaktanty a produkty. Reakce může mít vysokou bariéru (velké Ea) a přesto skončit na nižší energetické hladině, než na jaké začala (záporné ΔH, tedy exotermická reakce).

Proč výpočet potřebuje dvě teploty místo jedné? Jedno měření rychlosti poskytuje jednu rovnici se dvěma neznámými, A a Ea. Druhé měření při jiné teplotě umožní, aby se A v poměru vykrátilo, takže zůstane rovnice, z níž lze určit samotné Ea.

Jaké jednotky teploty lze zadat? Kelviny, stupně Celsia nebo stupně Fahrenheita. Kalkulačka před výpočtem Ea převede zvolenou jednotku na kelviny, protože Arrheniova rovnice vyžaduje absolutní teplotu.

Jak přesná je metoda dvou bodů? V běžném rozsahu teplot obvykle poskytne hodnotu s odchylkou do 5–10% od skutečné hodnoty, pokud se reakce řídí Arrheniovým vztahem. Měření při více než dvou teplotách a proložení přímky přes tato měření přesnost dále zvyšuje.

Implementace v kódu

Základní výpočet je jediný aritmetický řádek. Stejný vzorec je níže uveden v několika programovacích jazycích; každý z nich vrací Ea v 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

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 (10. vydání). Oxford University Press.
4. IUPAC. Compendium of Chemical Terminology (Zlatá kniha), "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