Constrained nuclear-electronic orbital second-order Moller-Plesset perturbation theory

arXiv:2602.10267v1 Announce Type: new
Abstract: A multicomponent second-order Mo ller-Plesset perturbation theory (MP2) method is derived and implemented within the constrained nuclear-electronic orbital (cNEO) framework from a multicomponent generalization of the Hylleraas functional. The cNEO-MP2 method includes electronic-nuclear and nuclear correlation in the calculation of vibrationally averaged molecular properties, and is the first post Hartree-Fock wavefunction-based cNEO method. Nuclear quantum effects like vibrational averaging, isotopic effects, and zero-point energy can be captured in a single calculation or geometry optimization with cNEO-MP2, eliminating the need to perform costly subsequent calculations to determine higher order force constants as required with many existing methods used to determine vibrational effects upon molecular properties. The cNEO-MP2 method is benchmarked on a test set of diatomic and small polyatomic molecules and ions. Herein, we present internuclear distances, bond angles, potential energy surfaces, and vibrational frequencies calculated with the cNEO-MP2 method to demonstrate that it correctly captures the effects of nuclear vibrational motion upon molecular properties.