A black-box, general purpose quadratic self-consistent field code with and without Cholesky Decomposition of the two-electron integrals
arXiv:2106.04836 · doi:10.1080/00268976.2021.1974590
Abstract
We present the implementation of a quadratically convergent Self-consistent field (QCSCF) algorithm based on an adaptive trust-radius optimization scheme for restricted open-shell Hartree-Fock (ROHF), restricted Hartree-Fock (RHF), and unrestricted Hartree-Fock (UHF) references. The algorithm can exploit Cholesky decomposition (CD) of the two-electron integrals to allow calculations on larger systems. The most important feature of the QCSCF code lies in its black-box nature -- probably the most important quality desired by a generic user. As shown for pilot applications, it does not require one to tune the self-consistent field (SCF) parameters (damping, Pulay's DIIS, and other similar techniques) in difficult-to-converge molecules. Also, it can be used to obtain a very thigh convergence with extended basis set - a situation often needed when computing high-order molecular properties - where the standard SCF algorithm starts to oscillate. Nevertheless, trouble may appear even with a QCSCF solver. In this respect, we discuss what can go wrong, focusing on the multiple UHF solutions of ortho-benzyne
References in corpus (1)
Cited by in corpus (6)
- Heptazine, Cyclazine, and Related Compounds: Chemically-Accurate Estimates of the Inverted Singlet-Triplet Gap
- Cholesky decomposition of complex two-electron integrals over GIAOs: Efficient MP2 computations for large molecules in strong magnetic fields
- Cholesky decomposition of two-electron integrals in quantum-chemical calculations with perturbative or finite magnetic fields using gauge-including atomic orbitals
- Economical Quasi-Newton Self Consistent Field Solver
- A trust-region augmented Hessian implementation for state-specific and state-averaged CASSCF wave functions
- A Reusable Library for Second-Order Orbital Optimization Using the Trust Region Method