Ratengleichung

Ratengleichungen sind Gleichungen, welche bei gekoppelten chemischen Reaktionen den zeitlichen Konzentrationsverlauf verschiedener chemischer Stoffe beschreiben oder bei An- und Abregungsprozessen den Besetzungsverlauf der Energieniveaus in Atomen oder Molekülen beschreiben.

Hierbei wird für jeden Stoff bzw. Zustand, den das System umfasst, eine Differentialgleichung aufgestellt. Diese Differentialgleichung beschreibt jeweils die Änderung der Konzentration eines beteiligten Stoffes – eben durch die Raten der Zu- und Abflüsse. Es müssen stets alle Beiträge zu den Zu- und Abflüssen betrachtet werden, die die Konzentration des betrachteten Stoffes ändern. Wird für jede Konzentration der an der Reaktion beteiligten Stoffe (Edukt, Produkt, Zwischenprodukte), eine Differentialgleichung aufgestellt, erhält man ein System von gekoppelten Differentialgleichungen. Dieses System gekoppelter Differentialgleichungen kann den zeitlichen Verlauf der betrachteten Ration vollständig beschreiben (Reaktionskinetik). Bei der Beschreibung von Reaktionskaskaden und vielen vorkommenden Stoffen kann die Anzahl der Terme in einer Differentialgleichung und die die Anzahl der Differentialgleichungen selbst schnell unübersichtlich werden. Aus praktischen Gründen werden daher oft von vornherein einige Übergänge ausgeschlossen und die Raten auf Null gesetzt.

Ratengleichung werden i. d. R. genutzt, um experimentell ermittelte Daten zu modellieren. Die Konzentrationsverläufe (Konzentrationsänderung in Abhängigkeit von der Zeit) werden i. d. R. experimentell ermittelt und als Eingaben verwendet. Die Ratenkoeffizienten und weitere Parameter werden numerisch ermittelt und als Ergebnis ausgegeben. Dabei beziehen sich die Konzentrationen nicht ausschließlich auf chemische Reaktionen, sondern auch auf Besetzungsdichten atomarer und molekularer Zustände. Der praktische Nutzen der Ratengleichungen besteht beispielsweise in

• Erkenntnissen über das grundlegenden Verständnis zu Reaktionen (z. B. aufgrund der aktiven Fläche oder der Mitwirkung einer nur kurzlebigen Spezies)
• Erkenntnissen zu Scale-up-Problematiken, da mit Messreihen im Labor auf einen industriellen Maßstab extrapoliert werden kann
• die quantitativen Vergleichbarkeit des Einflusses der im Experiment manipulierten Parameter (Katalysatoren, Gefäßgröße, Energiezufuhr in Form von Wärme, Strahlung, Druck, Quenscher …)
• Nachweis naturwissenschaftlicher Mechanismen, die den Wert eines Ratenkoeffizienten vorhersagen (z. B. Einsteinkoeffizienten, Kascha-Regel, Auswahlregeln)

Nicht naturwissenschaftliche Anwendungsgebiete sind z. B. Modellierungen von Mitarbeiterfluktuationen in den Ebenen eines Unternehmens.

Aus bekannten Raten können auch Vorhersagen und Optimierungssimulationen erstellt werden. Üblicherweise sind die Konzentrationen im Gleichgewicht und die Zeit zum Erreichen des Gleichgewichts von Interesse. Eine analoge Simulation eines Ratenprozesses ist das Stechheberexperiment.^[1]

Schwachpunkte bei der Anwendung von Ratengleichungen sind u. a.

• Anzahl der freien Parameter
• die Uneinigkeit über das zu verwendende System von Zuständen und deren Wechselwirkung
• die Darstellung von Messwerten, Fitwerten und den zugehörigen Fehlern

Die Änderungsrate der Konzentration ${\displaystyle c_{i}}$ einer Spezies ${\displaystyle i}$ ist die Summe der Änderungsraten der Konzentrationen, welche durch verschiedene Reaktionen hervorgerufen werden:

${\displaystyle {\frac {\mathrm {d} c_{i}}{\mathrm {d} t}}=\sum _{j=1}^{N_{\mathrm {R} }}{\frac {\mathrm {d} c_{i,r}}{\mathrm {d} t}}=\sum _{j=1}^{N_{\mathrm {R} }}\nu _{i,j}r_{j},}$,

wobei ${\displaystyle i}$ die betrachtete Spezies, ${\displaystyle c_{i}}$ ihre Konzentration und ${\displaystyle j}$ ein Index ist, der über alle auftretenden Reaktionen (also auch jeweils über die Hin- und Rückreaktion) läuft. ${\displaystyle r_{j}}$ ist die Reaktionsgeschwindigkeit der Reaktion ${\displaystyle j}$. Die Reaktionsgeschwindigkeit ${\displaystyle r_{j}}$ ist proportional zum Produkt der Eduktaktivitäten ${\displaystyle a_{j,l}^{|\nu _{j,l}|}}$ der Reaktion ${\displaystyle j}$, mit der Geschwindigkeitskonstante ${\displaystyle k_{j}}$ der Reaktion als Proportionalitätskonstante:

${\displaystyle r_{j}=k_{j}\prod _{l}a_{j,l}^{|\nu _{j,l}|}}$.

Dadurch ergibt sich die Ratengleichung als:

${\displaystyle {\frac {\mathrm {d} c_{i}}{\mathrm {d} t}}=\sum _{j=1}^{N_{\mathrm {R} }}\nu _{ij}k_{j}\prod _{k=1}^{N_{j}}a_{k}^{b_{kj}},}$

wobei

• ${\displaystyle a_{i}}$ die Aktivität der Spezies ${\displaystyle i\,}$ – häufig wird vereinfachend die Konzentration ${\displaystyle c_{i}}$ verwendet –,
• ${\displaystyle \nu _{ij}}$ die stöchiometrischen Koeffizienten der Spezies ${\displaystyle i}$ in der Reaktion ${\displaystyle j}$,
• ${\displaystyle b_{ij}}$ die Beträge der stöchiometrischen Koeffizienten, falls die Ratengleichung mit Aktivitäten aufgestellt ist, oder die partiellen Reaktionsordnungen der Spezies ${\displaystyle i}$ in der Reaktion ${\displaystyle j}$ (im Allgemeinen ungleich den stöchiometrischen Koeffizienten), falls die Ratengleichung mit Konzentrationen aufgestellt ist,
• ${\displaystyle k_{j}}$ die Ratenkoeffizienten (i. A. Geschwindigkeitskonstanten),
• ${\displaystyle N_{\mathrm {R} }}$ die Anzahl der Reaktionen und
• ${\displaystyle N_{j}}$ die Anzahl der beteiligten Stoffe in der Reaktion ${\displaystyle j}$ ist.

Bei den Ratengleichungen handelt es sich im Allgemeinen um ein System von gekoppelten, steifen, nichtlinearen Differentialgleichungen erster Ordnung, für die die Bausteinerhaltung gelten muss. Im stationären Fall ${\displaystyle (\mathrm {d} c_{i}/\mathrm {d} t=0)}$ ergibt sich das Massenwirkungsgesetz. Ratengleichungen können kompakt mithilfe der Stöchiometrischen Matrix dargestellt werden.

Durch Lösung der Differentialgleichungen erhält man den zeitlichen Verlauf der mittleren Konzentrationen. Um Realisierungen der Konzentrationen inklusive Fluktuationen zu erhalten, können stochastische Simulationen mit dem Gillespie-Algorithmus ausgeführt werden.

Die Ratengleichungen lassen sich für alle beteiligten Spezies herleiten, indem die Kontinuitätsgleichung mit Quell- beziehungsweise Senkterm (bzw. einer Bilanzgleichung) für die Teilchenkonzentrationen aufgestellt werden:

${\displaystyle {\frac {\partial c_{i}}{\partial t}}+\operatorname {div} {\vec {j}}_{i}=f_{i}(\{a_{i}\})}$,

wobei ${\displaystyle f_{i}(\{a_{i}\})}$ der Quellterm ist, welcher von den Aktivitäten ${\displaystyle a_{i}(\{c_{i}\})}$ abhängt. Diese Aktivitäten sind im Allgemeinen wiederum nichttrivial von allen Konzentrationen abhängig.

Da eine Gleichgewichtsreaktion immer eine Hinreaktion und eine Rückreaktion besitzt, existieren die Hin-Reaktionsrate ${\displaystyle k_{\text{+,r}}}$ und die Rück-Reaktionsrate ${\displaystyle k_{\text{-,r}}}$.

Der Quellterm ist durch eine Summe über alle Reaktionen gegeben:

${\displaystyle f_{i}(\{a_{i}\})=\sum _{r}|\nu _{i,r}|\left[-k_{\text{-,r}}\prod _{{\text{Reaktanden j}},r}a_{r,j}^{|\nu _{r,j}|}+k_{\text{+,r}}\prod _{{\text{Produkte j}},r}a_{r,j}^{|\nu _{r,j}|}\right].}$

Man beachte, dass die partielle Reaktionsordnung (der Exponent, mit dem die Konzentrationen eingehen) nur dann dem Betrag der stöchiometrischen Koeffizienten entspricht, wenn Aktivitäten verwendet werden. Werden statt Aktivitäten ebenfalls Konzentrationen im Quellterm verwendet und liegen Teilchenwechselwirkungen vor, so sind die Beträge der stöchiometrischen Koeffizienten mit den partiellen Reaktionsordnungen zu ersetzen. Die partielle Reaktionsordnung kann beliebige Werte annehmen (z. B. 0) und wird experimentell bestimmt.

• Im Gleichgewicht gibt es keine Teilchenströme (${\displaystyle {\vec {j}}_{i}={\vec {0}}\,\forall i}$) und die Teilchenkonzentrationen ändern sich nicht mehr zeitlich. Daher gilt im Gleichgewicht:

${\displaystyle 0=f_{i}(\{a_{i}\}).}$

Unter der Annahme, dass jede Reaktion (als Paar von Hin- und Rückreaktion) im Einzelnen ausgeglichen ist, erhält man für jede Reaktion das Massenwirkungsgesetz durch Umformung:

${\displaystyle K_{r}={\frac {k_{\text{-,r}}}{k_{\text{+,r}}}}={\frac {\prod _{{\text{Produkte j}},r}a_{r,j}^{|\nu _{r,j}|}}{\prod _{{\text{Reaktanden j}},r}a_{r,j}^{|\nu _{r,j}|}}}}$

• Ist das System im Nichtgleichgewicht, aber homogen, so treten keine Teilchenströme auf, jedoch ändern sich die Konzentrationen zeitlich, bis das Gleichgewicht erreicht ist:

${\displaystyle {\frac {\partial c_{i}}{\partial t}}=f_{i}(\{a_{i}\})}$

• Für den Fall, dass man ein inhomogenes System im Nichtgleichgewicht betrachtet, ist der Teilchenstrom ${\displaystyle {\vec {j}}_{i}\neq 0}$ und kann durch das erste Fick’sche Gesetz

${\displaystyle {\vec {j}}_{i}=-{\vec {\nabla }}(D_{i}c_{i}({\vec {r}}))+{\vec {j}}_{i}^{\text{excess}}}$

beschrieben werden (wobei der nicht-ideale Exzess-Term nur für nicht-ideale Systeme auftritt). Man erhält dann eine Reaktionsdiffusionsgleichung.

• Für den Fall, dass es zusätzlich Strömung im System gibt, ist Konvektion im Teilchenstrom zu berücksichtigen und man erhält die Konvektions-Diffusions-Gleichung.

Die in den Ratengleichungen auftretenden Reaktionsratenkoeffizienten können allgemein als beliebige Funktionen der jeweiligen, gegebenenfalls zeitabhängigen, Temperatur (siehe auch Plasmaphysik: Thermisches Gleichgewicht) betrachtet werden. Im Allgemeinen müssen Ratenkoeffizienten für chemische Prozesse der schweren Teilchen aus der Literatur entnommen werden ('Geschwindigkeitskonstante' einer chemischen Reaktion), die Ratenkoeffizienten für die elektronenstoßinduzierten Prozesse können mit Hilfe der Elektronenkinetik erhalten werden.

Grundlage für die kinetische Behandlung der Elektronen, sowohl zur Berechnung derartiger Ratenkoeffizienten, als auch elektronischer Transportprozesse (elektrische Leitfähigkeit) bildet die Boltzmann-Gleichung für die Elektronenenergieverteilung.

Zur Verdeutlichung wird die Wasserstoffoxidation herangezogen:

ein Teil dissoziiert

Die Ratengleichungen (Gl.1) für die fünf Spezies lauten:

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\mathrm {[H_{2}]} =-2k_{1}\mathrm {[H_{2}]^{2}[O_{2}]} }$

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\mathrm {[O_{2}]} =-k_{1}\mathrm {[H_{2}]^{2}[O_{2}]} }$

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\mathrm {[H_{2}O]} =+2k_{1}\mathrm {[H_{2}]^{2}[O_{2}]} -2k_{2}\mathrm {[H_{2}O]^{2}} }$

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\mathrm {[H_{3}O^{+}]} =+k_{2}\mathrm {[H_{2}O]^{2}} }$

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\mathrm {[OH^{-}]} =+k_{2}\mathrm {[H_{2}O]^{2}} }$

Die Konzentrationen der Spezies:

${\displaystyle c_{1}\equiv \mathrm {[H_{2}]} ,\;\;c_{2}\equiv \mathrm {[O_{2}]} ,\;\;c_{3}\equiv \mathrm {[H_{2}O]} ,\;\;c_{4}\equiv \mathrm {[H_{3}O^{+}]} ,\;\;c_{5}\equiv \mathrm {[OH^{-}]} }$

Oszillierende Reaktionen werden durch Ratengleichungen beschrieben. Für spezielle Modelle hierzu, siehe Oregonator und Brüsselator. Die numerische Lösung solcher Differentialgleichungssysteme liefert dann oszillierende chemische Konzentrationen^[2].

Die Wechselwirkung von Räuber- und Beutepopulationen wird durch die Lotka-Volterra-Gleichungen beschrieben.

Da es sich bei den Ratengleichungen um ein System von steifen Differentialgleichungen handelt, ist man gezwungen ein Verfahren mit einem möglichst großen Stabilitätsgebiet zu wählen, damit die Integrationsschritte nicht allzu klein werden. Am günstigsten sind A-stabile Verfahren.

Für die Ratengleichungen bedeutet 'steif', dass sich die Zeitkonstanten der verschiedenen Spezies sehr stark unterscheiden: Im Verhältnis zu anderen ändern sich einige Konzentrationen nur sehr langsam. Zwei Beispiele absolut steif-stabiler Integrationsverfahren sind die Implizite Trapez-Methode und die Implizite Euler-Methode, ebenso sind einige BDF-Verfahren (backward differentiation formula) geeignet.

Das Prinzip der Bausteinerhaltung liefert eine Möglichkeit, die Güte der numerischen Lösungen zu überprüfen, denn es gilt zu jedem Zeitpunkt:

${\displaystyle \sum _{i=1}^{N_{\mathrm {Sp} }}c_{i}\beta _{ik}=\gamma _{k}={\text{const.}},\quad \forall k=1,\dots ,N_{\mathrm {B} }}$

wobei

${\displaystyle N_{\mathrm {B} }\,}$    Minimale Anzahl der Bausteine,

${\displaystyle N_{\mathrm {Sp} }\,}$   Anzahl an den Reaktionen beteiligten Spezies.

Eine Spezies i, hier geschrieben als ${\displaystyle A_{i}}$ setzt sich dabei aus den Bausteinen ${\displaystyle B_{k}}$ folgendermaßen zusammen:

${\displaystyle A_{i}=\sum _{k=1}^{N_{\mathrm {B} }}\beta _{ik}B_{k}\quad \Rightarrow \quad \sum _{i=1}^{N_{\mathrm {Sp} }}\nu _{ij}\beta _{ik}={\underline {\underline {0}}}}$  .

in die Ratengleichung (Gl.1) eingesetzt und über alle Spezies summiert, liefert wegen ${\displaystyle \sum _{i=1}^{N_{\mathrm {Sp} }}{\dot {c}}_{i}\beta _{ik}=0}$ die oben genannte Bausteinerhaltung.

${\displaystyle {\begin{pmatrix}\mathrm {H_{2}} \\\mathrm {O_{2}} \\\mathrm {H_{2}O} \\\mathrm {H_{3}O^{+}} \\\mathrm {OH^{-}} \end{pmatrix}}=\underbrace {\begin{pmatrix}2&0&0\\0&2&2\\2&1&1\\3&1&0\\1&1&2\end{pmatrix}} _{\beta _{ik}}\cdot {\begin{pmatrix}\mathrm {H} \\\mathrm {O^{+}} \\e^{-}\end{pmatrix}}}$

• Mastergleichung

• W. Frie: Berechnung der Gaszusammensetzungen und der Materialfunktionen von ${\displaystyle SF_{6}}$. In: Zeitschrift für Physik, 201, 269, 1967; Springer-Verlag, Berlin / Heidelberg/ New York
• C. Schwab: Beiträge zur kinetischen Modellierung von teilweise ionisierten Nichtgleichgewichtsplasmen. Dissertation an der Fakultät für Physik der Eberhard Karls Universität Tübingen, 1989
• H. R. Schwarz: Numerische Mathematik. B. G. Teubner, Stuttgart 1986, ISBN 3-519-02960-X
• G. Wedler: Lehrbuch der Physikalischen Chemie, Wiley-VCH, 2004, ISBN 3-527-31066-5
• D. A. McQuarrie, J. D. Simon, J. Choi: Physical Chemistry: A Molecular Approach. University Science Books, 1997, ISBN 0-935702-99-7

1. ↑ 2.3 Modellexperiment - Stechheberversuch, auf w-hoelzel.de
2. ↑ 9.11: Oscillating Reactions. In: Chemistry Library. LibreTexts, 1. September 2020, abgerufen am 7. September 2020 (englisch).