Improved initial guess for minimum energy path calculations

Abstract

A method is presented for generating a good initial guess of a transition path between given initial and final states of a system without evaluation of the energy. An objective function surface is constructed using an interpolation of pairwise distances at each discretization point along the path and the nudged elastic band method then used to find an optimal path on this image dependent pair potential (IDPP) surface. This provides an initial path for the more computationally intensive calculations of the true minimum energy path using some method of choice for evaluating the energy and atomic forces, for example by ab initio or density functional theory. The optimal path on the IDPP surface is significantly closer to the true minimum energy path than a linear interpolation of the Cartesian coordinates and, therefore, reduces the number of iterations needed to reach convergence and averts divergence in the electronic structure calculations when atoms are brought too close to each other in the initial path. The method is illustrated with three examples: (1) rotation of a methyl group in an ethane molecule, (2) an exchange of atoms in an island on a crystal surface, and (3) an exchange of two Si-atoms in amorphous silicon. In all three cases, the computational effort in finding the minimum energy path with DFT was reduced by a factor ranging from 50 % to an order of magnitude by using an IDPP path as the initial path. The time required for parallel computations was reduced even more because of load imbalance when linear interpolation of Cartesian coordinates was used.

Cited by in corpus (41)

• QuantumATK: An integrated platform of electronic and atomic-scale modelling tools
• NWChem: Past, Present, and Future
• SIESTA: recent developments and applications
• GPAW: An open Python package for electronic-structure calculations
• Method for finding mechanism and activation energy of magnetic transitions, applied to skyrmion and antivortex annihilation
• Nudged elastic band calculations accelerated with Gaussian process regression
• ATK-ForceField: A New Generation Molecular Dynamics Software Package
• Predicting the Activity and Selectivity of Bimetallic Metal Catalysts for Ethanol Reforming using Machine Learning
• Gaussian Process Regression for Transition State Search
• Scaled and Dynamic Optimizations of Nudged Elastic Bands
• Minimum Energy Paths and Transition States by Curve Optimization
• Minimum energy path calculations with Gaussian process regression
• Analytical Gradients for Projection-Based Wavefunction-in-DFT Embedding
• Metadynamics sampling in atomic environment space for collecting training data for machine learning potentials
• Removing External Degrees of Freedom from Transition-State Search Methods using Quaternions
• Effects of hole self-trapping by polarons on transport and negative bias illumination stress in amorphous-IGZO
• Two-dimensional 5d multiferroic W3Cl8: breathing Kagome lattice and tunable magneto-optical Kerr effect
• Molecular Quantum Chemical Data Sets and Databases for Machine Learning Potentials
• Silver electrodes are highly selective for CO in CO${}_2$ electroreduction due to interplay between voltage dependent kinetics and thermodynamics
• Understanding Memristive Behavior: An Atomistic Study of the Influence of Grain Boundaries on Surface and Out-of-Plane Diffusion of Metallic Atoms
• Metatensor and metatomic: foundational libraries for interoperable atomistic machine learning
• Does Hessian Data Improve the Performance of Machine Learning Potentials?
• Towards targeted kinetic trapping of organic-inorganic interfaces: A computational case study
• Dynamic control of octahedral rotation in perovskites by defect engineering
• Minimum energy path for the nucleation of misfit dislocations in Ge/Si(001) heteroepitaxy
• Anticipating the Selectivity of Intramolecular Cyclization Reaction Pathways with Neural Network Potentials
• GPR_calculator: An On-the-Fly Surrogate Model to Accelerate Massive Nudged Elastic Band Calculations
• Breaking scaling relations with inverse catalysts: a machine learning exploration of trends in $\mathrm{C}{\mathrm{O}}_2$ hydrogenation energy barriers
• Structural Feature in Dynamical Processes Accelerated Transition State Calculations
• Implicit Neural Representations for Chemical Reaction Paths
• Leveraging transfer learning for accurate estimation of ionic migration barriers in solids
• Initial estimate for minimum energy pathways and transition states using velocities in internal coordinates
• Diffusion Mechanisms in Li${}_{0.5}$CoO${}_2$ -- A Computational Study
• Locating Ab Initio Transition States via Approximate Geodesics on Machine Learned Potential Energy Surfaces
• Unraveling Mn intercalation and diffusion in NbSe${}_2$ bilayers through DFTB simulations
• Structural Control of Atomic Silicon Wires
• Reproducible Orchestration of Best Practices for Reaction Path Optimization with the Nudged Elastic Band
• From AC-STEM Image to 3D Structure: A Systematic Analysis of Au55 nanocluster
• Efficient Chain-of-States Approach for Locating Transition State via Spherical Optimization
• Formation of methane and cyclohexane through the hydrogenation of toluene
• Indirect Mechanism of Au adatom Diffusion on the Si(100) Surface