Potentialenergieflächen und Reaktionspfade – Minima, Sattelpunkte, IRC, NEB und dynamische Verzweigungen

Titelbild: Eigene Illustration einer zweidimensionalen Potentiallandschaft mit Edukt- und Produktminimum, Übergangszustand und Reaktionspfad.

Eine chemische Reaktion ist eine Bewegung durch einen hochdimensionalen Raum atomarer Geometrien. Jeder Punkt dieses Raums entspricht einer bestimmten Kernanordnung, und jeder Kernanordnung ist im Rahmen der Born-Oppenheimer-Näherung eine elektronische Energie zugeordnet. Die Gesamtheit dieser Energien bildet eine Potentialenergiefläche.

Minima entsprechen stabilen Molekülen oder Zwischenprodukten. Sattelpunkte erster Ordnung markieren lokale Engpässe zwischen benachbarten Minima. Ein Reaktionspfad verbindet diese stationären Punkte. Doch eine statische Kurve auf der Potentialfläche ist noch keine vollständige Dynamik: Reale Moleküle besitzen Impuls, Schwingungsphase und Drehimpuls und können vom Minimum-Energy-Path abweichen.

Teil 5 entwickelt deshalb zwei Perspektiven gleichzeitig:

• die statische Topologie der Potentialenergiefläche mit Minima, Sattelpunkten, Gradienten und Hesse-Matrizen,
• die dynamische Bewegung von Trajektorien, einschließlich Verzweigungen, Valley Switching und post-transition-state bifurcations.

Einordnung: Teil 4 untersuchte einzelne Molekülstöße und ihre Winkelverteilungen. Teil 5 erklärt, auf welcher Energielandschaft diese Trajektorien verlaufen. Teil 6 behandelt anschließend Tunneleffekt und kinetische Isotopeneffekte.

Lernziele

Nach diesem Beitrag kannst du:

• eine Potentialenergiefläche in der Born-Oppenheimer-Näherung definieren,
• die Zahl innerer Freiheitsgrade eines Moleküls bestimmen,
• Gradient und Hesse-Matrix interpretieren,
• Minima und Sattelpunkte über Hesse-Eigenwerte klassifizieren,
• eine imaginäre Übergangszustandsfrequenz richtig deuten,
• lokale quadratische Potentiale verwenden,
• Reaktionskoordinate, Minimum-Energy-Path und IRC unterscheiden,
• die Grundidee von NEB-, Climbing-Image- und String-Verfahren erklären,
• Potentialenergie- und freie Energieflächen auseinanderhalten,
• den Committor als Test einer Reaktionskoordinate einordnen,
• Valley-Ridge-Inflection-Punkte und post-TS-Bifurkationen beschreiben,
• statische Pfade und reale dynamische Trajektorien unterscheiden.

1. Born-Oppenheimer-Potentialenergiefläche

Elektronen bewegen sich aufgrund ihrer viel kleineren Masse im Allgemeinen schneller als Atomkerne. In der Born-Oppenheimer-Näherung werden die Kernpositionen R zunächst festgehalten und die elektronische Schrödinger-Gleichung gelöst:

Ĥ_el(R)ψ_el(r;R)=E_el(R)ψ_el(r;R)

Zusammen mit der Kern-Kern-Abstoßung ergibt sich die potentielle Energie der Kerne:

V(R)=E_el(R)+V_NN(R)

Die Funktion V(R) ist eine adiabatische Potentialenergiefläche für einen bestimmten elektronischen Zustand.

2. Dimension der Energielandschaft

Ein Molekül mit N Kernen besitzt 3N kartesische Kernkoordinaten. Davon entfallen:

• drei auf die Gesamttranslation,
• drei auf die Gesamtrotation eines nichtlinearen Moleküls,
• zwei auf die Gesamtrotation eines linearen Moleküls.

Damit besitzt ein nichtlineares Molekül:

f=3N−6

innere Freiheitsgrade, ein lineares Molekül:

f=3N−5

Eine gewöhnliche zweidimensionale Konturkarte ist daher fast immer nur ein Ausschnitt oder eine Projektion einer wesentlich höherdimensionalen Fläche.

3. Wahl der Koordinaten

Potentialflächen können in verschiedenen Koordinatensystemen dargestellt werden:

• kartesische Koordinaten,
• Bindungslängen, Bindungswinkel und Torsionen,
• Normalschwingungskoordinaten,
• massengewichtete kartesische Koordinaten,
• kollektive Variablen wie Abstände oder Koordinationszahlen.

Die Form einer gezeichneten Fläche hängt von der Koordinatenwahl ab. Stationäre Punkte und Energiedifferenzen sind physikalisch invariant, lokale Krümmungen und Pfadlängen dagegen koordinatenabhängig.

4. Potentialenergie und freie Energie

Eine Potentialenergiefläche ordnet jeder Kerngeometrie eine Energie zu. Eine freie Energiefläche integriert dagegen nicht explizit dargestellte Freiheitsgrade statistisch aus.

A(ξ)=−k_BT ln P(ξ)+C

ξ ist eine kollektive Variable und P(ξ) ihre Gleichgewichtsverteilung. Eine freie Energiebarriere enthält daher sowohl energetische als auch entropische Beiträge.

In Lösung oder bei flexiblen Biomolekülen ist eine freie Energiefläche häufig aussagekräftiger als eine einzelne optimierte Gasphasen-Potentialfläche.

5. Gradient

Der Gradient enthält die ersten Ableitungen der Energie:

g_i(q)=∂V/∂q_i

Die Kraft ist das negative Gradientenfeld:

F_i=−g_i

Eine Geometrieoptimierung bewegt die Struktur idealerweise in eine Richtung, die V senkt, bis der Gradient verschwindet.

6. Stationärer Punkt

Ein stationärer Punkt erfüllt:

∇V(q_s)=0

Diese Bedingung allein unterscheidet noch nicht zwischen Minimum, Maximum und Sattelpunkt. Dafür wird die lokale Krümmung benötigt.

7. Hesse-Matrix

Die Hesse-Matrix besteht aus zweiten Ableitungen:

H_ij=∂²V/(∂q_i∂q_j)

Sie ist für glatte Potentiale symmetrisch. Ihre Eigenvektoren geben lokale Hauptkrümmungsrichtungen an; die Eigenwerte λ_k geben Stärke und Vorzeichen der Krümmung an:

H e_k=λ_ke_k

8. Massengewichtete Hesse-Matrix

Für Schwingungsfrequenzen wird die Hesse-Matrix massengewichtet:

F=M^−1/2HM^−1/2

Die Eigenwerte der massengewichteten Matrix stehen mit den quadratischen Kreisfrequenzen in Beziehung:

λ_k=ω_k²

Positive λ erzeugen reelle Frequenzen. Negative λ erzeugen eine imaginäre Frequenz und zeigen eine instabile lokale Richtung.

9. Minimum

Ein nichtentartetes lokales Minimum besitzt:

• verschwindenden Gradienten,
• positive Hesse-Eigenwerte in allen inneren Freiheitsgraden.

Es entspricht einer stabilen Gleichgewichtsstruktur, einem Konformer, Reaktanden-, Produkt- oder Zwischenproduktminimum.

10. Übergangszustand

Ein Übergangszustand im üblichen Sinn ist ein Sattelpunkt erster Ordnung:

• ∇V=0,
• genau ein negativer Hesse-Eigenwert,
• alle übrigen inneren Eigenwerte positiv.

Die Eigenrichtung des negativen Eigenwerts zeigt lokal die Bindungsänderung entlang der Reaktionskoordinate.

11. Imaginäre Frequenz

Da λ=ω² negativ ist, wird ω formal imaginär. Quantenchemieprogramme geben dann beispielsweise eine Frequenz als −800 cm⁻¹ oder 800i cm⁻¹ aus.

Diese Zahl ist keine beobachtbare stabile Schwingung. Sie beschreibt die lokale Instabilität des Sattelpunkts. Eine plausible Übergangszustandsmode sollte die erwartete Bindungsbildung, Bindungsspaltung oder Umlagerung zeigen.

12. Sattelpunkte höherer Ordnung

Ein Sattelpunkt zweiter Ordnung besitzt zwei negative Eigenwerte, ein Sattelpunkt dritter Ordnung drei und so weiter.

Solche Punkte können für die globale Topologie wichtig sein, sind aber gewöhnlich nicht die elementaren Engpässe einer einzelnen Reaktion. Optimierungsverfahren können versehentlich auf ihnen enden.

13. Null- und weiche Moden

Bei einer nicht massengewichteten kartesischen Analyse treten für ein isoliertes Molekül Translationen und Rotationen als theoretische Nullmoden auf. Numerisch sind sie oft kleine positive oder negative Werte.

Sehr kleine innere Eigenwerte können außerdem auf:

• weiche Torsionen,
• nahezu freie Rotation,
• flache Potentialtäler,
• unzureichende numerische Konvergenz

hinweisen.

14. Lokale Taylor-Entwicklung

In der Nähe eines stationären Punkts q_s gilt:

V(q_s+Δq)≈V(q_s)+1/2 Δq^THΔq

Der lineare Term verschwindet, weil ∇V=0. In den Hesse-Eigenkoordinaten:

ΔV≈1/2 Σ_kλ_kQ_k²

Am Minimum steigt die Energie in jeder Richtung. Am Übergangszustand fällt sie entlang genau einer Richtung und steigt in allen anderen.

Abbildung 1. Gradient und Hesse-Matrix klassifizieren die lokale Topologie: Minima besitzen nur positive innere Krümmungen, Übergangszustände genau eine negative.

15. Reaktionskoordinate

Eine Reaktionskoordinate ist eine reduzierte Variable, die den Fortschritt von Edukten zu Produkten beschreibt. Sie kann eine einzelne Bindungslänge sein, ist aber häufig eine Kombination vieler Kernbewegungen.

Eine gute Reaktionskoordinate trennt reaktive und nichtreaktive Konfigurationen möglichst eindeutig. Eine schlechte Koordinate kann scheinbare Barrieren oder verborgene Zwischenzustände erzeugen.

16. Minimum-Energy-Path

Ein Minimum-Energy-Path – MEP – ist ein Pfad, der lokal entlang des Talbodens der Potentialenergiefläche verläuft und Reaktanden-, Übergangszustands- und Produktregion verbindet.

Auf einem idealen MEP besitzt der Gradient keine Komponente senkrecht zum Pfad:

g_⊥=0

Der Gradient zeigt nur entlang der Pfadtangente.

17. Intrinsische Reaktionskoordinate

Die intrinsische Reaktionskoordinate – IRC – ist der massengewichtete steilste Abstiegsweg vom Übergangszustand zu benachbarten Minima.

In massengewichteten Koordinaten q lautet eine idealisierte Weggleichung:

dq/ds=−g(q)/|g(q)|

Für die andere Seite des Übergangszustands wird das Vorzeichen umgekehrt.

18. Start der IRC

Direkt am Sattelpunkt ist g=0, sodass die Differentialgleichung dort nicht unmittelbar gestartet werden kann. Man verschiebt die Geometrie zunächst um einen kleinen Betrag entlang der positiven beziehungsweise negativen instabilen Eigenmode.

Von diesen beiden Startpunkten wird der steilste Abstieg numerisch integriert.

19. Zweck einer IRC-Rechnung

Eine IRC-Rechnung beantwortet vor allem:

• Welche beiden Minima verbindet der gefundene Sattelpunkt?
• Entspricht die imaginäre Mode der erwarteten Reaktion?
• Gibt es entlang des Abstiegs unerwartete Strukturänderungen?

Sie ist damit eine zentrale Validierung eines berechneten Übergangszustands.

20. Grenzen der IRC

Die IRC ist ein geometrischer Pfad ohne Impulsinformation. Sie sagt nicht voraus:

• welche Trajektorien bei gegebener Temperatur tatsächlich genommen werden,
• wie Energie zwischen Moden verteilt wird,
• wie oft Rekreuzung auftritt,
• welches Produkt bei einer dynamischen Bifurkation entsteht.

Ein Reaktionsmechanismus ist daher nicht vollständig durch eine IRC-Kurve definiert.

21. Reaktionspfadlänge

In massengewichteten Koordinaten wird ein Bogenelement häufig als:

ds²=Σ_idq_i²

definiert. Durch die Massengewichtung erhalten Bewegungen schwerer und leichter Kerne eine dynamisch sinnvollere Skalierung als in gewöhnlichen kartesischen Koordinaten.

Die absolute Pfadlänge ist dennoch kein direkt beobachtbarer Stoffparameter. Sie hängt von Koordinatenwahl und Parametrisierung ab.

22. Reaktionspfadkrümmung

Die normierte Tangente des Pfads ist:

t(s)=dq/ds

Die Krümmung kann durch:

κ_Pfad(s)=|dt/ds|

charakterisiert werden. Große Pfadkrümmung zeigt, dass sich die Richtung der Kernbewegung rasch ändert. Dadurch kann Energie aus der Reaktionskoordinate in orthogonale Schwingungsmoden übertragen werden.

23. Reaktionspfad-Hamiltonoperator

In einer lokalen Beschreibung wird die Bewegung in eine Pfadkoordinate s und orthogonale Moden Q_k zerlegt. Die Frequenzen dieser Moden können sich entlang des Pfads verändern.

Kopplungsterme zwischen Pfadkrümmung und orthogonalen Schwingungen erklären, warum ein Molekül selbst bei Bewegung nahe dem MEP Energie zwischen Freiheitsgraden austauscht.

24. Valley und Ridge

Ein Valley ist eine Region, die quer zum Pfad positive Krümmung besitzt. Kleine seitliche Auslenkungen werden zurück zum Talboden gedrückt.

Eine Ridge besitzt in mindestens einer seitlichen Richtung negative Krümmung. Auslenkungen wachsen und können Trajektorien in verschiedene Täler lenken.

25. Geometrieoptimierung von Minima

Minima werden durch Algorithmen gesucht, die Energie und Gradient verwenden. Verbreitete Methoden sind:

• steilstes Gefälle,
• konjugierte Gradienten,
• Newton-Raphson,
• Quasi-Newton-Verfahren wie BFGS,
• rationale Funktionsoptimierung.

Newton-Raphson verwendet lokal:

Δq=−H⁻¹g

In der Nähe eines gutartigen Minimums konvergiert es schnell. Weit davon entfernt können negative oder sehr kleine Hesse-Eigenwerte problematisch sein.

26. Übergangszustandssuche

Ein Übergangszustand ist kein Energiemaximum in allen Richtungen. Eine gewöhnliche Maximierung würde deshalb nicht funktionieren. Gesucht wird ein Maximum entlang einer ausgewählten Richtung und ein Minimum in allen orthogonalen Richtungen.

Eigenvector-Following invertiert die Optimierungsrichtung entlang der vermuteten instabilen Mode, während alle übrigen Richtungen minimiert werden.

27. Berny- und RFO-Verfahren

Quantenchemische Programme verwenden häufig interne Koordinaten, approximierte Hesse-Matrizen und rational-function optimization. Die Methode stabilisiert Schritte, wenn die quadratische Näherung allein zu große Bewegungen vorschlagen würde.

Eine gute Übergangszustandsstartgeometrie bleibt entscheidend. Ohne chemische Vorinformation kann die Suche zu einem unerwünschten Sattelpunkt konvergieren.

28. Dimer-Methode

Die Dimer-Methode benötigt keinen vollständigen Hesse-Matrixaufbau. Zwei nahe Bilder – der „Dimer“ – rotieren, um die Richtung kleinster Krümmung zu finden. Anschließend wird entlang dieser Richtung bergauf und senkrecht dazu bergab optimiert.

Sie ist besonders nützlich für große Systeme, bei denen die Hesse-Matrix teuer ist.

29. Eigenmode und chemische Interpretation

Eine formal korrekte imaginäre Frequenz genügt nicht. Die animierte Eigenmode sollte chemisch plausibel sein:

• Bindungsbruch oder -bildung,
• Protonentransfer,
• Torsion oder Ringöffnung,
• koordinierte Mehrzentrenbewegung.

Eine rein translatorische oder rotationsartige imaginäre Mode deutet häufig auf schlechte Konvergenz oder numerische Artefakte hin.

30. Nudged Elastic Band

Die Nudged Elastic Band – NEB – sucht einen MEP zwischen bekannten Anfangs- und Endminima. Zwischen beiden Strukturen wird eine Kette von Bildern angelegt.

Jedes Bild erfährt:

• die wahre Potentialkraft senkrecht zum Pfad,
• eine künstliche Federkraft entlang des Pfads.

Die „Nudging“-Projektion verhindert, dass die Federkraft den Pfad seitlich von der Potentialtalsohle wegzieht.

31. NEB-Kraft

Schematisch:

F_i^NEB=−∇V(R_i)_⊥+F_i^Feder_∥

Die Federkomponente hält die Bilder entlang des Pfads verteilt. Ohne sie würden sich die meisten Bilder in den Minima sammeln und die Barrierenregion schlecht auflösen.

32. Climbing-Image NEB

Bei Climbing-Image NEB wird das energiereichste Bild besonders behandelt. Die Federkraft wird entfernt und die Potentialkraft entlang der Pfadtangente umgekehrt.

Das Bild steigt dadurch zum Sattelpunkt erster Ordnung, während es senkrecht zum Pfad minimiert wird.

33. Zahl und Verteilung der Bilder

Zu wenige Bilder können einen gekrümmten Pfad abschneiden oder Zwischenprodukte übersehen. Zu viele Bilder erhöhen den Rechenaufwand und können numerische Steifigkeit erzeugen.

Adaptive oder energiegewichtete Bildverteilungen verbessern die Auflösung in Barrierenregionen.

34. String-Methode

Die String-Methode beschreibt den Reaktionspfad ebenfalls durch diskrete Bilder. Nach einem Relaxationsschritt senkrecht zum Pfad wird der String regelmäßig neu parametrisiert, sodass die Bilder gleichmäßig oder nach einer gewünschten Metrik verteilt bleiben.

Im Gegensatz zur klassischen NEB sind keine expliziten Federn erforderlich.

35. Growing-String-Methode

Bei Growing-String-Verfahren wachsen zwei Pfadenden von den Minima aufeinander zu oder ein Pfad wächst von einem bekannten Minimum in Richtung eines vermuteten Produkts.

Dies kann effizienter sein, wenn eine lineare Interpolation stark verzerrte oder chemisch unsinnige Zwischenbilder erzeugt.

Abbildung 2. IRC folgt vom bekannten Sattelpunkt bergab. NEB und String suchen dagegen einen gesamten Pfad zwischen bekannten Minima.

36. Lineare Interpolation ist kein Reaktionspfad

Eine lineare kartesische Interpolation zwischen Edukt und Produkt kann Atome unphysikalisch übereinander schieben, Bindungen gleichzeitig verzerren und unrealistisch hohe Energien erzeugen.

Sie ist lediglich eine Startschätzung für ein Pfadverfahren. Interne Koordinaten oder bildbasierte Geometrieinterpolation sind häufig besser.

37. Mehrere Reaktionspfade

Zwischen denselben Minima können mehrere Übergangszustände und Pfade existieren. Der niedrigste Potentialbarrierenweg muss nicht bei jeder Temperatur dominieren, weil Entropie, Tunneln und dynamische Faktoren unterschiedlich sein können.

Für konkurrierende Wege addieren sich die Teilraten:

k_ges=Σ_jk_j

38. Potentialbarriere und freie Barriere

Ein elektronischer Sattelpunkt liefert zunächst eine elektronische Energiebarriere. Für einen Vergleich mit Geschwindigkeitskonstanten müssen unter anderem berücksichtigt werden:

• Nullpunktsenergien,
• thermische Enthalpiekorrekturen,
• Entropie,
• Standardzustand,
• Solvatation,
• Transmissionskoeffizient.

Der kinetisch relevante Engpass ist daher ein Maximum der freien Energie, nicht zwangsläufig der höchste Punkt einer reinen Potentialkurve.

39. Potential of Mean Force

Für eine kollektive Variable ξ ist das Potential of Mean Force:

A(ξ)=−k_BT ln P(ξ)+C

Es mittelt alle übrigen Freiheitsgrade bei festem ξ. Die mittlere Kraft ist:

dA/dξ=−⟨F_ξ⟩_ξ

Ein PMF ist temperatur- und ensembleabhängig.

40. Umbrella Sampling

Hohe freie Energiebarrieren werden in einer gewöhnlichen Gleichgewichtssimulation selten besucht. Beim Umbrella Sampling wird deshalb ein Bias-Potential verwendet:

U_Bias(ξ)=1/2 k(ξ−ξ₀)²

Mehrere überlappende Fenster decken den Reaktionsbereich ab. Anschließend wird der Bias rechnerisch entfernt und die freie Energie rekonstruiert.

41. Metadynamik

Metadynamik fügt entlang gewählter kollektiver Variablen wiederholt kleine Bias-Potentiale hinzu. Bereits besuchte Regionen werden dadurch energetisch angehoben und das System entkommt aus freien Energieminima.

Das Ergebnis hängt kritisch davon ab, ob die gewählten kollektiven Variablen die langsamen Prozesse tatsächlich beschreiben.

42. Committor

Der Committor p_B(x) ist die Wahrscheinlichkeit, dass eine Trajektorie aus dem Phasenraumpunkt x zuerst Produktbecken B erreicht, bevor sie zu Eduktbecken A zurückkehrt.

• p_B≈0: klar auf Eduktseite,
• p_B≈1: klar auf Produktseite,
• p_B≈1/2: dynamische Trennregion.

Der Committor gilt als besonders strenger Test einer Reaktionskoordinate.

43. Konfigurations- und Phasenraum

Eine Potentialenergiefläche lebt im Konfigurationsraum der Kernkoordinaten. Dynamik lebt im Phasenraum aus Koordinaten und Impulsen:

x=(q,p)

Zwei Trajektorien können am selben geometrischen Punkt starten, aber aufgrund verschiedener Impulse in unterschiedliche Produkte laufen. Eine rein geometrische Reaktionskoordinate kann diese Fälle nicht unterscheiden.

44. Transition State als Trennfläche im Phasenraum

Der geometrische Sattelpunkt ist ein wichtiger Anker. Die dynamisch optimale Trennfläche ist jedoch eine Struktur im Phasenraum, die reaktive und nichtreaktive Bewegung möglichst ohne Rekreuzung trennt.

Für mehrdimensionale Systeme wird diese Idee mit normally hyperbolic invariant manifolds und zugehörigen Trennflächen formalisiert.

45. Valley-Ridge-Inflection-Punkt

An einem Valley-Ridge-Inflection-Punkt – VRI – ändert sich die transversale Krümmung einer Landschaft von talartig zu kammartig.

Ein lokales Kriterium ist:

• ein Hesse-Eigenwert wird null,
• der Gradient ist orthogonal zur zugehörigen Nullkrümmungs-Eigenrichtung,
• der Punkt ist kein stationärer Punkt.

Ein VRI kann den Übergang von einem einzigen Reaktionstal zu einer verzweigenden Landschaft markieren.

46. VRI ist nicht gleich Übergangszustand

Ein Übergangszustand erfüllt ∇V=0 und besitzt eine negative Krümmung. Ein VRI liegt im Allgemeinen bei ∇V≠0 und besitzt eine seitliche Nullkrümmung.

Beide topologischen Merkmale können auf demselben Reaktionsweg auftreten, sind aber mathematisch und physikalisch verschieden.

47. Post-transition-state bifurcation

Nach Überquerung eines gemeinsamen Übergangszustands kann sich die Landschaft in zwei oder mehr Produkttäler teilen, ohne dass zwischen ihnen ein weiterer Sattelpunkt liegt.

Dann lässt sich die Produktselektivität nicht durch zwei unabhängige Eyring-Barrieren beschreiben. Die Verzweigung wird durch Trajektoriendynamik bestimmt.

48. Valley Switching

Eine Trajektorie kann zunächst einem Tal folgen, aufgrund von Impuls und Pfadkrümmung aber in ein benachbartes Tal wechseln. Dies wird als Valley Switching bezeichnet.

Der Wechsel kann auftreten, obwohl der Minimum-Energy-Path zu einem anderen Produkt führt.

49. Dynamic Matching

Dynamic Matching beschreibt die Tendenz, dass die Impulsrichtung aus dem Übergangszustand bevorzugt mit einer bestimmten nachfolgenden Bewegung gekoppelt ist.

Eine im Übergangszustand aufgebaute Schwingungsphase kann dadurch die Produktselektivität beeinflussen, obwohl beide Produkte energetisch ähnlich zugänglich sind.

50. Reaktionspfadverzweigung und Selektivität

Bei einer rein statistischen Verteilung würde man Produktverhältnisse aus Zustandszahlen oder freien Barrieren ableiten. Bei einer dynamischen Bifurkation müssen Ensembles von Trajektorien gestartet werden.

Y_j=N_j/N_gesamt

Die Unsicherheit einer binomialen Verzweigung kann näherungsweise mit:

σ_Y≈√[Y(1−Y)/N]

abgeschätzt werden.

51. Direkte Dynamik

Bei Direct Dynamics wird die elektronische Energie und ihr Gradient während jeder Trajektorie „on the fly“ berechnet. Eine vorher global angepasste analytische Potentialfläche ist nicht erforderlich.

Vorteile sind Flexibilität und konsistente Elektronenstruktur. Nachteile sind hohe Kosten und begrenzte Zahl beziehungsweise Dauer der Trajektorien.

52. Ab-initio-Molekulardynamik

Ab-initio-Molekulardynamik integriert klassische Kernbewegung mit elektronischen Kräften aus Dichtefunktionaltheorie oder Wellenfunktionsmethoden.

Born-Oppenheimer-Molekulardynamik konvergiert die Elektronenstruktur in jedem Schritt. Car-Parrinello-artige Ansätze koppeln zusätzliche elektronische Freiheitsgrade dynamisch an.

53. Quasiklassische Trajektorien

Quasiklassische Trajektorien initialisieren quantisierte Rotations- und Schwingungsenergien, entwickeln die Kerne danach aber klassisch.

Probleme sind:

• Nullpunktsenergieleckage,
• fehlende Interferenz,
• kein echtes Tunneln,
• unscharfe Zuordnung klassischer Endenergien zu Quantenzuständen.

54. Tunneln und Corner Cutting

Eine Tunnelbahn muss nicht dem klassischen MEP folgen. In massengewichteten Koordinaten kann sie eine kürzere Strecke durch eine etwas höhere Region nehmen. Dieses „Corner Cutting“ ist besonders bei leichten Atomen wichtig.

Die optimale Tunnelbahn minimiert eine Wirkung, nicht die klassische Potentialenergie an jedem Punkt.

55. Mehrere elektronische Flächen

Photochemische und spinverbotene Prozesse können nicht auf einer einzigen adiabatischen Fläche beschrieben werden. Elektronische Zustände können sich annähern, vermeiden oder tatsächlich schneiden.

Nichtadiabatische Kopplung erlaubt Übergänge zwischen Flächen. Konische Durchschneidungen bilden dabei hochdimensionale Trichter für sehr schnelle strahlungslose Relaxation.

56. Typische Fehler

Fehler	Korrektur
eine zweidimensionale Grafik als vollständige PES behandeln	Sie ist fast immer nur Projektion oder Schnitt eines 3N−6-dimensionalen Raums.
jeden Punkt mit ∇V=0 als Minimum bezeichnen	Die Hesse-Eigenwerte bestimmen Minimum oder Sattelpunkt.
negative Frequenz als physikalische stabile Schwingung interpretieren	Sie kennzeichnet eine lokale Instabilität am Sattelpunkt.
eine kleine imaginäre Frequenz automatisch als Übergangszustand akzeptieren	Translations-, Rotations- und numerische Artefakte müssen ausgeschlossen werden.
IRC und reale Trajektorie gleichsetzen	IRC besitzt keine Impuls- oder Phaseninformation.
lineare Interpolation als MEP ansehen	Sie ist nur eine mögliche Startkette für Pfadoptimierung.
NEB-Federkräfte vollständig auf die Bilder anwenden	Nur die parallele Federkomponente gehört in die genudgte Kraft.
elektronische Barriere direkt mit Δ‡G vergleichen	Nullpunktsenergie, Wärme, Entropie, Standardzustand und Solvatation fehlen.
eine freie Energiefläche als temperaturunabhängig behandeln	Sie hängt von T, Ensemble und gewählten kollektiven Variablen ab.
pB=1/2 nur aus einer Geometrie ableiten	Der Committor ist eine dynamische Wahrscheinlichkeit aus Trajektorien.
VRI und Sattelpunkt gleichsetzen	Ein VRI ist gewöhnlich nichtstationär und besitzt eine seitliche Nullkrümmung.
Produktverhältnisse nach einer Bifurkation nur aus Barrieren bestimmen	Impuls, Schwingungsphase und Trajektorienensemble können dominieren.

Abbildung 3. Nach einem gemeinsamen Übergangszustand kann sich eine Landschaft ohne weiteren Sattelpunkt verzweigen. Impuls und Schwingungsphase bestimmen dann, welches Produkttal erreicht wird.

57. Übungsaufgaben

Aufgabe 1: Dimension einer Potentialenergiefläche

Bestimme die Zahl der inneren Freiheitsgrade für:

1. H₂O als nichtlineares dreiatomiges Molekül,
2. CO₂ als lineares dreiatomiges Molekül,
3. CH₄ als nichtlineares fünfatomiges Molekül.

Aufgabe 2: Klassifikation stationärer Punkte

Drei massengewichtete Hesse-Matrizen besitzen – nach Entfernung der Translations- und Rotationsmoden – folgende Eigenwertmuster in willkürlichen Einheiten:

Punkt	Eigenwerte
A	2,0; 5,0; 8,0
B	−3,0; 2,0; 6,0
C	−4,0; −1,0; 7,0

Klassifiziere A, B und C.

Aufgabe 3: Quadratische Potentialnäherung

In zwei Hesse-Eigenkoordinaten gelten λ₁=4,0 und λ₂=9,0 in konsistenten Einheiten. Für eine Verschiebung Q₁=0,20 und Q₂=−0,10 berechne:

ΔV=1/2(λ₁Q₁²+λ₂Q₂²)

Aufgabe 4: Ein IRC-Eulerschritt

In massengewichteten Koordinaten befindet sich ein IRC-Punkt bei q=(1,20; −0,50). Der Gradient ist g=(3,00; 4,00). Führe mit:

q_neu=q−Δs g/|g|

einen Abstiegsschritt mit Δs=0,100 aus.

Aufgabe 5: Barrierenunterschied und Rate

Zwei konkurrierende Reaktionspfade besitzen bei 298,15 K freie Aktivierungsenergien, die sich um 6,00 kJ mol⁻¹ unterscheiden. Transmissionskoeffizienten und Vorfaktoren seien gleich. Berechne das Geschwindigkeitsverhältnis des schnelleren zum langsameren Weg.

Aufgabe 6: NEB-Federkraft

Ein NEB-Bild liegt entlang der Pfadtangente 0,80 Å vom vorherigen und 1,20 Å vom nächsten Bild entfernt. Die Federkonstante beträgt 5,00 eV Å⁻². Berechne die parallele Federkraft mit:

F_Feder,∥=k_s(d_nächster−d_vorheriger)

Gib Betrag und Richtung an.

Aufgabe 7: Freie Energie aus Wahrscheinlichkeiten

Bei 298,15 K werden für zwei Werte einer kollektiven Variable die Gleichgewichtswahrscheinlichkeiten P₁=0,0500 und P₂=0,200 gemessen. Berechne A₂−A₁ aus:

A₂−A₁=−RT ln(P₂/P₁)

Aufgabe 8: Committor-Schätzung

Von 100 Trajektorien, die mit identischer Geometrie und zufällig gezogenen thermischen Impulsen gestartet werden, erreichen 47 zuerst das Produktbecken B. Schätze p_B und den binomialen Standardfehler:

σ≈√[p_B(1−p_B)/N]

Ist die Konfiguration mit einem idealen p_B=0,5 vereinbar?

Aufgabe 9: Dynamische Produktverzweigung

In einem Trajektorienensemble führen 360 von 600 Trajektorien zu Produkt A und 240 zu Produkt B. Berechne:

1. die beiden Produktausbeuten,
2. das Verhältnis A:B,
3. den binomialen Standardfehler der Ausbeute von A,
4. den scheinbaren freien Barrierenunterschied bei 298,15 K, falls man das Verhältnis fälschlich wie ein Verhältnis zweier Eyring-Raten interpretieren würde.

Aufgabe 10: Nullpunktsenergiekorrigierte Barriere

Eine vereinfachte Modellreaktion besitzt eine elektronische Barriere von 50,0 kJ mol⁻¹. Die stabilen Reaktantenschwingungen seien 500, 1000 und 1500 cm⁻¹. Der Übergangszustand besitzt zwei stabile Frequenzen von 600 und 1100 cm⁻¹ sowie eine imaginäre Mode, die nicht in die Nullpunktsenergie eingeht.

Berechne die Nullpunktsenergiekorrektur und die korrigierte Barriere. Verwende:

1 cm⁻¹=0,01196266 kJ mol⁻¹

58. Vollständige Lösungen

Lösung 1

H₂O: nichtlinear, N=3:

f=3N−6=9−6=3

CO₂: linear, N=3:

f=3N−5=9−5=4

CH₄: nichtlinear, N=5:

f=3N−6=15−6=9

Die vier Schwingungsfreiheitsgrade von CO₂ umfassen zwei entartete Biegungsrichtungen, die als zwei Freiheitsgrade gezählt werden.

Lösung 2

A: Alle Eigenwerte sind positiv. A ist ein lokales Minimum.

B: Genau ein Eigenwert ist negativ. B ist ein Sattelpunkt erster Ordnung und damit ein Kandidat für einen Übergangszustand.

C: Zwei Eigenwerte sind negativ. C ist ein Sattelpunkt zweiter Ordnung.

Die Klassifikation setzt voraus, dass Translations- und Rotationsmoden bereits entfernt wurden und keine numerischen Nullmoden vorliegen.

Lösung 3

ΔV=1/2[4,0·(0,20)²+9,0·(−0,10)²]

ΔV=1/2[4,0·0,0400+9,0·0,0100]

ΔV=1/2(0,160+0,090)=0,125

Das Vorzeichen von Q₂ verschwindet in der quadratischen Näherung. Beide positiven Eigenwerte erhöhen die Energie.

Lösung 4

Betrag des Gradienten:

|g|=√(3,00²+4,00²)=5,00

Normierter Gradient:

g/|g|=(0,600; 0,800)

Schritt:

q_neu=(1,20;−0,50)−0,100(0,600;0,800)

q_neu=(1,14;−0,58)

Der Schritt führt in Richtung des negativen Gradienten und senkt daher lokal die Energie.

Lösung 5

Bei gleichen Vorfaktoren:

k_schnell/k_langsam=exp(ΔΔG^‡/(RT))

=exp[6000/(8,314462618·298,15)]

k_schnell/k_langsam≈11,25

Eine Absenkung der freien Barriere um nur 6,00 kJ mol⁻¹ erhöht die Rate bei Raumtemperatur um mehr als eine Größenordnung.

Lösung 6

F_Feder,∥=5,00(1,20−0,80)

F_Feder,∥=2,00 eV Å⁻¹

Der positive Wert zeigt entlang der gewählten Tangentenorientierung zum weiter entfernten nächsten Bild. Die Federkraft versucht, die ungleichen Bildabstände auszugleichen.

Lösung 7

A₂−A₁=−RT ln(0,200/0,0500)

=−RT ln4

A₂−A₁≈−3,437 kJ mol⁻¹

Zustand 2 liegt in der freien Energie niedriger, weil er viermal wahrscheinlicher beobachtet wird.

Lösung 8

p_B=47/100=0,470

σ≈√[0,470·0,530/100]

σ≈0,0499

Damit:

p_B≈0,470±0,050

Der ideale Wert 0,5 liegt klar innerhalb einer Standardabweichung. Die Konfiguration ist daher mit einer dynamischen Trennregion vereinbar, obwohl ein größeres Ensemble die Aussage verbessern würde.

Lösung 9

Ausbeuten:

Y_A=360/600=0,600

Y_B=240/600=0,400

Verhältnis:

A:B=360:240=1,50:1

Standardfehler von Y_A:

σ_A≈√[0,600·0,400/600]=0,0200

Würde man das Verhältnis unzulässig als Verhältnis zweier unabhängiger Eyring-Raten deuten:

ΔG_A^‡−ΔG_B^‡=−RT ln(1,50)

≈−1,005 kJ mol⁻¹

Man würde also scheinbar eine um etwa 1,01 kJ mol⁻¹ niedrigere Barriere für A ableiten. Bei einer post-TS-Bifurkation ist diese Interpretation jedoch nicht gerechtfertigt, weil die Verzweigung nicht aus zwei unabhängigen Gleichgewichts-Übergangszuständen entstehen muss.

Lösung 10

Reaktanten-Nullpunktsenergie:

ZPE_R=1/2(500+1000+1500)·0,01196266

ZPE_R≈17,944 kJ mol⁻¹

Übergangszustands-Nullpunktsenergie; die imaginäre Mode wird ausgeschlossen:

ZPE^‡=1/2(600+1100)·0,01196266

ZPE^‡≈10,169 kJ mol⁻¹

Korrektur:

ΔZPE^‡=ZPE^‡−ZPE_R

ΔZPE^‡≈−7,776 kJ mol⁻¹

Korrigierte Barriere:

E₀^‡=50,0−7,776≈42,22 kJ mol⁻¹

Das vereinfachte Beispiel zeigt, dass Nullpunktsenergien eine Barriere deutlich verändern können. In realen Rechnungen müssen alle stabilen Moden konsistent einbezogen werden.

59. Zusammenfassung

• Eine Potentialenergiefläche ordnet jeder Kerngeometrie eine elektronische Energie zu.
• Nichtlineare Moleküle besitzen 3N−6, lineare 3N−5 innere Freiheitsgrade.
• Eine zweidimensionale Energiekarte ist nur ein Schnitt oder eine Projektion.
• Der Gradient ist die negative Kraft bis auf das Vorzeichen.
• Stationäre Punkte erfüllen ∇V=0.
• Die Hesse-Matrix beschreibt lokale Krümmungen.
• Minima besitzen nur positive innere Hesse-Eigenwerte.
• Ein Übergangszustand besitzt genau einen negativen Eigenwert.
• Die zugehörige imaginäre Frequenz beschreibt die instabile Reaktionsrichtung.
• Die quadratische Näherung ist nur lokal gültig.
• Ein MEP verläuft entlang eines Potentialtals.
• Die IRC ist der massengewichtete steilste Abstieg vom Sattelpunkt.
• Eine IRC validiert, welche Minima ein Übergangszustand verbindet.
• IRC und reale dynamische Trajektorie sind nicht identisch.
• Pfadkrümmung koppelt Reaktionsbewegung und orthogonale Moden.
• Eigenvector-Following und Dimer-Methoden suchen Sattelpunkte.
• NEB verwendet eine Bildkette mit senkrechter Potential- und paralleler Federkraft.
• Climbing-Image NEB optimiert das höchste Bild zum Sattelpunkt.
• String-Methoden relaxieren und reparametrisieren einen Pfad.
• Potential- und freie Energiebarrieren müssen unterschieden werden.
• Der Committor prüft, ob eine Konfiguration dynamisch zur Trennregion gehört.
• Ein VRI markiert eine Änderung seitlicher Krümmung und ist kein Übergangszustand.
• Post-TS-Bifurkationen können mehrere Produkte ohne weiteren Sattelpunkt erzeugen.
• Impuls und Schwingungsphase können die Produktselektivität steuern.
• Direkte Dynamik berechnet Kräfte während der Trajektorie.
• Tunnelpfade können den klassischen MEP durch Corner Cutting verlassen.
• Nichtadiabatische Reaktionen benötigen mehrere elektronische Flächen.

60. Häufig gestellte Fragen

Ist eine Potentialenergiefläche wirklich eine Fläche?

Nur bei zwei inneren Koordinaten. Für reale Moleküle ist sie meist ein hochdimensionales Hyperflächenobjekt. Die Bezeichnung „Fläche“ wird trotzdem allgemein verwendet.

Warum reicht der Gradient nicht zur Klassifikation?

Gradient null bedeutet nur stationär. Erst die Hesse-Eigenwerte zeigen, ob die Energie in den verschiedenen Richtungen steigt oder fällt.

Warum besitzt ein Übergangszustand eine imaginäre Frequenz?

Die lokale Krümmung entlang der Reaktionsrichtung ist negativ. Daher ist ω² negativ und die harmonische Frequenz formal imaginär.

Ist die IRC der wahrscheinlichste Reaktionsweg?

Nicht unbedingt. Sie ist ein massengewichteter steilster Abstiegsweg ohne Impulsinformation. Reale Trajektorien können aufgrund ihrer Dynamik abweichen.

Was unterscheidet IRC und NEB?

IRC startet an einem bekannten Sattelpunkt und folgt bergab. NEB startet mit bekannten Endminima und optimiert eine ganze Bildkette zwischen ihnen.

Warum braucht NEB künstliche Federn?

Ohne Kopplung würden Bilder in energetisch günstige Minima rutschen. Die parallele Federkraft hält eine nutzbare Verteilung entlang des Pfads aufrecht.

Was ist eine freie Energiefläche?

Sie beschreibt die logarithmische Wahrscheinlichkeit einer kollektiven Variable und enthält die statistisch integrierten Beiträge aller übrigen Freiheitsgrade.

Was bedeutet p_B=0,5?

Aus diesem Phasenraumpunkt erreichen gleich viele thermisch gestartete Trajektorien zuerst Produkt B wie Edukt A. Er liegt damit in einer dynamischen Trennregion.

Kann eine Reaktion nach einem Übergangszustand zwei Produkte bilden?

Ja. Bei einer post-transition-state bifurcation verzweigt sich die Landschaft hinter einem gemeinsamen Sattelpunkt, ohne dass ein zweiter Übergangszustand nötig ist.

Warum kann ein Produkt entstehen, das nicht auf der IRC liegt?

Die IRC enthält keine Impulse. Reale Trajektorien können durch Pfadkrümmung, Schwingungsphase und dynamische Kopplung in ein anderes Tal wechseln.

61. Ausblick auf Teil 6

Teil 6 behandelt Tunneleffekt und kinetische Isotopeneffekte. Wigner- und Eckart-Korrekturen, semiklassische Tunnelwahrscheinlichkeiten, primäre und sekundäre H/D-Isotopeneffekte, Temperaturabhängigkeit und Protonentransfer werden mit den Potentialbarrieren dieses Teils verbunden.

Literatur und weiterführende Quellen

1. P. Atkins, J. de Paula & J. Keeler: Atkins’ Physical Chemistry. Oxford University Press.
2. W. H. Miller, N. C. Handy & J. E. Adams: Reaction Path Hamiltonian for Polyatomic Molecules. Journal of Chemical Physics.
3. K. Fukui: The Path of Chemical Reactions – the IRC Approach. Accounts of Chemical Research.
4. G. Henkelman, B. P. Uberuaga & H. Jónsson: A Climbing Image Nudged Elastic Band Method. Journal of Chemical Physics.
5. B. K. Carpenter: Determination of Organic Reaction Mechanisms. Wiley.

Didaktischer Hinweis: Stationäre Punkte und Reaktionspfade hängen von Elektronenstrukturmethode, Basissatz, Solvatationsmodell und Koordinatenwahl ab. Eine belastbare Mechanismusanalyse kombiniert Frequenzen, IRC oder Pfadverfahren, freie Energien, Konformationssuche und bei dynamisch verzweigenden Systemen Trajektorienensembles.

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.