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

Minimum Energy Paths and Transition States by Curve Optimization

Locating Transition States by Variational Reaction Path Optimization with an Energy-Derivative-Free Objective Function

Minima Hopping Guided Path Search: An Efficient Method for Finding Complex Chemical Reaction Pathways

Entropic Path Sampling: Computational Protocol to Evaluate Entropic Profile along a Reaction Path

Optimization-based String Method for Finding Minimum Energy Path

Integrated Reaction Path Processing from Sampled Structure Sequences

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

Identifying a Feasible Transition Pathway between Two Conformational States for a Protein

Minimal Curvature Trajectories: Riemannian Geometry Concepts for Model Reduction in Chemical Kinetics

Accelerated Relaxation Engines for Optimizing to Minimum Energy Path

Diffusion Methods for Generating Transition Paths

Finding Reaction Pathways with Optimal Atomic Index Mappings

Efficient Chain-of-States Approach for Locating Transition State via Spherical Optimization.

Constrained Broyden Dimer Method with Bias Potential for Exploring Potential Energy Surface of Multistep Reaction Process.

Adaptive Learning Reaction Pathways Search

Modified String Method for Finding Minimum Energy Path

Exploration of Reaction Pathways and Chemical Transformation Networks

Transition Pathways in Complex Systems: Reaction Coordinates, Isocommittor Surfaces, and Transition Tubes

Physical Prior Mean Function-Driven Gaussian Processes Search for Minimum-Energy Reaction Paths with a Climbing-Image Nudged Elastic Band: A General Method for Gas-Phase, Interfacial, and Bulk-Phase Reactions

A Flat-bottom Elastic Network Model for Generating Improved Plausible Reaction Paths

Computationally Efficient Characterization of Potential Energy Surfaces Based on Fingerprint Distances