Analytical optimization of orbital exponents in Gaussian-type functions for molecular systems based on MCSCF and MP2 levels of fully variational molecular orbital method

We have analyzed the basis function series in molecular systems by optimization of orbital exponents in Gaussian-type functions (GTFs) including the electron correlation effects with multiconfiguration self-consistent field (MCSCF) and M?ller?CPlesset second-order perturbation (MP2) methods. First, we have derived and implemented the gradient formulas of MCSCF and MP2 energies with respect to GTF exponent, as well as GTF center and nuclear geometry, based on the fully variational molecular orbital (FVMO) method. Second, we have applied these electron-correlated FVMO methods to H₂, LiH, and hydrocarbon (CH₄, C₂H₆, C₂H₄, and C₂H₂) molecules. We have clearly demonstrated that the optimized exponent values with electron-correlated methods are different from those with simple Hartree?CFock method, since adequate basis functions for adequate virtual orbitals are indispensable to describe the accurate wave function and geometry for electron-correlated calculations.