Αριθμομηχανή Ενέργειας Ενεργοποίησης (Εξίσωση Arrhenius)

Τι κάνει αυτό το εργαλείο

Οι αντιδράσεις επιταχύνονται όταν τις θερμαίνεις. Αυτό το εργαλείο μετατρέπει αυτή την καθημερινή παρατήρηση σε ένα νούμερο: μέτρα πόσο γρήγορα σχηματίζεται μια αντίδραση σε δύο θερμοκρασίες, και θα σου πει την ενέργεια ενεργοποίησης της αντίδρασης (Ea) — το μέγεθος του ενεργειακού «λόφου» που πρέπει να ανεβούν τα μόρια για να αντιδράσουν.

Χρειάζεσαι τέσσερις τιμές: δύο σταθερές ταχύτητας (k₁, k₂) και τις δύο θερμοκρασίες στις οποίες μετρήθηκαν (T₁, T₂). Εισάγετε θερμοκρασίες σε Κέλβιν, Κελσίου ή Φαρενάιτ — το εργαλείο θα κάνει τη μετατροπή για σας. Το αποτέλεσμα εμφανίζεται σε kJ/mol, τη σταθερή μονάδα στη χημεία.

Ενέργεια ενεργοποίησης σε απλά λόγια

Φανταστείτε μια μπάλα που κάθεται σε μια κοιλάδα. Για να την κυλήσετε στην επόμενη κοιλάδα, πρέπει πρώτα να την ώσετε πάνω από μια κορυφογραμμή. Η ενέργεια ενεργοποίησης είναι το ύψος αυτής της κορυφογραμμής. Ακόμη και οι αντιδράσεις που απελευθερώνουν ενέργεια συνολικά χρειάζονται αυτή την αρχική ώθηση για να σπάσουν τους αρχικούς δεσμούς.

Η θέρμανση μιας αντίδρασης δεν χαμηλώνει τη κορυφογραμμή — δίνει σε περισσότερα μόρια αρκετή ενέργεια για να τη διασχίσουν, έτσι περισσότερα από αυτά αντιδρούν κάθε δευτερόλεπτο. Γι' αυτό μια αύξηση 10 °C μπορεί να επιταχύνει αισθητά τα πράγματα.

Ο τύπος

Η εξίσωση Arrhenius συνδέει μια σταθερά ταχύτητας 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). Γεωμετρικά, είναι η κλίση μιας ευθείας γραμμής σε ένα διάγραμμα Arrhenius του 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. Η ενέργεια ενεργοποίησης είναι ιδιότητα του μονοπατιού, όχι της θερμοκρασίας — η θέρμανση αλλάζει πόσα μόρια ξεπερνούν το εμπόδιο, όχι το ύψος του εμποδίου.

Πέρα από τη μέθοδο δύο σημείων

Ο τύπος δύο σημείων υποθέτει τέλειο συμπεριφορά Arrhenius ανάμεσα στα σημεία σας. Όταν δεν είναι γραμμική:

• Καμπύλη μη-Arrhenius εμφανίζεται με αποδιάταξη ενζύμων, κβαντικό σήραγγα σε χαμηλή θερμοκρασία, ή αντιδράσεις περιορισμένες από διάχυση. Ένα καμπύλο διάγραμμα ln k έναντι 1/T είναι το σημάδι — τότε η τιμή δύο σημείων είναι ένα αποτελεσματικό μέσο όρο, όχι ένα σταθερό εμπόδιο.
• Μια αρνητική Ea είναι αδύνατη για ένα στοιχειώδες βήμα. Αν υπολογίσετε ένα, ο μηχανισμός έχει πολλά βήματα με μια θερμοκρασία-ευαίσθητη προ-ισορροπία.
• Για ακρίβεια, η θεωρία της κατάστασης μετάβασης χρησιμοποιεί την εξίσωση Eyring $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, χωρίζοντας το εμπόδιο σε συνεισφορές ενθαλπίας και εντροπίας. Για καλύτερη πειραματική ακρίβεια, μετρήστε το k σε 4–5 θερμοκρασίες και προσαρμόστε την κλίση του διαγράμματος Arrhenius χρησιμοποιώντας γραμμική παλινδρόμηση αντί να χρησιμοποιήσετε μόνο δύο σημεία.

Υλοποιήσεις κώδικα

Ο πυρήνας είναι μια γραμμή αριθμητικής. Ταξινόμηση 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% σε ένα μέτρια εύρος θερμοκρασίας, υποθέτοντας γραμμική συμπεριφορά Arrhenius. Για περισσότερη ακρίβεια, προσαρμόστε 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