We present the results of an ab initio study of the interaction of electronically excited NO(A (2)Sigma(+)) with rare gas (Rg) atoms. The bound states of each NO(A)-Rg species are determined from potential energy surfaces calculated at the RCCSD(T) level of theory. Making use of the NO(X (2)Pi)-Rg vibrational wavefunctions, we then simulate electronic spectra. For NO-Kr and NO-Xe we obtain good qualitative agreement with the previously published experimental spectra. For NO-Ar, the shallowness of the surface gives rise to agreement that is less satisfactory, but a global scaling provides better qualitative agreement. The assignment of the spectra is far from straightforward and is only possible with guidance from the calculated energies and wavefunctions of the energy levels of the complex. Previous assignments are discussed in the light of this conclusion.

Intermolecular potentials for the three lowest multiplet states (singlet, triplet and quintet) of the O2(3Sigma)-O2(3Sigma) dimer have been investigated in detail by means of high level ab initio calculations. The methods used include MRCI, ACPF, CASPT2, using different active spaces and basis sets. The results for the quintet state are compared with benchmark CCSD(T) calculations. As expected, the former methods do not account accurately for dispersion interactions, although the CASPT2 method performs better than the CI based ones. On the other hand, it is shown that highly correlated methods are necessary to accurately describe the splittings among the multiplet states. We propose to obtain singlet and triplet interaction potentials by combining CCSD(T) quintet potentials and multiconfigurational singlet-quintet and triplet-quintet splittings, respectively. The calculated splittings are quite stable regarding the method employed, except for the well region of the singlet and triplet states within the rectangular configuration, which corresponds to the absolute minima of these multiplet states. Nevertheless, we have been able to assess adequate upper and lower bounds to the interaction potential for this particular region.

Halogen bonds have received a great deal of attention in recent years. Their properties, sometimes paralleled with those of hydrogen bonds, have not yet been fully understood. In this work, we investigate the nature of the intermolecular interactions between Cl(2) and Br(2) with water. Our analysis of several features of MP2/aug-cc-pVDZ-optimized stable clusters with different number and arrangement of water molecules shows that two different kinds of halogen-water coordination patterns are involved in the stability and properties found for these systems: halogen bonds (X-X...O) and halogen-hydrogen interactions, (X-X...H-O-H). Both types of interactions result in a large polarization of the halogen molecule, which leads to important cooperative effects on these structures. Although the general structural aspects of these clusters can be understood in terms of dipole-quadrupole forces at long range, where it is the dominant term, the SAPT analysis shows that factors such as polarization of pi densities and dispersion become increasingly important close to equilibrium. In particular, we show that the halogen-hydrogen interactions are weaker than halogen-oxygen interactions mainly due to the electrostatic and dispersion forces. We also calculate vibrational and electronic shifts that should be helpful for the interpretation of experimental results and for investigating the microsolvation phenomena for halogens in an aqueous environment.

Completely ab initio global potential energy surfaces (PESs) for the singlet and triplet spin multiplicities of rigid O(2)((3)Σ(g)(-))+O(2)((3)Σ(g)(-)) are reported for the first time. They have been obtained by combining an accurate restricted coupled cluster theory with singles, doubles, and perturbative triple excitations [RCCSD(T)] quintet potential [Bartolomei et al., J. Chem. Phys. 128, 214304 (2008)] with complete active space second order perturbation theory (CASPT2) or, alternatively, multireference configuration interaction (MRCI) calculations of the singlet-quintet and triplet-quintet splittings. Spherical harmonic expansions, containing a large number of terms due to the high anisotropy of the interaction, have been built from the ab initio data. The radial coefficients of these expansions are matched at long range distances with analytical functions based on recent ab initio calculations of the electric properties of the monomers [M. Bartolomei, E. Carmona-Novillo, M. I. Hernández, J. Campos-Martínez, and R. Hernández-Lamoneda, J. Comput. Chem. (2010) (in press)]. The singlet and triplet PESs obtained from either RCCSD(T)-CASPT2 or RCCSD(T)-MRCI calculations are quite similar, although quantitative differences appear in specific terms of the expansion. CASPT2 calculations are the ones giving rise to larger splittings and more attractive interactions, particularly in the region of the absolute minima (in the rectangular D(2h) geometry). The new singlet, triplet, and quintet PESs are tested against second virial coefficient B(T) data and, their spherically averaged components, against integral cross sections measured with rotationally hot effusive beams. Both types of multiconfigurational approaches provide quite similar results, which, in turn, are in good agreement with the measurements. It is found that discrepancies with the experiments could be removed if the PESs were slightly more attractive. In this regard, the most attractive RCCSD(T)-CASPT2 PESs perform slightly better than the RCCSD(T)-MRCI counterpart.

A new potential energy surface ͑PES͒ for the quintet state of rigid O 2 ͑ 3 ⌺ g − ͒ +O 2 ͑ 3 ⌺ g − ͒ has been obtained using restricted coupled-cluster theory with singles, doubles, and perturbative triple excitations ͓RCCSD͑T͔͒. A large number of relative orientations of the monomers ͑65͒ and intermolecular distances ͑17͒ have been considered. A spherical harmonic expansion of the interaction potential has been built from the ab initio data. It involves 29 terms, as a consequence of the large anisotropy of the interaction. The spherically averaged term agrees quite well with the one obtained from analysis of total integral cross sections. The absolute minimum of the PES corresponds to the crossed ͑D 2d ͒ structure ͑X shape͒ with an intermolecular distance of 6.224 bohrs and a well depth of 16.27 meV. Interestingly, the PES presents another ͑local͒ minimum close in energy ͑15.66 meV͒ at 6.50 bohrs and within a planar skewed geometry ͑S shape͒. We find that the origin of this second structure is due to the orientational dependence of the spin-exchange interactions which break the spin degeneracy and leads to three distinct intermolecular PESs with singlet, triplet, and quintet multiplicities. The lowest vibrational bound states of the O 2 -O 2 dimer have been obtained and it is found that they reflect the above mentioned topological features of the PES: The first allowed bound state for the 16 O isotope has an X structure but the next state is just 0.12 meV higher in energy and exhibits an S shape.