Mint´ azatk´ epz˝ od´ es k´ emiai reakci´ o-diff´ uzi´orendszerekben Szalai Istv´ an E¨otv¨ os Lor´ and Tudom´ anyegyetem Budapest Szalai Istv´ an ELTE and Tudom´ 2011 1 / 28
Mintazatkepzodes kemiai reakcio-diffuzio rendszerekben
Szalai Istvan
Eotvos Lorand TudomanyegyetemBudapest
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 1 / 28
Az eloadas vazlata
1 Bevezetes
2 Kıserleti eszkoz: nyitott gelreaktorok
3 Terbeli bistabilitas
4 Oszcillacio terben es idoben
5 Turing mintazatok
6 Labirintus mintazatok
7 Lokalizalt mintazatok
8 Osszefoglalas
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 2 / 28
Bevezetes
Bevezetes
Motivacio
,,egy elettelen, de kelloen osszetett kemiaju bolygon nagyon valoszınu,hogy az elet spontan modon kialakul, es egyre osszetettebb esbonyolultabb formakba rendezi onmagat. Csak azt kell megerteni,hogy. . . milyen fizikai torvenyek tettek az elet fele vezeto elso lepestnemcsak lehetsegesse, de elkerulhetetlenne is.”(Ian Stewart: A vegtelen megszelıdıtese, Helikon, 2008)
Milyen kemiai reakciok es milyen korulmenyek szuksegesek komplexviselkedesmodok kialakulasahoz?
Mennyire kulonlegesek, ritkak az ilyen kemiai reakciok?
Hogyan szabalyozhato a dinamikai viselkedes?
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 3 / 28
Bevezetes
Bevezetes
Motivacio
,,egy elettelen, de kelloen osszetett kemiaju bolygon nagyon valoszınu,hogy az elet spontan modon kialakul, es egyre osszetettebb esbonyolultabb formakba rendezi onmagat. Csak azt kell megerteni,hogy. . . milyen fizikai torvenyek tettek az elet fele vezeto elso lepestnemcsak lehetsegesse, de elkerulhetetlenne is.”(Ian Stewart: A vegtelen megszelıdıtese, Helikon, 2008)
Milyen kemiai reakciok es milyen korulmenyek szuksegesek komplexviselkedesmodok kialakulasahoz?
Mennyire kulonlegesek, ritkak az ilyen kemiai reakciok?
Hogyan szabalyozhato a dinamikai viselkedes?
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 3 / 28
Bevezetes
Bevezetes
Reakcio-diffuzio rendszerek
∂c
∂t= f(c, ki ) + D∇2c
f(c, ki ) jellemzoen a koncentraciok polinomialis fuggvenye,D diagonalis matrix
aktivator: ∂fact/∂cact > 0inhibitor: ∂fact/∂cinh < 0idoskalak: τact ∼ 1
kact, τinh ∼ 1
kinh
tavolsagskalak: `act ∼√τactDact , `inh ∼
√τinhDinh
Instabilitasok: pl. Hopf-bifurkacio (τact � τinh),Turing-bifurkacio (`act � `inh). . .
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 4 / 28
Bevezetes
Bevezetes
Reakcio-diffuzio rendszerek
∂c
∂t= f(c, ki ) + D∇2c
f(c, ki ) jellemzoen a koncentraciok polinomialis fuggvenye,D diagonalis matrixaktivator: ∂fact/∂cact > 0inhibitor: ∂fact/∂cinh < 0
idoskalak: τact ∼ 1kact
, τinh ∼ 1kinh
tavolsagskalak: `act ∼√τactDact , `inh ∼
√τinhDinh
Instabilitasok: pl. Hopf-bifurkacio (τact � τinh),Turing-bifurkacio (`act � `inh). . .
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 4 / 28
Bevezetes
Bevezetes
Reakcio-diffuzio rendszerek
∂c
∂t= f(c, ki ) + D∇2c
f(c, ki ) jellemzoen a koncentraciok polinomialis fuggvenye,D diagonalis matrixaktivator: ∂fact/∂cact > 0inhibitor: ∂fact/∂cinh < 0idoskalak: τact ∼ 1
kact, τinh ∼ 1
kinh
tavolsagskalak: `act ∼√τactDact , `inh ∼
√τinhDinh
Instabilitasok: pl. Hopf-bifurkacio (τact � τinh),Turing-bifurkacio (`act � `inh). . .
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 4 / 28
Bevezetes
Bevezetes
Reakcio-diffuzio rendszerek
∂c
∂t= f(c, ki ) + D∇2c
f(c, ki ) jellemzoen a koncentraciok polinomialis fuggvenye,D diagonalis matrixaktivator: ∂fact/∂cact > 0inhibitor: ∂fact/∂cinh < 0idoskalak: τact ∼ 1
kact, τinh ∼ 1
kinh
tavolsagskalak: `act ∼√τactDact , `inh ∼
√τinhDinh
Instabilitasok: pl. Hopf-bifurkacio (τact � τinh),Turing-bifurkacio (`act � `inh). . .
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 4 / 28
Bevezetes
Bevezetes
Autokatalızis
A + mX → nX (n>m)
dcX
dt= (n −m)kcAc
mX
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 5 / 28
Bevezetes
Bevezetes
,,Zart” rendszerben is (tranziens) komplex viselkedest mutato reakciok
Pelda: Belouszov-Zsabotyinszkıj reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 6 / 28
Bevezetes
Bevezetes
,,Zart” rendszerben is (tranziens) komplex viselkedest mutato reakciok
Pelda: Belouszov-Zsabotyinszkıj reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 6 / 28
Bevezetes
Bevezetes
Kevert, nyitott rendszerben komplex viselkedest mutato reakciok
dx
dt= −(x3 − µ0x + λ)− ky
µ0 > 0
dy
dt= (x − y)/τ
J. Boissonade and P. De Kepper
J. Chem. Phys. 1980, 84,
501-506
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 7 / 28
Bevezetes
Bevezetes
Hidrogen-peroxid–szulfit reakcio
SO2−3 + H+ HSO−
3
Autokatalızis: H2O2 + HSO−3 + H+ → SO2−
4 + H2O + 2H+
Inhibicio: H2O2 + 2[Fe(CN)]4−6 + 2H+ → 2[Fe(CN)]3−6 + 2H2O
Rabai,G; Kustin, K; Epstein, I. R. JACS (1989)
CSTR dinamika
Bistabilitas
dccstr
dt= f(ccstr )+
1
τ(c0−ccstr )
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 8 / 28
Bevezetes
Bevezetes
Hidrogen-peroxid–szulfit reakcio
SO2−3 + H+ HSO−
3
Autokatalızis: H2O2 + HSO−3 + H+ → SO2−
4 + H2O + 2H+
Inhibicio: H2O2 + 2[Fe(CN)]4−6 + 2H+ → 2[Fe(CN)]3−6 + 2H2ORabai,G; Kustin, K; Epstein, I. R. JACS (1989)
CSTR dinamika
Oszcillacio
dccstr
dt= f(ccstr )+
1
τ(c0−ccstr )
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 8 / 28
Bevezetes
Bevezetes
Hidrogen-peroxid–szulfit reakcio
SO2−3 + H+ HSO−
3
Autokatalızis: H2O2 + HSO−3 + H+ → SO2−
4 + H2O + 2H+
Inhibicio: H2O2 + 2[Fe(CN)]4−6 + 2H+ → 2[Fe(CN)]3−6 + 2H2ORabai,G; Kustin, K; Epstein, I. R. JACS (1989)
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 9 / 28
Bevezetes
Kevert, nyitott reaktorokban oszcillaciot mutato reakciok
F. Sagues and I. R. Epstein Dalton Trans., 2003, 1201-1217
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 10 / 28
Bevezetes
Reakcio-diffuzio mintazatok
Klorit
Klorit–jodid–malonsav (CIMA)Klordioxid-jod-malonsav (CDIMA) reakcio
Bromat
Belouszov-Zsabotyinszkıj (BZ) reakcio
Jodat
Jodat-szulfit-ferrocianid (FIS) reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 11 / 28
Bevezetes
Reakcio-diffuzio mintazatok
Klorit
Klorit–jodid–malonsav (CIMA)Klordioxid-jod-malonsav (CDIMA) reakcio
Bromat
Belouszov-Zsabotyinszkıj (BZ) reakcio
Jodat
Jodat-szulfit-ferrocianid (FIS) reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 11 / 28
Bevezetes
Reakcio-diffuzio mintazatok
Klorit
Klorit–jodid–malonsav (CIMA)Klordioxid-jod-malonsav (CDIMA) reakcio
Bromat
Belouszov-Zsabotyinszkıj (BZ) reakcio
Jodat
Jodat-szulfit-ferrocianid (FIS) reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 11 / 28
Bevezetes
Mintazatkepzodes
Reakcio-diffuzio rendszerek
∂c
∂t= f(c, ki ) + D∇2c
Instabilitasok: kinetikai vagy diffuziovezerelt
A homogen stacionarius allapot infinitezimalis perturbaciojanakeredmenyekent kialakulo globalis mintaztatok (pl. szuperkritikusHopf, Turing instabilitas)
A homogen stacionarius allapot veges, egy kuszoberteket meghaladoperturbaciojanak eredmenyekent kialakulo mintaztatok. Ilyen modonlokalizalt mintazatok is letre johetnek. (pl. szubkritikus instabilitasok,gerjesztheto es bistabil rendszerek)
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 12 / 28
Kıserleti eszkoz: nyitott gelreaktorok
Nyitott gelreaktorok
CSTR:
dcr
dt= f(cr , k1, k2 . . . kn) +
1
tr(c0 − cr ) + ρV
D
W
[∂c
∂x
]x=0
(1)
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 13 / 28
Kıserleti eszkoz: nyitott gelreaktorok
Nyitott gelreaktorok
GEL:∂c
∂t= f(c) + D∇2c (2)
A CSTR-ban kialakulo koncentraciok Dirichlet hatarfeltetelkent, a masikoldalon a fal Neumann hatarfeltetelkent szolgalnak.
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 14 / 28
Terbeli bistabilitas
Terbeli bistabilitas
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 15 / 28
Oszcillacio terben es idoben
Oszcillacio terben es idoben
τact � τinh
Jodat-szulfit-tiokarbamid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 16 / 28
Oszcillacio terben es idoben
Oszcillacio terben es idoben
τact � τinh
Jodat-szulfit-tiokarbamid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 17 / 28
Oszcillacio terben es idoben
Oszcillacio terben es idoben
τact � τinh
Jodat-szulfit-ferrocianid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 18 / 28
Oszcillacio terben es idoben
Oszcillacio terben es idoben
τact � τinh
Jodat-szulfit-ferrocianid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 19 / 28
Oszcillacio terben es idoben
Oszcillacio terben es idoben: szubkritikus instabilitas(?)
τact � τinh
Hidrogen-peroxid-szulfit-ferrocianid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 20 / 28
Oszcillacio terben es idoben
Diffuziovezerelt instabilitasok
Amennyiben egy pozitıv es negatıv visszacsatolasokat tartalmazoreakcio-diffuzio rendszerben Dact � Dinh, akkor diffuziovezereltinstabilitasok johetnek letre (Turing Philos. Trans. R. Soc. London(1952)).
Megvalosıthato makromolekularis komplexkepzo segıtsegevel:
X + S SX
∂cX
∂t= σ(fX + DX∇2cX )
σ =1
1 + KeqcS
DS = DSX ≈ 0, cS � cX
(Lengyel, Epstein PNAS (1992))
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 21 / 28
Turing mintazatok
Turing mintazatok
`act � `inh
Jodat-szulfit-tiokarbamid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 22 / 28
Turing mintazatok
Turing mintazatok
`act � `inh
Jodat-szulfit-tiokarbamid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 23 / 28
Labirintus mintazatok
Labirintus mintazatok: frontok kozotti kolcsonhatasok
`act � `inh
Jodat-szulfit-ferrocianid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 24 / 28
Lokalizalt mintazatok
Lokalizalt mintaztatok: szubkritikus instabilitas(?)
`act � `inh
Hidrogen-peroxid-szulfit-ferrocianid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 25 / 28
Lokalizalt mintazatok
Lokalizalt mintaztatok: szubkritikus instabilitas(?)
`act � `inh
Hidrogen-peroxid-szulfit-ferrocianid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 26 / 28
Lokalizalt mintazatok
Lokalizalt mintaztatok: szubkritikus instabilitas(?)
`act � `inh
Hidrogen-peroxid-szulfit-ferrocianid reakcio
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 27 / 28
Osszefoglalas
Tanulsag
A kemiailag kulonbozo reakcio-diffuzio rendszerek vizsgalata lehetoeget adarra, hogy a dinamikai viselkedesmodok szeles valasztekat vizsgalhassuk.
CikkekJ. Boissonade, P. De Kepper, F. Gauffre, I. Szalai: ,,Spatial Bistability: A source ofComplex Dymamics. From Spatio-Temporal Reaction-Diffusion Patterns toChemomechanical Structures”Chaos, 2006, 16, 037110.
Szalai, I.; De Kepper, P.: ,,Pattern formation of the ferrocyanide-iodate-sulfite reaction:The control of space scale separation”Chaos, 2008, 18, 026105
Horvath, J.; Szalai, I.; De Kepper, P.: ,,An Experimental Design Method Leading toChemical Turing Patterns”Science, 2009, 324, 772
Horvath, J.; Szalai, I.; De Kepper, P.: ,,Pattern formation in the Thiourea-Iodate-Sulfitesystem: spatial bistability, waves, and stationary patterns”
Physica D, 2009, 239, 776
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 28 / 28
Osszefoglalas
Tanulsag
A kemiailag kulonbozo reakcio-diffuzio rendszerek vizsgalata lehetoeget adarra, hogy a dinamikai viselkedesmodok szeles valasztekat vizsgalhassuk.
CikkekJ. Boissonade, P. De Kepper, F. Gauffre, I. Szalai: ,,Spatial Bistability: A source ofComplex Dymamics. From Spatio-Temporal Reaction-Diffusion Patterns toChemomechanical Structures”Chaos, 2006, 16, 037110.
Szalai, I.; De Kepper, P.: ,,Pattern formation of the ferrocyanide-iodate-sulfite reaction:The control of space scale separation”Chaos, 2008, 18, 026105
Horvath, J.; Szalai, I.; De Kepper, P.: ,,An Experimental Design Method Leading toChemical Turing Patterns”Science, 2009, 324, 772
Horvath, J.; Szalai, I.; De Kepper, P.: ,,Pattern formation in the Thiourea-Iodate-Sulfitesystem: spatial bistability, waves, and stationary patterns”
Physica D, 2009, 239, 776
Szalai Istvan ELTE ( Eotvos Lorand Tudomanyegyetem Budapest )2011 28 / 28