🌀 Reakcja Biełousowa-Żabotyńskiego

Oregonator to uproszczony trójzmienny model ODE reakcji BZ, zaproponowany przez Fielda i Noyesa w 1974 roku, sprowadzający złożony mechanizm FKN (Field–Körös–Noyes) do trzech kluczowych gatunków: HBrO₂ (X), utlenionego katalizatora (Y) i kwasu bromomalonowego (Z). Równania dX/dt = s(qY + X(1−X) − XY)/ε, dY/dt = (−qY − XY + fZ)/ε′, dZ/dt = X − Z opisują pętle sprzężenia dodatniego (X) i ujemnego (Y, Z). Mimo swojej prostoty Oregonator odtwarza ilościowo okres oscylacji i prędkość fali z rzeczywistych eksperymentów BZ.

W używanym tutaj automacie komórkowym Greenberga-Hastingsa próg określa, ilu wzbudzonych sąsiadów Moore'a (z 8) potrzebuje komórka w spoczynku, aby się zapalić. Niski próg (1 sąsiad) oznacza, że nawet pojedynczy wzbudzony sąsiad wyzwala zapłon — fale rozprzestrzeniają się szybko, dążą do wypełnienia całej domeny, a spirale wirują szybko z krótką długością fali. Wysoki próg (5–6) oznacza, że tylko gęsto wzbudzone obszary wyzwalają sąsiadów, co daje szeroko rozstawione, wolno poruszające się koncentryczne pierścienie, które łatwo wygasają. Rzeczywiste ośrodki BZ mają efektywny próg pobudliwości regulowany przez zmianę stężenia bromianu lub kwasu.

Tak — pobudliwe fale spiralne występują powszechnie w biologii. Migotanie serca (zarówno przedsionków, jak i komór) obejmuje obracające się fale spiralne aktywacji elektrycznej, które zaburzają normalne, okresowe bicie serca; defibrylacja przerywa te spirale poprzez jednoczesną depolaryzację całej tkanki. Ameba społeczna Dictyostelium discoideum wykorzystuje spiralne fale cAMP do koordynowania agregacji podczas głodu. Korowa fala depresji rozprzestrzeniającej się — neurologiczna podstawa aury migrenowej — rozchodzi się jako powolna fala depolaryzacji z prędkością ~3 mm/min. Reakcja BZ była wykorzystywana jako fizyczny model do badania wszystkich tych zjawisk.

Okres refrakcji określa, jak długo komórka musi czekać po zapłonie, zanim będzie mogła zostać ponownie wzbudzona. Dłuższy okres refrakcji (więcej kroków) zwiększa długość przestrzenną spiral — komórki za czołem fali potrzebują więcej czasu na regenerację, więc kolejny grzbiet fali musi przebyć większą odległość, zanim napotka zregenerowany obszar. Bardzo długie okresy refrakcji mogą całkowicie tłumić formowanie się spirali, pozostawiając jedynie izolowaną emisję fal docelowych z obszarów rozrusznikowych. Ten parametr odpowiada bezpośrednio względnemu okresowi refrakcji w elektrofizjologii serca, który reguluje podatność na arytmię.

Tak — układy reakcyjno-dyfuzyjne BZ zostały zaproponowane i zademonstrowane jako niekonwencjonalne podłoża obliczeniowe. Adamatzky i współpracownicy wykorzystali fale BZ w ośrodkach żelowych do rozwiązywania problemów labiryntowych (fala znajduje najkrótszą drogę), implementacji bramek logicznych (kolizja fal → AND/OR), a nawet do przetwarzania obrazu. Wewnętrzny paralelizm układu — każda lokalizacja przestrzenna ewoluuje jednocześnie — oferuje potencjalne korzyści dla pewnych masowo równoległych obliczeń. Bardziej praktycznie, obliczenia reakcyjno-dyfuzyjne pomagają zrozumieć obliczenia morfogenetyczne w biologii rozwoju.

Częstość spontanicznego zapłonu kontroluje, jak często komórki w spoczynku zapalają się bez wystarczającej liczby wzbudzonych sąsiadów — reprezentując szum termiczny lub heterogeniczne miejsca nukleacji w rzeczywistym ośrodku BZ. Przy zerowej częstości spontanicznej siatka ostatecznie osiągnie ciche stan stały, chyba że zostanie ręcznie zasiana. Niskie częstości spontaniczne podtrzymują wzór, okazjonalnie tworząc nowe spirale, gdy istniejące zanikają. Wysokie częstości spontaniczne powodują gęsty, chaotyczny zapłon, który niszczy zorganizowane spirale i tworzy nieuporządkowane fale docelowe — analogicznie do eksperymentu BZ prowadzonego w wysokiej temperaturze, gdzie dominują fluktuacje termiczne.

Praca Alana Turinga z 1952 roku „Chemiczne podstawy morfogenezy” przewidywała, że dwa dyfundujące gatunki chemiczne — aktywator i inhibitor — przy czym inhibitor dyfunduje szybciej niż aktywator, mogą spontanicznie złamać symetrię przestrzenną i generować okresowe wzory (plamy, paski, labirynty). Reakcja BZ demonstruje tę klasę niestabilności Turinga w trzech wymiarach: fale spiralne i docelowe to czasoprzestrzenne wzory Turinga. Ten sam matematyczny szkielet wyjaśnia wzory umaszczenia zwierząt (plamy lamparta, paski zebry), rozstaw palców w rozwijających się kończynach oraz rozmieszczenie fotoreceptorów w oku.

Anihilacja fal przy zderzeniu jest cechą charakterystyczną ośrodków pobudliwych i wynika z okresu refrakcji. Gdy dwa czoła fal się spotykają, komórki w punkcie kolizji zostały właśnie wzbudzone przez oba nadchodzące czoła; zamiast generować nowe wzbudzenie, komórki te wchodzą w stan refrakcji. Za każdym czołem znajduje się strefa refrakcji, która nie może propagować wzbudzenia w przeciwnym kierunku, więc żadna z fal nie kontynuuje ruchu poza punktem kolizji — anihilują się. Jest to zasadniczo inne od ośrodków oscylacyjnych, w których dwie fale przechodzą przez siebie, oraz od fal liniowych, w których zachodzi interferencja konstruktywna/destruktywna.