Abstract
The density matrix form of Hartree-Fock perturbation theory is developed for the case in which the basis functions themselves are perturbation-dependent. Energy expressions are derived, through second order, for both single and double perturbations.
The theory is applied in the calculation of electric dipole polarizabilities and hyperpolarizabilities for atoms (He, Be) and molecules (H2, LiH), with excellent results. The high computational efficiency of the method is discussed and possibilities of further development are outlined.
Notes
This research was partly supported by the Institute of Low Temperatures and Structure Research of the Polish Academy of Sciences.