The diatomic neutral oxides and their ions, MO(0,+/-), M = Sc, Ti, Cr, and Mn, have been studied through multireference configuration interaction and coupled-cluster methods. With the purpose to paint a more comprehensive and detailed picture on these not so easily tamed systems, we have constructed complete potential energy curves for a large number of states of all MO(0,+/-)'s reporting structural and spectroscopic properties. Our overall results are in very good agreement with the, in general limited, experimental data. The always difficult to be pinpointed "nature of the chemical bond" becomes more recondite for these highly open ionic-covalent species. We have tried to give some answers as to the bonding interactions using simple valence-bond-Lewis diagrams in conjunction with Mulliken populations and the symmetry of the in situ atoms. It is our belief that, particularly for this kind of molecule, molecular orbital concepts are of limited help for a consistent rationalization of the bond formation.

We report the first optimum geometries and harmonic vibrational frequencies for the ring pentamer and several water hexamer (prism, cage, cyclic and two book) at the coupled-cluster including single, double, and full perturbative triple excitations (CCSD(T))/aug-cc-pVDZ level of theory. All five examined hexamer isomer minima previously reported by Møller-Plesset perturbation theory (MP2) are also minima on the CCSD(T) potential energy surface (PES). In addition, all CCSD(T) minimum energy structures for the n = 2-6 cluster isomers are quite close to the ones previously obtained by MP2 on the respective PESs, as confirmed by a modified Procrustes analysis that quantifies the difference between any two cluster geometries. The CCSD(T) results confirm the cooperative effect of the homodromic ring networks (systematic contraction of the nearest-neighbor (nn) intermolecular separations with cluster size) previously reported by MP2, albeit with O-O distances shorter by ~0.02 Å, indicating that MP2 overcorrects this effect. The harmonic frequencies at the minimum geometries were obtained by the double differentiation of the CCSD(T) energy using an efficient scheme based on internal coordinates that reduces the number of required single point energy evaluations by ~15% when compared to the corresponding double differentiation using Cartesian coordinates. Negligible differences between MP2 and CCSD(T) frequencies are found for the librational modes, while uniform increases of ~15 and ~25 cm(-1) are observed for the bending and "free" OH harmonic frequencies. The largest differences between CCSD(T) and MP2 are observed for the harmonic hydrogen bonded frequencies, for which the former produces larger absolute values than the latter. Their CCSD(T) redshifts from the monomer values (Δω) are smaller than the MP2 ones, due to the fact that CCSD(T) produces shorter elongations (ΔR) of the respective hydrogen bonded OH lengths from the monomer value with respect to MP2. Both the MP2 and CCSD(T) results for the hydrogen bonded frequencies were found to closely follow the relation -Δω = s · ΔR, with a rate of s = 20.2 cm(-1)/0.001 Å for hydrogen bonded frequencies with IR intensities >400 km/mol. The CCSD(T) harmonic frequencies, when corrected using the MP2 anharmonicities obtained from second order vibrational perturbation theory, produce anharmonic CCSD(T) estimates that are within <60 cm(-1) from the measured infrared (IR) active bands of the n = 2-6 clusters. Furthermore, the CCSD(T) harmonic redshifts (with respect to the monomer) trace the measured ones quite accurately. The energetic order between the various hexamer isomers on the PES (prism has the lowest energy) previously reported at MP2 was found to be preserved at the CCSD(T) level, whereas the inclusion of anharmonic corrections further stabilizes the cage among the hexamer isomers.

We establish a new estimate for the binding energy between two benzene molecules in the parallel-displaced (PD) conformation by systematically converging (i) the intra- and intermolecular geometry at the minimum, (ii) the expansion of the orbital basis set, and (iii) the level of electron correlation. The calculations were performed at the second-order Møller-Plesset perturbation (MP2) and the coupled cluster including singles, doubles, and a perturbative estimate of triples replacement [CCSD(T)] levels of electronic structure theory. At both levels of theory, by including results corrected for basis set superposition error (BSSE), we have estimated the complete basis set (CBS) limit by employing the family of Dunning's correlation-consistent polarized valence basis sets. The largest MP2 calculation was performed with the cc-pV6Z basis set (2772 basis functions), whereas the largest CCSD(T) calculation was with the cc-pV5Z basis set (1752 basis functions). The cluster geometries were optimized with basis sets up to quadruple-ζ quality, observing that both its intra- and intermolecular parts have practically converged with the triple-ζ quality sets. The use of converged geometries was found to play an important role for obtaining accurate estimates for the CBS limits. Our results demonstrate that the binding energies with the families of the plain (cc-pVnZ) and augmented (aug-cc-pVnZ) sets converge [within <0.01 kcal/mol for MP2 and <0.15 kcal/mol for CCSD(T)] to the same CBS limit. In addition, the average of the uncorrected and BSSE-corrected binding energies was found to converge to the same CBS limit much faster than either of the two constituents (uncorrected or BSSE-corrected binding energies). Due to the fact that the family of augmented basis sets (especially for the larger sets) causes serious linear dependency problems, the plain basis sets (for which no linear dependencies were found) are deemed as a more efficient and straightforward path for obtaining an accurate CBS limit. We considered extrapolations of the uncorrected (ΔE) and BSSE-corrected (ΔEcp) binding energies, their average value (ΔEave), as well as the average of the latter over the plain and augmented sets (ΔẼave) with the cardinal number of the basis set n. Our best estimate of the CCSD(T)/CBS limit for the π-π binding energy in the PD benzene dimer is De = -2.65 ± 0.02 kcal/mol. The best CCSD(T)/cc-pV5Z calculated value is -2.62 kcal/mol, just 0.03 kcal/mol away from the CBS limit. For comparison, the MP2/CBS limit estimate is -5.00 ± 0.01 kcal/mol, demonstrating a 90% overbinding with respect to CCSD(T). The spin-component-scaled (SCS) MP2 variant was found to closely reproduce the CCSD(T) results for each basis set, while scaled opposite spin (SOS) MP2 yielded results that are too low when compared to CCSD(T).