Hayır, olmaz.
Bu postada başlıktaki soruya verdiğimiz olumsuz yanıtı bir karşı örnekle gerekçelendireceğiz. Biyokimyada, bazı moleküllerin farklı izomerlerinin birbirlerine dönüşme kinetiği, oyun şebekesi (play network) dediğimiz yandaki döngüsel sisteme benzer. Şebekede yer alan üç reaksiyonun hız sabitleri sırasıyla α,β,γ ile ilgili reaksiyon okunun üzerinde gösteriliyor. Kimyasal kinetiğin en temel varsayımından, kütle aksiyon kanunundan başlayarak bu şebekede gerçekleşen reaksiyonlarda yer alan maddelerin konsantrasyonları için üç adet hareket denklemi yazacağız.
˙a(t)=−αa(t)+γc(t)˙b(t)=αa(t)−βb(t)˙c(t)=βb(t)−γc(t)
Adet olduğu üzere, bu sistem için denge durumunu soruşturmakla işe başlayacağız. Dinamik bir sistemi oluşturan durum değişkenlerinin hepsinin zamana göre türevinin sıfır olduğu noktalara denge noktaları denir. Bu tanımı kullandığımızda ˙a=˙b=˙c=0 denklemlerinin ortak çözümü, bize denge konsantrasyonlarını verecektir. Basit bir alıştırma ile bu denge konsantrasyonlarını hesaplayabiliriz. ad=βγαβ+αγ+βγmbd=αγαβ+αγ+βγmcd=αβαβ+αγ+βγm
Tanım: (Tersinirlik ve zayıf tersinirlik) Bir kimyasal reaksiyon şebekesinde yer alan bütün tepkimeler tersinir (reversible) ise, o zaman o şebekeye tersinir şebeke denir. Bir kimyasal reaksiyon şebekesinde yer alan reaksiyon oklarının her iki tarafında yer alan reaktant ve ürünlere topluca kompleks denir. Bir kimyasal reaksiyon şebekesinde yer alan herhangi bir C1→C2 reaksiyonu için, C2 kompleksi ile başlayıp C1 kompleksi ile biten bir yol (yani reaksiyon zinciri) bulunabiliyorsa, o zaman o şebekeye zayıf tersinir şebeke (weakly reversible) denir.Bütün tersinir şebekelerin aynı zamanda zayıf tersinir olduğu çok barizdir. Burada çalıştığımız oyun şebekesi tersinir değil. Çünkü, örneğin A→B reaksiyonu var ama B→A yok. Öte yandan B→C→A kanalıyla B maddesinden A maddesine ulaşmak mümkün olduğundan çalıştığımız oyun şebekesi zayıf tersinirdir.
Şebekenin dinamiğini çözmeden önce birimsiz niceliklere geçeceğiz. t=:τ/α, A:=a/m, B:=b/m, C:=c/m, g:=γ/α ve h:=β/α tanımlayalım. A+B+C=1 olduğunu gözleyiniz. Dahası kütlenin korunumu gereğince C=1−A−B yazabiliriz. Bu, ˙c için yazılan diferansiyel denklemi fuzuli (redundant) kılar. İndirgenmiş birimlerde çalışmamız gereken dinamik sistem aşağıdaki gibidir. A′(τ)=−(1+g)A(τ)−gB(τ)+g (1)B′(τ)=A(τ)−hB(τ) (2)
(4) nolu denklem ikinci dereceden, sabit katsayılı, homojen olmayan, lineer bir adi diferansiyel denklemdir ve adi diferansiyel denklemlerin teorisinde çözüm yöntemi vardır. İlkin tekabül eden karakteristik denklemi çözeceğiz. λ2+(1+g+h)λ+(g+h+gh)=0 (5)
Ap:=gh/(g+h+gh) ifadesinin (4) nolu denklemi sağladığını gözleyiniz. (Basitçe yerine koymanız yeterlidir.) İndirgenmiş birimlerde A malzemesinin denge konsantrasyonunu da veren bu ifadeye kısmi çözüm diyeceğiz. Kısmi çözüm aynı zamanda Ap=ad/m denklemini de sağladığından, aslında indirgenmiş birimlerde A maddesinin denge konsantrasyonudur. Artık aradığımız çözümü nihayet verebiliriz. A(τ)=Ap+c1eλ1τ+c2eλ2τ (6a)
Durum I. Δ>0 için (6a) nolu denklemi kullanarak limτ→∞=Ap olduğunu gösterebiliriz. Dahası yine aynı denklemin türevini aldıktan sonra limτ→∞A′(τ)=0 olduğunu da gösterebiliriz. Bu bize (1) nolu denklem kanalıyla limτ→∞B(τ)=g/(g+h+gh) sonucunu verir. İndirgenmiş birimlerde bu, B maddesinin denge konsantrasyonundan başka bir şey değildir. Kütle korunumu ile C maddesinin de τ→∞ limitinde, dengeye geldiği gösterilir. Sistemin dengeye varma süresi, indirgenmiş birimlerde |λ2|−1 ile kestirilebilir. (Neden |λ1|−1 değil?) Birimli niceliklerde bu değer T∼(2/α)/(1+g+h−√Δ) kadardır.
Durum II. Δ=0 için λ1=λ2=:Λ=−(1+g+h)/2 olacaktır. Bu durumda çözümü A(τ)=Ap+d1eΛτ+d2τeΛτ (11)
Durum III. Δ<0 için özdeğerlerin gerçel olmadığını gözleyiniz. Bu durumu çalışmak için öncelikle δ:=√|Δ| tanımıyla işe başlıyoruz. O zaman özdeğerler λ1,2=Λ±iδ2 ile verilecektir. Burada Λ bir önceki paragrafta tanımlandığı gibidir. Bu tanımlarla ve trigonometrik fonksiyonların (2isin(θ)=eiθ−e−iθ ve 2cos(θ)=eiθ+e−iθ gibi) bazı özelliklerini kullanarak A maddesinin konsantrasyonu aşağıdaki gibi verilir. A(τ)=Ap+eΛτ(A(0)−Ap)cos(δτ2)−eΛτ2δ((1+g)A(0)+gB(0)−g+Λ(A(0)−Ap))sin(δτ2)
İşaret. Λ<0 olduğundan, burada da limτ→∞A(τ)=Ap olur, yani tepkime dengeye üstel hızda gelir. Ama (12) nolu denklem şimdiye kadar hiç karşılaşmadığımız bir davranışa,
salınımlara (oscillations) sahiptir. Biyokimyada ve kimya mühendisliğinde salınım yapan tepkimeler önemli bir yer tutar. Ne yazık ki üstel terimin hızla sıfıra gitmesinden ötürü, oyun şebekesinde salınımları uzun süre gözlemek mümkün değildir. Bu meyanda Deficiency-0 teoreminin
pek çok kimyasal sistemde salınımları gözlemenin imkansız olduğunu söylediğini kaydedelim.
İşaret. Birimsiz niceliklerde tepkime zamanı |Λ|−1=2/(1+g+h) ile kestirilebilir. Öte yandan salınımların periyotu tam olarak 4πδ=4π√2(g+h)−(g−h)2−1
Hiç yorum yok:
Yorum Gönder