Gradient of molecular Hartree–Fock–Bogoliubov energy with a linear combination of atomic orbital quasiparticle wave functions

Abstract

The analytical gradient for the atomic-orbital-based Hartree-Fock-Bogoliubov (HFB) energy functional, the modified form of which was proposed by Staroverov and Scuseria to account for the static electron correlation [J. Chem. Phys. 117, 11107 (2002)], is derived. Interestingly, the Pulay force for the HFB energy is expressed with the same formula as that for the Hartree-Fock method. The efficiency of the present HFB energy gradient is demonstrated in the geometry optimizations of conjugated and biradical systems. The geometries optimized by using the HFB method with the appropriate factor ζ, which controls the degree of static correlation included, are found to show good agreement with those obtained by using a complete active-space self-consistent field method, although they are significantly dependent on ζ.