Калкулатор активне енергије (Аренијусова једначина)

Шта овај калкулатор ради

Реакције се убрзавају када их загријете. Овај алат претвара то свакодневно посматрање у број: измерите брзину реакције на две температуре, и добијете активациону енергију реакције Ea — величину енергетског "брега" кроз који морају пролазити молекули да би реаговали.

Требају вам четири вредности: две константе брзине (k₁, k₂) и две температуре на којима су измерене (T₁, T₂). Унесите температуре у Келвинима, Целзијусима или Фаренхајтима — калкулатор ће их конвертовати за вас. Резултат долази у kJ/mol, стандардној јединици у хемији.

Активациона енергија у простим речима

Замислите лопту која лежи у долини. Да би је скатили у следећу долину, прво мораме да је гурнемо преко гребена. Активациона енергија је висина тог гребена. Чак и реакције које ослобађају енергију укупно и даље требају овај почетни гурам да прекину почетне везе.

Загревање реакције не мањи гребен — то даје више молекула довољно енергије да га пређу, па више њих реагује у секунди. Због тога пораст од 10 °C може приметно убрзати ствари.

Формула

Аренијусова једначина повезује константу брзине k са температуром:

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

где је A пред-експоненцијални фактор, R = 8.314 J/(mol·K) је гасна константа, и T је апсолутна температура у Келвинима.

Измерите k на две температуре и непозната A се поништава када узмете однос:

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

Ово је тачно оно што калкулатор вреднује (затим дели са 1000 да пријави kJ/mol). Геометријски, то је нагиб праве линије на Аренијусовом графику од ln k наспрам 1/T: нагиб је једнак −Ea/R. Алат црта тај график од ваше две тачке.

Прорађен пример

Са k₁ = 0.0025 s⁻¹ на 300 K и k₂ = 0.035 s⁻¹ на 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}$

Добијање поузданог одговора

• Користите Келвине унутра. Целзијус или Фаренхајт су у реду за куцање — алат ће их претворити — али математика функционише само на апсолутној температури. Заборављање на ово је најчешћа грешка.
• Раширите своје температуре. Пропуст од 20–50 °C даје стабилне резултате. Тачке које су превише близу (однос k₂/k₁ близу 1) раздувавају грешку, јер узимате логаритам броја близу 1.
• Одржавајте јединице конзистентним. Пошто користите однос k₂/k₁, јединице за k се поништавају — s⁻¹, M⁻¹s⁻¹, било шта — све док обе константе не користе исту јединицу и не долазе из исте кинетичке анализе.

Читање резултата

Катализатори и ензими раде достављајући нижи путању барајера. На пример, каталаза опада барајер декомпозиције водиковог пероксида са ~75 kJ/mol на отприлике 8 kJ/mol. Активациона енергија је svojство путање, не температуре — загревање мења колико молекула прелази барајер, не висину барајера.

Изван метода две тачке

Формула две тачке претпоставља савршено Аренијусово понашање између ваших тачака. Када то није линеарно:

• Не-Аренијусова кривина се појављује са денатурирањем ензима, квантним тунелирањем на ниској температури или дифузијом ограниченим реакцијама. Крива ln k наспрам 1/T граф је знак — затим је вредност две тачке ефективна просечна, не фиксни барајер.
• Негативна Ea је немогућа за елементарни корак. Ако је израчунате, механизам има више корака са температурно осетљивом пред-равнотежом.
• За строгост, теорија преходног стања користи Eyringову једначину $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, раздвајајући барајер на енталпију и ентропију доприносе. За најбољу експерименталну тачност, измерите k на 4–5 температура и уклопите нагиб Аренијусовог графика линеарном регресијом уместо коришћења само две тачке.

Имплементације кода

Језгро је једна линија аритметике. Наредбени 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

Честа питања

Да ли температура мења активациону енергију? Не. Ea је фиксна за дату путању. Подизање температуре повећава део молекула који прелазе барајер, због чега брзина расте.

Како се Ea разликује од ΔH? Ea је висина брега кроз који морате пролазити да бисте реаговали; ΔH је нето промена висине од реагенса до производа. Реакција може имати висок барајер (велик Ea) а да заврши нижа него што је почела (негативан ΔH, егзотермно).

Зашто требам две температуре? Једна мерења даје једну једначину са две непознате (Ea и A). Друга мерења омогућава A да се поништи па можете решити за Ea директно.

Колико је тачна метода две тачке? Типично унутар 5–10% преко умереног температурног опсега, под претпоставком линеарног Аренијусовог понашања. За већу тачност, уклопите ln k наспрам 1/T преко 4–5 температура.

Литература

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