An energy decomposition analysis for second-order Møller-Plesset perturbation theory based on absolutely localized molecular orbitals.

Abstract

An energy decomposition analysis (EDA) of intermolecular interactions is proposed for second-order Møller-Plesset perturbation theory (MP2) based on absolutely localized molecular orbitals (ALMOs), as an extension to a previous ALMO-based EDA for self-consistent field methods. It decomposes the canonical MP2 binding energy by dividing the double excitations that contribute to the MP2 wave function into classes based on how the excitations involve different molecules. The MP2 contribution to the binding energy is decomposed into four components: frozen interaction, polarization, charge transfer, and dispersion. Charge transfer is defined by excitations that change the number of electrons on a molecule, dispersion by intermolecular excitations that do not transfer charge, and polarization and frozen interactions by intra-molecular excitations. The final two are separated by evaluations of the frozen, isolated wave functions in the presence of the other molecules, with adjustments for orbital response. Unlike previous EDAs for electron correlation methods, this one includes components for the electrostatics, which is vital as adjustment to the electrostatic behavior of the system is in some cases the dominant effect of the treatment of electron correlation. The proposed EDA is then applied to a variety of different systems to demonstrate that all proposed components behave correctly. This includes systems with one molecule and an external electric perturbation to test the separation between polarization and frozen interactions and various bimolecular systems in the equilibrium range and beyond to test the rest of the EDA. We find that it performs well on these tests. We then apply the EDA to a halogen bonded system to investigate the nature of the halogen bond.