Steering Orbital Optimization out of Local Minima and Saddle Points Toward Lower Energy

Abstract

The general procedure underlying Hartree-Fock and Kohn-Sham density functional theory calculations consists in optimizing orbitals for a self-consistent solution of the Roothaan-Hall equations in an iterative process. It is often ignored that multiple self-consistent solutions can exist, several of which may correspond to minima of the energy functional. In addition to the difficulty sometimes encountered to converge the calculation to a self-consistent solution, one must ensure that the correct self-consistent solution was found, typically the one with the lowest electronic energy. Convergence to an unwanted solution is in general not trivial to detect and will deliver incorrect energy and molecular properties, and accordingly a misleading description of chemical reactivity. Wrong conclusions based on incorrect self-consistent field convergence are particularly cumbersome in automated calculations met in high-throughput virtual screening, structure optimizations, ab initio molecular dynamics, and in real-time explorations of chemical reactivity, where the vast amount of data can hardly be manually inspected. Here, we introduce a fast and automated approach to detect and cure incorrect orbital convergence, which is especially suited for electronic structure calculations on sequences of molecular structures. Our approach consists of a randomized perturbation of the converged electron density (matrix) intended to push orbital convergence to solutions that correspond to another stationary point (of potentially lower electronic energy) in the variational parameter space of an electronic wave function approximation.

arXiv comment12 pages, 12 figures

References in corpus (11)

• Big Data meets Quantum Chemistry Approximations: The $\Delta$-Machine Learning Approach
• Spin in Density-Functional Theory
• Heuristics-Guided Exploration of Reaction Mechanisms
• Projected Quasiparticle Theory for Molecular Electronic Structure
• Exploring the Quantum Speed Limit with Computer Games
• Interactive Chemical Reactivity Exploration
• Guaranteed convergence of the Kohn-Sham equations
• Real-time Quantum Chemistry
• Real-time feedback from iterative electronic structure calculations
• Accelerating Wave Function Convergence in Interactive Quantum Chemical Reactivity Studies
• Molecular propensity as a driver for explorative reactivity studies

Cited by in corpus (16)

• Exploration of Reaction Pathways and Chemical Transformation Networks
• The Exploration of Chemical Reaction Networks
• Variational density functional calculations of excited states via direct optimization
• Chemoton 2.0: Autonomous Exploration of Chemical Reaction Networks
• Calculation of Ligand Dissociation Energies in Large Transition-Metal Complexes
• Minimum Energy Paths and Transition States by Curve Optimization
• Excited states, symmetry breaking, and unphysical solutions in state-specific CASSCF theory
• Direct Energy Minimization Based on Exponential Transformation in Density Functional Calculations of Finite and Extended Systems
• SCINE -- Software for Chemical Interaction Networks
• Second-order self-consistent field algorithms: from classical to quantum nuclei
• Integrated Reaction Path Processing from Sampled Structure Sequences
• Holomorphic Density Functional Theory
• Weighted Active Space Protocol for Multireference Machine-Learned Potentials
• Heron: Visualizing and Controlling Chemical Reaction Explorations and Networks
• Direct minimization on the complex Stiefel manifold in Kohn-Sham density functional theory for finite and extended systems
• Variational Density Functional Calculations of Excited States: Conical Intersection and Avoided Crossing in Ethylene Bond Twisting