Force-field-based reactive

Force-Jield-based reactive MD Empirical force fields can also be used to study chemical reactions. Starting from diatomics in molecules (DIM) [47, 48], global... 

Quantum mechanical methods can now be applied to systems with up to 1000 atoms 87 this capacity is not only from advances in computer technology but also from improvements in algorithms. Recent developments in reactive classical force fields promise to allow the study of significantly larger systems.88 Many approximations can also be made to yield a variety of methods, each of which can address a range of questions based on the inherent accuracy of the method chosen. We now discuss a range of quantum mechanical-based methods that one can use to answer specific questions regarding shock-induced detonation conditions. 

There are various methods for the prediction of stereoselectivities, both in racemate separation (thermodynamics, usually computed by force field methods) and enantioselective catalysis (reactivity, usually computed by quantum mechanics)18. Promising recent developments, primarily based on force field and statistical methods, but also involving QM modeling, are based on stereocartography, the computation of the chirality content and the evaluation of chirophores.151 154... 

The results of molecular dynamics simulations are critically dependent on the quality of the force field. In principle, any potential that accurately describes SWNT should have terms that reflect the curvature of the tube. However, curvature effects are largely ignored in force fields published in the literature.20-23 Some potentials, such as those based on reactive bond-order formalism, correctly account for the local curvature in describing the covalent bonding interactions.24-26 But the adsorption... 

Although a valence-type force field of the type illustrated by Eq. [1] is most suitable for modeling molecular systems, the electronegativity equalization approach to treating polarization can be coupled equally well to other types of potentials. Streitz and Mintmire used an EE-based model in conjunction with an embedded atom method (EAM) potential to treat polarization effects in bulk metals and oxides. The resulting ES + EAM model has been parameterized for aluminum and titanium oxides, and has been used to study both charge-transfer effects and reactivity at interfaces. 

In 2003, Han et al. reported a theoretical study of the self-capping of zigzag n, 0) SWCNTs using classical MD techniques based on a semi-empirical reactive bond order [8] (Tersoff) force field at 3000 Kelvin [9]. They found the importance of... 

ReaxFF describes the total energy of an atomistic system with three main terms i) covalent (bonds, angles, torsions, etc.), ii) electrostatics with environment-dependent charges, and iii) van der Waals interactions. Covalent interactions are based on the concept of partial bond orders that are calculated solely from atomic positions (no pre-determined connectivities). Once the bond order between every pair of atoms is known, bond energies, angles, and torsions are determined. The second key concept in reactive force fields (also used in the... 

A step forward along the route to the correct modelling of the spectroscopy and photochemical reactivity of photoreactive proteins is represented by the implementation of a Quantum Mechanics/Molecular Mechanics (QM/MM) computational strategy based on a suitable QM part coupled with a protein force field such as AMBER [34] (or CHARMM [35]). Very recently a CASPT2//CASSCF/AMBER method for rhodopsin has been implemented in our laboratory [36,37] within the QM/MM hnk-atom scheme [38]. Special care has been taken in the parametrization of the protonated Schiff base linkage region that describes the dehcate border region between the MM (the protein)... 

The AIMD method, based on the Carr and ParrineUo approach [127], has also been applied in the study of electrochemistry [128]. Reactive Force Field approaches are now being used to study the ionomer/water/catalyst interfaces during an electrochemical reaction [129]. Neurock et al. developed a detailed first-principles approach that employs a double-reference method to simulate the influence of the electrochemical potential on the chemistry at the metal/solution interface [130]. hi this method the aqueous solution metal interface and the interfacial potential drop are explicitly treated. However the choice of an appropriate water surface structure is critical for establishing the appropriate electrochemical behavior at the atomistic scale. This method has been applied to smdy some electrochemical steps involved in the ORR and methanol oxidation on Pt (e.g. [131, 132]). 

Besides electrode surface structure, the nature of the electrode and solvent also affect the evolution of the SEI. An interesting and detailed simulation of formation and growth of SEI on Li metal surface in EC, DMC, and EC mixed with DMC electrolyte was done by Kim et al., using reactive force field (ReaxFF) molecular simulations [61]. The SEI film was found to grow faster in EC-based electrolyte compared to DMC, generating thicker SEI film, and EC mixed with DMC electrolyte came in between, as shown in Fig. 5.18. This simulation result agrees with the... 

Reactive molecular dynamics (ReaxMD) is based on the Reax force field (ReaxFF) parameterized by fitting to a training set of QM data. Compared to first-principles simulations, with the ReaxFF model it is possible to speed up the calculation by several orders of magnitude. However, due to the enormous complexity of the underlying mathematical expressions, ReaxFF is -10 to 100 times more expensive computationally than simple (nonreactive) FFs. 

Mainly in the field of materials science various types of potentials have been developed based on the concept of the bond order. " Like for reactive force fields also for the application of these potentials a specification of the atomic positions is sufficient. Although many of these potentials like the Tersoff potential, the Stillinger-Weber potential, the Breimer potential and many others have been introduced already one or two decades ago, they are still frequently used in materials simulations, in particular for semiconductors. For metallic systems the embedded atom method (EAM) and the modified embedded atom method (MEAM) introduced by Baskes and coworkers are widely distributed. 

Why then bother about much more expensive QM-based models One reason is that MM may only lead to accurate results for molecules of the same type used for the optimization and validation of the force field, i.e. extrapolation is seen to be dangerous if not impossible [9], This also extends to transition states and shortlived, unstable intermediates and therefore to chemical reactivity. Since electrons are not considered explicitly in MM, electronic effects related to structural distortions, specific stabilities and spectroscopy cannot be modeled by MM. However, in all other areas, there is no good reason for not using a well-optimized and validated MM model. Also, there are MM-based approaches to deal with most of the deficiencies listed above [9,20-28]. In the last decade, there have been a number of approaches, which have, based on simple rules [29], valence bond theory [30-33] and ligand-field theory [20-23], allowed the simplification of the force-field optimization and validation procedures and/or inclusion of electronic effects in MM models. 

---