Chemia i kinetyka chemiczna — reakcja-dyfuzja, spalanie i równowagi kwasowo-zasadowe

Chemia jako układ dynamiczny

Każda reakcja chemiczna jest układem dynamicznym: stężenia zmieniają się w czasie zgodnie z równaniami różniczkowymi, równowagi pojawiają się jako punkty stałe, a nieliniowe sprzężenie między reakcjami może wytwarzać oscylacje, fale i samoorganizujące się wzory przestrzenne. Ten sam aparat matematyczny, który opisuje wahadło — portrety fazowe, analizę stabilności, diagramy bifurkacji — stosuje się wprost do kinetyki chemicznej.

To sprawia, że chemia wyjątkowo dobrze nadaje się do symulacji. W przeciwieństwie do żywego eksperymentu laboratoryjnego symulacja pozwala zamrozić czas, cofnąć go, przybliżyć pojedynczy front reakcji i natychmiast zmienić stałe szybkości, by zobaczyć, co dzieje się na granicy niestabilności. Sześć symulacji w tym spotlightcie obejmuje zakres od reakcji trwających mikrosekundy po cykle biogeochemiczne rozciągnięte na stulecia.

Warstwa 1: Samoorganizacja — wzory Turinga

Symulacja Reakcja-Dyfuzja

W 1952 roku Alan Turing opublikował artykuł, który nie miał nic wspólnego z komputerami. Pokazał, że układ dwóch substancji chemicznych — aktywatora, który wspomaga własną produkcję, i inhibitora, który ją tłumi — w połączeniu z różnymi szybkościami dyfuzji może samoistnie złamać symetrię przestrzenną i utworzyć stabilne wzory: plamki, paski i labirynty. Te „wzory Turinga" tłumaczą pasy zebry, plamki lamparta i grzbiety linii papilarnych.

Równania różniczkowe cząstkowe (ciągłe):
  ∂U/∂t = D_U · ∇²U  -  U·V²  +  f·(1-U)
  ∂V/∂t = D_V · ∇²V  +  U·V²  -  (f+k)·V

U = stężenie aktywatora  [0, 1]
V = stężenie inhibitora  [0, 1]
D_U, D_V = współczynniki dyfuzji  (D_U > D_V dla powstawania wzorów)
f = szybkość zasilania (uzupełnianie U)
k = szybkość usuwania (usuwanie V)

Reżimy wzorów (przestrzeń parametrów f, k):
  f=0.035, k=0.065 → plamki (lampart)
  f=0.060, k=0.062 → paski (zebra)
  f=0.025, k=0.060 → poruszające się plamki (koral)
  f=0.014, k=0.054 → labirynt (odcisk palca)
  f=0.039, k=0.058 → robaki + plamki

Implementacja WebGL:
  Bufory ramek ping-pong (A↔B w każdej klatce)
  Fragment shader: Laplasjan przez 9-punktowy szablon
  Format tekstury: RG32F (U w R, V w G)
  Rozdzielczość: 512×512 → 262 143 komórki aktualizowane na klatkę

Symulacja Reakcja-Dyfuzja implementuje układ Gray-Scotta na siatce WebGL 512×512 aktualizowanej w czasie rzeczywistym. Przesuwaj się po przestrzeni parametrów (f, k) i obserwuj, jak wzór się przekształca — od izolowanych plamek przez labiryntowe kanały po fale wędrujące — bez restartowania symulacji. Maluj kursorem początkowe zarodki i patrz, jak wzór rośnie spod twojego dotyku.

Warstwa 2: Szybka kinetyka — spalanie

Symulacja Spalanie

Spalanie to najbardziej energetycznie doniosła chemia na Ziemi — napędza każdy samolot, statek i większość sieci elektroenergetycznych. To też jedna z najbardziej złożonych reakcji: prawdziwe płomienie węglowodorowe obejmują setki jednoczesnych reakcji, rozgałęzienie łańcuchowe, wolne rodniki i sprzężenie zwrotne cieplne. Symulacja skupia się na kluczowej fenomenologii: zapłonie, propagacji frontu płomienia i warunkach rozróżniających deflagrację od detonacji.

Prawo szybkości Arrheniusa:
  k(T) = A · exp(-E_a / R·T)
  A   = czynnik przedwykładniczy [mol⁻¹·m³·s⁻¹]
  E_a = energia aktywacji [J/mol]
  R   = 8,314 J/(mol·K)
  Podwojenie T z 300K do 600K może zwiększyć k nawet 10⁶ razy

Rozgałęzienie łańcuchowe H₂/O₂ (uproszczone):
  H + O₂ → OH + O    (rozgałęzienie łańcuchowe, E_a=70 kJ/mol)
  O + H₂ → OH + H    (rozgałęzienie łańcuchowe)
  H + O₂ + M → HO₂ + M  (zakończenie łańcucha przy wysokim P)

Granice wybuchowości:
  Pierwsza (niskie P):   zakończenie na ścianach > rozgałęzienie → brak wybuchu
  Druga (~1 atm): rozgałęzienie > zakończenie  → wybuch
  Trzecia (wysokie P):  ucieczka cieplna → wybuch

Deflagracja a detonacja:
  Deflagracja: płomień porusza się z prędkością 0,5–5 m/s, poddźwiękowo
  Detonacja:   płomień sprzężony z falą uderzeniową, 1500–3000 m/s, naddźwiękowo
  Prędkość Chapmana-Jougueta: D_CJ = √(2(γ²-1)·q)  (q = wydzielone ciepło)

Symulacja Spalanie modeluje 2D przepływ reaktywny za pomocą uproszczonego jednostopniowego mechanizmu Arrheniusa. Ustaw stosunek równoważnikowy paliwa do powietrza, temperaturę początkową i ciśnienie, a następnie zapal iskrą. Obserwuj, jak front płomienia przyspiesza, gdy ciepło wydzielane w strefie reakcji podgrzewa niespaloną mieszankę. Zmień stosunek równoważnikowy poniżej granicy ubogiej i zobacz, jak płomień gaśnie; idź w kierunku bogatej mieszanki i zobacz niepełne spalanie oraz powstawanie prekursorów sadzy.

Warstwa 3: Chemia równowag — reakcje kwasowo-zasadowe

Symulacja Kwasy-Zasady

Chemia kwasowo-zasadowa leży u podstaw biochemii (centra aktywne enzymów, buforowanie krwi), nauk o środowisku (zakwaszenie oceanów, kwaśne deszcze) i procesów przemysłowych (miareczkowania, kontrola pH w reaktorach). Matematyka równowagi jest elegancka: garstka stałych równowagi w pełni określa pH i specjację nawet złożonych układów wieloskładnikowych.

Równowaga wody:
  Kw = [H⁺][OH⁻] = 1×10⁻¹⁴  (w 25°C)
  pH + pOH = 14

Równowaga słabego kwasu:
  HA ⇌ H⁺ + A⁻
  Ka = [H⁺][A⁻] / [HA]       pKa = -log₁₀(Ka)

Równanie Hendersona-Hasselbalcha:
  pH = pKa + log₁₀([A⁻]/[HA])
  Zakres buforowy: pKa ± 1  (stosunek 10:1 do 1:10)

Krzywa miareczkowania (mocny kwas kontra mocna zasada):
  Przed pkt. rów.:  pH = -log[H⁺]_nadmiar
  W pkt. rów.:      pH = 7 (czysta woda)
  Po pkt. rów.:     pH = 14 + log[OH⁻]_nadmiar

Miareczkowanie słabego kwasu:
  Pkt. półrów.:    pH = pKa  (środek strefy buforowej)
  Pkt. równoważnikowy:  pH > 7  (hydroliza sprzężonej zasady)

Przykłady:
  Kwas octowy:       pKa = 4,76  (ocet)
  Kwas węglowy:      pKa = 6,35  (krew, CO₂ oceaniczny)
  Amon:              pKa = 9,25  (środki czyszczące)
  Kwas fosforowy:    pKa = 2,15, 7,20, 12,35  (bufor biologiczny)

Symulacja Kwasy-Zasady renderuje interaktywne krzywe miareczkowania dla dowolnej kombinacji mocnych/słabych kwasów i zasad. Dodawaj wirtualne krople titranta suwakiem i obserwuj, jak krzywa pH gwałtownie skacze w punkcie równoważnikowym. Włącz nakładkę wskaźnikową — lakmus, fenoloftaleinę, błękit bromotymolowy — i zobacz, jaki zakres barw obejmuje każdy z nich. Wykres pojemności buforowej pokazuje dokładnie, gdzie roztwór najskuteczniej opiera się zmianie pH.

Warstwa 4: Przejście fazowe — wzrost kryształów

Symulacja Wzrost Kryształów

Gdy przesycony roztwór lub przechłodzony stop zaczyna krystalizować, kształt rosnącego kryształu określa rywalizacja między anizotropią energii powierzchniowej a transportem ograniczonym dyfuzją. Ta sama matematyka rządzi dendrytami płatków śniegu, wzrostem ziaren metalu w stopach oraz strukturami mineralnymi w formacjach geologicznych.

Agregacja ograniczona dyfuzją (DLA):
  Losowy błądzący uwalniany z granicy
  Przykleja się przy pierwszym kontakcie z rosnącym klastrem
  Wymiar fraktalny D ≈ 1,71 w 2D

Model pola fazowego (Kobayashi 1993):
  ∂φ/∂t = ε²·∇²φ + φ(1-φ)(φ - 1/2 + m)
  ∂T/∂t = ∇²T + ∂φ/∂t      (sprzężenie z ciepłem utajonym)

  φ ∈ [0,1]: parametr porządku (ciecz=0, ciało stałe=1)
  ε = parametr grubości granicy fazowej
  m = siła napędowa = (T_topnienia - T) / T_topnienia

Anizotropia (ramiona dendrytyczne):
  ε(θ) = ε̄ · (1 + δ·cos(N·θ))
  N=4 → kryształ kubiczny (6-krotna → heksagonalna)
  δ=0,05: słaba anizotropia → zwarty kryształ
  δ=0,15: silna anizotropia → wyraźne dendryty

Kryterium płatka śniegu:
  N=6, δ≈0,05, D=0,5 → heksagonalne dendryty płatka śniegu
  Rozstaw ramion ∝ √(D·t)  (ograniczone dyfuzją)

Symulacja Wzrost Kryształów pozwala wybrać między DLA (agregacją losowych błądzeń) a modelem pola fazowego. Obejrzyj, jak DLA tworzy charakterystyczny fraktalny kształt kalafiora, a następnie przełącz na pole fazowe i dostrój siłę anizotropii, aby wyhodować dendryty przypominające płatki śniegu. Zmiana szybkości chłodzenia przesuwa morfologię od zwartych ziaren po długie, wiotkie ramiona — dokładnie to zachowanie, które metalurdzy kontrolują, projektując mikrostrukturę odlewów.

Warstwa 5: Dynamika leków — farmakokinetyka

Symulacja Farmakokinetyka

Po połknięciu tabletki lek przechodzi precyzyjnie mierzalną drogę: wchłanianie z jelit do osocza, dystrybucję do przedziałów ustrojowych, metabolizm (głównie w wątrobie) i wydalanie (głównie przez nerki). Ten proces ADME jest rządzony równaniami różniczkowymi pierwszego rzędu, które musi znać każdy farmakolog i projektant leków.

Model jednokompartmentowy (bolus dożylny):
  dC/dt = -k_e · C
  C(t) = C₀ · e^{-k_e·t}
  Okres półtrwania: t₁/₂ = ln(2)/k_e = 0,693/k_e

Model dwukompartmentowy:
  dC_p/dt = -(k_e + k_pt)·C_p + k_tp·C_t + F·D(t)/V_p
  dC_t/dt =   k_pt·C_p - k_tp·C_t

  C_p = stężenie w osoczu
  C_t = stężenie w tkance
  k_pt = dystrybucja do tkanki, k_tp = redystrybucja
  F = biodostępność (doustnie: F < 1)

Kluczowe parametry:
  V_d = objętość dystrybucji = dawka / C₀
  CL = klirens = k_e · V_d  [L/h]
  AUC = ∫₀^∞ C(t)dt = dawka·F / CL
  Stan stacjonarny (dawki wielokrotne):
    C_ss = F·dawka / (CL·τ)     τ = odstęp między dawkami

Okno terapeutyczne:
  MEC (min. skuteczne) < C_ss < MTC (min. toksyczne)
  Zaprojektuj schemat dawkowania tak, by pozostać w oknie

Symulacja Farmakokinetyka implementuje model dwukompartmentowy z konfigurowalną dawką początkową, biodostępnością, klirensem i objętością dystrybucji. Wybierz schemat dawkowania — pojedynczy bolus dożylny, powtarzane dawki doustne lub wlew — i obserwuj, jak zmienia się krzywa stężenia w osoczu. Nakładka okna terapeutycznego od razu pokazuje, czy wybrany schemat utrzymuje pacjenta w strefie skutecznej, czy przekracza próg toksyczności.

Warstwa 6: Chemia planetarna — cykl węglowy

Symulacja Cykl Węglowy

Globalny cykl węglowy to chemia w skali planetarnej: atomy węgla przepływają między atmosferą, oceanem, biosferą lądową i zbiornikami geologicznymi w skalach czasowych od sekund (fotosynteza) po miliony lat (formowanie się skał węglanowych). Matematyka tych przepływów — w istocie ogromny układ równań różniczkowych zwyczajnych pierwszego rzędu ze sprzężeniami nieliniowymi — decyduje o stężeniu CO₂ w atmosferze, a więc i o globalnej temperaturze.

Cztery główne zbiorniki (GtC):
  Atmosfera:  ~870  (rośnie ~5 GtC/rok)
  Biosfera lądowa: ~2600  (NPP ≈ 120 GtC/rok)
  Powierzchnia oceanu:  ~900  (wymiana ~90 GtC/rok)
  Głęboki ocean: ~37 000

Powiązanie z chemią oceanu:
  CO₂(g) ⇌ CO₂(aq) → H₂CO₃ → H⁺ + HCO₃⁻ → H⁺ + CO₃²⁻
  pCO₂ = α · [CO₂(aq)]     Prawo Henry'ego
  pH = -log[H⁺]           Zakwaszenie oceanu: ΔpH = -0,1 od 1750 roku

Model boxowy (uproszczony):
  dC_atm/dt = E_paliwa_kopalne + E_zmiana_użytkowania - F_ląd - F_ocean
  dC_ocean/dt = F_ocean - F_węglan - F_eksport
  E_paliwa_kopalne ≈ 10 GtC/rok  (lata 2020.)
  Frakcja pozostająca w atmosferze: ~44% emisji

Sprzężenie zwrotne temperatury:
  ΔT = λ · ΔF = λ · 3,7 · ln(CO₂/CO₂_0)  (wymuszenie w W/m²)
  λ ≈ 0,8 K/(W/m²)   (równowagowa wrażliwość klimatyczna)

Symulacja Cykl Węglowy animuje przepływy węgla między zbiornikami atmosfery, lądu i oceanu w czasie rzeczywistym. Zwiększ emisje paliw kopalnych, wywołaj wylesianie albo pobudź nawożenie żelazem oceanów i obserwuj, jak w odpowiedzi rośnie CO₂ i temperatura. Obok biegnie wykres zakwaszenia oceanu, pokazujący, jak chemia węglanowa ogranicza zdolność buforową oceanu na przestrzeni dekad.

Pełna kolekcja chemiczna

Połączenia między kolekcjami

Chemia łączy się niemal z każdą inną kategorią na platformie. Prawo szybkości Arrheniusa rządzące rozgałęzieniem łańcuchowym w spalaniu to ten sam czynnik wykładniczy co w symulacji Maxwell-Boltzmann, która tłumaczy, dlaczego gorące cząsteczki gazu reagują szybciej — obie to czynniki Boltzmanna. Mechanizm formowania wzorów reakcja-dyfuzja powraca w symulacji Automaty Komórkowe — reguła Gray-Scotta i reguła CA typu aktywator-inhibitor mają tę samą strukturę matematyczną. Farmakokinetyka jest bezpośrednim zastosowaniem metod równań różniczkowych omówionych w Learning #16 (równania różniczkowe w biologii), a model boxowy cyklu węglowego wykorzystuje dokładnie to samo podejście kompartmentowe. Wzrost kryształów metodą DLA łączy się z geometrią fraktalną w Spotlight #15 (Chaos) — obie generują struktury o niecałkowitym wymiarze fraktalnym.

Algorytmy i metody w tej kolekcji