Schnelle Reaktionsmethoden und Relaxationskinetik – Stopped-Flow, Flash-Photolyse, Temperatur- und Drucksprung, Eigenmoden und globale Analyse

Titelbild: Eigene Illustration schneller Transientensignale, exponentieller Relaxation und zentraler Beziehungen der linearen Relaxationskinetik.

Viele chemische Elementarschritte laufen in Mikrosekunden, Nanosekunden oder noch schneller ab. Eine gewöhnliche Probenahme ist dann zu langsam. Statt die Reaktion in einzelnen Zeitpunkten anzuhalten, muss sie abrupt gestartet oder ein Gleichgewicht sehr schnell gestört werden. Anschließend wird die Rückkehr oder Produktbildung mit hoher Zeitauflösung verfolgt.

Die wichtigsten Strategien sind:

• schnelles Mischen, etwa mit Stopped-Flow, Continuous-Flow oder Quench-Flow,
• Pulsanregung, etwa durch Flash-Photolyse oder Laser-Pump-Probe,
• Relaxationsmethoden, etwa Temperatur-, Druck-, Feld- oder Konzentrationssprung.

Die Relaxationskinetik besitzt einen besonderen Vorteil: Sie untersucht die Rückkehr eines Systems zum Gleichgewicht nach einer kleinen Störung. Für ein reversibles Zweizustandssystem:

A ⇌ B

relaxiert eine kleine Abweichung exponentiell:

Δx(t)=Δx₀e^−t/τ

mit:

τ⁻¹=k₁+k₋₁

Damit liefert eine Relaxationsmessung direkt die Summe der Vorwärts- und Rückwärtskonstanten. Zusammen mit der Gleichgewichtskonstante können beide Einzelkonstanten bestimmt werden.

Einordnung: Teil 7 behandelte diffusionskontrollierte Reaktionen und Lösungsmittelreibung. Teil 8 zeigt, wie schnelle Prozesse experimentell ausgelöst, beobachtet und mit linearen sowie nichtlinearen kinetischen Modellen ausgewertet werden. Teil 9 behandelt anschließend photochemische Kinetik und angeregte Zustände.

Lernziele

Nach diesem Beitrag kannst du:

• direkte Initiations- und Relaxationsmethoden unterscheiden,
• Aufbau, Dead Time und Messbereich eines Stopped-Flow-Experiments erklären,
• Flash-Photolyse, Laser-T-Jump und Drucksprung einordnen,
• pseudo-erste Ordnung gezielt verwenden,
• die Relaxationszeit eines reversiblen Zweizustandssystems herleiten,
• kleine Störungen mit einer Jacobi-Matrix linearisieren,
• Relaxationseigenwerte und Eigenmoden interpretieren,
• mehrfache Exponentialanteile von mechanistischen Spezies unterscheiden,
• Instrumentenantwort und Totzeit in ein kinetisches Modell einbeziehen,
• globale spektral-kinetische Anpassungen beurteilen,
• Identifizierbarkeit, Residuen und Parameterkorrelationen prüfen,
• zehn typische Rechnungen der schnellen Kinetik durchführen.

1. Chemische Zeitskalen

Chemische Dynamik umfasst einen enormen Zeitbereich:

Zeitskala	typische Prozesse	geeignete Methoden
Femtosekunden bis Pikosekunden	elektronische Relaxation, Schwingungsbewegung, ultrafast Protonen- und Elektronentransfer	Pump-Probe, transient absorption, Fluoreszenz-Upconversion
Nanosekunden bis Mikrosekunden	Radikale, angeregte Zustände, schnelle Konformationsänderungen	Flash-Photolyse, Laser-T-Jump, zeitaufgelöste Fluoreszenz
Millisekunden bis Sekunden	Enzymreaktionen, Ligandenbindung, Ionenreaktionen, Faltungsprozesse	Stopped-Flow, Drucksprung, Leitfähigkeits- und Absorptionskinetik
Sekunden bis Stunden	langsame Umsetzungen, Alterung, Aggregation	konventionelle Zeitreihen und Probenahme

Die Zeitauflösung der Methode muss deutlich kürzer sein als die schnellste relevante kinetische Komponente.

2. Direkte Initiation und Relaxation

Bei einer direkten Initiationsmethode wird eine neue Population erzeugt:

• Reaktanden werden gemischt,
• ein Zwischenprodukt wird photochemisch erzeugt,
• ein gebundener Ligand wird freigesetzt,
• ein Elektron oder Proton wird injiziert.

Bei einer Relaxationsmethode befindet sich das System zunächst im Gleichgewicht. Ein äußerer Parameter wird rasch verändert, wodurch das alte Gleichgewicht nicht mehr zum neuen Wert von Temperatur, Druck oder Feld passt. Das System relaxiert zum neuen Gleichgewicht.

3. Kleine und große Störungen

Eine kleine Störung erlaubt eine Linearisierung der nichtlinearen Geschwindigkeitsgleichungen um den Gleichgewichtspunkt. Die Relaxation wird dann durch Exponentialmoden beschrieben.

Eine große Störung kann:

• stark nichtlineare Konzentrationsabhängigkeiten,
• Mechanismuswechsel,
• Temperatur- oder Druckabhängigkeit der Raten während des Sprungs,
• nicht-exponentielle Rückkehr

erzeugen. Dann muss das vollständige kinetische Modell numerisch integriert werden.

4. Stopped-Flow-Prinzip

Beim Stopped-Flow werden zwei oder mehr Lösungen aus Spritzen schnell durch einen Mischer in eine Beobachtungszelle gedrückt. Ein Stoppspritzenmechanismus beendet den Fluss abrupt. Ab diesem Zeitpunkt wird das Signal als Funktion der Zeit gemessen.

Die Reaktion startet bereits im Mischer. Deshalb ist die früheste beobachtbare Zeit nicht null, sondern die Summe aus Misch- und Transportzeit.

5. Komponenten eines Stopped-Flow-Geräts

• Präzisionsspritzen oder pneumatische Antriebe,
• Mikromischer mit kleinem Totvolumen,
• Beobachtungszelle,
• Stoppspritze oder mechanischer Anschlag,
• optische, elektrische oder kalorimetrische Detektion,
• schnelle Datenerfassung und Triggerung.

Die Mischqualität hängt von Strömungsgeometrie, Viskosität, Dichte und Flussrate ab.

6. Dead Time

Die Dead Time ist die Zeit zwischen dem ersten wirksamen Kontakt der Lösungen und dem frühesten zuverlässigen Messpunkt. Eine grobe Abschätzung lautet:

t_dead≈V_tot/Q

V_tot ist das Volumen zwischen Mischregion und Beobachtung, Q der Volumenstrom.

Liegt die Reaktionszeit in der Größenordnung der Dead Time, ist bereits ein erheblicher Umsatz erfolgt, bevor die Messung beginnt.

7. Mischzeit und chemische Zeit

Die Mischzeit ist nicht nur eine instrumentelle Zahl. Wenn eine Reaktion während des Mischens läuft, überlagern sich räumliche Konzentrationsgradienten und Chemie.

Ein brauchbares Experiment verlangt im Allgemeinen:

t_mix≪τ_chem

Andernfalls muss der Mischprozess modelliert oder eine schnellere Methode gewählt werden.

8. Continuous-Flow

Beim Continuous-Flow werden Reaktanden kontinuierlich gemischt und entlang eines Strömungskanals beobachtet. Die Position x wird über die Strömungsgeschwindigkeit u in eine Reaktionszeit übersetzt:

t=x/u

Sehr kurze Zeiten sind erreichbar, allerdings bei hohem Probenverbrauch.

9. Quench-Flow

Rapid Quench-Flow startet die Reaktion durch Mischen und stoppt sie nach einer definierten Zeit chemisch, thermisch oder durch Lösungsmittelwechsel. Die Produkte werden anschließend analytisch bestimmt.

Die Methode eignet sich besonders, wenn kein kontinuierliches optisches Signal verfügbar ist.

10. Pseudo-erste Ordnung

Für eine bimolekulare Reaktion:

A+B→P

mit:

v=k₂[A][B]

wird B im großen Überschuss eingesetzt. Dann bleibt [B] näherungsweise konstant:

k_obs=k₂[B]₀

d[A]/dt=−k_obs[A]

Die Auswertung wird exponentiell, und mehrere [B]₀-Werte erlauben die Bestimmung von k₂ aus der Steigung von k_obs gegen [B]₀.

11. Reaktion mit zusätzlichem unimolekularem Weg

Besitzt A zusätzlich einen unimolekularen Zerfall mit k₁:

k_obs=k₂[B]₀+k₁

Die Gerade besitzt dann einen von null verschiedenen Achsenabschnitt.

12. Flash-Photolyse

Bei der Flash-Photolyse erzeugt ein kurzer Lichtpuls eine nichtthermische Population:

A+hν→A*

A* kann angeregter Zustand, Radikal, Ion, Isomer oder Photoprodukt sein. Ein verzögerter Analysepuls oder kontinuierliches Messlicht verfolgt das zeitabhängige Spektrum.

13. Transiente Absorption

Die Änderung der Absorbanz lautet:

ΔA(λ,t)=A(λ,t)−A_{vor Puls}(λ)

Positive Signale können neue Absorber anzeigen; negative Signale können Grundzustandsbleaching oder stimulierte Emission enthalten.

Die spektrale Zuordnung erfordert daher mehr als die Betrachtung einer einzigen Wellenlänge.

14. Pump-Probe-Prinzip

Ein Pump-Puls startet die Dynamik. Ein zeitlich verzögerter Probe-Puls misst den Zustand. Durch Variation der Verzögerung wird die Transiente rekonstruiert.

Die Zeitauflösung wird durch Pulsdauer, Synchronisation, Dispersion und Instrumentenantwort begrenzt.

15. Temperatur-Sprung

Ein Temperatur-Sprung verschiebt Gleichgewichtskonstanten und Geschwindigkeitskonstanten abrupt. Das System besitzt unmittelbar nach dem Sprung noch die alte Zusammensetzung, relaxiert aber unter den neuen Raten zum neuen Gleichgewicht.

Historische T-Jump-Geräte nutzten elektrische Entladung in leitfähigen Lösungen. Moderne Laser-T-Jump-Experimente absorbieren einen kurzen IR-Puls im Lösungsmittel.

16. Laser-T-Jump

Ein Laser-T-Jump kann die Temperatur innerhalb von Nanosekunden erhöhen. Die Temperaturamplitude wird über Lösungsmittelabsorption, Pulsenergie, Wärmeleitung und Zellgeometrie bestimmt.

Geeignete Anwendungen sind:

• Protein- und RNA-Faltung,
• Ligandenbindung,
• Protonentransfer,
• Konformationsgleichgewichte,
• supramolekulare Assoziation.

17. Drucksprung

Ein Drucksprung verschiebt Gleichgewichte, wenn Reaktand und Produkt unterschiedliche partielle Molvolumina besitzen. Für eine Gleichgewichtskonstante gilt:

(∂lnK/∂p)_T=−ΔV°/(RT)

Die Relaxationszeit nach dem Drucksprung liefert kinetische Information, während die Amplitude mit dem Gleichgewichtsvolumen verknüpft ist.

18. Feldsprung

Ein elektrischer Feldsprung verändert die Orientierung polarer Moleküle, Ionenverteilungen oder Dissoziationsgleichgewichte. Nach Abschalten oder Ändern des Feldes relaxiert die Polarisation oder Konzentration.

Dielektrische Relaxation und elektrische Doppeltschichtprozesse sind typische Anwendungen.

19. Konzentrationssprung und Caged Compounds

Eine Konzentration kann durch schnelles Mischen oder photochemische Freisetzung verändert werden. Caged Compounds enthalten eine inaktive, photolabile Schutzgruppe. Ein Lichtpuls setzt beispielsweise ATP, Ca²⁺, Protonen oder Neurotransmitter frei.

Der photochemische Freisetzungsschritt muss schneller als die nachfolgende Reaktion und quantitativ charakterisiert sein.

20. Messgrößen

Schnelle Kinetik kann über viele Observablen verfolgt werden:

• Absorbanz und transiente Absorption,
• Fluoreszenzintensität und -lebensdauer,
• Fluoreszenzanisotropie,
• Leitfähigkeit und Impedanz,
• Brechungsindex und Lichtstreuung,
• IR-, Raman- und EPR-Signale,
• Temperatur, Druck oder Wärmefluss.

Die gemessene Observable ist nicht automatisch identisch mit der Konzentration einer einzelnen Spezies.

Abbildung 1. Schnelles Mischen startet eine Reaktion durch Kontakt der Reaktanden. Sprung- und Pulsmethoden verändern dagegen abrupt Lichtpopulation, Temperatur, Druck, Feld oder Konzentration.

21. Beer-Lambert-Signal

Für nichtwechselwirkende absorbierende Spezies gilt:

A(λ,t)=lΣ_iε_i(λ)c_i(t)

Das Signal ist eine lineare Kombination der Konzentrationen. Ein einzelner exponentieller Konzentrationsverlauf kann je nach Extinktionskoeffizienten als steigendes, fallendes oder nahezu verschwindendes optisches Signal erscheinen.

22. Reversibles Zweizustandssystem

Für:

A ⇌ B

mit Vorwärtskonstante k₁ und Rückwärtskonstante k₋₁ gelten:

d[A]/dt=−k₁[A]+k₋₁[B]

d[B]/dt=k₁[A]−k₋₁[B]

Die Gesamtmenge c_tot=[A]+[B] bleibt erhalten.

23. Gleichgewicht

Im Gleichgewicht ist der Nettostrom null:

k₁[A]_eq=k₋₁[B]_eq

Daraus folgt:

K=[B]_eq/[A]_eq=k₁/k₋₁

Mit c_tot:

[A]_eq=c_tot/(1+K)

[B]_eq=Kc_tot/(1+K)

24. Relaxation nach einer kleinen Störung

Nach einem kleinen Sprung wird geschrieben:

[A](t)=[A]_eq+δA(t)

[B](t)=[B]_eq+δB(t)

Wegen der Stoffmengenerhaltung gilt δB=−δA. Einsetzen liefert:

dδA/dt=−(k₁+k₋₁)δA

25. Relaxationszeit

Die Lösung lautet:

δA(t)=δA(0)e^−t/τ

mit:

τ⁻¹=k₁+k₋₁

Die Relaxation misst also nicht nur die Vorwärts- oder Rückwärtsrate, sondern ihre Summe.

26. Bestimmung beider Einzelkonstanten

Aus:

K=k₁/k₋₁

und:

τ⁻¹=k₁+k₋₁

folgt:

k₁=[K/(1+K)]τ⁻¹

k₋₁=[1/(1+K)]τ⁻¹

Eine Gleichgewichts- und eine Relaxationsmessung reichen im idealen Zweizustandsfall aus.

27. Relaxationsamplitude

Die Relaxationszeit beschreibt die Geschwindigkeit, die Amplitude dagegen die Verschiebung des Gleichgewichts und die Empfindlichkeit der Messgröße.

Eine kleine Amplitude kann entstehen, obwohl die Reaktion schnell abläuft, wenn:

• der Sprung das Gleichgewicht nur schwach verschiebt,
• die beobachteten Spezies ähnliche Spektren besitzen,
• die Konzentration gering ist,
• mehrere Signalanteile einander aufheben.

28. Kleine-Störungs-Näherung

Für ein allgemeines kinetisches System:

dc/dt=f(c)

wird um den Gleichgewichtspunkt c_eq entwickelt:

c=c_eq+δc

dδc/dt≈Jδc

J ist die Jacobi-Matrix:

J_ij=∂f_i/∂c_j|_eq

29. Relaxationseigenwerte

Die Eigenwertgleichung lautet:

Jv_i=λ_iv_i

Für ein stabiles Gleichgewicht besitzen die relevanten Eigenwerte negative Realteile. Jede Mode relaxiert wie:

a_i(t)=a_i(0)e^λ_it

Bei reellen negativen Eigenwerten ist:

τ_i=−1/λ_i

30. Relaxationseigenmoden

Ein Eigenvektor beschreibt eine bestimmte gekoppelte Konzentrationsänderung. Eine Mode kann beispielsweise A und B gegensinnig verändern, während eine andere B und C koppelt.

Die Eigenmoden sind keine einzelnen chemischen Spezies. Sie sind mathematische Kombinationen von Konzentrationsabweichungen.

31. Erhaltungsgesetze und Nullmoden

Stoffmengenerhaltung kann einen formalen Eigenwert null erzeugen, wenn alle Konzentrationen ohne Eliminierung einer abhängigen Variable betrachtet werden. Diese Mode beschreibt keine Relaxation, sondern die erhaltene Gesamtmenge.

Für die kinetische Auswertung sollte das Modell auf unabhängige Variablen reduziert oder die Nullmode korrekt behandelt werden.

32. Komplexe Eigenwerte

Komplex konjugierte Eigenwerte erzeugen gedämpfte Oszillationen:

δc(t)∝e^αtcos(βt+φ)

Solche Moden sind in gewöhnlichen detailliert balancierten chemischen Gleichgewichtssystemen selten, können aber in offenen Reaktionsnetzwerken, Rückkopplungssystemen und chemischen Oszillatoren auftreten.

33. Mehrere Exponentialanteile

Ein Signal kann als Summe geschrieben werden:

S(t)=S_∞+Σ_iA_ie^−t/τ_i

Mehrere Exponentialterme können aus:

• mehreren Relaxationseigenmoden,
• parallelen Populationen,
• räumlicher Heterogenität,
• Instrumentenantwort,
• Verteilungen von Raten

entstehen.

34. Amplituden sind beobachtungsabhängig

Dieselben kinetischen Eigenwerte erscheinen an verschiedenen Wellenlängen, aber ihre Amplituden können Vorzeichen und Größe wechseln. Ein Zwischenprodukt kann bei einer Wellenlänge stark, bei einer anderen nahezu unsichtbar sein.

Deshalb ist ein einzelner Zeittrace häufig mechanistisch mehrdeutig.

Abbildung 2. Im Zweizustandssystem ist die Relaxationsrate die Summe von Vorwärts- und Rückwärtskonstante. In größeren Netzwerken liefern die Eigenwerte der Jacobi-Matrix mehrere gekoppelte Zeitkonstanten.

35. Sequentielle Reaktion

Für:

A → B → C

besitzt B typischerweise einen Anstieg und anschließenden Abfall. Bei Konstanten k₁ und k₂ lautet für [B](0)=0:

[B](t)=[A]₀k₁[e^−k₁t−e^−k₂t]/(k₂−k₁)

Ein reiner monoexponentieller Fit kann einen solchen Zwischenproduktverlauf nicht angemessen beschreiben.

36. Parallele Reaktionen

Für:

A → P₁    mit k₁

A → P₂    mit k₂

zerfällt A mit:

[A](t)=[A]₀e^{−(k₁+k₂)t}

Der Zerfall von A ist monoexponentiell, obwohl zwei Produktwege vorliegen. Das Produktverhältnis liefert zusätzliche Information:

[P₁]_∞/[P₂]_∞=k₁/k₂

37. Große Sprünge und numerische Integration

Bei großen Störungen ist die Jacobi-Matrix nicht konstant. Dann werden die vollständigen Differentialgleichungen integriert:

dc/dt=f(c,k)

Das berechnete Signal wird mit dem Experiment verglichen. Steife Systeme benötigen geeignete implizite Integratoren.

38. Globale kinetische Anpassung

Bei einer globalen Anpassung werden viele Zeittraces oder Wellenlängen gemeinsam ausgewertet. Die kinetischen Parameter gelten für alle Datensätze, während spektrale Amplituden oder Anfangsbedingungen variieren dürfen.

Für spektrale Daten:

S(λ,t)=Σ_ic_i(t)ε_i(λ)

Die gemeinsame Struktur erhöht die statistische Information und kann mechanistische Spezies trennen.

39. Lifetime Analysis

Bei einer modellfreien oder phänomenologischen Lifetime Analysis werden Daten als Summe exponentieller Zeitkonstanten beschrieben, ohne diese direkt einem Mechanismus zuzuordnen.

Das Ergebnis liefert decay-associated spectra oder amplitudenassoziierte Spektren. Diese Spektren sind mathematische Komponenten, nicht zwangsläufig Spektren reiner Zwischenprodukte.

40. Target Analysis

Target Analysis verwendet ein konkretes Reaktionsnetzwerk. Die Konzentrationsprofile werden aus den Geschwindigkeitsgleichungen berechnet, und den Spezies werden Spektren zugeordnet.

Sie ist mechanistisch aussagekräftiger, aber auch stärker von Modellannahmen abhängig.

41. Singularwertzerlegung

Eine Datenmatrix D aus Wellenlängen und Zeiten kann zerlegt werden:

D=UΣV^T

Große Singulärwerte zeigen dominante unabhängige Datenkomponenten. SVD hilft bei:

• Rauschunterdrückung,
• Abschätzung der Datenrangzahl,
• Erkennung spektraler und zeitlicher Strukturen,
• Initialisierung globaler Fits.

Die Zahl signifikanter Singulärwerte ist jedoch nicht automatisch gleich der Zahl chemischer Spezies.

42. Instrumentenantwort

Ein reales Instrument misst nicht das ideale Signal S_true(t), sondern dessen Faltung mit der Instrument Response Function – IRF:

S_mess(t)=∫IRF(t−t′)S_true(t′)dt′

Bei ultrakurzen Pulsen umfasst die IRF Pump- und Probe-Pulsdauer, Timing-Jitter und Dispersion. Beim Stopped-Flow umfasst sie Misch- und Transportverteilung.

43. Faltung statt naive Korrektur

Eine robuste Auswertung faltet das kinetische Modell mit der gemessenen oder parametrisierten IRF und passt das gefaltete Signal direkt an die Daten an.

Eine numerische Entfaltung kann Rauschen stark verstärken und ist häufig instabil.

44. Dead-Time-Verlust

Für einen Zerfall erster Ordnung ist der verbleibende Anteil am ersten Messpunkt:

f_rest=e^{−k t_dead}

Ist kt_dead groß, geht der schnelle Anteil vor der Messung verloren. Ein Fit der sichtbaren Restkurve kann die ursprüngliche Amplitude nicht ohne Zusatzinformation rekonstruieren.

45. Zeitnullpunkt

Der experimentelle Zeitnullpunkt kann unsicher sein. Besonders bei kurzen Zeitkonstanten korreliert eine Verschiebung t₀ stark mit Rate und Amplitude.

t₀ sollte durch unabhängige Kalibration oder als gemeinsamer globaler Parameter behandelt werden.

46. Basislinie und Drift

Langsame Instrumentendrift, Photobleaching oder Zellbewegung können Basislinien verändern. Ein Fit kann enthalten:

S(t)=S_kin(t)+b₀+b₁t

Zu flexible Basislinien können jedoch echte langsame Kinetik absorbieren.

47. Rauschmodelle

Least-Squares-Anpassung setzt häufig unabhängiges, normalverteiltes Rauschen mit konstanter Varianz voraus. Reale Daten können besitzen:

• zählstatistisches Rauschen,
• wellenlängenabhängige Varianz,
• zeitliche Korrelation,
• Ausreißer,
• Laserpuls-zu-Puls-Schwankungen.

Gewichtung und Likelihood sollten zum Messprozess passen.

48. Residuenanalyse

Residuen sind:

r_i=y_i−ŷ_i

Ein gutes Modell erzeugt Residuen ohne systematische Zeit-, Wellenlängen- oder Amplitudenstruktur. Typische Warnsignale sind:

• abwechselnde positive und negative Bögen,
• ein nichtmodellierter schneller Peak,
• langsame Drift,
• wellenlängenspezifische Muster,
• Varianz, die mit der Signalhöhe wächst.

49. Parameterkorrelation

Zwei Parameter sind stark korreliert, wenn verschiedene Kombinationen nahezu dasselbe Signal erzeugen. Beispiele:

• Amplitude und Konzentration,
• Dead Time und schnelle Rate,
• zwei ähnliche Zeitkonstanten,
• Vorwärts- und Rückwärtsrate ohne Gleichgewichtsinformation.

Ein kleiner Fitfehler garantiert daher keine präzisen Einzelparameter.

50. Strukturelle Identifizierbarkeit

Ein Modell ist strukturell identifizierbar, wenn ideale rauschfreie Daten die Parameter eindeutig bestimmen könnten.

Ein nichtidentifizierbares Modell besitzt Parameterkombinationen, die exakt dieselbe Observable erzeugen. Mehr Datenpunkte desselben Typs lösen dieses Problem nicht.

51. Praktische Identifizierbarkeit

Ein strukturell identifizierbarer Parameter kann praktisch schlecht bestimmt sein, wenn:

• das Rauschen groß ist,
• der Zeitbereich ungeeignet ist,
• die Komponente kleine Amplitude besitzt,
• Zeitkonstanten zu ähnlich sind,
• die IRF zu breit ist.

52. Konfidenzintervalle

Unsicherheiten können bestimmt werden durch:

• Kovarianzmatrix in lokal linearer Näherung,
• Profil-Likelihood,
• Bootstrap-Resampling,
• Bayes'sche Posteriorverteilungen,
• Monte-Carlo-Simulation synthetischer Datensätze.

Bei stark korrelierten oder begrenzten Parametern sind symmetrische Standardfehler oft unzureichend.

53. Modellvergleich

Ein komplexeres Modell passt Daten fast immer besser. Kriterien wie AIC oder BIC bestrafen zusätzliche Parameter:

AIC=2p−2lnL

BIC=p ln n−2lnL

Statistische Kriterien ersetzen keine chemische Plausibilität, Massenbilanz oder unabhängige Spektralinformation.

54. Temperaturserien

Raten über mehrere Temperaturen können gemeinsam mit Arrhenius- oder Eyring-Beziehungen angepasst werden:

k(T)=A e^−E_a/(RT)

k(T)=(k_BT/h)e^ΔS‡/Re^{−ΔH‡/(RT)}

Eine globale Temperaturauswertung reduziert die Zahl unabhängiger Parameter und testet die Konsistenz eines Mechanismus.

55. Druckserien

Die Druckabhängigkeit einer Geschwindigkeitskonstante ist mit dem Aktivierungsvolumen verknüpft:

(∂lnk/∂p)_T=−ΔV^‡/(RT)

Ein negatives ΔV^‡ führt zu einer Beschleunigung mit Druck. Drucksprung und statische Druckserien ergänzen einander.

56. Rapid Freeze-Quench

Rapid Freeze-Quench mischt Reaktanden, lässt sie für eine definierte kurze Zeit reagieren und friert die Probe rasch ein. Die Zwischenzustände werden anschließend beispielsweise mit EPR, Mößbauer-, Raman- oder Röntgenspektroskopie untersucht.

Die Zeitauflösung wird durch Mischen, Flugstrecke und Gefrierprozess bestimmt.

57. Typische Artefakte schneller Methoden

Methode	typische Artefakte	Kontrolle
Stopped-Flow	unvollständiges Mischen, Blasen, Schlieren, Zellverschmutzung, Dead-Time-Verlust	Farbstoff-Mischtest, Leerexperimente, Fluss- und Viskositätsvariation
Flash-Photolyse	Mehrphotonenprozesse, Photobleaching, Streulicht, Pulsenergievariation	Pulsenergie-Reihe, Wiederholungsrate, Spektralkontrolle
Laser-T-Jump	räumlich inhomogene Erwärmung, akustische Wellen, Brechungsindexsignal	Temperaturkalibration, Referenzprobe, Pulsenergievariation
Drucksprung	mechanische Schwingung, Kompressionserwärmung, Kavitation	druckunempfindliche Referenz und Temperaturkontrolle
Leitfähigkeit	Elektrodenpolarisation und Joule-Heizung	Frequenzvariation, geringe Messleistung und Blindprobe

58. Häufige Auswertungsfehler

• beliebig viele Exponentialterme ohne mechanistische Prüfung,
• Ignorieren der Instrumentenantwort,
• Fit ab t=0, obwohl t₀ unbekannt ist,
• ungewichtete Fits bei stark heteroskedastischem Rauschen,
• Überinterpretation decay-associated spectra als Reinspektren,
• fehlende Prüfung von Parameterkorrelation und Identifizierbarkeit,
• nur visuelle Beurteilung ohne Residuenanalyse,
• fehlende Wiederholungsmessungen und Kalibrationen.

Abbildung 3. Das gemessene Signal ist eine instrumentell verbreiterte Observable. Globale Analyse nutzt gemeinsame kinetische Parameter, während Residuen, SVD und Unsicherheiten die Modellqualität prüfen.

59. Übungsaufgaben

Aufgabe 1: Reversible Zweizustandsrelaxation

Für:

A ⇌ B

gelten k₁=120 s⁻¹ und k₋₁=30,0 s⁻¹. Berechne:

1. die Relaxationszeit τ,
2. die Gleichgewichtskonstante K=[B]/[A],
3. die Gleichgewichtsanteile von A und B.

Aufgabe 2: Exponentielle Rückkehr

Nach einem Sprung beträgt die anfängliche Konzentrationsabweichung Δc₀=0,0800 mmol L⁻¹. Die Relaxationszeit ist τ=2,50 ms. Berechne Δc nach 2,50 ms, 5,00 ms und 10,0 ms.

Aufgabe 3: Signal nach einem Temperatur-Sprung

Die Absorbanz eines Systems beträgt unmittelbar nach einem T-Sprung noch A(0)=0,200. Das neue Gleichgewicht besitzt A_∞=0,260. Die Relaxationszeit ist 80,0 μs. Berechne die Absorbanz nach 100 μs mit:

A(t)=A_∞+[A(0)−A_∞]e^−t/τ

Aufgabe 4: Dead Time eines Stopped-Flow-Geräts

Zwischen Mischer und Beobachtungsvolumen befinden sich 12,0 μL. Der gesamte Volumenstrom beträgt 6,00 mL s⁻¹. Schätze die Dead Time.

Aufgabe 5: Pseudo-erste Ordnung

A reagiert mit B nach zweiter Ordnung mit k₂=2,50×10⁶ L mol⁻¹ s⁻¹. B liegt mit 0,0200 mol L⁻¹ im großen Überschuss vor. Zusätzlich zerfällt A unimolekular mit k₁=100 s⁻¹.

Berechne k_obs und die Halbwertszeit von A.

Aufgabe 6: Flash-Photolyse

Eine transiente Spezies erzeugt unmittelbar nach dem Laserpuls eine Absorbanzänderung ΔA=0,240. Ihr Extinktionskoeffizient beträgt ε=1,20×10⁴ L mol⁻¹ cm⁻¹, die Schichtlänge l=1,00 cm.

Berechne die erzeugte Konzentration. Die Spezies zerfällt anschließend erster Ordnung mit k=3,50×10³ s⁻¹. Berechne ihre Halbwertszeit.

Aufgabe 7: Biexponentielles Signal

Ein normiertes Signal lautet, wenn t in Millisekunden eingesetzt wird:

S(t)=0,700e^{−t/(2,00 ms)}+0,300e^{−t/(20,0 ms)}

Berechne S nach 5,00 ms. Welcher Anteil dominiert zu diesem Zeitpunkt?

Aufgabe 8: Eigenmoden eines gekoppelten Systems

Für die Abweichungen δA und δB gilt:

d/dt [δA, δB]^T = J[δA,δB]^T

mit:

J = [ [−40, 10], [40, −15] ] s⁻¹

Bestimme die beiden Eigenwerte und Relaxationszeiten.

Aufgabe 9: Apparente Aktivierungsenergie

Eine Relaxationsrate beträgt bei 290 K 1,00×10³ s⁻¹ und bei 310 K 2,50×10³ s⁻¹. Berechne mit einer Zwei-Punkt-Arrhenius-Auswertung:

E_a=R ln(k₂/k₁)/(1/T₁−1/T₂)

Aufgabe 10: Signalverlust während der Dead Time

Eine Reaktion zerfällt erster Ordnung mit k=400 s⁻¹. Das Stopped-Flow-Gerät besitzt t_dead=1,50 ms. Welcher Anteil der anfänglichen Spezies ist am ersten beobachtbaren Zeitpunkt noch vorhanden, und welcher Anteil ist bereits umgesetzt?

60. Vollständige Lösungen

Lösung 1

Die Relaxationsrate ist:

τ⁻¹=k₁+k₋₁=120+30,0=150 s⁻¹

τ=1/150=6,667×10⁻³ s=6,67 ms

Die Gleichgewichtskonstante:

K=k₁/k₋₁=120/30,0=4,00

Die Anteile bei erhaltener Gesamtmenge sind:

x_A=1/(1+K)=0,200

x_B=K/(1+K)=0,800

Im Gleichgewicht liegen somit 20 Prozent als A und 80 Prozent als B vor.

Lösung 2

Δc(t)=0,0800e^{−t/(2,50 ms)} mmol L⁻¹

Nach einer Relaxationszeit:

Δc(2,50 ms)=0,0800e⁻¹=0,02943 mmol L⁻¹

Nach zwei Relaxationszeiten:

Δc(5,00 ms)=0,0800e⁻²=0,01083 mmol L⁻¹

Nach vier Relaxationszeiten:

Δc(10,0 ms)=0,0800e⁻⁴=0,001465 mmol L⁻¹

Nach vier Zeitkonstanten sind nur noch etwa 1,83 Prozent der Anfangsabweichung vorhanden.

Lösung 3

A(100 μs)=0,260+(0,200−0,260)e^−100/80

A(100 μs)=0,260−0,0600e^−1,25

A(100 μs)≈0,24281

Das Signal hat sich dem neuen Gleichgewicht bereits deutlich angenähert, aber etwa 28,7 Prozent der anfänglichen Signalabweichung verbleiben.

Lösung 4

t_dead≈V/Q

t_dead=12,0 μL/(6000 μL s⁻¹)

t_dead=2,00×10⁻³ s=2,00 ms

Diese Abschätzung berücksichtigt noch keine endliche Mischzeit oder Verteilung der Verweilzeiten.

Lösung 5

Die pseudo-erste Ordnung ergibt:

k_obs=k₂[B]₀+k₁

k_obs=(2,50×10⁶)(0,0200)+100

k_obs=5,01×10⁴ s⁻¹

Halbwertszeit:

t_1/2=ln2/k_obs

t_1/2≈1,384×10⁻⁵ s=13,84 μs

Diese Reaktion ist für ein gewöhnliches Millisekunden-Stopped-Flow zu schnell und verlangt eine schnellere Misch- oder Sprungmethode.

Lösung 6

Beer-Lambert:

c=ΔA/(εl)

c=0,240/[(1,20×10⁴)(1,00)]

c=2,00×10⁻⁵ mol L⁻¹=20,0 μmol L⁻¹

Halbwertszeit:

t_1/2=ln2/(3,50×10³)

t_1/2=1,980×10⁻⁴ s=0,198 ms

Lösung 7

S(5 ms)=0,700e^−5/2+0,300e^−5/20

S(5 ms)=0,700e^−2,5+0,300e^−0,25

S(5 ms)≈0,05746+0,23364=0,29110

Obwohl die schnelle Komponente anfangs die größere Amplitude besitzt, ist sie nach 5 ms weitgehend abgeklungen. Zu diesem Zeitpunkt dominiert die langsame 20-ms-Komponente.

Lösung 8

Die charakteristische Gleichung lautet:

det(J−λI)=0

(−40−λ)(−15−λ)−400=0

λ²+55λ+200=0

Damit:

λ₁≈−51,085 s⁻¹

λ₂≈−3,915 s⁻¹

Relaxationszeiten:

τ₁=−1/λ₁≈0,01958 s=19,6 ms

τ₂=−1/λ₂≈0,2554 s

Das System besitzt eine schnelle und eine etwa dreizehnmal langsamere gekoppelte Relaxationsmode.

Lösung 9

E_a=R ln(2,50)/(1/290−1/310)

E_a≈3,4245×10⁴ J mol⁻¹

E_a≈34,25 kJ mol⁻¹

Bei einer Relaxationsrate kann dies eine apparente Aktivierungsenergie der Kombination mehrerer Vorwärts- und Rückwärtsprozesse sein. Sie muss nicht der Barriere eines einzelnen Elementarschritts entsprechen.

Lösung 10

Verbleibender Anteil:

f=e^−kt_dead

f=e^{−400·0,00150}=e^−0,600

f≈0,5488

Es sind am ersten Messpunkt noch 54,88 Prozent vorhanden. Umgesetzt wurden:

1−f≈0,4512

Also sind bereits 45,12 Prozent der Anfangspopulation während der Dead Time verschwunden.

61. Zusammenfassung

• Schnelle Reaktionen werden durch abruptes Mischen, Pulsanregung oder Gleichgewichtssprünge untersucht.
• Stopped-Flow deckt typischerweise Millisekunden bis Sekunden ab.
• Dead Time umfasst Misch- und Transportverzögerung.
• Continuous-Flow übersetzt räumliche Position in Reaktionszeit.
• Quench-Flow friert die Zusammensetzung nach einer definierten Zeit chemisch oder thermisch ein.
• Pseudo-erste Ordnung vereinfacht bimolekulare Kinetik bei großem Überschuss eines Reaktanden.
• Flash-Photolyse erzeugt transiente angeregte Zustände, Radikale oder Ionen.
• Pump-Probe rekonstruiert Dynamik aus variabler Pulsverzögerung.
• Laser-T-Jump verschiebt ein Gleichgewicht innerhalb sehr kurzer Zeit.
• Drucksprung koppelt Gleichgewicht und Kinetik an Volumenänderungen.
• Relaxation untersucht die Rückkehr zum Gleichgewicht.
• Für A⇌B gilt τ⁻¹=k₁+k₋₁.
• Gleichgewichtskonstante und Relaxationszeit bestimmen beide Einzelraten.
• Kleine Störungen werden mit der Jacobi-Matrix linearisiert.
• Ihre Eigenwerte liefern Relaxationsraten, die Eigenvektoren gekoppelte Konzentrationsmoden.
• Eine Relaxationsmode ist nicht mit einer einzelnen Spezies gleichzusetzen.
• Mehrere Exponentialterme können unterschiedliche mechanistische oder instrumentelle Ursachen besitzen.
• Globale Analyse verwendet gemeinsame kinetische Parameter für viele Datensätze.
• Lifetime Analysis ist phänomenologisch, Target Analysis mechanistisch.
• SVD schätzt den effektiven Datenrang, aber nicht automatisch die Spezieszahl.
• Das gemessene Signal ist die Faltung des wahren Signals mit der Instrumentenantwort.
• Dead Time und Zeitnullpunkt korrelieren besonders mit schnellen Raten.
• Residuen müssen strukturfrei und mit dem Rauschmodell vereinbar sein.
• Parameterkorrelation und Identifizierbarkeit begrenzen mechanistische Aussagen.
• AIC und BIC unterstützen Modellvergleich, ersetzen aber keine Chemie.
• Temperatur- und Druckserien liefern Aktivierungsparameter.
• Rapid Freeze-Quench konserviert kurzlebige Zwischenzustände für spätere Spektroskopie.

62. Häufig gestellte Fragen

Wie schnell ist Stopped-Flow?

Moderne Geräte besitzen häufig Dead Times im Bereich von etwa einer Millisekunde oder darunter. Der genaue Wert hängt von Mischer, Flussrate, Viskosität und Beobachtungsgeometrie ab.

Warum misst eine Relaxation die Summe von Vorwärts- und Rückwärtsrate?

Eine Abweichung vom Gleichgewicht wird gleichzeitig durch den Vorwärts- und Rückwärtsprozess abgebaut. Beide tragen zur Wiederherstellung des Gleichgewichts bei.

Ist jede exponentielle Kurve eine Reaktion erster Ordnung?

Nein. Ein exponentieller Relaxationsmodus kann aus einem gekoppelten linearen Netzwerk entstehen. Auch parallele Reaktionen können einen monoexponentiellen Reaktandenzerfall erzeugen.

Was ist der Unterschied zwischen Dead Time und Zeitauflösung?

Dead Time ist die früheste beobachtbare Zeit nach Initiation. Zeitauflösung beschreibt zusätzlich, wie stark schnelle Änderungen durch Misch- oder Pulsbreite verschmiert werden.

Warum reicht ein Fit mit mehreren Exponentialtermen nicht?

Viele mathematische Zerlegungen können ähnlich gut passen. Ohne globale Spektren, Konzentrationsvariation oder mechanistische Nebeninformation bleiben die Komponenten häufig mehrdeutig.

Was zeigt die Singularwertzerlegung?

Sie zeigt, wie viele unabhängige mathematische Komponenten oberhalb des Rauschens in der Datenmatrix vorhanden sind. Diese Zahl ist nicht zwingend identisch mit der Zahl chemischer Spezies.

Warum wird das Modell mit der IRF gefaltet?

Das Instrument verbreitert das ideale Signal. Direkte Modellfaltung ist stabiler als eine nachträgliche Entfaltung verrauschter Daten.

Was bedeutet ein strukturierter Residuenverlauf?

Das Modell, die Instrumentenantwort, die Basislinie oder das Rauschmodell beschreibt einen systematischen Datenanteil nicht korrekt.

Kann man aus einer Relaxationszeit alle Einzelraten bestimmen?

Im Allgemeinen nicht. Beim idealen Zweizustandssystem liefert τ nur die Summe. Zusätzliche Gleichgewichts- oder Konzentrationsinformation ist erforderlich.

Wann ist eine große Störung problematisch?

Wenn die lineare Näherung nicht mehr gilt. Dann ändern sich effektive Raten während der Relaxation, und die vollständigen nichtlinearen Differentialgleichungen müssen angepasst werden.

63. Ausblick auf Teil 9

Teil 9 behandelt photochemische Kinetik und angeregte Zustände. Photonenabsorption, Quantenausbeuten, Konkurrenz von Fluoreszenz, innerer Konversion und Intersystem Crossing, Stern–Volmer-Kinetik, Energie- und Elektronentransfer sowie photostationäre Zustände werden kinetisch entwickelt.

Literatur und weiterführende Quellen

1. P. Atkins, J. de Paula & J. Keeler: Atkins’ Physical Chemistry. Oxford University Press.
2. M. Eigen & L. De Maeyer: Relaxation Methods. In: Technique of Organic Chemistry.
3. J. F. K. Ewing: Fast Reactions in Solution. Wiley.
4. I. H. Segel: Enzyme Kinetics. Wiley.
5. I. H. van Stokkum, D. S. Larsen & R. van Grondelle: Global and Target Analysis of Time-Resolved Spectra. Biochimica et Biophysica Acta.

Didaktischer Hinweis: Relaxationszeiten, exponentielle Komponenten und globale Fits sind modellabhängig. Eine belastbare schnelle Kinetik kombiniert instrumentelle Kalibration, unabhängige Konzentrations- und Gleichgewichtsinformation, mehrere Messbedingungen, physikalisch sinnvolle Rauschmodelle und eine explizite Prüfung der Parameteridentifizierbarkeit.

Bild- und Nutzungsnachweis: Das Titelbild und die Abbildungen 1–3 wurden eigens für diesen Blogartikel mit selbst erstelltem Programmcode erzeugt. Es wurden keine Grafiken aus Lehrbüchern oder anderen Publikationen übernommen. Die Abbildungen sind für die Verwendung in diesem monetarisierbaren Blog vorgesehen.