Titelbild: Eigene Illustration einer quantenmechanischen Wellenfunktion an einer Barriere sowie der unterschiedlichen Tunnelwahrscheinlichkeit leichter und schwerer Isotope.
Quantenmechanische Teilchen können eine klassische Energiebarriere durchdringen, obwohl ihre Energie unterhalb des Barrierenmaximums liegt. Dieser Tunneleffekt ist für Elektronen besonders ausgeprägt, kann aber auch bei leichten Atomkernen chemische Reaktionsraten wesentlich verändern. Protonen- und Wasserstofftransferreaktionen sind die wichtigsten Beispiele.
Die Tunnelwahrscheinlichkeit hängt exponentiell von drei Größen ab:
- der Masse des übertragenen Teilchens,
- der Breite der Barriere,
- dem Energieabstand zwischen Teilchenenergie und Barriere.
Deshalb reagieren H- und D-substituierte Moleküle häufig deutlich unterschiedlich. Deuterium ist ungefähr doppelt so schwer wie Protium, besitzt niedrigere Schwingungsfrequenzen und Nullpunktsenergien und tunnelt wesentlich schlechter. Der Quotient:
wird als kinetischer H/D-Isotopeneffekt bezeichnet.
Lernziele
- klassisch erlaubte und verbotene Bereiche unterscheiden,
- die WKB-Tunnelwahrscheinlichkeit interpretieren und berechnen,
- den Einfluss von Masse, Energie und Barrierenbreite erklären,
- Rechteck-, Parabel- und Eckart-Barrieren unterscheiden,
- Wigner-Korrektur und Crossover-Temperatur berechnen,
- Nullpunktsenergie und Frequenzverschiebung bei Isotopensubstitution bestimmen,
- primäre, sekundäre und Lösungsmittel-Isotopeneffekte unterscheiden,
- klassische und tunnelbedingte Beiträge zu einem KIE trennen,
- Temperaturabhängigkeiten von Isotopeneffekten auswerten,
- experimentelle Hinweise auf Protonentunneln kritisch beurteilen,
- Corner Cutting, Instanton und Proton-Coupled Electron Transfer einordnen.
1. Klassische Barrierenüberquerung
In der klassischen Mechanik kann ein Teilchen eine eindimensionale Barriere V(x) nur überqueren, wenn seine Gesamtenergie E mindestens so groß wie das Barrierenmaximum ist:
Ist E kleiner, existieren Umkehrpunkte x1 und x2, an denen:
Zwischen diesen Punkten wäre die klassische kinetische Energie negativ. Der Bereich ist klassisch verboten.
2. Quantenmechanische Wellenfunktion
Ein quantenmechanisches Teilchen wird durch eine Wellenfunktion beschrieben. An einer Barriere wird die einlaufende Welle teilweise reflektiert und teilweise transmittiert.
Selbst für E<Vmax verschwindet die Wellenfunktion im verbotenen Bereich nicht abrupt. Sie fällt exponentiell ab und kann auf der anderen Seite der Barriere eine endliche Amplitude besitzen.
3. Transmissionswahrscheinlichkeit
Die Transmissionswahrscheinlichkeit T ist das Verhältnis von transmittiertem zu einfallendem Wahrscheinlichkeitsstrom:
Für eine eindimensionale stationäre Streuung gilt außerdem:
wenn keine Absorption oder zusätzliche Kanäle vorliegen.
4. Exponentieller Abfall im verbotenen Bereich
Für eine lokal konstante Barriere mit V>E lautet die Schrödinger-Gleichung:
mit:
Die Lösungen sind exponentiell wachsend oder fallend. Physikalisch relevant ist innerhalb einer ausgedehnten Barriere die abfallende Komponente:
5. WKB-Näherung
Die Wentzel-Kramers-Brillouin-Näherung – WKB – beschreibt eine langsam veränderliche Barriere semiklassisch. Zwischen den klassischen Umkehrpunkten lautet die Tunnelwahrscheinlichkeit in führender Näherung:
Der Exponent ist eine Wirkung in Einheiten von ℏ. Er misst die „Schwierigkeit“ des Durchtritts durch die verbotene Region.
6. Rechteckbarriere
Für eine Rechteckbarriere der Höhe V0 und Breite a wird das WKB-Integral einfach:
Diese tiefe Tunnelgrenze ist besonders anschaulich, aber reale chemische Barrieren sind glatt und meist asymmetrisch.
7. Einfluss der Teilchenmasse
Die Masse erscheint als Quadratwurzel im Exponenten:
Verdopplung der Masse halbiert die Tunnelwahrscheinlichkeit daher nicht, sondern kann sie um mehrere Größenordnungen verringern. Genau deshalb unterscheiden sich H- und D-Transfer so deutlich.
8. Einfluss der Barrierenbreite
Eine kleine Vergrößerung des Donor-Akzeptor-Abstands kann die Rate stark senken. In Enzymen, Wasserstoffbrücken und Protonentransferketten ist daher nicht nur die Barrierenhöhe, sondern auch die geometrische Kompression der Reaktionsstelle wichtig.
9. Einfluss der Teilchenenergie
Mit wachsendem E wird V−E kleiner. Die verbotene Region wird energetisch weniger und häufig auch räumlich schmaler. T steigt deshalb stark an.
Thermisch aktiviertes Tunneln kombiniert beide Mechanismen: Das System wird zunächst auf ein höheres Schwingungsniveau angeregt und tunnelt von dort durch den verbleibenden Teil der Barriere.
10. Barrierenform
Zwei Barrieren mit gleicher Höhe können sehr unterschiedliche Tunnelwahrscheinlichkeiten besitzen. Entscheidend ist das gesamte Integral über √(V−E).
Eine schmale hohe Barriere kann leichter durchtunnelbar sein als eine breite niedrigere. Aus einer einzigen Aktivierungsenergie lässt sich daher keine zuverlässige Tunnelrate ableiten.
11. Grenzen der WKB-Näherung
WKB funktioniert am besten, wenn sich die lokale de-Broglie-Wellenlänge nur langsam ändert. Schwierigkeiten entstehen:
- direkt an Umkehrpunkten,
- bei sehr dünnen abrupten Barrieren,
- nahe dem Barrierenmaximum,
- bei starken Mehrdimensionalitätseffekten,
- bei Resonanzen und Interferenz.
Verbindungsformeln oder exaktere Barrierenmodelle korrigieren diese Bereiche.
Abbildung 1. Die WKB-Wirkung wächst mit Teilchenmasse, Barrierenbreite und Energieabstand zur Barriere. Deuterium wird daher wesentlich stärker unterdrückt als Protium.
12. Invertierte Parabel
Nahe einem glatten Barrierenmaximum kann das Potential als invertierte Parabel angenähert werden:
ωb ist der Betrag der imaginären Übergangszustandsfrequenz entlang der Reaktionskoordinate. Eine große |ωb| entspricht einer schmalen, stark gekrümmten Barriere.
13. Transmission durch eine Parabelbarriere
Für die ideale invertierte Parabel ergibt sich eine Fermi-artige Transmission:
Unterhalb der Barriere beschreibt sie Tunneln, oberhalb berücksichtigt sie weiterhin quantenmechanische Reflexion.
14. Wigner-Tunnelkorrektur
Die einfachste Korrektur einer klassischen Übergangszustandstheorie ist die Wigner-Näherung:
Mit der imaginären Wellenzahl ṽ‡:
Sie ist eine Hochtemperaturentwicklung und eignet sich nur für schwaches Tunneln. Große Korrekturfaktoren aus dieser Formel sind selbst ein Warnsignal, dass die Näherung außerhalb ihres Gültigkeitsbereichs verwendet wird.
15. Eckart-Barriere
Die Eckart-Barriere ist eine glatte analytische Funktion, die asymmetrische Reaktionen beschreiben kann. Ihre Parameter werden typischerweise aus:
- Vorwärtsbarriere,
- Rückwärtsbarriere oder Reaktionsenergie,
- imaginärer Übergangszustandsfrequenz
bestimmt. Die Transmission kann analytisch berechnet und über die thermische Energieverteilung gemittelt werden.
16. Vorteil der Eckart-Korrektur
Im Gegensatz zur Wigner-Korrektur verwendet die Eckart-Näherung die gesamte modellierte Barrierenform. Sie kann daher tiefes Tunneln und Asymmetrie besser erfassen.
Ihre Genauigkeit bleibt jedoch davon abhängig, wie gut eine eindimensionale Eckart-Barriere den mehrdimensionalen Reaktionsweg repräsentiert.
17. Bell-Korrekturen
Unter dem Begriff Bell-Tunnelkorrektur werden verschiedene thermisch gemittelte eindimensionale Modelle zusammengefasst. Sie verbinden klassische Überbarrierenreaktion und quantenmechanische Transmission.
Historisch waren sie wichtig, weil sie aus wenigen kinetischen Parametern handhabbare Korrekturen lieferten. Moderne Rechnungen verwenden häufig Eckart-, variationale oder Instanton-Methoden.
18. Crossover-Temperatur
Für eine parabolische Barriere wird eine charakteristische Temperatur definiert:
oder mit der imaginären Wellenzahl:
Oberhalb Tc ist Tunneln meist eine Korrektur der klassischen Barrierenüberquerung. Unterhalb Tc wird eine ausgedehnte Tunnelbahn – ein Instanton – relevant.
19. Instanton-Theorie
Instanton-Theorie beschreibt Tunneln durch eine periodische klassische Bahn in imaginärer Zeit. Diese Bahn minimiert die euklidische Wirkung und liefert den dominanten semiklassischen Tunnelpfad.
Sie ist mehrdimensional und kann daher vom klassischen Minimum-Energy-Path abweichen.
20. Corner Cutting
Der klassische Reaktionspfad minimiert lokal die Potentialenergie. Der Tunnelpfad minimiert dagegen eine Wirkung, die sowohl Barrierenhöhe als auch Weglänge und Masse enthält.
Er kann deshalb eine geometrische Ecke abschneiden: eine etwas höhere, aber deutlich schmalere Region kann günstiger sein als der längere Weg entlang des MEP.
21. Warum Isotope chemisch ähnlich, kinetisch aber verschieden sind
Isotope eines Elements besitzen dieselbe Kernladung und nahezu dieselbe elektronische Struktur. Im Born-Oppenheimer-Potential V(R) ist die Elektronenenergie daher im Wesentlichen gleich.
Die Kernmassen unterscheiden sich jedoch. Dadurch ändern sich:
- Schwingungsfrequenzen,
- Nullpunktsenergien,
- Rotationskonstanten,
- thermische Zustandssummen,
- klassische Geschwindigkeiten,
- Tunnelwahrscheinlichkeiten.
Isotopeneffekte sind deshalb empfindliche Sonden der Kernbewegung im geschwindigkeitsbestimmenden Bereich einer Reaktion.
22. Frequenzverschiebung bei Isotopensubstitution
Für einen harmonischen Oszillator:
Bei unveränderter Kraftkonstante gilt:
Für eine weitgehend lokalisierte X–H-Schwingung ist die X–D-Frequenz ungefähr um den Faktor 1/√2 kleiner, sofern X wesentlich schwerer als H ist.
23. Nullpunktsenergie
Der harmonische Grundzustand besitzt selbst bei T=0 eine Energie:
Da H-Schwingungen höhere Frequenzen als D-Schwingungen besitzen, liegt ihre Nullpunktsenergie höher.
Für eine Aktivierungsbarriere zählt die Differenz der Nullpunktsenergien zwischen Übergangszustand und Reaktand:
24. Ursprung eines normalen primären H/D-Effekts
Wird eine X–H-Bindung im Übergangszustand stark geschwächt oder gebrochen, sinkt ihre stabile Schwingungsfrequenz. Die hohe H-Nullpunktsenergie des Reaktanten wird dabei stärker abgesenkt als die D-Nullpunktsenergie.
Dadurch ist die nullpunktsenergiekorrigierte Aktivierungsbarriere für H häufig kleiner:
und damit:
Dies ist ein normaler kinetischer Isotopeneffekt.
25. Vollständige KIE-Beziehung
Aus der Eyring-Gleichung folgt für zwei Isotope:
Der Isotopeneffekt enthält damit zwei konzeptionelle Beiträge:
- einen freien Energieunterschied der Aktivierungsbarrieren, einschließlich Nullpunktsenergie und Zustandssummen,
- das Verhältnis der Transmissionskoeffizienten, in dem Tunneln besonders stark erscheinen kann.
26. Primärer kinetischer Isotopeneffekt
Ein primärer KIE liegt vor, wenn die isotopisch markierte Bindung im geschwindigkeitsbestimmenden Schritt gebrochen oder gebildet wird.
Typische Beispiele:
- Protonentransfer,
- Hydridtransfer,
- C–H-Bindungsspaltung,
- H-Atom-Transfer.
Primäre H/D-Effekte sind häufig groß, aber ihre genaue Größe hängt von Übergangszustandsstruktur, Temperatur, Tunneln und Mechanismus ab.
27. Sekundärer kinetischer Isotopeneffekt
Bei einem sekundären KIE bleibt die isotopierte Bindung erhalten. Der Isotopenaustausch liegt jedoch in der Nähe des Reaktionszentrums und beeinflusst Schwingungen oder Hyperkonjugation.
Häufige Kategorien sind:
- α-sekundärer Effekt am Reaktionszentrum,
- β-sekundärer Effekt am benachbarten Atom,
- Effekte durch Hybridisierungsänderungen,
- Effekte durch Konformations- und Hyperkonjugationsänderungen.
Sekundäre Effekte sind meist kleiner als primäre, können aber mechanistisch sehr aussagekräftig sein.
28. Normaler und inverser Isotopeneffekt
Ein normaler Effekt besitzt kH/kD>1. Ein inverser Effekt besitzt:
Inverse Effekte können auftreten, wenn die isotopensensitive Bindung im Übergangszustand steifer wird, wenn sich Gleichgewichte vor dem geschwindigkeitsbestimmenden Schritt verschieben oder wenn mehrere Beiträge gegeneinander wirken.
29. Gleichgewichts- und kinetischer Isotopeneffekt
Ein Gleichgewichtsisotopeneffekt beschreibt die unterschiedliche Verteilung von Isotopen zwischen zwei Gleichgewichtszuständen:
Ein kinetischer Isotopeneffekt beschreibt das Verhältnis von Geschwindigkeitskonstanten. In mehrstufigen Mechanismen kann ein beobachteter KIE aus Vor-Gleichgewicht und intrinsischem KIE zusammengesetzt sein.
30. Intrinsischer und beobachteter KIE
Der intrinsische KIE gehört zum eigentlichen isotopensensitiven Elementarschritt. Der beobachtete KIE kann kleiner sein, wenn:
- ein anderer Schritt teilweise geschwindigkeitsbestimmend ist,
- Vor-Gleichgewichte den Effekt maskieren,
- Rückreaktionen auftreten,
- mehrere parallele Wege beitragen.
Diese Abschwächung wird häufig als kinetic complexity oder masking bezeichnet.
Abbildung 2. H und D bewegen sich auf nahezu derselben elektronischen Potentialfläche, besitzen aber unterschiedliche Schwingungsfrequenzen, Nullpunktsenergien und Tunneltransmissionen.
31. Temperaturabhängigkeit des KIE
Werden Aktivierungsparameter verwendet, gilt näherungsweise:
oder in Eyring-Form:
Ein großer positiver Unterschied ΔHD‡−ΔHH‡ lässt den KIE bei sinkender Temperatur wachsen.
32. Arrhenius-Vorfaktoren als Tunnelhinweis
Bei einer rein klassischen einfachen ZPE-Deutung liegen AH/AD häufig in einem begrenzten Bereich. Ungewöhnlich kleine oder große Vorfaktorverhältnisse wurden historisch als Tunnelhinweis verwendet.
Diese Diagnose ist nicht eindeutig, weil Vorfaktoren auch durch:
- Mechanismuswechsel,
- Entropieänderungen,
- gekoppelte Gleichgewichte,
- Temperaturabhängigkeit der Aktivierungswärmekapazität
beeinflusst werden.
33. Nichtlineare Arrhenius-Darstellung
Ein Tunnelbeitrag kann bei niedriger Temperatur relativ wichtiger werden. Dann kann ln k gegen 1/T gekrümmt sein.
Aber auch folgende Prozesse erzeugen Krümmung:
- mehrere Reaktionswege,
- Konformationswechsel,
- Phasen- oder Lösungsmitteländerungen,
- temperaturabhängige Wärmekapazitäten,
- Fall-off-Kinetik.
Krümmung ist daher ein Hinweis, kein Beweis.
34. Große H/D-Effekte
Ein sehr großer primärer H/D-KIE kann mit Protonentunneln vereinbar sein. Die häufig zitierte „klassische Obergrenze“ ist jedoch keine universelle harte Zahl. Sie hängt von Schwingungsfrequenzen, Übergangszustandsstruktur und Temperatur ab.
Ein belastbarer Nachweis vergleicht experimentelle Daten mit einer konsistenten quantenstatistischen oder dynamischen Rechnung.
35. Isotopenabhängigkeit der Tunnelwahrscheinlichkeit
Im WKB-Exponenten steht √m. Für dieselbe Barriere gilt daher näherungsweise:
wenn Deuterium näherungsweise die doppelte Masse besitzt. Schon ein moderater H-Tunnelexponent kann dadurch einen sehr großen H/D-Unterschied erzeugen.
36. Protonentransfer
Bei einem Protonentransfer bewegt sich H zwischen Donor D und Akzeptor A:
Die effektive Barriere hängt stark ab von:
- Donor-Akzeptor-Abstand,
- Protonenaffinitäten,
- Solvatation,
- Orientierung,
- Kopplung an schwere Atombewegungen.
37. Donor-Akzeptor-Kompression
Eine Bewegung schwerer Atome kann den Donor-Akzeptor-Abstand verkürzen und dadurch die Protonenbarriere schmaler machen. Diese promoting vibration oder gating coordinate kann Tunneln stark erleichtern.
Das Proton tunnelt dann nicht auf einer starren Barriere. Die Barriere fluktuiert mit dem molekularen und lösungsmittelbedingten Umfeld.
38. Adiabatischer und nichtadiabatischer Protonentransfer
Bei adiabatischem Protonentransfer bleibt die elektronische Wellenfunktion der Kernbewegung angepasst. Die Protonendynamik kann auf einer effektiven adiabatischen Fläche behandelt werden.
Bei nichtadiabatischem Protonentransfer ist die Kopplung zwischen elektronischen oder protonischen Zuständen klein. Die Rate enthält dann einen elektronischen beziehungsweise vibronischen Kopplungsfaktor und Franck-Condon-Überlappungen.
39. Proton-Coupled Electron Transfer
Beim protonengekoppelten Elektronentransfer – PCET – ändern sich Protonen- und Elektronenzustand gekoppelt. Proton und Elektron können:
- konzertiert übertragen werden,
- schrittweise über Protonentransfer und Elektronentransfer reagieren,
- an unterschiedliche Donor-Akzeptor-Paare gekoppelt sein.
PCET ist zentral in Enzymkatalyse, Photosynthese, Brennstoffzellen und Redoxchemie.
40. H-Atom- und Hydridtransfer
Ein H-Atom-Transfer bewegt formal H•, also Proton und Elektron zusammen. Ein Hydridtransfer bewegt H−, formal zwei Elektronen und ein Proton.
Die übertragene effektive Masse und die gekoppelte Elektronenumordnung unterscheiden sich. Ein beobachteter H/D-KIE allein identifiziert daher nicht eindeutig, ob Proton, H-Atom oder Hydrid übertragen wird.
41. Lösungsmittel-Isotopeneffekt
Beim kinetischen Lösungsmittelisotopeneffekt wird die Reaktion in H2O und D2O verglichen:
Der Effekt kann Protonentransfer aus dem Lösungsmittel anzeigen, enthält aber auch Änderungen von:
- Wasserstoffbrücken,
- pKa-Werten,
- Viskosität,
- Solvatation,
- Vor-Gleichgewichten.
42. Proton Inventory
Bei einem Proton Inventory wird die Rate in Mischungen aus H2O und D2O gemessen. Die Abhängigkeit von der Deuteriumfraktion kann Hinweise auf die Zahl und Art austauschbarer Protonen im geschwindigkeitsbestimmenden Bereich geben.
Die Interpretation verwendet Fractionation Factors und ist empfindlich gegenüber gekoppelten Gleichgewichten.
43. Sekundäre α-Isotopeneffekte
Ein α-sekundärer KIE tritt auf, wenn das isotopierte Atom direkt am Reaktionszentrum sitzt, seine C–H- oder C–D-Bindung aber nicht gebrochen wird.
Hybridisierungsänderungen, beispielsweise sp³→sp², verändern Biege- und Streckfrequenzen und damit die Nullpunktsenergie.
44. Sekundäre β-Isotopeneffekte
β-Effekte entstehen häufig durch Hyperkonjugation. Eine C–H- oder C–D-Bindung am Nachbaratom koppelt unterschiedlich stark an ein entstehendes leeres oder teilweise besetztes Orbital.
Sie werden oft zur Untersuchung von Carbokationen, Radikalen und Eliminierungsmechanismen eingesetzt.
45. Schweratom-Isotopeneffekte
Auch 13C, 15N, 18O, 34S und andere Isotope erzeugen KIEs. Wegen der kleineren relativen Massenänderung sind diese Effekte meist viel kleiner als H/D-Effekte.
Präzise Messungen natürlicher Isotopenhäufigkeiten können dennoch zeigen, welche Bindungen im Übergangszustand verändert werden.
46. Tunneln schwerer Atome
Tunneln ist nicht auf Protonen beschränkt. Bei sehr niedriger Temperatur, schmalen Barrieren und geeigneter Geometrie können auch C-, N- oder O-Atome tunneln.
Die große Masse unterdrückt die Wahrscheinlichkeit stark, aber kryogene Matrixexperimente und spezielle Umlagerungen können messbare Schweratom-Tunnelprozesse zeigen.
47. Enzymtunneln
Enzyme können Protonen-, H-Atom- oder Hydridtransfer mit großen H/D-Effekten katalysieren. Moderne Deutungen vermeiden die zu einfache Aussage, das Enzym „erzeuge“ Tunneln.
Vielmehr kann die Proteinstruktur:
- Donor und Akzeptor vororganisieren,
- den Abstand komprimieren,
- elektrostatische Barrieren senken,
- geeignete Konformationen selektieren,
- schwere Atombewegungen mit dem Transfer koppeln.
48. Dynamische Gating-Modelle
In Gating-Modellen wird zunächst eine geeignete Donor-Akzeptor-Geometrie thermisch erreicht. Aus dieser Konfiguration tunnelt das leichte Teilchen.
Die beobachtete Rate kann daher sowohl von einer klassischen Konformationsbarriere als auch von einer quantenmechanischen Transferwahrscheinlichkeit abhängen.
49. Instanton und Ringpolymer
Pfadintegralmethoden stellen ein Quantenteilchen als zyklische Kette klassischer Repliken – Beads – dar. Unterhalb Tc lokalisiert sich die dominante Ringpolymerbahn entlang eines Instantons.
Ring-Polymer-Instanton-Theorie und verwandte Methoden erlauben mehrdimensionale Tunnelraten für komplexe Moleküle.
50. Kinetic Quantum Sieving
Poren und Molekülgitter können Isotope aufgrund unterschiedlicher Nullpunktsenergie und Tunnelmobilität trennen. Dieser Effekt wird als kinetic quantum sieving bezeichnet.
Er ist besonders bei H2/D2 und sehr engen Poren bei niedriger Temperatur relevant.
51. Messstrategien
Eine überzeugende Untersuchung eines kinetischen Isotopeneffekts umfasst möglichst:
- identische Reaktionsbedingungen für beide Isotopologe,
- mehrere Temperaturen,
- Kontrolle von Gleichgewichts- und Lösungsmittelisotopeneffekten,
- Prüfung auf Mechanismuswechsel,
- vollständige Fehleranalyse,
- Berechnung von ZPE- und Tunnelbeiträgen.
52. Konkurrenz paralleler Wege
Besitzt eine Reaktion einen H-sensitiven und einen H-unempfindlichen Weg:
dann wird der beobachtete KIE durch den unempfindlichen Weg in Richtung eins verdünnt.
53. Commitment und Maskierung
In Enzymkinetik kann ein Substrat nach Bindung bereits mit hoher Wahrscheinlichkeit zu Produkt weiterreagieren. Vorwärts- und Rückwärts-Commitments verändern dann den beobachteten Isotopeneffekt.
Der intrinsische Effekt des chemischen Schritts kann wesentlich größer sein als der Effekt auf kcat oder kcat/KM.
54. Rechnerische Bestimmung von KIEs
Eine klassische harmonische TST-Rechnung benötigt für jedes Isotopolog:
- dieselbe elektronische Potentialfläche,
- isotopenspezifische Frequenzen und Zustandssummen,
- Nullpunktsenergien,
- konsistente Standardzustände,
- gegebenenfalls Tunnelkorrekturen.
Da die elektronische Energie bei Born-Oppenheimer-Näherung gleich bleibt, entstehen Unterschiede hauptsächlich aus Kernbewegung und Masse.
55. Typische Fehler
| Fehler | Korrektur |
|---|---|
| Tunneln als klassisches „Ausleihen“ von Energie erklären | Tunneln folgt aus der Wellenfunktion und einer endlichen Transmission durch den verbotenen Bereich. |
| nur die Barrierenhöhe betrachten | Breite, Form, Masse und Ausgangsenergie bestimmen den WKB-Exponenten. |
| Wigner-Korrektur bei sehr großem Faktor unkritisch verwenden | Sie ist nur eine Hochtemperaturnäherung für schwaches Tunneln. |
| die imaginäre Frequenz als stabile Schwingung in die ZPE einsetzen | Die instabile Reaktionsmode wird nicht als gebundener Oszillator gezählt. |
| jeden H/D-KIE ausschließlich als Tunnelbeweis interpretieren | Nullpunktsenergie, Gleichgewichte und Mechanismuskomplexität können große Effekte erzeugen. |
| eine universelle klassische Obergrenze für kH/kD annehmen | Die klassische Erwartung hängt von Frequenzen, Temperatur und Übergangszustand ab. |
| beobachteten und intrinsischen KIE gleichsetzen | Vor-Gleichgewichte, Commitments und parallele Wege können maskieren. |
| Lösungsmittelisotopeneffekt direkt als einzelnes übertragenes Proton zählen | Solvatation, pKa, Viskosität und mehrere Protonen können beitragen. |
| inverse Effekte als Messfehler behandeln | Sie können aus steiferer Bindung, Vor-Gleichgewicht oder sekundären Effekten entstehen. |
| IRC und Tunnelpfad gleichsetzen | Ein Instanton kann durch Corner Cutting deutlich vom MEP abweichen. |
| Arrhenius-Krümmung allein als Tunnelbeweis ansehen | Mehrere klassische Ursachen erzeugen ebenfalls Krümmung. |
| Isotopensubstitution als Änderung der elektronischen PES behandeln | In Born-Oppenheimer-Näherung bleibt die PES nahezu gleich; Kernbewegung und Zustandssummen ändern sich. |
Abbildung 3. Tunnelmodelle besitzen unterschiedliche Gültigkeitsbereiche. Ein belastbarer Nachweis kombiniert Isotopeneffekt, Temperaturabhängigkeit, Barrierenform und mechanistische Kontrolle.
56. Übungsaufgaben
Aufgabe 1: WKB-Tunneln von H und D
Ein H- beziehungsweise D-Teilchen tunnelt durch eine rechteckige Barriere. Der Energieabstand beträgt V0−E=0,150 eV, die Breite a=0,700 Å. Verwende:
mit mH=1,007825 u und mD=2,014102 u. Berechne TH, TD und TH/TD.
Aufgabe 2: Empfindlichkeit gegenüber der Barrierenbreite
Für ein Proton sei V0−E=0,200 eV. Berechne die WKB-Wahrscheinlichkeit für a=0,700 Å und a=1,000 Å. Um welchen Faktor sinkt T bei der Verbreiterung?
Aufgabe 3: Wigner-Korrektur
Ein Übergangszustand besitzt eine imaginäre Wellenzahl vom Betrag 1000 cm−1. Berechne bei 300 K:
Aufgabe 4: Crossover-Temperatur
Berechne für |ṽ‡|=1200 cm−1:
Ordne ein, ob 298 K oberhalb oder unterhalb dieser Temperatur liegt.
Aufgabe 5: Nullpunktsenergie-Modell eines primären KIE
Eine reaktive X–H-Streckschwingung besitzt im Reaktanten 3000 cm−1 und im Übergangszustand eine stabile Restfrequenz von 1200 cm−1. Für D skaliere beide Frequenzen mit 1/√2. Berechne:
- die D-Frequenzen,
- den ZPE-Beitrag zur Aktivierungsenergie für H und D,
- ΔΔE‡=ΔED‡−ΔEH‡,
- den daraus folgenden klassischen KIE bei 298,15 K.
Verwende 1 cm−1=0,01196266 kJ mol−1.
Aufgabe 6: KIE aus freier Aktivierungsenergie
Die freie Aktivierungsbarriere für D liegt bei 298,15 K um 4,50 kJ mol−1 über der H-Barriere. Transmissionskoeffizienten seien gleich. Berechne kH/kD.
Aufgabe 7: Temperaturabhängigkeit eines KIE
Für die Isotopendifferenzen seien:
Berechne den KIE bei 250 K, 298,15 K und 350 K mit:
Aufgabe 8: Kinetischer Lösungsmittelisotopeneffekt
Eine Reaktion besitzt in H2O die Konstante 0,180 s−1 und in D2O 0,0300 s−1. Berechne den KSIE. Nenne mindestens drei mögliche Beiträge, die vor einer direkten Zuordnung zu einem einzelnen Protonentransfer geprüft werden müssen.
Aufgabe 9: Klassischer und tunnelbedingter Beitrag
Eine harmonische TST-Rechnung ohne Tunneln ergibt kH/kD=3,57. Eine Tunnelrechnung liefert κH/κD=5,00. Berechne den gesamten KIE.
Aufgabe 10: Arrhenius-Parameter eines Isotopeneffekts
Für zwei Isotopologe gilt:
Berechne kH/kD bei 298,15 K aus:
57. Vollständige Lösungen
Für H ergibt sich der WKB-Exponent:
Damit:
Für D:
Verhältnis:
Obwohl D nur ungefähr doppelt so schwer ist, ist seine Tunnelwahrscheinlichkeit hier um fast zwei Größenordnungen kleiner.
Für a=0,700 Å:
Für a=1,000 Å:
Abnahmefaktor:
Eine Verbreiterung um nur 0,300 Å senkt die Wahrscheinlichkeit in diesem Modell um mehr als den Faktor 360.
Der dimensionslose Parameter ist:
Damit:
Die Wigner-Näherung würde die klassische TST-Rate also ungefähr um den Faktor 1,96 erhöhen. Da die Korrektur bereits groß ist, sollte ein besseres Barrierenmodell geprüft werden.
298 K liegt oberhalb von Tc. In dieser groben Einordnung ist Tunneln bei 298 K eher eine Korrektur zur thermischen Barrierenüberquerung. Die Nähe zu Tc zeigt jedoch, dass eine einfache Hochtemperaturkorrektur bereits ungenau werden kann.
1. D-Frequenzen:
2. ZPE-Beitrag zur Barriere:
3. Isotopendifferenz:
4. Klassischer KIE:
Dieses vereinfachte Ergebnis entsteht allein aus einer isotopenabhängigen Nullpunktsenergiekorrektur. Zusätzliche Tunnelbeiträge können den Effekt weiter vergrößern.
Die H-Reaktion ist bei gleichen Transmissionskoeffizienten ungefähr sechsmal schneller.
Bei 250 K:
Bei 298,15 K:
Bei 350 K:
Der Effekt wächst beim Abkühlen, weil die positive Aktivierungsenthalpiedifferenz stärker durch RT gewichtet wird.
Vor einer direkten Zuordnung zu einem einzelnen Protonentransfer müssen unter anderem geprüft werden:
- Verschiebung von Säure-Base-Gleichgewichten und pKa-Werten,
- unterschiedliche Solvatation und Wasserstoffbrückenstruktur,
- Viskositäts- und Diffusionseffekte,
- mehrere austauschbare Protonen,
- Mechanismus- oder geschwindigkeitsbestimmender-Schritt-Wechsel.
Die vollständige Beziehung ist:
Mit den ungerundeten Ausgangswerten ergäbe sich ungefähr 17,84. Der Tunnelbeitrag verstärkt den bereits vorhandenen Nullpunktsenergieeffekt erheblich.
Die größere Aktivierungsenergie des D-Isotopologs überkompensiert das kleine Vorfaktorverhältnis. Ein ungewöhnliches AH/AD kann ein Tunnelhinweis sein, beweist Tunneln aber nicht allein.
58. Zusammenfassung
- Quantenmechanische Teilchen besitzen eine endliche Transmission durch klassisch verbotene Bereiche.
- Im verbotenen Bereich fällt die Wellenfunktion exponentiell ab.
- Die WKB-Wahrscheinlichkeit enthält das Integral über √[2m(V−E)]/ℏ.
- Teilchenmasse, Barrierenbreite und Barrierenform wirken exponentiell auf die Rate.
- H tunnelt wesentlich leichter als D.
- Eine invertierte Parabel beschreibt die lokale Krümmung am Übergangszustand.
- Die Wigner-Korrektur ist nur für schwaches Tunneln und hohe Temperaturen geeignet.
- Die Eckart-Barriere berücksichtigt Asymmetrie und eine vollständige eindimensionale Barrierenform.
- Die Crossover-Temperatur trennt grob Hochtemperaturkorrektur und tiefe Instantonregion.
- Instanton-Tunnelpfade können den klassischen MEP durch Corner Cutting verlassen.
- Isotopensubstitution ändert Schwingungsfrequenzen und Nullpunktsenergien.
- Ein primärer KIE betrifft eine gebrochene oder gebildete isotopierte Bindung.
- Sekundäre KIEs entstehen ohne Bruch der isotopierten Bindung.
- Normale KIEs besitzen kH/kD>1, inverse Werte kleiner als eins.
- Der beobachtete KIE kann durch Vor-Gleichgewichte, Commitments und Parallelwege maskiert sein.
- Die vollständige KIE-Beziehung enthält freie Barrieren und Transmissionskoeffizienten.
- Temperaturabhängige KIEs liefern Information über Aktivierungsenthalpie, Entropie und Tunneln.
- Große KIEs und Arrhenius-Krümmung sind Hinweise, aber keine alleinigen Beweise.
- Protonentransfer ist stark vom Donor-Akzeptor-Abstand abhängig.
- Schwere Atombewegungen können die Protonenbarriere dynamisch komprimieren.
- PCET koppelt Protonen- und Elektronentransfer.
- Lösungsmittelisotopeneffekte enthalten mehr als nur die Masse eines übertragenen Protons.
- Pfadintegral- und Instantonmethoden behandeln mehrdimensionales Tunneln.
- Eine belastbare Analyse kombiniert Temperaturreihen, Isotopologe, Mechanismuskontrolle und Theorie.
59. Häufig gestellte Fragen
Woher nimmt ein tunnelndes Teilchen die fehlende Energie?
Es „leiht“ keine Energie. Der quantenmechanische Zustand besitzt eine endliche Wellenfunktionsamplitude auf der anderen Seite der Barriere. Die Gesamtenergie bleibt erhalten.
Warum tunnelt H besser als D?
Der WKB-Exponent wächst mit der Quadratwurzel der Masse. Die größere D-Masse führt daher zu einem wesentlich stärkeren exponentiellen Abfall.
Ist jeder große H/D-KIE ein Beweis für Tunneln?
Nein. Nullpunktsenergien, Vor-Gleichgewichte und Mechanismuskomplexität können große Effekte erzeugen. Ein Tunnelbeitrag sollte quantitativ modelliert werden.
Warum ist die Barrierenbreite so wichtig?
Sie steht linear im WKB-Exponenten. Bereits Zehntel Ångström können Raten um Größenordnungen verändern.
Was bedeutet die imaginäre Übergangszustandsfrequenz?
Sie misst die negative lokale Krümmung entlang der instabilen Reaktionsrichtung und dient als Maß für die Barrierenkrümmung.
Wann ist die Wigner-Korrektur ungeeignet?
Bei niedrigen Temperaturen, sehr großen imaginären Frequenzen, starkem Tunneln oder stark asymmetrischen beziehungsweise mehrdimensionalen Barrieren.
Was ist der Unterschied zwischen primärem und sekundärem KIE?
Beim primären Effekt wird die isotopierte Bindung im relevanten Schritt gebrochen oder gebildet. Beim sekundären bleibt sie erhalten, beeinflusst aber den Übergangszustand.
Kann ein kinetischer Isotopeneffekt invers sein?
Ja. Eine steifere isotopensensitive Bindung im Übergangszustand, Vor-Gleichgewichte oder konkurrierende Schwingungsbeiträge können kH/kD<1 erzeugen.
Was ist Corner Cutting?
Der Tunnelpfad nimmt eine kürzere Route durch eine etwas höhere Potentialregion, statt dem längeren klassischen Minimum-Energy-Path zu folgen.
Warum genügt eine einzelne Messtemperatur nicht?
Die Temperaturabhängigkeit hilft, freie Energie-, Entropie- und Tunnelbeiträge zu trennen und Mechanismuswechsel oder Krümmung zu erkennen.
60. Ausblick auf Teil 7
Teil 7 behandelt diffusionskontrollierte Reaktionen und Lösungsmittelreibung. Smoluchowski- und Collins-Kimball-Modelle, Kramers-Barrierenüberquerung, Viskosität, Käfigeffekt, Encounter Complexes und Reaktionen in Flüssigkeiten werden mit den dynamischen Tunnel- und Übergangszustandskonzepten verbunden.
Literatur und weiterführende Quellen
- P. Atkins, J. de Paula & J. Keeler: Atkins’ Physical Chemistry. Oxford University Press.
- R. P. Bell: The Tunnel Effect in Chemistry. Chapman and Hall.
- D. G. Truhlar, B. C. Garrett & S. J. Klippenstein: Current Status of Transition-State Theory. Journal of Physical Chemistry.
- A. Kohen & H.-H. Limbach, Hrsg.: Isotope Effects in Chemistry and Biology. CRC Press.
- J. Meisner & J. Kästner: Atom Tunneling in Chemistry. Angewandte Chemie International Edition.
Didaktischer Hinweis: Eindimensionale Rechteck-, Parabel- und Eckart-Barrieren sind kontrollierte Modelle. Reale Tunnelreaktionen sind häufig mehrdimensional, an schwere Atombewegungen und Lösungsmittel gekoppelt und können mehrere elektronische oder protonische Zustände umfassen.
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