Samstag, 1. August 2026

Diffusionskontrollierte Reaktionen und Lösungsmittelreibung – Smoluchowski, Collins–Kimball, Käfigeffekt und Kramers-Theorie




Titelbild: Eigene Illustration Brownscher Molekülbewegung, eines Encounter-Komplexes und zentraler Beziehungen der Diffusions- und Reibungskinetik.

In einer Flüssigkeit müssen zwei Reaktionspartner einander zunächst durch Brownsche Bewegung finden, bevor sie chemisch reagieren können. Für sehr schnelle Elementarschritte ist deshalb nicht mehr die Aktivierungsbarriere allein geschwindigkeitsbestimmend. Der molekulare Transport durch das Lösungsmittel setzt eine obere Grenze für bimolekulare Geschwindigkeitskonstanten.

Die einfachste Theorie dieses Grenzfalls stammt von Smoluchowski. Zwei kugelförmige Reaktanden diffundieren relativ zueinander. Sobald ihr Abstand den Reaktionsradius R erreicht, werden sie vollständig absorbiert und bilden Produkt. Daraus folgt:

kdiff=4πR(DA+DB)NA

In realen Systemen reagiert jedoch nicht jede Begegnung. Orientierung, Solvatation, Spin, elektronische Kopplung oder eine verbleibende chemische Barriere können die Reaktion am Kontakt begrenzen. Die Collins–Kimball-Theorie kombiniert deshalb Diffusion und intrinsische Oberflächenreaktivität:

1/kobs=1/kdiff+1/kchem

Außerdem kann das Lösungsmittel die Partner nach einer Begegnung in einem gemeinsamen Käfig festhalten. Wiederholte Nahkontakte, geminate Rekombination und Käfigflucht werden dann Teil des Mechanismus. Bei intramolekularer Barrierenüberquerung wirkt dieselbe Flüssigkeit als Reibungsbad; Kramers-Theorie beschreibt, wie die Reaktionsrate von der Reibung abhängt.

Einordnung: Teil 6 behandelte quantenmechanisches Tunneln und kinetische Isotopeneffekte. Teil 7 ergänzt diese mikroskopische Barrierenphysik um Stofftransport, Lösungsmittelkäfig und Reibung. Teil 8 behandelt anschließend schnelle Reaktionsmethoden und Relaxationskinetik.

Lernziele

Nach diesem Beitrag kannst du:
  • Brownsche Bewegung und Diffusion mikroskopisch erklären,
  • Einstein- und Stokes–Einstein-Beziehung anwenden,
  • relative Diffusion zweier Reaktionspartner bestimmen,
  • den Smoluchowski-Grenzwert für bimolekulare Reaktionen berechnen,
  • Encounter-Komplex, Käfigeffekt und geminate Rekombination unterscheiden,
  • die Collins–Kimball-Gleichung interpretieren,
  • diffusions-, aktivierungs- und gemischt kontrollierte Regime erkennen,
  • intermolekulare Potentiale in der Debye–Smoluchowski-Theorie einordnen,
  • Viskositäts- und Orientierungseffekte beurteilen,
  • den Kramers-Turnover und die hohe-Reibungs-Grenze erklären,
  • Fluoreszenzquenching und Elektrodenkinetik als Diffusionsprobleme analysieren,
  • Grenzen klassischer Diffusionstheorien in komplexen Medien benennen.

1. Warum Flüssigkeitskinetik besonders ist

In der Gasphase bewegen sich Moleküle zwischen Kollisionen annähernd geradlinig. In einer Flüssigkeit erfolgen dagegen fortlaufend Stöße mit Lösungsmittelmolekülen. Die Bewegung eines gelösten Teilchens ist unregelmäßig, stark gedämpft und auf langen Zeitskalen diffusiv.

Eine schnelle bimolekulare Reaktion in Lösung besteht deshalb mindestens aus zwei Schritten:

A+B ⇌ [A···B] → P

Der erste Schritt beschreibt die diffusiv erzeugte Begegnung. Der zweite beschreibt die chemische Umwandlung innerhalb des Encounter-Komplexes.

2. Brownsche Bewegung

Brownsche Bewegung entsteht aus der Summe zahlreicher ungeordneter Stöße des Lösungsmittels mit einem gelösten Teilchen. Auf sehr kurzen Zeitskalen besitzt das Teilchen Trägheit; auf längeren Zeitskalen ist seine Bahn ein Zufallsweg.

Für isotrope dreidimensionale Diffusion gilt für die mittlere quadratische Verschiebung:

⟨r²(t)⟩=6Dt

D ist der Diffusionskoeffizient mit der Einheit m² s−1.

3. Diffusionsgleichung

Die Konzentrationsverteilung c(r,t) eines frei diffundierenden Stoffes erfüllt Ficks zweites Gesetz:

∂c/∂t=D∇²c

Es beschreibt, wie Konzentrationsunterschiede mit der Zeit geglättet werden. Der zugehörige Fluss lautet nach Ficks erstem Gesetz:

J=−D∇c

Der Fluss zeigt in Richtung abnehmender Konzentration.

4. Langevin-Gleichung

Eine mikroskopische Beschreibung kombiniert Trägheit, Reibung, systematische Kraft und zufällige Lösungsmittelstöße:

m dv/dt=−ζv+F(r)+R(t)

ζ ist der Reibungskoeffizient. R(t) ist eine stochastische Kraft mit Mittelwert null. Reibung und Fluktuation sind nicht unabhängig; das Fluktuations-Dissipations-Theorem verknüpft sie so, dass thermisches Gleichgewicht entsteht.

5. Einstein-Beziehung

Für die translatorische Diffusion gilt:

D=kBT/ζ

Stärkere Reibung senkt D, höhere Temperatur erhöht D. Diese Beziehung ist eine Form des Fluktuations-Dissipations-Theorems.

6. Stokes-Reibung

Für eine starre Kugel mit Radius a in einer kontinuierlichen Flüssigkeit mit dynamischer Viskosität η lautet bei haftender Randbedingung:

ζ=6πηa

Einsetzen in die Einstein-Beziehung liefert die Stokes–Einstein-Gleichung:

D=kBT/(6πηa)

7. Physikalische Bedeutung der Stokes–Einstein-Gleichung

  • D steigt linear mit T.
  • D fällt invers mit der Viskosität.
  • D fällt invers mit dem hydrodynamischen Radius.

Die Gleichung ist für Teilchen geeignet, die gegenüber den Lösungsmittelmolekülen hinreichend groß sind. Für kleine Moleküle, strukturierte Lösungsmittel und starke Solvatation ist a ein effektiver hydrodynamischer Radius.

8. Slip- und Stick-Randbedingungen

Die klassische Stokes-Reibung setzt voraus, dass die Flüssigkeit an der Teilchenoberfläche haftet. Bei vollständigem Gleiten wäre der Reibungskoeffizient kleiner. Molekulare Systeme liegen häufig zwischen diesen Grenzfällen.

Die Wahl der hydrodynamischen Randbedingung beeinflusst berechnete Diffusionskoeffizienten und Rotationsdiffusion.

9. Grenzen von Stokes–Einstein

Abweichungen treten auf bei:

  • sehr kleinen gelösten Molekülen,
  • hochviskosen oder unterkühlten Flüssigkeiten,
  • Wasserstoffbrückennetzwerken,
  • Polymerlösungen und Gläsern,
  • stark anisotropen Molekülen,
  • räumlich heterogenen Medien.

In solchen Fällen kann eine fraktionelle Beziehung D∝η−α mit α<1 beobachtet werden.

10. Relative Diffusion

Für die Begegnung zweier unabhängig diffundierender Teilchen ist die Relativkoordinate r=rA−rB entscheidend. Ihr Diffusionskoeffizient ist:

Drel=DA+DB

Beide Bewegungen tragen also dazu bei, dass sich der Abstand ändert.

11. Reaktionsradius

Der Reaktionsradius R ist der Abstand, bei dem eine Begegnung als reaktiv behandelt wird. Für harte Kugeln kann er näherungsweise die Summe der Radien sein:

R≈aA+aB

Bei langreichweitigen Wechselwirkungen, Elektronentransfer oder anisotropen Reaktionen ist R eine effektive Größe und nicht zwingend ein direkter Kontaktabstand.

12. Smoluchowski-Modell

Man fixiert A im Ursprung und lässt B mit Drel diffundieren. Im stationären Zustand erfüllt die radialsymmetrische Konzentration c(r):

0=Drel(1/r²)d/dr[r² dc/dr]

Randbedingungen:

c(R)=0
c(r→∞)=c

c(R)=0 bedeutet vollständige Absorption: Jeder Kontakt reagiert sofort.

13. Stationäre Konzentrationsverteilung

Die Lösung lautet:

c(r)=c(1−R/r)

In der Nähe der absorbierenden Oberfläche ist die Konzentration abgesenkt. Dieses Konzentrationsloch treibt den diffusiven Zustrom.

14. Smoluchowski-Fluss

Der radiale Fluss bei r=R ist:

J(R)=−Drel dc/dr|R

Multiplikation mit der Kugeloberfläche 4πR² ergibt die Zahl der Begegnungen pro Zeit:

Rate=4πR Drelc

Daraus folgt der mikroskopische Volumenquerschnitt pro Zeit:

kpair=4πR Drel

15. Molarer Diffusionsgrenzwert

Für eine molare Geschwindigkeitskonstante wird mit der Avogadro-Konstante multipliziert. Bei SI-Einheiten ergibt sich zunächst m³ mol−1 s−1:

kdiff=4πR DrelNA

Zur Umrechnung in L mol−1 s−1 wird zusätzlich mit 1000 multipliziert.

Typische Werte kleiner Moleküle in niedrigviskosen Flüssigkeiten liegen ungefähr bei 109 bis 1010 L mol−1 s−1.



Abbildung 1. Zwei unabhängig diffundierende Partner können als ein Relativteilchen mit Drel=DA+DB behandelt werden. Eine vollständig absorbierende Kugeloberfläche liefert den Smoluchowski-Grenzwert.

16. Zeitabhängige Diffusionsrate

Direkt nach dem Start einer Reaktion ist noch kein stationäres Konzentrationsprofil aufgebaut. Der zeitabhängige Smoluchowski-Koeffizient besitzt für eine absorbierende Kugel näherungsweise:

k(t)=4πR Drel[1+R/√(πDrelt)]

Bei sehr kurzen Zeiten ist der Anfangsfluss größer, weil sich Reaktanden bereits nahe der Kontaktoberfläche befinden können. Für t→∞ nähert sich k(t) dem stationären Grenzwert.

17. Encounter-Komplex

Ein Encounter-Komplex [A···B] ist kein notwendigerweise gebundenes Molekül. Er bezeichnet eine räumlich nahe Paaranordnung, die durch das Lösungsmittel mehrfach zusammenstoßen kann.

Innerhalb des Encounter-Komplexes können auftreten:

  • chemische Reaktion,
  • Rotation und Neuorientierung,
  • Energie- oder Elektronentransfer,
  • Käfigflucht,
  • Rückkehr zu erneutem Kontakt.

18. Käfigeffekt

Das Lösungsmittel bildet um zwei frisch erzeugte Fragmente einen transienten Käfig. Statt sich sofort weit zu trennen, stoßen sie wiederholt miteinander und mit Nachbarmolekülen zusammen.

Der Käfigeffekt erhöht die Wahrscheinlichkeit geminater Rekombination und kann Produktverteilungen verändern.

19. Geminate Rekombination

Geminate Partner stammen aus demselben Vorläufer oder demselben primären Ereignis. Beispiele sind:

  • Photodissoziation eines Moleküls,
  • Radikalpaarbildung,
  • Ladungstrennung,
  • Bindungshomolyse.

Ihre frühe Rekombination ist keine zufällige Begegnung zweier unabhängiger Teilchen, sondern eine korrelierte Rückreaktion im gemeinsamen Lösungsmittelkäfig.

20. Käfigflucht

Die Partner können den Käfig verlassen und in das Volumen diffundieren. Ein einfaches Konkurrenzmodell lautet:

[A···B] → P    mit kchem
[A···B] → A+B    mit kesc

Die Reaktionswahrscheinlichkeit pro Begegnung ist:

Preakt=kchem/(kchem+kesc)

Die mittlere Lebensdauer des Käfigzustands beträgt:

τcage=1/(kchem+kesc)

21. Teilreaktive Kontaktfläche

Die Smoluchowski-Randbedingung c(R)=0 setzt voraus, dass jeder Kontakt reagiert. Reale Reaktionen benötigen möglicherweise eine bestimmte Orientierung oder müssen am Kontakt noch eine Aktivierungsbarriere überwinden.

Collins und Kimball ersetzten die vollständig absorbierende Oberfläche durch eine Strahlungsrandbedingung:

Drel∂c/∂r|Rsc(R)

κs besitzt die Einheit einer Geschwindigkeit und beschreibt die intrinsische Reaktivität der Kontaktoberfläche.

22. Collins–Kimball-Geschwindigkeitskonstante

Für eine kugelsymmetrische Kontaktreaktion ergibt sich:

kpair=4πR Drel/[1+Drel/(κsR)]

Definiert man eine intrinsische chemische Konstante kchem, kann dieselbe Physik als Widerstandsaddition geschrieben werden:

1/kobs=1/kdiff+1/kchem

Diffusion und chemische Umwandlung wirken wie zwei aufeinanderfolgende kinetische Widerstände.

23. Diffusionskontrolliertes Regime

Für kchem≫kdiff gilt:

kobs≈kdiff

Fast jede Begegnung führt zum Produkt. Eine weitere Absenkung der chemischen Barriere erhöht die beobachtete Rate kaum noch, weil der Transport begrenzt.

24. Aktivierungskontrolliertes Regime

Für kchem≪kdiff gilt:

kobs≈kchem

Die Partner begegnen sich häufig, reagieren aber nur selten. Die Rate wird hauptsächlich durch Orientierung, elektronische Kopplung oder eine lokale Aktivierungsbarriere bestimmt.

25. Gemischte Kontrolle

Wenn kchem und kdiff ähnlich groß sind, beeinflussen beide Prozesse die beobachtete Konstante. Eine praktische dimensionslose Kennzahl ist:

Da=kchem/kdiff
  • Da≫1: Diffusionskontrolle,
  • Da≪1: chemische Kontrolle,
  • Da≈1: gemischte Kontrolle.

26. Wiederholte Begegnungen

Nach einer nichtreaktiven Trennung am Kontakt sind die Partner häufig noch räumlich korreliert. Sie können erneut zusammenstoßen, bevor sie vollständig in das Volumen entkommen.

Eine einzelne „Begegnung“ kann daher aus vielen mikroskopischen Kontaktversuchen bestehen. Die effektive Reaktionswahrscheinlichkeit eines Encounter-Komplexes ist größer als die Wahrscheinlichkeit eines einzelnen Kontakts.

27. Orientierungsgating

Viele Moleküle reagieren nur, wenn bestimmte reaktive Gruppen zueinander ausgerichtet sind. Neben der Translationsdiffusion wird dann Rotationsdiffusion wichtig.

Ein einfaches Modell zerlegt den Prozess:

A+B ⇌ [A···B]falsch ⇌ [A···B]reaktiv → P

Große Proteine oder anisotrope Moleküle können eine hohe Kontaktfrequenz, aber eine kleine orientierungsgewichtete Reaktionsrate besitzen.

28. Sterischer Faktor in Lösung

Ein effektiver sterischer Faktor p kann den Smoluchowski-Grenzwert reduzieren:

k≈p kdiff

Diese Form ist nur eine grobe Näherung. Rotation, wiederholte Kontakte und anisotrope Potentiale koppeln Orientierung und Translation dynamisch.



Abbildung 2. Eine Begegnung kann zur Reaktion oder zur Käfigflucht führen. Die Collins–Kimball-Gleichung verbindet den Diffusionsgrenzwert mit der intrinsischen chemischen Reaktivität.

29. Zwischenmolekulare Potentiale

Reaktionspartner diffundieren nicht immer frei. Elektrostatische, van-der-Waals- oder hydrophobe Wechselwirkungen erzeugen ein Potential of Mean Force U(r), das den Zustrom verändert.

Eine attraktive Wechselwirkung erhöht die Kontaktwahrscheinlichkeit, eine repulsive verringert sie.

30. Debye–Smoluchowski-Gleichung

Für ein radialsymmetrisches Potential U(r) lautet der stationäre diffusionskontrollierte Paar-Koeffizient:

kpair=4πDrel/[∫ReU(r)/(kBT)dr/r²]

Die molare Konstante folgt wieder durch Multiplikation mit NA und gegebenenfalls dem Literfaktor 1000.

Für U(r)=0 ist das Integral 1/R, und die Smoluchowski-Form 4πRDrel wird zurückgewonnen.

31. Elektrostatische Anziehung und Abstoßung

Für zwei punktförmige Ionen:

U(r)=zAzBe²/(4πε0εrr)

Entgegengesetzte Ladungen besitzen U<0 und werden zum Kontakt fokussiert. Gleichnamige Ladungen werden abgestoßen.

Die Stärke hängt stark von der relativen Permittivität εr des Lösungsmittels ab.

32. Bjerrum-Länge

Die Bjerrum-Länge ist der Abstand, bei dem die Coulomb-Energie zweier Elementarladungen dem thermischen Energiebetrag kBT entspricht:

lB=e²/(4πε0εrkBT)

In Wasser bei Raumtemperatur beträgt sie ungefähr 0,7 nm. In weniger polaren Lösungsmitteln ist sie größer, und Ionenwechselwirkungen reichen weiter.

33. Ionenstärke und Abschirmung

Elektrolytionen schirmen Coulomb-Wechselwirkungen ab. In der Debye-Hückel-Näherung fällt das Potential ungefähr exponentiell mit der Debye-Länge.

Steigende Ionenstärke kann daher:

  • die Anziehung entgegengesetzter Reaktanden schwächen,
  • die Abstoßung gleichnamiger Reaktanden verringern,
  • diffusionskontrollierte Ionenraten verändern.

34. Hydrophobe und spezifische Wechselwirkungen

In wässriger Lösung können hydrophobe Gruppen einen attraktiven Potential-of-Mean-Force-Beitrag besitzen. Wasserstoffbrücken oder Metallkoordination können orientierungsspezifische Encounter-Komplexe stabilisieren.

Eine stärkere Begegnungsbindung erhöht jedoch nicht automatisch die Produktbildung. Ein zu tiefes nichtreaktives Minimum kann die Partner in einer falschen Orientierung festhalten.

35. Käfigeffekt bei Radikalen

Nach photochemischer oder thermischer Homolyse entstehen zwei Radikale zunächst benachbart. Der Lösungsmittelkäfig fördert:

  • geminate Rekombination,
  • Disproportionierung,
  • Spinselektivität,
  • mehrfache Rückstöße.

Eine höhere Viskosität kann die Käfigflucht verlangsamen und die geminate Ausbeute erhöhen.

36. Radikalpaar-Spinchemie

Radikalpaare können als Singulett- oder Triplettpaar entstehen. Rekombination ist häufig nur aus einem bestimmten Spinzustand erlaubt.

Hyperfeinwechselwirkungen und Magnetfelder können Singulett-Triplett-Umwandlung beeinflussen. Diffusion und Käfiglebensdauer bestimmen, wie viel Zeit für diese Spindynamik verfügbar ist.

37. Fluoreszenzquenching

Ein angeregtes Fluorophor F* kann bei Begegnung mit einem Quencher Q deaktiviert werden:

F*+Q → F+Q

Für dynamisches Quenching gilt die Stern–Volmer-Beziehung:

I0/I=τ0/τ=1+kqτ0[Q]

kq ist eine bimolekulare Quenching-Konstante. Liegt sie nahe 1010 L mol−1 s−1, ist Diffusionskontrolle plausibel.

38. Dynamisches und statisches Quenching

Merkmaldynamisches Quenchingstatisches Quenching
UrsprungKollision mit angeregtem Fluorophornichtfluoreszierender Grundzustandskomplex
Fluoreszenzintensitätsinktsinkt
LebensdauersinktLebensdauer des verbleibenden freien Fluorophors bleibt idealerweise unverändert
Temperaturwirkungoft stärker bei höherer DiffusionKomplex kann bei höherer Temperatur instabiler werden

39. Sauerstoffquenching

Molekularer Sauerstoff quencht viele Triplett- und Singulettzustände effizient. Die beobachtete Rate kann durch Sauerstoffdiffusion, Spinwechselwirkung und elektronische Energien bestimmt werden.

In hochviskosen Medien oder Polymeren ist der Sauerstofftransport oft langsamer und räumlich heterogen.

40. Diffusionsgrenze in Enzymreaktionen

Für eine einfache Enzym-Substrat-Reaktion:

E+S ⇌ ES → E+P

kcat/KM nähert sich bei sehr effizienter Chemie der diffusiven Assoziationsgrenze. Elektrostatische Lenkung und orientierte Bindung können den einfachen Kugelwert verändern.

Ein hoher Wert allein beweist nicht, dass jeder geometrische Kontakt erfolgreich ist; Proteine besitzen komplexe Bindungstrichter.

41. Elektrochemischer Stofftransport

An einer Elektrode müssen elektroaktive Spezies zur Oberfläche diffundieren. Bei sehr schneller Elektronenübertragung wird der Strom durch Stofftransport begrenzt.

Für eine planare Elektrode wächst die Diffusionsschicht mit:

δ(t)∼√(Dt)

Der Cottrell-Strom nach einem Potentialsprung fällt deshalb wie t−1/2.

42. Mikroelektroden

An kleinen Elektroden entsteht radiale statt rein planarer Diffusion. Dadurch kann sich ein stationärer Stoffstrom einstellen, analog zum Smoluchowski-Zustrom zu einer Kugel.

Mikroelektroden besitzen deshalb hohe massentransportbezogene Stromdichten und schnelle Zeitantwort.

43. Reaktionen auf Membranen

In Membranen diffundieren Moleküle näherungsweise zweidimensional. Zweidimensionale Diffusion besitzt andere Langzeit- und Größenabhängigkeiten als dreidimensionale Diffusion.

Zusätzlich sind Orientierung, Domänen, Membranviskosität und laterale Heterogenität wichtig.

44. Makromolekulare Überfüllung

In Zellen und konzentrierten Polymerlösungen ist der Raum durch viele Objekte belegt. Crowding kann:

  • Diffusion verlangsamen,
  • Encounter-Komplexe durch ausgeschlossenen Raum begünstigen,
  • Konformationen verschieben,
  • anomale Diffusion erzeugen.

Eine einzelne makroskopische Viskosität beschreibt solche Medien häufig unzureichend.

45. Anomale Diffusion

Bei normaler Diffusion gilt ⟨r²⟩∝t. Allgemeiner kann beobachtet werden:

⟨r²(t)⟩∝tα
  • α=1: normale Diffusion,
  • α<1: Subdiffusion,
  • α>1: Superdiffusion.

Subdiffusion kann aus Bindung, Hindernissen, viskoelastischer Erinnerung oder heterogenen Diffusionskoeffizienten entstehen.

46. Makroviskosität und Mikroviskosität

Die mit einem Viskosimeter gemessene Makroviskosität muss nicht der lokalen Reibung entsprechen, die ein molekularer Reaktionskoordinat oder Fluoreszenzrotor erfährt.

Mikroviskosität hängt von Sondengröße, Zeitskala, lokalen Wechselwirkungen und struktureller Heterogenität ab.

47. Kramers-Problem

Kramers betrachtete die Flucht eines Teilchens aus einem Potentialminimum über eine Barriere unter dem Einfluss eines Wärmebades. Eine eindimensionale Langevin-Gleichung lautet:

m q̈=−dV/dq−mγq̇+R(t)

γ ist eine Reibungsfrequenz. Das Bad liefert gleichzeitig zufällige Aktivierung und dissipiert Bewegungsenergie.

48. Sehr kleine Reibung

Bei extrem kleiner Reibung tauscht das System nur langsam Energie mit dem Bad aus. Ein Teilchen kann viele Schwingungen im Reaktandental ausführen, bevor es genügend Energie für die Barriere erhält.

Die Rate wächst in diesem Energiediffusionsregime näherungsweise mit γ.

49. Mittlere Reibung

Mit steigender Reibung wird Energie schneller thermalisiert, und die Rate erreicht ein Maximum. Dieser Bereich bildet den Kramers-Turnover.

Die maximale Rate bleibt unter der idealen Übergangszustandsrate, weil Rekreuzung und Dissipation berücksichtigt werden.

50. Hohe Reibung

Bei starker Reibung wird die Bewegung über die Barriere räumlich diffusionsartig. Das Teilchen wird am Vorankommen gehindert. Die Rate sinkt ungefähr invers mit γ.

Für eine parabolische Barriere ist ein gebräuchlicher Transmissionsfaktor:

κK=[√((γ/2)²+ωb²)−γ/2]/ωb

und:

k=κKkTST

51. Hohe-Reibungs-Grenze

Für γ≫ωb wird:

κK≈ωb

und damit näherungsweise k∝1/γ. Wird γ proportional zur Viskosität angenommen, folgt k∝1/η.

52. Viskositätsabhängigkeit

Experimentell wird oft eine fraktionelle Beziehung gefunden:

k∝η−α

mit 0<α<1. Gründe sind:

  • interne Molekülreibung,
  • unvollständige Kopplung an die makroskopische Viskosität,
  • mehrdimensionale Reaktionskoordinaten,
  • viskoelastische Lösungsmittelantwort,
  • Konformationsgating.

53. Generalisierte Langevin-Gleichung

Reale Lösungsmittel besitzen Erinnerung. Die Reibung hängt dann von der vergangenen Geschwindigkeit ab:

m q̈(t)=−dV/dq−∫0tΓ(t−t′)q̇(t′)dt′+R(t)

Γ(t) ist ein Gedächtniskern. Das zugehörige farbige Rauschen erfüllt eine verallgemeinerte Fluktuations-Dissipations-Beziehung.

54. Grote–Hynes-Theorie

Grote–Hynes-Theorie erweitert Kramers auf frequenzabhängige Reibung. Der Transmissionskoeffizient lautet:

κGHrb

λr ist die positive reaktive Frequenz, die aus der Laplace-transformierten Reibung bestimmt wird. Schnelle und langsame Lösungsmittelmoden koppeln dadurch unterschiedlich an die Barrierenbewegung.

55. Lösungsmittelkoordinate

Bei Elektronen- und Protonentransfer kann das Lösungsmittel eine kollektive Koordinate besitzen. Die Reaktion erfolgt erst, wenn Polarisation und Geometrie eine geeignete Energielücke erzeugen.

Dann kann die Rate sowohl durch die eigentliche quantenmechanische Kopplung als auch durch die diffusive Bewegung der Lösungsmittelkoordinate bestimmt sein.

56. Temperaturabhängigkeit diffusionskontrollierter Raten

Mit Stokes–Einstein und Smoluchowski gilt näherungsweise:

kdiff∝T/η(T)

Die Viskosität besitzt häufig eine starke Temperaturabhängigkeit. Deshalb kann eine diffusionskontrollierte Rate scheinbar Arrhenius-artig sein, obwohl keine gewöhnliche chemische Aktivierungsbarriere dominiert.

57. Simulation diffusionskontrollierter Reaktionen

Verbreitete Verfahren sind:

  • Brownsche Dynamik,
  • Langevin-Molekulardynamik,
  • reaktive Randbedingungen,
  • milestoning und Markov-State-Modelle,
  • partikelbasierte Reaktions-Diffusionssimulationen.

Für Proteine müssen elektrostatische Lenkung, Formanisotropie und orientierungsabhängige Kontaktbedingungen berücksichtigt werden.

58. Typische Fehler

FehlerKorrektur
jede schnelle bimolekulare Rate als diffusionskontrolliert bezeichnenDer gemessene Wert muss mit einem realistischen kdiff verglichen werden.
Drel mit nur einem Diffusionskoeffizienten ansetzenFür unabhängige Partner gilt Drel=DA+DB.
den Literfaktor bei der Smoluchowski-Konstante vergessenm³ mol−1 s−1 wird durch Multiplikation mit 1000 zu L mol−1 s−1.
den Reaktionsradius als exakte geometrische Molekülgröße behandelnR ist oft ein effektiver Kontakt- oder Elektronentransferabstand.
c(R)=0 für jede Reaktion verwendenTeilreaktive Kontakte benötigen Collins–Kimball- oder orientierungsabhängige Bedingungen.
Encounter-Komplex und stabiles Zwischenprodukt gleichsetzenEin Encounter-Komplex ist häufig nur eine transiente räumliche Korrelation.
geminate Rekombination als gewöhnliche bimolekulare Rückreaktion behandelnGeminate Partner sind räumlich und häufig spinmäßig korreliert.
Makroviskosität ohne Prüfung als lokale Reibung einsetzenMikroviskosität und frequenzabhängige Reibung können abweichen.
Kramers-Theorie nur als k∝1/η darstellenBei kleiner Reibung steigt die Rate zunächst mit der Reibung; es gibt einen Turnover.
Stern–Volmer-Linearität automatisch als dynamisches Quenching beweisenLebensdauermessungen und statische Komplexbildung müssen geprüft werden.
elektrostatische Effekte ignorierenPotential of Mean Force und Ionenstärke können den Diffusionsfluss stark verändern.
anomale Diffusion mit einem konstanten D beschreibenBei ⟨r²⟩∝tα ist die normale Smoluchowski-Theorie nicht direkt anwendbar.


Abbildung 3. Bei kleiner Reibung begrenzt der Energieaustausch mit dem Bad, bei großer Reibung die räumliche Bewegung über die Barriere. Dazwischen liegt das Kramers-Maximum.

59. Übungsaufgaben

Aufgabe 1: Stokes–Einstein-Diffusion

Ein kugelförmiges Molekül besitzt den hydrodynamischen Radius a=0,500 nm. Es diffundiert bei 298,15 K in Wasser mit η=0,890 mPa s. Berechne:

D=kBT/(6πηa)

Aufgabe 2: Smoluchowski-Grenzwert

Zwei Reaktanden besitzen:

DA=8,00×10−10 m² s−1
DB=6,00×10−10 m² s−1
R=0,800 nm

Berechne kdiff in m³ mol−1 s−1 und L mol−1 s−1.

Aufgabe 3: Diffusionslänge

Ein Molekül besitzt D=1,00×10−9 m² s−1. Welche Wurzel der mittleren quadratischen Verschiebung erreicht es in t=1,00 μs?

rrms=√(6Dt)

Aufgabe 4: Collins–Kimball-Kinetik

Für eine bimolekulare Reaktion gelten:

kdiff=8,00×109 L mol−1 s−1
kchem=2,00×109 L mol−1 s−1

Berechne kobs aus der Widerstandsaddition und ordne das Kontrollregime ein.

Aufgabe 5: Reaktion oder Käfigflucht

Ein Encounter-Komplex reagiert mit kchem=3,00×108 s−1 und entkommt mit kesc=7,00×108 s−1. Berechne Reaktionswahrscheinlichkeit, Fluchtwahrscheinlichkeit und mittlere Käfiglebensdauer.

Aufgabe 6: Kramers-Transmissionskoeffizient

Für eine Barrierenmode seien ωb=1,00×1013 s−1 und γ=5,00×1013 s−1. Berechne:

κK=[√((γ/2)²+ωb²)−γ/2]/ωb

Wie groß ist k, wenn kTST=2,00×106 s−1?

Aufgabe 7: Coulomb-Anziehung am Kontakt

Ein einfach positiv und ein einfach negativ geladenes Ion befinden sich in Wasser bei 298,15 K im Abstand R=0,500 nm. Verwende εr=78,5. Berechne:

  1. U(R) in kJ mol−1,
  2. U/(kBT),
  3. den lokalen Boltzmann-Verstärkungsfaktor exp[−U/(kBT)].

Aufgabe 8: Verdopplung der Viskosität

Eine vollständig diffusionskontrollierte Reaktion folgt Stokes–Einstein und Smoluchowski. Die Viskosität wird bei gleicher Temperatur von 0,890 auf 1,780 mPa s verdoppelt. Wie ändern sich D und kdiff im idealen Modell?

Aufgabe 9: Dynamisches Fluoreszenzquenching

Ein Fluorophor besitzt τ0=4,00 ns. Die bimolekulare Quenching-Konstante beträgt kq=1,00×1010 L mol−1 s−1, die Quencherkonzentration 0,0500 mol L−1. Berechne I0/I und τ/τ0.

Aufgabe 10: Diffusive Begegnungszeitskala

Zwei Partner sind zunächst ungefähr 100 nm voneinander entfernt und besitzen Drel=2,00×10−9 m² s−1. Schätze mit:

t≈r²/(6Drel)

die Diffusionszeitskala zum Überbrücken dieser Entfernung.

60. Vollständige Lösungen

Lösung 1

Einsetzen in die Stokes–Einstein-Gleichung:

D=(1,380649×10−23·298,15)/[6π·0,890×10−3·0,500×10−9]
D≈4,91×10−10 m² s−1

Der Wert liegt in der typischen Größenordnung kleiner gelöster Moleküle in Wasser.

Lösung 2

Relative Diffusion:

Drel=8,00×10−10+6,00×10−10
Drel=1,40×10−9 m² s−1

Smoluchowski-Konstante:

kdiff=4π(0,800×10−9)(1,40×10−9)(6,02214076×1023)
kdiff≈8,476×106 m³ mol−1 s−1

Umrechnung in Liter:

kdiff≈8,476×109 L mol−1 s−1

Dies ist ein typischer Diffusionsgrenzwert in einer niedrigviskosen Flüssigkeit.

Lösung 3
rrms=√[6·1,00×10−9·1,00×10−6]
rrms=7,746×10−8 m
rrms≈77,5 nm

Die Größe ist ein Ensemble-Mittel der quadratischen Verschiebung, keine geradlinig zurückgelegte Strecke.

Lösung 4
1/kobs=1/(8,00×109)+1/(2,00×109)
kobs=1,60×109 L mol−1 s−1

Die chemische Kontaktreaktion ist viermal langsamer als der Diffusionsgrenzwert. Das System ist überwiegend aktivierungskontrolliert, besitzt aber einen merklichen Diffusionsbeitrag.

Die Damköhler-Zahl ist:

Da=kchem/kdiff=0,250
Lösung 5

Gesamtverlustrate:

kges=3,00×108+7,00×108=1,00×109 s−1

Reaktionswahrscheinlichkeit:

Preakt=3,00×108/1,00×109=0,300

Fluchtwahrscheinlichkeit:

Pesc=0,700

Lebensdauer:

τcage=1/(1,00×109)=1,00 ns

Pro gebildetem Encounter-Komplex reagieren im Modell 30 Prozent, während 70 Prozent entkommen.

Lösung 6
κK=[√((2,50×1013)²+(1,00×1013)²)−2,50×1013]/(1,00×1013)
κK≈0,1926

Damit:

k=0,1926·2,00×106
k≈3,85×105 s−1

Die Lösungsmittelreibung reduziert die ideale TST-Rate in diesem Modell auf etwa 19 Prozent.

Lösung 7

Coulomb-Energie eines einfach geladenen Ionenpaars:

U=−e²/(4πε0εrR)
U≈−5,877×10−21 J pro Paar

Molar:

U≈−3,540 kJ mol−1

Thermisch normiert:

U/(kBT)≈−1,428

Lokaler Boltzmannfaktor:

exp[−U/(kBT)]=exp(1,428)≈4,17

Die einfache Kontaktbetrachtung sagt eine etwa vierfach erhöhte lokale Paarwahrscheinlichkeit voraus. Eine vollständige Ratenberechnung benötigt das gesamte radiale Potential und die Debye–Smoluchowski-Integration.

Lösung 8

Bei gleicher Temperatur und Teilchengröße gilt:

D∝1/η

Verdopplung von η halbiert daher beide Diffusionskoeffizienten. Somit wird auch:

Drel,neu=Drel,alt/2

Da kdiff∝Drel:

kdiff,neu=kdiff,alt/2

Das Ergebnis gilt nur im idealen Stokes–Einstein- und vollständig diffusionskontrollierten Grenzfall.

Lösung 9
I0/I=1+kqτ0[Q]
=1+(1,00×1010)(4,00×10−9)(0,0500)
I0/I=3,00

Für rein dynamisches Quenching gilt dasselbe Verhältnis für die Lebensdauern:

τ0/τ=3,00
τ/τ0=1/3≈0,333

Die beobachtete Lebensdauer beträgt daher ungefähr 1,33 ns.

Lösung 10
t≈(100×10−9)²/[6·2,00×10−9]
t≈8,33×10−7 s
t≈0,833 μs

Dies ist eine grobe Diffusionszeitskala. Eine exakte mittlere Erstpassagezeit hängt von Geometrie, Randbedingungen und Reaktionsradius ab.

61. Zusammenfassung

  • Brownsche Bewegung entsteht aus ungeordneten Lösungsmittelstößen.
  • Für normale dreidimensionale Diffusion gilt ⟨r²⟩=6Dt.
  • Die Einstein-Beziehung lautet D=kBT/ζ.
  • Für eine Kugel mit haftender Randbedingung gilt ζ=6πηa.
  • Stokes–Einstein liefert D=kBT/(6πηa).
  • Für zwei unabhängige Reaktanden ist Drel=DA+DB.
  • Der Reaktionsradius ist häufig eine effektive Kontaktgröße.
  • Smoluchowski nimmt eine vollständig absorbierende Kontaktfläche an.
  • Der molare Grenzwert ist kdiff=4πRDrelNA.
  • Typische Diffusionsgrenzen kleiner Moleküle liegen bei 109 bis 1010 L mol−1 s−1.
  • Ein Encounter-Komplex ist eine transiente Nahkontaktpopulation.
  • Der Käfigeffekt erzeugt wiederholte Kontakte und geminate Rekombination.
  • Collins–Kimball erlaubt eine nur teilweise reaktive Kontaktfläche.
  • Diffusions- und chemischer Widerstand addieren sich reziprok.
  • Orientierung und Rotationsdiffusion können die Rate weiter begrenzen.
  • Debye–Smoluchowski berücksichtigt radiale intermolekulare Potentiale.
  • Elektrostatische Anziehung kann den Kontaktfluss erhöhen, Abstoßung senkt ihn.
  • Fluoreszenzquenching kann eine direkt messbare diffusionskontrollierte Rate liefern.
  • Dynamisches Quenching verkürzt Intensität und Lebensdauer.
  • Kramers-Theorie beschreibt Barrierenüberquerung in einem Wärmebad.
  • Bei sehr kleiner Reibung begrenzt Energiediffusion.
  • Bei hoher Reibung begrenzt räumliche Diffusion, und k fällt ungefähr wie 1/η.
  • Zwischen beiden Grenzfällen liegt der Kramers-Turnover.
  • Generalisierte Langevin- und Grote–Hynes-Theorien berücksichtigen Lösungsmittelerinnerung.
  • Crowding, Membranen und heterogene Medien können anomale Diffusion erzeugen.
  • Makroviskosität und lokale molekulare Reibung sind nicht immer identisch.

62. Häufig gestellte Fragen

Was bedeutet „diffusionskontrolliert“?

Die chemische Umwandlung am Kontakt ist so schnell, dass hauptsächlich die Zeit bis zur Begegnung der Reaktanden die beobachtete Rate bestimmt.

Ist 1010 L mol−1 s−1 eine universelle Obergrenze?

Nein. Der Wert hängt von Diffusionskoeffizienten, Reaktionsradius, Temperatur, Viskosität und intermolekularen Potentialen ab. Er ist nur eine typische Größenordnung.

Warum werden die Diffusionskoeffizienten addiert?

Beide Moleküle bewegen sich unabhängig. Die Varianzen ihrer relativen Verschiebung addieren sich, daher Drel=DA+DB.

Was ist der Unterschied zwischen Encounter-Komplex und Zwischenprodukt?

Ein Encounter-Komplex ist eine räumlich korrelierte Nahkontaktpopulation. Ein Zwischenprodukt entspricht gewöhnlich einem lokalen Minimum mit eigener chemischer Struktur.

Warum reagiert nicht jede Begegnung?

Falsche Orientierung, Solvathülle, Spin, geringe elektronische Kopplung oder eine verbleibende Aktivierungsbarriere können die Reaktion verhindern.

Was bedeutet der Käfigeffekt?

Das Lösungsmittel hält frisch erzeugte Partner vorübergehend nahe beieinander und ermöglicht wiederholte Kontakte oder geminate Rekombination.

Warum kann höhere Viskosität eine Rekombination erhöhen?

Sie verlangsamt die Käfigflucht. Dadurch bleibt mehr Zeit für erneute Kontakte und Rückreaktion.

Warum steigt die Kramers-Rate bei sehr kleiner Reibung zunächst an?

Ohne ausreichende Kopplung an das Wärmebad erhält das System nur langsam die notwendige Aktivierungsenergie. Etwas stärkere Reibung beschleunigt diesen Energieaustausch.

Warum sinkt die Rate bei großer Reibung?

Die Reaktionskoordinate kann sich nur langsam über die Barriere bewegen und rekreuzt häufiger. Die räumliche Diffusion entlang der Koordinate wird begrenzend.

Wie erkennt man dynamisches Fluoreszenzquenching?

Intensität und Fluoreszenzlebensdauer sinken im gleichen Stern–Volmer-Verhältnis. Bei rein statischem Quenching sinkt die Intensität, die Lebensdauer der verbleibenden freien Fluorophore jedoch nicht.

63. Ausblick auf Teil 8

Teil 8 behandelt schnelle Reaktionsmethoden und Relaxationskinetik. Stopped-Flow, Flash-Photolyse, Temperatur- und Drucksprung, Relaxation nach kleinen Störungen, gekoppelte Eigenmoden und globale kinetische Anpassung werden systematisch entwickelt.

Literatur und weiterführende Quellen

  1. P. Atkins, J. de Paula & J. Keeler: Atkins’ Physical Chemistry. Oxford University Press.
  2. M. von Smoluchowski: Versuch einer mathematischen Theorie der Koagulationskinetik kolloider Lösungen. Zeitschrift für Physikalische Chemie.
  3. F. C. Collins & G. E. Kimball: Diffusion-Controlled Reaction Rates. Journal of Colloid Science.
  4. H. A. Kramers: Brownian Motion in a Field of Force and the Diffusion Model of Chemical Reactions. Physica.
  5. N. Agmon & A. Szabo: Theory of Reversible Diffusion-Influenced Reactions. Journal of Chemical Physics.

Didaktischer Hinweis: Smoluchowski-, Collins–Kimball- und Kramers-Modelle sind idealisierte Grenzbeschreibungen. Reale Flüssigkeitsreaktionen können anisotrope Molekülformen, orientierungsabhängige Reaktivität, nichtmarkovsche Reibung, elektrostatische Lenkung, anomale Diffusion und gekoppelte chemische Schritte enthalten.

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.

Keine Kommentare:

Kommentar veröffentlichen