Finding reaction pathways with optimal atomic index mappings
arXiv:1906.06077 · doi:10.1103/PhysRevLett.123.206102
Abstract
Finding complex reaction and transformation pathways, involving many intermediate states, is in general not possible on the DFT level with existing simulation methods due to the very large number of required energy and force evaluations. This is due to a large extent to the fact that for complex reactions, it is not possible to determine which atom in the educt is mapped onto which atom in the product. Trying out all possible atomic index mappings is not feasible because of the factorial increase in the number of possible mappings. By using a penalty function that is invariant under index permutations, we can bias the potential energy surface in such a way that it obtains the characteristics of a structure seeker whose global minimum is the product. By performing a Minima Hopping based global optimization on this biased potential energy surface we can rapidly find intermediate states that lead into the global minimum. Based on this information we can then extract the full reaction pathway. We first demonstrate for a benchmark system, namely LJ38 that our method allows to preferentially find intermediate states that are relevant for the lowest energy reaction pathway and that we, therefore, need a much smaller number of intermediate states than previous methods to find the lowest energy reaction pathway. Finally, we apply the method to two real systems, C60 and C20H20 and show that the found reaction pathway contains valuable information on how the system can be synthesized.
References in corpus (14)
- Escaping free-energy minima
- Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms
- Daubechies wavelets as a basis set for density functional pseudopotential calculations
- The double-funnel energy landscape of the 38-atom Lennard-Jones cluster
- Exploration of Reaction Pathways and Chemical Transformation Networks
- Phase changes in 38 atom Lennard-Jones clusters. I: A parallel tempering study in the canonical ensemble
- Evolution of the Potential Energy Surface with Size for Lennard-Jones Clusters
- Metrics for measuring distances in configuration spaces
- Norm-conserving pseudopotentials with chemical accuracy compared to all-electron calculations
- Predicting polymorphism in molecular crystals using orientational entropy
- A Bell-Evans-Polanyi principle for molecular dynamics trajectories and its implications for global optimization
- Efficient moves for global geometry optimization methods and their application to binary systems
- Minima Hopping Guided Path Search: An Efficient Method for Finding Complex Chemical Reaction Pathways
- Stabilized Quasi-Newton Optimization of Noisy Potential Energy Surfaces