ReaxFF (for “reactive force field”) is a bond order-based
force field developed by Adri van Duin,
William A. Goddard, III, and co-workers at the
California Institute of Technology
The California Institute of Technology (branded as Caltech) is a private research university in Pasadena, California, United States. The university is responsible for many modern scientific advancements and is among a small group of institutes ...
. One of its applications is
molecular dynamics
Molecular dynamics (MD) is a computer simulation method for analyzing the Motion (physics), physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a fixed period of time, giving a view of the dynamics ( ...
simulations. Whereas traditional force fields are unable to model chemical reactions because
of the requirement of breaking and forming bonds (a force field's functional form depends on having all bonds defined
explicitly), ReaxFF eschews explicit bonds in favor of
bond orders, which allows for continuous bond formation/breaking. ReaxFF aims to be as general as possible and has been parameterized and tested for hydrocarbon reactions, alkoxysilane gelation, transition-metal-catalyzed nanotube formation, and many advanced material applications such as Li ion batteries, TiO
2, polymers, and high-energy materials.
To be able to deal with bond breaking and formation whilst having only 1 single atom type for each element, ReaxFF is a fairly complex force field with many parameters. Therefore an extensive training set is necessary covering the relevant chemical phase space, including bond and angle stretches, activation and reaction energies, equation of state, surface energies, and much more. Usually, but not necessarily, the training data is generated with electronic structure methods. In practice, often
DFT calculations are used as a pragmatic approach, especially since more accurate functionals are available.
For the parameterization of such a complex force field,
global optimization
Global optimization is a branch of operations research, applied mathematics, and numerical analysis that attempts to find the global minimum or maximum of a function or a set of functions on a given set. It is usually described as a minimization ...
techniques offer the best chance to get a parameter set that most closely describes the training data.
References
*
*
*
*
*
*
*
*
*{{Cite journal , doi = 10.1021/jp504138r, title = Parallel Optimization of a Reactive Force Field for Polycondensation of Alkoxysilanes, journal = The Journal of Physical Chemistry B, volume = 118, issue = 37, pages = 10966–10978, year = 2014, last1 = Deetz , first1 = J. D. , last2 = Faller , first2 = R. , pmid = 25153668
External links
Adri van Duin's WebsiteReaxFF in the Amsterdam Modeling SuiteReaxFF in PuReMD (Purdue Reactive MD) suite
Force fields (chemistry)