Abstract:In reaction path optimization, such as the calculation of a minimum energy path (MEP) between given reactant and product configurations of atoms, it is advantageous to start with an initial guess where close proximity of atoms is avoided and bonds are not unnecessarily broken only to be reformed later. When the configurations of the atoms are described with Cartesian coordinates, a linear interpolation between the endpoints can be problematic, and a better option is provided by the so-called image dependent pair potential (IDPP) approach where interpolated pairwise distances are generated to form an objective function that can be used to generate an improved initial path. When started with a linear interpolation, this method can, however, still lead to unnecessary bond breaking in, for example, reactions where a molecular subgroup undergoes significant rotation. In the method presented here, this problem is addressed by constructing the path gradually, introducing images sequentially starting from the vicinity of the endpoints while the distance between images in the central region is larger. The distribution of images is controlled by systematically scaling the tightness of springs acting between the images until the desired number of images is obtained and they are evenly spaced. This procedure generates an initial path on the IDPP surface, a task that requires insignificant computational effort, as no evaluation of the energy of the system is needed. The calculation of the MEP, typically using electronic structure calculations, is then subsequently carried out in a way that makes efficient use of parallel computing with the nudged elastic band method. Several examples of reactions are given where the linear interpolation IDPP (LI-IDPP) method yields problematic paths with unnecessary bond breaking in some of the intermediate images, while the sequential IDPP (S-IDPP) method yields pat

Potential Energy Interpolation with Target-Customized Weighting Coordinates: Application to Excited-State Dynamics of Photoactive Yellow Protein Chromophore in Water

Energy Transition Rate At Peptide-Bond Using Stochastic Averaging Method

Interpolating many-body wave functions for accelerated molecular dynamics on the near-exact electronic surface

A combined on-the-fly/interpolation procedure for evaluating energy values needed in molecular simulations

Advanced Potential Energy Surfaces for Molecular Simulation

Data Quality, Data Sampling and Data Fitting: A Tutorial Guide for Constructing Full-Dimensional Accurate Potential Energy Surfaces (PESs) of Molecules and Reactions

Potential energy surface interpolation with neural networks for instanton rate calculations

Electronic and Steric Tuning of Excited State Proton Transfer in Water

Interpolation and extrapolation of global potential energy surfaces for polyatomic systems by Gaussian processes with composite kernels

An overview of nonadiabatic dynamics simulations methods, with focus on the direct approach versus the fitting of potential energy surfaces

Improved description of environment and vibronic effects with electrostatically embedded ML potentials

Efficient Sampling over Rough Energy Landscapes with High Barriers: A Combination of Metadynamics with Integrated Tempering Sampling.

Refining Potential Energy Surface through Dynamical Properties via Differentiable Molecular Simulation

Application of Interpolated Potential Energy Surfaces to Quantum Reactive Scattering

Combined Self-Consistent-Field and Spin-Flip Tamm-Dancoff Density Functional Approach to Potential Energy Surfaces for Photochemistry

Improved Initialization of Optimal Path Calculations Using Sequential Traversal over the Image Dependent Pair Potential Surface

Instanton rate constant calculations using interpolated potential energy surfaces in non-redundant, rotationally and translationally invariant coordinates

Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation

Transfer-Learned Potential Energy Surfaces: Towards Microsecond-Scale Molecular Dynamics Simulations in the Gas Phase at CCSD(T) Quality

Multistate, polarizable QM/MM embedding scheme based on the direct reaction field method: Solvatochromic shifts, analytical gradients and optimization of conical intersections in solution

Quasiclassical Trajectory Simulation as a Protocol to Build Locally Accurate Machine Learning Potentials