Máy tính Năng lượng Kích hoạt (Phương trình Arrhenius)

Máy tính này làm gì

Phản ứng tăng tốc độ khi bạn làm nóng chúng. Công cụ này biến quan sát hàng ngày đó thành một con số: đo tốc độ phản ứng ở hai nhiệt độ, và nó sẽ cho bạn biết năng lượng kích hoạt (Ea) của phản ứng — kích thước của "ngọn đồi" năng lượng mà các phân tử phải trèo để phản ứng.

Bạn cần bốn giá trị: hai hằng số tốc độ (k₁, k₂) và hai nhiệt độ được đo (T₁, T₂). Nhập nhiệt độ bằng Kelvin, Celsius hoặc Fahrenheit — máy tính sẽ chuyển đổi cho bạn. Kết quả tính bằng kJ/mol, đơn vị chuẩn trong hóa học.

Năng lượng kích hoạt theo cách đơn giản

Hãy tưởng tượng một quả bóng nằm trong một thung lũng. Để lăn nó vào thung lũng tiếp theo, trước tiên bạn phải đẩy nó lên một gọc. Năng lượng kích hoạt là chiều cao của gọc đó. Ngay cả những phản ứng giải phóng năng lượng cũng cần lực đẩy ban đầu này để phá vỡ các liên kết ban đầu.

Làm nóng một phản ứng không hạ thấp gọc — nó cung cấp cho nhiều phân tử hơn đủ năng lượng để vượt qua, vì vậy nhiều phân tử hơn phản ứng mỗi giây. Đó là lý do tại sao sự tăng 10 °C có thể tăng tốc độ đáng chú ý.

Công thức

Phương trình Arrhenius liên kết hằng số tốc độ k với nhiệt độ:

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

trong đó A là thừa số trước hàm số mũ, R = 8.314 J/(mol·K) là hằng số khí, và T là nhiệt độ tuyệt đối tính bằng Kelvin.

Đo k ở hai nhiệt độ và A chưa biết sẽ triệt tiêu khi bạn lấy tỷ lệ:

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

Đây chính xác là những gì máy tính đánh giá (sau đó chia cho 1000 để báo cáo kJ/mol). Về mặt hình học, nó là độ dốc của một đường thẳng trên biểu đồ Arrhenius của ln k so với 1/T: độ dốc bằng −Ea/R. Công cụ vẽ biểu đồ đó từ hai điểm của bạn.

Ví dụ đã giải

Với k₁ = 0.0025 s⁻¹ ở 300 K và 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}$

Nhận câu trả lời đáng tin cậy

• Sử dụng Kelvin nội bộ. Celsius hoặc Fahrenheit có thể nhập — công cụ chuyển đổi — nhưng toán học chỉ hoạt động với nhiệt độ tuyệt đối. Quên điều này là lỗi phổ biến nhất.
• Trải rộng nhiệt độ của bạn. Khoảng cách 20–50 °C cho kết quả ổn định. Các điểm quá gần nhau (tỷ lệ k₂/k₁ gần 1) làm cho sai số tăng vọt, vì bạn đang lấy logarit của một số gần 1.
• Giữ đơn vị nhất quán. Vì bạn sử dụng tỷ lệ k₂/k₁, các đơn vị của k triệt tiêu — s⁻¹, M⁻¹s⁻¹, bất kỳ điều gì — miễn là cả hai hằng số sử dụng cùng đơn vị và đến từ cùng phân tích động học.

Đọc kết quả

Các chất xúc tác và enzyme hoạt động bằng cách cung cấp một con đường rào cản thấp hơn. Catalase, chẳng hạn, làm giảm rào cản phân hủy peroxide hydro từ ~75 kJ/mol xuống khoảng 8 kJ/mol. Năng lượng kích hoạt là một thuộc tính của con đường, không phải nhiệt độ — làm nóng thay đổi bao nhiêu phân tử vượt qua rào cản, không phải chiều cao rào cản.

Vượt ra ngoài phương pháp hai điểm

Công thức hai điểm giả định hành vi Arrhenius hoàn hảo giữa các điểm của bạn. Khi nó không tuyến tính:

• Độ cong không phải Arrhenius xuất hiện với biến tính enzyme, đường hầm lượng tử ở nhiệt độ thấp, hoặc phản ứng giới hạn khuếch tán. Một biểu đồ ln k so với 1/T cong là dấu hiệu — khi đó giá trị hai điểm là trung bình hiệu quả, không phải rào cản cố định.
• Một Ea âm là không thể có đối với một bước sơ cấp. Nếu bạn tính toán một, cơ chế có nhiều bước với cân bằng trước nhạy cảm với nhiệt độ.
• Để chặt chẽ, lý thuyết trạng thái chuyển tiếp sử dụng phương trình Eyring $k = \frac{k_B T}{h}\,e^{-\Delta G^{\ddagger}/RT}$, chia rào cản thành đóng góp enthalpic và entropy. Để có độ chính xác thí nghiệm tốt nhất, đo k ở 4–5 nhiệt độ và điều chỉnh độ dốc của biểu đồ Arrhenius bằng hồi quy tuyến tính thay vì chỉ sử dụng hai điểm.

Triển khai mã

Cốt lõi là một dòng số học. Sắp xếp 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

Câu hỏi thường gặp

Có phải nhiệt độ thay đổi năng lượng kích hoạt không? Không. Ea được cố định cho một con đường nhất định. Tăng nhiệt độ làm tăng phần của các phân tử vượt qua rào cản, đó là lý do tại sao tốc độ tăng lên.

Ea khác với ΔH như thế nào? Ea là chiều cao của ngọn đồi bạn phải leo lên để phản ứng; ΔH là sự thay đổi ròng về độ cao từ chất phản ứng đến sản phẩm. Một phản ứng có thể có rào cản cao (Ea lớn) nhưng kết thúc thấp hơn khi nó bắt đầu (ΔH âm, phát tán).

Tại sao tôi cần hai nhiệt độ? Một phép đo cung cấp một phương trình với hai ẩn số (Ea và A). Một phép đo thứ hai cho phép A triệt tiêu để bạn có thể giải Ea trực tiếp.

Phương pháp hai điểm chính xác đến mức nào? Thông thường trong vòng 5–10% trên phạm vi nhiệt độ vừa phải, giả định hành vi Arrhenius tuyến tính. Để có độ chính xác cao hơn, điều chỉnh ln k so với 1/T trên 4–5 nhiệt độ.

Tài liệu tham khảo

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