ReaxFF (für englisch reactive force-field, Reaktions-Kraftfeld) ist ein Kraftfeld zur Simulation chemischer Reaktionen, welches auf der Bindungsordnung basiert. Es wurde von Adri van Duin, William A. Goddard, III et al. am California Institute of Technology entwickelt. Kraftfelder herkömmlicher Molekulardynamik-Simulationen (z. B. mit AMBER, CHARMM) sind nicht in der Lage, chemische Reaktionen zu modellieren, da die verwendeten harmonischen Potentiale bzw. festen Bindungslisten den Bruch und die Bildung chemischer Bindungen nicht beschreiben können. ReaxFF besitzt keine festen Bindungslisten. Stattdessen werden Potentiale verwendet, die von der Bindungsordnung abhängig sind. Dies gewährleistet, dass die Potentialbeiträge für verschwindende Bindungsordnung (große interatomare Abstände) gegen Null konvergieren.

In der ursprünglichen Variante von ReaxFF wurde eine Parametrisierung für Kohlenwasserstoffe entwickelt.[1] Mittlerweile existiert eine Vielzahl von Varianten, die sich durch zusätzliche bzw. modifizierte Potentiale, sowie Parametersätze unterscheiden.

Version 2001

Bearbeiten

Die 2001 veröffentlichte Version wurde speziell für die Simulation von Kohlenwasserstoffen entwickelt. Die Systemenergie setzt sich aus folgenden Termen zusammen:

 .

Bis auf die Terme   und   hängen alle Terme von einer oder mehrerer Bindungsordnungen ab. Die unkorrigierte Bindungsordnung zwischen den Atomen i und j ist definiert durch:

 

mit den Parametern   und  . Die unkorrigierte Bindungsordnung   ist somit die Summe aus den Bindungsordnungen der Sigma- ( ), Pi-( ), und Doppelpibindungsordnungen ( ). Nach dieser Definition kann die Bindungsordnung maximal den Wert 3 annehmen, wodurch auch Doppel- und Dreifachbindungen beschrieben werden können.

Um eine Überkoordination von Atomen in einem Molekül zu vermeiden, werden die Bindungsordnungen einer Korrektur unterzogen. Diese Korrektur bewirkt, dass langreichweitige Bindungen mit niedriger Bindungsordnung, deren beteiligte Atome bereits Gesamtbindungsordnungen nahe deren Valenz besitzen, eine Bindungsordnung von Null bekommen.

Der Term   entspricht den Bindungsenergien und ist durch

 ,

mit den Parametern   und   gegeben. Trotz Bindungsordnungskorrektur kann es in Molekülen zur Über- bzw. Unterkoordination kommen. Im Falle der Überkoordination liefert der Term

 

einen positiven Energiebeitrag. Darin sind   und   Parameter. Die Überkoordination des Atoms i ergibt sich durch

 .
  ist die Valenz des Atoms i (z. B. 4 für Kohlenstoff).

Der Term   gibt einen zusätzlichen Energiebeitrag durch die Resonanz der  -Elektronen zwischen benachbarten unterkoordinierten Atomen. In einem Molekül oder Festkörper besitzen Atomtripel ijk einen optimalen Bindungswinkel (Gleichgewichtswinkel), welcher von den Elementen und der chemischen Umgebung abhängt (z. B. beträgt der H-C-H Winkel in CH4 109.7°). Eine Abweichung von diesem Gleichgewichtswinkel führt zu einer Verschlechterung der Gesamtenergie. In ReaxFF wird dies durch den Term   realisiert. Im Gegensatz zu klassischen Kraftfeldern hängt dieser selbst von der chemischen Umgebung (Summer der  -Bindungsordnungen) ab.

Literatur

Bearbeiten
  1. A. C. T. van Duin, S. Dasgupta, F. Lorant, W. A. Goddard III, „ReaxFF: A Reactive Force Field for Hydrocarbons“, The Journal of Physical Chemistry A, 105, 9396–9409 (2001); doi:10.1021/jp004368u.
Bearbeiten