Abstract:The importance of Localized Molecular Orbitals in correlation treatments beyond mean-field calculation and in the illustration of chemical bonding can hardly be overstated. However, generation of orthonormal localized occupied MOs is significantly more straightforward than obtaining orthonormal localized virtual MOs. Orthonormal molecular orbitals allow facile use of highly efficient group theoretical methods for calculation of Hamiltonian matrix elements in multireference configuration interaction calculations and in quasidegenerate perturbation treatments, such as the GVVPT2 Theory. Moreover, localized MOs can elucidate qualitative understanding of bonding in molecules, in addition to high accuracy quantitative descriptions. We adopt the powers of the fourth moment cost function introduced by Jørgensen and coworkers. Because the fourth moment cost functions are prone to having multiple negative Hessian eigenvalues when starting from easily available canonical MOs, standard optimization algorithms can fail to obtain the orbitals of the virtual or partially occupied spaces. To overcome this drawback, we applied a Trust Region algorithm on an orthonormal Riemannian manifold with an approximate retraction from the tangent space built into first and second derivatives of the cost function. Moreover, the Riemannian Trust Region outer iterations were coupled to truncated conjugate gradient inner loops, which avoided any costly solutions of simultaneous linear equations or eigenvector/eigenvalue solutions. Numerical examples are provided on model systems, including the high connectivity H10 set in 1-, 2- and 3-dimensional arrangements, and on a chemically realistic description of cyclobutadiene, and the propargyl radical. In addition to demonstrating the algorithm on occupied and virtual blocks of orbitals, the method is also shown to work on the active space at the MCSCF level of theory.

The Use of Trust Regions in Kohn-Sham Total EnergyMinimization

A Trust Region Direct Constrained Minimization Algorithm for the Kohn-Sham Equation.

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

A Constrained Optimization Algorithm for Total Energy Minimization in Electronic Structure Calculations

On the Convergence of the Self-Consistent Field Iteration in Kohn--Sham Density Functional Theory

Minimizing the Kohn–Sham Total Energy for Periodic Systems

On the Analysis of the Discretized Kohn-Sham Density Functional Theory

Globally Convergent Composite-Step Trust-Region Framework for Model-Based Real-Time Optimization

A Story of Three Levels of Sophistication in SCF/KS-DFT Orbital Optimization Procedures

Riemannian Trust Region Method for Minimization of Fourth Central Moment for Localized Molecular Orbitals

Adaptive Regularized Self-Consistent Field Iteration with Exact Hessian for Electronic Structure Calculation.

Towards Translational Invariance of Total Energy with Finite Element Methods for Kohn-Sham Equation

Direct Unconstrained Optimization of Molecular Orbital Coefficients in Density Functional Theory

Riemannian Trust Region Methods for SC$^1$ Minimization

A trust-region augmented Hessian implementation for state-specific and state-averaged CASSCF wave functions

Exploiting the Hessian for a Better Convergence of the SCF-RDMFT Procedure

Computing the Self-Consistent Field in Kohn-Sham Density Functional Theory

Adaptive Local Basis Set for Kohn-Sham Density Functional Theory in a Discontinuous Galerkin Framework I: Total Energy Calculation

By-passing the Kohn-Sham equations with machine learning

D4FT: A Deep Learning Approach to Kohn-Sham Density Functional Theory

NeuralSCF: Neural network self-consistent fields for density functional theory