Abstract:Conspectus The development of more efficient and more accurate ways to represent reactive potential energy surfaces is a requirement for extending the simulation of large systems to more complex systems, longer-time dynamical processes, and more complete statistical mechanical sampling. One way to treat large systems is by direct dynamics fragment methods. Another way is by fitting system-specific analytic potential energy functions with methods adapted to large systems. Here we consider both approaches. First we consider three fragment methods that allow a given monomer to appear in more than one fragment. The first two approaches are the electrostatically embedded many-body (EE-MB) expansion and the electrostatically embedded many-body expansion of the correlation energy (EE-MB-CE), which we have shown to yield quite accurate results even when one restricts the calculations to include only electrostatically embedded dimers. The third fragment method is the electrostatically embedded molecular tailoring approach (EE-MTA), which is more flexible than EE-MB and EE-MB-CE. We show that electrostatic embedding greatly improves the accuracy of these approaches compared with the original unembedded approaches. Quantum mechanical fragment methods share with combined quantum mechanical/molecular mechanical (QM/MM) methods the need to treat a quantum mechanical fragment in the presence of the rest of the system, which is especially challenging for those parts of the rest of the system that are close to the boundary of the quantum mechanical fragment. This is a delicate matter even for fragments that are not covalently bonded to the rest of the system, but it becomes even more difficult when the boundary of the quantum mechanical fragment cuts a bond. We have developed a suite of methods for more realistically treating interactions across such boundaries. These methods include redistributing and balancing the external partial atomic charges and the use of tuned fluorine atoms for capping dangling bonds, and we have shown that they can greatly improve the accuracy. Finally we present a new approach that goes beyond QM/MM by combining the convenience of molecular mechanics with the accuracy of fitting a potential function to electronic structure calculations on a specific system. To make the latter practical for systems with a large number of degrees of freedom, we developed a method to interpolate between local internal-coordinate fits to the potential energy. A key issue for the application to large systems is that rather than assigning the atoms or monomers to fragments, we assign the internal coordinates to reaction, secondary, and tertiary sets. Thus, we make a partition in coordinate space rather than atom space. Fits to the local dependence of the potential energy on tertiary coordinates are arrayed along a preselected reaction coordinate at a sequence of geometries called anchor points; the potential energy function is called an anchor points reactive potential. Electrostatically embedded fragment methods and the anchor points reactive potential, because they are based on treating an entire system by quantum mechanical electronic structure methods but are affordable for large and complex systems, have the potential to open new areas for accurate simulations where combined QM/MM methods are inadequate.

Anchor Points Reactive Potential for Bond-Breaking Reactions

A brand new reactive potential made molecular dynamics on chemical reaction possible

From Ab Initio Potential Energy Surfaces to State-Resolved Reactivities: X + H2O ↔ HX + OH [X = F, Cl, and O(3P)] Reactions

Implementing Reactivity in Molecular Dynamics Simulations with Harmonic Force Fields

Methods for Calculating Intermolecular Potential Energy Surfaces

A Global Potential Energy Surface for the H2 + OH ↔ H2O + H Reaction Using Neural Networks

Global Analytical Potential Energy Surfaces for HO2(X̃A″2) Based on High-Level Ab Initio Calculations

Efficient estimation of bond dissociation energies of organic compounds

Permutation Invariant Potential Energy Surfaces for Polyatomic Reactions Using Atomistic Neural Networks

Near Equivalence of Polarizability and Bond Flux for Describing Covalent Bond Rearrangements

Accurate Potential Energy Surfaces for Hydrogen Abstraction Reactions: A Benchmark Study on the Xyg3 Doubly Hybrid Density Functional

An Accurate Full-Dimensional H4O Potential Energy Surface and Dynamics of an Exchange Reaction.

A New Ab Initio Based Global HOOH(13A″) Potential Energy Surface for the O(3P) + H2O(X1A1) ↔ OH(X2Π) + OH(X2Π) Reaction

Full-dimensional Global Potential Energy Surfaces Describing Abstraction and Exchange for the H + H2S Reaction.

A second-generation reactive empirical bond order (REBO) potential energy expression for hydrocarbons

Validating Experiments for the Reaction H2 + NH2- by Dynamical Calculations on an Accurate Full-Dimensional Potential Energy Surface.

Near Equivalence of Polarizability and Bond Order Flux Metrics for Describing Covalent Bond Rearrangements

Training Transferable Interatomic Neural Network Potentials for Reactive Chemistry: Improved Chemical Space Sampling

Newab Initiopotential Energy Surfaces for Cl(2P3/2, 2P1/2)+H2 Reaction

Potential Energy Surface and Quantum Dynamics Calculation of SH2− (2A′) Based on Ab Initio Scaled External Correlation Correction

Quantum mechanical fragment methods based on partitioning atoms or partitioning coordinates.