Огляд №39 – Фізична хімія та динаміка реакцій: перехідні стани, теорія Маркуса, спектроскопія та фемтохімія

Фізична хімія займає простір на перетині фізики й хімії: вона застосовує строгий кількісний апарат фізики для розуміння структури, динаміки та перетворень хімічних систем. Ця галузь охоплює все — від рівноважної термодинаміки та електрохімії до молекулярної квантової механіки та надшвидкої лазерної спектроскопії. Цей огляд охоплює шість центральних тем, кожна з математичним підґрунтям, що робить фізичну хімію дійсно передбачуваною, а не лише описовою.

1. Теорія перехідного стану

Класична теорія зіткнень дає оцінки передекспоненційного множника, але не враховує геометрію поверхні потенціальної енергії (ППЕ). Теорія перехідного стану (ТПС, 1935) Генрі Ейрінга розглядає конфігурацію в сідловій точці (сам «перехідний стан», або активований комплекс) як таку, що перебуває у квазітермодинамічній рівновазі з реагентами, а потім обчислює потік через бар'єр за допомогою статистичної механіки.

Закон швидкості Арреніуса:
  k(T) = A ⋅ e^(−E_a/RT)
  A = передекспоненційний (частотний) множник [с−¹ або М−¹с−¹]
  E_a = енергія активації [кДж/моль]

Рівняння Ейрінга (перехідного стану):
  k(T) = (k_B T / h) ⋅ (Κ)† ⋅ exp(−ΔG‡/RT)
  k_B T/h ≈ 6,25 × 10^12 с−¹ при 298 К
  ΔG‡ = ΔH‡ − TΔS‡   (вільна енергія активації)
  κ = коефіцієнт проходження (≤1; враховує повторні переходи)

Графік Ейрінга:
  ln(k/T) = −ΔH‡/R ⋅ (1/T) + ln(k_B/h) + ΔS‡/R
  Нахил = −ΔH‡/R,  перетин дає ΔS‡

Постулат Хеммонда:
  Екзотермічна: ПС подібний до реагенту → ΔE‡ ≈ малий
  Ендотермічна: ПС подібний до продукту → ΔE‡ ≈ ΔH_rxn
  Еванс-Полані: ΔE‡ = αΔH_rxn + β (лінійне співвідношення вільних енергій)

Поправки на квантове тунелювання:
  Поправка Вігнера: κ_tun ≈ 1 + (1/24)(hν‡/k_BT)²
  Поправка Белла (ширша парабола): складніша; суттєва для перенесення протона/Н
  Методи DFT/континуального інтеграла: використовуються для ферментативних реакцій та тунелювання Н
          

2. Теорія перенесення електрона Маркуса

Рудольф Маркус (Нобелівська премія 1992) розробив кількісну теорію реакцій перенесення електрона (ПЕ) у розчині, яка передбачила дивовижну «інвертовану область»: понад певний поріг рушійної сили збільшення термодинамічної рушійної сили зменшує швидкість. Це нелогічне на перший погляд передбачення, підтверджене у 1980-х роках, є центральним для фотосинтезу, сонячних елементів та транспорту заряду в органічній електроніці.

Перенесення електрона зовнішньої сфери:
  D + A → D•+ + A•−   (без розриву зв'язку, реорганізація розчинника)

Константа швидкості Маркуса:
  k_ET = (2π/ℏ) |H_DA|² (1/√(4πλk_BT)) exp(−(ΔG° + λ)²/(4λk_BT))

  H_DA = матричний елемент електронного зв'язку (перекриття донор-акцептор)
  λ   = повна енергія реорганізації = λ_inner + λ_outer
  ΔG° = стандартна зміна вільної енергії реакції ПЕ

Енергії реорганізації:
  λ_inner (внутрішньомолекулярна, зміна довжини зв'язку):
    λ_i = (1/2)(f_D f_A / (f_D+f_A)) ⋅ (Δq_eq)²
  λ_outer (поляризація розчинника):
    λ_o = (e²/2)(1/r_D + 1/r_A − 1/r_DA)(1/ε_op − 1/ε_s)

Три режими (|H_DA| мале → неадіабатичний):
  Нормальний:    |ΔG°| < λ  → бар'єр = (ΔG°+λ)²/4λ > 0
  Безбар'єрний:  ΔG° = −λ → максимальна швидкість (нульовий бар'єр)
  Інвертований:  |ΔG°| > λ  → швидкість ЗНИЖУЄТЬСЯ при подальшому зростанні |ΔG°|
    Підтверджено Міллером, Калькаттера, Клоссом (1984) для внутрішньомолекулярного ПЕ

Фотосинтетичне ПЕ в реакційному центрі фотосистеми II:
  P680* → феофітин → QA → QB  (послідовне ПЕ протягом ~200 пс)
  Інвертована область запобігає зворотному ПЕ на кожному кроці: квантова ефективність ~99%

Напівпровідникові сонячні елементи:
  Комірка Гретцеля (DSSC): збуджений стан барвника вприскує e− у зону провідності TiO&sub2;
  Маркусівське λ спроєктоване для швидкого прямого ПЕ, повільного зворотного ПЕ через ΔG° в інвертованій області
          

3. Хімічна термодинаміка

Рівноважний стан хімічної системи визначається мінімізацією вільної енергії Гіббса G за сталої температури й тиску. Хімічний потенціал μ_i кожного компонента — це ворота між термодинамічними таблицями та константами рівноваги, принципом Ле Шательє й фазовими діаграмами.

Вільна енергія Гіббса:
  G = H − TS = U + pV − TS
  dG = −S dT + V dp + Σ_i μ_i dn_i
  Рівновага при сталих T,p: dG = 0

Хімічний потенціал:
  μ_i = (∂G/∂n_i)_{T,p,n_j}
  Ідеальна суміш: μ_i = μ_i°(T) + RT ln x_i   (мольна частка x_i)
  Реальний розчин: μ_i = μ_i°(T) + RT ln a_i   (активність a_i = γ_i x_i)
  Летючість (реальні гази): μ = μ° + RT ln(f/p°)

Константа рівноваги:
  K(T) = exp(−Δ_r G°/RT) = ∏_i a_i^(ν_i)
  Ізохора Вант-Гоффа: d ln K / dT = Δ_r H°/(RT²)

Рівняння Клапейрона (нахил межі фаз):
  dP/dT = ΔH/(TΔV)
  Тверде тіло-рідина: великий dP/dT (малий ΔV); 1-3 МПа/К для більшості речовин
  Клаузіус-Клапейрон (пара-рідина, наближення ідеального газу):
    d ln P/dT = Δ_vap H / RT²
    P(T) = P_0 exp(−Δ_vap H/R ⋅ (1/T − 1/T_0))

Правило фаз Гіббса:
  F = C − P + 2    (F = ступені свободи, C = компоненти, P = фази)
  Бінарна двофазна (напр. рідина+пара): F=1 (ізотермічна ↔ ізобарична еквівалентність)
  Потрійна точка: F=0 (унікальні T,P для чистої речовини)
          

4. Спектроскопія

Спектроскопія — головний експериментальний інструмент фізичної хімії. Закон Бера-Ламберта кількісно описує УФ-видиме поглинання; обертальна й коливальна спектроскопія (ІЧ та КР) зондує молекулярну структуру; ЯМР розрізняє хімічні оточення та мережі спін-спінових взаємодій. Кожен метод пов'язаний із квантовою механікою через правила відбору та вирази для рівнів енергії.

Закон Бера-Ламберта:
  A = εlc = −log(I/I_0)
  ε = молярний коефіцієнт поглинання [М−¹см−¹],  l = довжина шляху [см]
  Лінійна залежність від концентрації для ідеальних розведених розчинів

Обертальна енергія жорсткого ротора:
  E_J = hcBJ(J+1),  B = h/(8π²Ic)  [см−¹]
  I = μr_e²,  μ = зведена маса
  Правило відбору: ΔJ = ±1 (активний у мікрохвильовому діапазоні → потребує постійного дипольного моменту)
  Обертальний інтервал: Δν = 2B(J+1)

Коливальна енергія гармонічного осцилятора:
  E_v = hν_e(v + 1/2),  v = 0,1,2,...
  ν_e = (1/2πc)√(k/μ)  [см−¹],  k = силова константа
  Правило відбору ІЧ: Δv = ±1, потребує dμ/dq ≠ 0 (зміна дипольного моменту)
  Правило відбору КР: потребує dα/dq ≠ 0 (зміна поляризовності)
  Правило взаємного виключення: активна ІЧ ↔ неактивна КР для центросиметричних молекул

Ангармонічність:
  E_v = hν_e(v+1/2) − hν_eχ_e(v+1/2)²
  D_e = hν_e/(4χ_e)  (енергія дисоціації з потенціалу Морзе)
  Обертони: Δv = ±2,±3 дозволені (слабкі)

Імпульсний Фур'є-ЯМР:
  Ларморівська частота: ω_L = γ B_0  (для протона γ/2π = 267,5 МГц/Тл)
  Хімічний зсув δ (ppm): відображає екранування електронами
  Спін-спінова взаємодія J (Гц): через-зв'язкове перекриття орбіталей
  ССІ (спад вільної індукції) → перетворення Фур'є → спектр хімічних зсувів
  2D COSY (зв'язані спіни), NOESY (просторова близькість <5Å): визначення структури
          

5. Квантова хімія

Квантова хімія розв'язує електронне рівняння Шредінгера, щоб отримати молекулярні енергії, структури та властивості. Наближення Борна-Оппенгеймера розділяє електронний і ядерний рух; теорія молекулярних орбіталей (МО) будує хвильові функції з атомних базисних функцій; а теорія функціонала густини (DFT) переформульовує задачу в термінах густини електронів, а не багатоелектронних хвильових функцій.

Наближення Борна-Оппенгеймера:
  Ψ_total ≈ Ψ_electronic(r;R) ċ Ψ_nuclear(R)
  Електронне рівняння Шредінгера: H_el Ψ_el = E_el(R) Ψ_el   (параметричне за ядерними координатами R)
  ППЕ = E_el(R) = поверхня потенціальної енергії, на якій рухаються ядра

Теорія ЛКАО-МО (молекула = лінійна комбінація атомних орбіталей):
  ψ = Σ_μ c_μ φ_μ
  Секулярне рівняння: det(H − ES) = 0
  H_μν = ⟨φ_μ|H|φ_ν⟩,  S_μν = ⟨φ_μ|φ_ν⟩

Метод Хюккеля (лише π-електрони, S_μν = δ_μν):
  H_μμ = α,  H_μν{ сусідні } = β (<0)
  Бензен C_6: E_k = α + 2βcos(2πk/6), k=0,1,2,3,4,5
  Енергія делокалізації = 2|β| на π-зв'язок відносно еталону етилену

Хартрі-Фок (ХФ):
  Наближення Ψ одним детермінантом Слейтера
  Оператор Фока: F(φ_i) = hφ_i + Σ_j(2J_j − K_j)φ_i
  Кореляційна енергія: E_corr = E_exact − E_HF  (відсутня в ХФ; ~1 еВ/зв'язок)
  MP2, CCSD, CCSD(T): систематичні покращення

Теорія функціонала густини (DFT):
  Теореми Хоенберга-Кона (1964): E однозначно визначається через ρ(r)
  Рівняння Кона-Шема: допоміжна невзаємодіюча система
    [−ℏ²∇²/(2m) + V_ext + V_H + V_xc] ψ_i = ε_i ψ_i
  ρ(r) = Σ_i|ψ_i(r)|²
  Обмінно-кореляційний функціонал E_xc[ρ]: наближується (LDA, GGA, гібридний B3LYP)
  Обчислювальна вартість: O(N³) проти O(N^5) для MP2; робочий кінь обчислювальної хімії
          

6. Фемтохімія

У 1999 році Ахмед Зевейл отримав Нобелівську премію з хімії за розробку фемтосекундної спектроскопії — використання лазерних імпульсів тривалістю 10⁻¹⁵ секунд для «зйомки фільму» хімічних реакцій у реальному часі. Фемтохімія показала, що перехідні стани — не просто математичні конструкції, а спостережувані перехідні структури з часом життя 100–500 фс.

Експеримент накачки-зондування:
  Імпульс накачки: збуджує молекулу до збудженого електронного стану при t=0
  Імпульс зондування (затриманий на Δt): виявляє продукт або проміжну сполуку за LIF або поглинанням
  Часова роздільна здатність обмежена тривалістю імпульсу Δt ~ 50–200 фс

Принцип Франка-Кондона:
  Електронний перехід: ядра не рухаються під час електронного збудження (швидко)
  Фактор ФК: S_vv' = |⟨χ_v'|χ_v⟩|²  (перекриття ядерних хвильових функцій)
  Великий фактор ФК для переходів, що зберігають ядерну геометрію

Динаміка хвильового пакета на збудженій ППЕ:
  Ψ(R,t) = Σ_v c_v χ_v(R) exp(−iE_vt/ℏ)  (когерентна суперпозиція)
  Хвильовий пакет коливається з періодом T_vib ~ 200–1000 фс
  Когерентні коливання видно як модуляції сигналу зондування залежно від Δt

Конічні перетини:
  Точка виродження між двома електронними ППЕ (S1 та S0 для етилену)
  Надшвидка безвипромінювальна релаксація (10−100 фс) через конічний перетин
  Ключова для фотостабільності основ ДНК (розпад тиміну: ~1 пс через S2→CI→S0)

XFEL (рентгенівський лазер на вільних електронах) – з 2010-х донині:
  LCLS (SLAC), European XFEL, SACLA: рентгенівські імпульси менше 100 фс
  Серійна кристалографія: молекулярні фільми динаміки білків
  Динаміка сольватованих електронів (кристал LiF), дифракція нанокристалів
  Δt_XRD ≈ 25 фс → структурна динаміка довжин зв'язків під час реакції
          

Спробуйте ці симуляції