Ця симуляція моделює послідовну реакцію A → B → C, відстежуючи, як концентрації реагенту A, проміжного продукту B та продукту C змінюються протягом 120-секундного інтервалу. Вона інтегрує зв'язані рівняння швидкості dA/dt = −k₁A + k₋₁B, dB/dt = k₁A − k₋₁B − k₂B та dC/dt = k₂B за допомогою простої явної схеми Ейлера. Константи швидкості випливають із рівняння Арреніуса, k = A·e−Ea/RT, де R = 8.314 Дж/моль·K.
Повзунок температури (273–700 K) і два повзунки енергії активації (Ea₁ та Ea₂, 10–120 кДж/моль) задають константи швидкості k₁ та k₂. Перемикач каталізатора знижує обидві енергії активації на 30%, а перемикач оборотності додає зворотну реакцію B → A. Подібна послідовна кінетика керує всім: від промислового синтезу та метаболізму ліків до хімії атмосфери, де швидкоплинний проміжний продукт накопичується, а потім розпадається.
Що насправді показує ця симуляція?
Вона будує графік концентрацій трьох речовин у послідовній реакції: A перетворюється на B, який потім перетворюється на C. Реагент A (червоний) спадає, проміжний продукт B (жовтий) зростає до піку, а потім спадає, а продукт C (зелений) накопичується. Ви спостерігаєте всю криву, а живий часовий маркер рухається вздовж неї в міру перебігу реакції.
Яке рівняння Арреніуса тут використовується?
Рівняння Арреніуса k = A помножене на e в степені мінус Ea поділене на RT пов'язує константу швидкості з температурою та енергією активації. Тут A — фіксований передекспоненційний множник 1 помножене на 10 у восьмому степені на секунду, R — газова стала 8.314 джоуля на моль на кельвін, T — температура в кельвінах, а Ea — енергія активації. Симуляція обчислює його окремо для кожного етапу.
Що роблять повзунки температури та енергії активації?
Температура (273 до 700 кельвінів) і дві енергії активації Ea1 та Ea2 (10 до 120 кілоджоулів на моль) безпосередньо подаються в рівняння Арреніуса, щоб задати константи швидкості k1 та k2. Підвищення температури або зниження енергії активації збільшує константи швидкості, прискорюючи реакцію. Панель показує отримані k1, k2 та період напіврозпаду A.
У послідовній реакції B утворюється з A, але також витрачається на утворення C. На початку A багато, тому B накопичується швидше, ніж розпадається, даючи зростаючу криву. Щойно A виснажується, утворення сповільнюється, а споживання триває, тож B проходить через максимум, а потім спадає. Висота й час цього піку залежать від співвідношення k1 до k2.
Прапорець каталізатора множить обидві енергії активації на 0.7, знижуючи їх на 30%. Оскільки енергія активації стоїть в експоненті рівняння Арреніуса, навіть таке помірне зниження може суттєво підвищити константи швидкості та різко прискорити реакцію. Важливо, що каталізатор змінює кінетику, а не термодинаміку, тож він не змінює загальної зміни вільної енергії чи кінцевих продуктів.
Увімкнення перемикача оборотності додає зворотну реакцію B перетворюється на A з власною константою швидкості. Це означає, що A більше не витрачається повністю, бо частина B постійно повертається до A. Система прямує до динамічної рівноваги, а не до повного перетворення, тож криві A і B виходять на ненульові значення, замість того щоб A спадало до нуля.
Вона інтегрує три зв'язані звичайні диференціальні рівняння прямим (явним) методом Ейлера, використовуючи малий фіксований крок за часом 0.02 секунди для попереднього обчислення повної кривої та 20 підкроків на кадр анімації для живого запуску. Концентрації обмежуються знизу нулем, щоб уникнути малих від'ємних значень через округлення. Метод Ейлера простий і швидкий, але має лише перший порядок точності.
Показаний період напіврозпаду — це час, за який A спадає до половини свого початкового значення, обчислений як натуральний логарифм 2 поділений на k1, приблизно 0.693 поділене на k1. Це стандартний період напіврозпаду першого порядку, і він точний лише для простої втрати A на першому етапі; зворотна реакція та подальші етапи означають, що реальний розпад може дещо відрізнятися, коли вони активні.
Схема реакції та температурна залежність за Арреніусом фізично коректні, але передекспоненційний множник і концентрації — це ілюстративні значення за замовчуванням, а не дані для конкретної реальної реакції. Спрощене інтегрування за Ейлером і фіксований множник каталізатора — це навчальні ідеалізації. Вона достовірно відтворює якісну поведінку послідовної кінетики першого порядку, а не калібрований промисловий процес.
Послідовна кінетика A до B до C трапляється повсюди в хімії та біології: фармацевтичні препарати метаболізуються через проміжні сполуки, полімери утворюються поетапно, а атмосферні забруднювачі розкладаються через реакційноздатні проміжні продукти. Розуміння того, коли проміжний продукт сягає піку, допомагає хімікам розрахувати час реакції, щоб максимізувати вихід B, перш ніж він перетвориться на C, — поширена мета в синтезі.