Fdez Galván I, Weßling D, Lindh R
J Chem Theory Comput 21 (24) 12674-12685 [2025-12-23; online 2025-12-07]
We present enhancements to the S-GEK/RVO method for self-consistent field (SCF) orbital optimization, aimed at improving computational efficiency and robustness. Building on a gradient-enhanced Kriging surrogate model and restricted-variance optimization, we introduce three key modifications: (i) a cost-effective subspace expansion using r-GDIIS or BFGS displacement predictions, (ii) a systematic undershoot mitigation strategy in flat energy regions, and (iii) rigorous coordinate and gradient transformations consistent with the exponential parametrization of orbital rotations. Benchmarking across an extensive set of molecular systems─including organic molecules, radicals, and transition-metal complexes─demonstrates that the new S-GEK/RVO variants consistently outperform the default (in OpenMolcas) r-GDIIS method in iteration count, convergence reliability, and wall time. These improvements make S-GEK/RVO a competitive alternative for SCF optimization and suggest broader applicability to other orbital optimization and localization problems.
PubMed 41355174
DOI 10.1021/acs.jctc.5c01714
Crossref 10.1021/acs.jctc.5c01714