The ammonia dimer equilibrium dissociation energy: convergence to the basis set limit at the correlated level

Abstract

The low-energy region of the intermolecular potential energy hypersurface (PES) of the ammonia dimer was studied at the level of second-order Moller-Plesset perturbation theory (MP2) using a very large basis set. Individual minima were located on the PES employing the counterpoise (CP) correction to account for the basis set superposition error (BSSE). Apart from these canonical MP2 calculations local MP2 (LMP2) calculations were performed. For the latter the BSSE at the correlated level is inherently absent by virtue of the local truncation of the virtual space. Results from canonical and local MP2 calculations are compared and the reliability of the LMP2 method for intermolecular complexes and clusters is discussed. The canonical MP2 calculations predicted five minimum structures, the four most stable ones lying energetically very close. For these four structures single point MP2 energy calculations with a further extended basis set (1024 functions for the ammonia dimer) were performed. The equilibrium dissociation energies so obtained are close to the one-particle basis set limit, as illustrated by a remaining BSSE of less than 0.2 kJ mol⁻¹. The geometry optimizations at the LMP2 level, using the three most stable canonical MP2 structures as initial geometries, all collapsed to a single minimum corresponding to an asymmetric structural arrangement. A canonical MP2 single point calculation, at that geometry, revealed that the LMP2 minimum structure is virtually as stable as the lowest minima on the canonical MP2 PES. Based on these calculations the global minimum of the ammonia dimer was assigned to a part of the PES represented by an asymmetric structure with an equilibrium dissociation energy of 13.5±0.3 kJ mol⁻¹