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

Шта је калкулатор енергије активације?

Калкулатор енергије активације одређује енергију активације (Ea) хемијске реакције на основу две константе брзине измерене на две различите температуре. Он примењује Аренијусову једначину, стандардну формулу која повезује брзину реакције с температуром. Корисник уноси константу брзине и температуру за свако од два мерења, а алат израчунава Ea у килоџулима по молу (kJ/mol), уз приказ вредности у kcal/mol и електронволтима по молекулу.

Шта је енергија активације?

Свака хемијска реакција мора да савлада енергетску баријеру пре него што реактанти могу да пређу у производе. Молекулима је потребна минимална количина енергије да би раскинули старе везе и почели да образују нове. Та минимална енергија назива се енергија активације и означава се са Ea.

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

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

Аренијусова једначина

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

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

Овде A означава предекспоненцијални фактор (приближно, учесталост судара молекула у одговарајућој оријентацији), R моларна гасна константа, а T апсолутна температура у келвинима. Од 2019. ревизије SI-ја, R има тачну вредност: 8,31446261815324 J/(mol·K). Често се наводи заокружена вредност 8,314. Калкулатор користи тачну вредност.

Једно мерење k не може да раздвоји A од Ea, јер једна једначина има две непознате. Мерење k на две температуре решава овај проблем, јер се A поништава када се једна једначина подели другом.

Формула за енергију активације

Дељењем Аренијусове једначине за T₂ истом једначином за T₁, а затим узимањем природног логаритма обе стране, добија се:

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

Издвајање Ea:

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

То је једначина коју калкулатор израчунава. Температуре се могу унети у келвинима, Целзијусовим или Фаренхајтовим степенима; алат их најпре претвара у келвине, јер формула важи само за апсолутну температуру. Резултат се дели са 1000 да би се изразио у kJ/mol, а исти број се приказује и у kcal/mol и електронволтима по молекулу.

Када се две мерене вредности нанесу као ln k у односу на 1/T, оне леже на правој линији чији је нагиб једнак −Ea/R. Такав график се назива Аренијусовим дијаграмом, а калкулатор га исцртава на основу унетих тачака.

Решени пример

Размотримо реакцију чија је константа брзине k₁ = 0,0025 s⁻¹ при T₁ = 300 K, а k₂ = 0,035 s⁻¹ при T₂ = 350 K.

1. Однос константи брзине: 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

Исти резултат износи приближно 11,01 kcal/mol, односно око 0,478 eV за сваки молекул који реагује.

Како тумачити резултат

Катализатори и ензими убрзавају реакције тако што нуде пут с нижом баријером. Каталаза, на пример, снижава баријеру за разлагање водоник-пероксида са око 75 kJ/mol на приближно 8 kJ/mol. Енергија активације односи се на реакциони пут, а не на температуру. Загревање мења број молекула који савладају баријеру, али не и саму висину баријере.

Ограничења методе са две тачке

Формула са две тачке претпоставља да реакција тачно следи Аренијусову једначину између две измерене тачке. Та претпоставка може да престане да важи:

• Закривљени графици. Денатурација ензима, квантно тунеловање на ниским температурама и реакције ограничене дифузијом могу да закриве линију ln k у односу на 1/T. Закривљени графикон значи да је Ea добијена из две тачке само просечна вредност за тај температурни опсег, а не фиксна баријера.
• Негативна Ea. Појединачни корак реакције не може имати негативну енергију активације. Негативан резултат обично значи да се реакција одвија у више корака, при чему је у механизму скривена равнотежа зависна од температуре.
• Већа прецизност. Мерење k на четири или пет температура и уклапање праве кроз све њих, уместо коришћења само две тачке, даје поузданији нагиб и тачнију вредност Ea.

За строжији приступ, према теорији прелазног стања, константа брзине дата је Ајринговом једначином, $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, која раздваја енергетску баријеру на засебне доприносе енталпије и ентропије.

Како добити тачан резултат

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

Најчешће постављана питања

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

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

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

Које температурне јединице се могу унети? Келвини, Целзијусови или Фаренхајтови степени. Калкулатор претвара изабрану јединицу у келвине пре израчунавања Ea, јер Аренијусова једначина захтева апсолутну температуру.

Колико је тачна метода са две тачке? Она обично даје резултат унутар 5–10% од стварне вредности у умереном температурном опсегу, под условом да реакција следи Аренијусово понашање. Мерење на више од две температуре и уклапање линије додатно побољшавају тачност.

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

Основни прорачун је један ред аритметике. Иста формула је у наставку приказана на више језика, при чему сваки враћа Ea у 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

Литература

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. издање). Oxford University Press.
4. IUPAC. Compendium of Chemical Terminology (Златна књига), "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