Global ab initio potential energy surface for the O2((3)Σ(g)(-)) + N2((1)Σ(g)(-)) interaction. Applications to the collisional, spectroscopic, and thermodynamic properties of the complex.

A detailed characterization of the interaction between the most abundant molecules in air is important for the understanding of a variety of phenomena in atmospherical science. A completely ab initio global potential energy surface (PES) for the O2((3)Σg(­)) + N2((1)Σg(+)) interaction is reported for the first time. It has been obtained with the symmetry-adapted perturbation theory utilizing a density functional description of monomers [SAPT(DFT)] extended to treat the interaction involving high-spin open-shell complexes. The computed interaction energies of the complex are in a good agreement with those obtained by using the spin-restricted coupled cluster methodology with singles, doubles, and noniterative triple excitations [RCCSD(T)]. A spherical harmonics expansion of the interaction potential containing a large number of terms due to the anisotropy of the interaction has been built from the ab initio data. The expansion coefficients, which are functions of the intermolecular distance, are matched in the long-range with the analytical functions based on the recent ab initio calculations of the electric properties of the monomers [M. Bartolomei et al. J. Comput. Chem. 2011, 32, 279]. The PES is tested against the second virial coefficient B(T) data and the integral cross sections measured with rotationally hot effusive beams, leading in both cases to a very good agreement. The lowest lying states of the complex have been computed and relevant spectroscopic features of the interacting complex are reported. A comparison with a previous experimentally derived PES is also provided.

Ab initio rovibrational structure of the lowest singlet state of O2-O2.

Rovibrational bound states of the O(2)((3)Σ(g)(-), v = 0)-O(2)((3)Σ(g)(-), v = 0) dimer in its singlet electronic state have been obtained by solving the time-independent Schrödinger equation for the nuclear degrees of freedom. We have employed two different ab initio potential energy surfaces, based on high level multiconfigurational methods, which are expected to give upper and lower bounds for the real values of the interaction. Results are compared with spectroscopy experiments as well as with calculations using other semi ab initio and empirical interaction potentials. For the two ab initio potentials studied here, the ground vibrational state has a rectangular geometry and behaves as a semi-rigid molecule. The associated rotational constant is found in very good agreement with high resolution spectra. However, the computed dissociation energy and the frequency of the torsion mode are larger than previous experimental determinations, and possible reasons for these discrepancies are discussed. On the other hand, we have computed the splitting between the rovibrational states of the singlet and triplet electronic states and have found a fair agreement with measurements of the dimer spectra in a solid rare gas host.

Long-range interaction for dimers of atmospheric interest: dispersion, induction and electrostatic contributions for O(2)-O(2), N(2)-N(2) and O(2)-N(2).

Electric multipole moments, static dipole polarizabilities, and dynamic dipole, quadrupole, and mixed dipole-octupole polarizabilities of molecular oxygen and nitrogen in their ground electronic states have been obtained by means of high level multiconfigurational ab initio calculations. From these properties, we have obtained electrostatic, dispersion, and induction coefficients for the long-range interactions of the O(2)-O(2) , N(2)-N(2) , and O(2)-N(2) dimers. Our data is a comprehensive and consistent set that for N(2)-N(2) shows a very good agreement with previous accurate calculations, whereas for quantities involving open-shell O(2) represents a considerable improvement over previous estimations. Moreover, the long-range interaction is analyzed and compared for the different interacting partners. It is found that the C(8) dispersion interaction plays a nonnegligible role and that the induction component is only important for a detailed description of the highest order anisotropy terms in the spherical harmonics expansion of the long-range potential. It is also found that the total long-range interaction is quite similar in O(2)-O(2) and O(2)-N(2) , and that differences with N(2)-N(2) are mainly because of the important role of the electrostatic interaction in that dimer. Comparison with high level supermolecular calculations indicates that the present long-range potentials are accurate for intermolecular distances larger than about 15 bohr.

Global ab initio potential energy surfaces for the O2(3Σg-)+O2(3Σg-) interaction.

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.

Accurate ab initio intermolecular potential energy surface for the quintet state of the O2(3Sigma(g)-)-O2(3Sigma(g)-) dimer.

A new potential energy surface (PES) for the quintet state of rigid O(2)((3)Sigma(g)(-)) + O(2)((3)Sigma(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.

Penetration Barrier of Water through Graphynes' Pores: First-Principles Predictions and Force Field Optimization.

Graphynes are novel two-dimensional carbon-based materials that have been proposed as molecular filters, especially for water purification technologies. We carry out first-principles electronic structure calculations at the MP2C level of theory to assess the interaction between water and graphyne, graphdiyne, and graphtriyne pores. The computed penetration barriers suggest that water transport is unfeasible through graphyne while being unimpeded for graphtriyne. For graphdiyne, with a pore size almost matching that of water, a low barrier is found that in turn disappears if an active hydrogen bond with an additional water molecule on the opposite side of the opening is considered. Thus, in contrast with previous determinations, our results do not exclude graphdiyne as a promising membrane for water filtration. In fact, present calculations lead to water permeation probabilities that are 2 orders of magnitude larger than estimations based on common force fields. A new pair potential for the water-carbon noncovalent component of the interaction is proposed for molecular dynamics simulations involving graphdiyne and water.

The intermolecular potentials of the O2-O2 dimer: a detailed ab initio study of the energy splittings for the three lowest multiplet states.

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.

On the dynamics of the H+ +D2(v=0,j=0)-->HD+D + reaction: a comparison between theory and experiment.

The H+ +D2(v=0,j=0)-->HD+D + reaction has been theoretically investigated by means of a time independent exact quantum mechanical approach, a quantum wave packet calculation within an adiabatic centrifugal sudden approximation, a statistical quantum model, and a quasiclassical trajectory calculation. Besides reaction probabilities as a function of collision energy at different values of the total angular momentum, J, special emphasis has been made at two specific collision energies, 0.1 and 0.524 eV. The occurrence of distinctive dynamical behavior at these two energies is analyzed in some detail. An extensive comparison with previous experimental measurements on the Rydberg H atom with D2 molecules has been carried out at the higher collision energy. In particular, the present theoretical results have been employed to perform simulations of the experimental kinetic energy spectra.

Intermolecular interaction potentials for the Ar-C2H2, Kr-C2H2, and Xe-C2H2 weakly bound complexes: information from molecular beam scattering, pressure broadening coefficients, and rovibrational spectroscopy.

Integral cross sections and pressure broadening coefficients have been measured for the acetylene-krypton complex, by molecular beam scattering and by high resolution IR spectroscopy, respectively. A new potential energy surface (PES) is proposed to describe structure and dynamical properties of this prototypical weakly bound complex. The PES has been parametrized exploiting a novel atom-bond pairwise additive scheme and has been fitted to the experimental data. A similar PES has been obtained for the acetylene-xenon system by a proper scaling of the interaction parameters of the krypton case, based on empirical considerations. These PESs together with that recently proposed by the same authors [J. Phys. Chem. 109, 8471 (2005)] for the acetylene-argon case have been employed for close coupling calculations of the pressure broadening cross sections and for a characterization of the rovibrational structure of the complexes.

Quasiclassical trajectory study of reactive and dissociative processes in H2+H2: comparison with quantum-mechanical calculations.

Quasiclassical trajectory calculations have been carried out for H(2)(v(1)=high)+H(2)(v(2)=low) collisions within a three degrees of freedom model where five different geometries of the colliding complex were considered. Within this approach, probabilities for different competitive processes are studied: four center reaction, collision induced dissociation, reactive dissociation, and three-body complex formation. The purpose is to compare in detail with equivalent quantum-mechanical wave packet calculations [Bartolomei et al., J. Chem. Phys 122, 064305 (2005)], especially the behavior of the probabilities near reaction thresholds. Quasiclassical calculations compare quite well with the quantum-mechanical ones for collision induced dissociation as well as for the four center reaction, although quantum effects become very important near thresholds, particularly for lower v(1)'s and for the four center process. Less quantitative agreement is found for reactive dissociation and three-body complex formation. It is found that most quantum effects are due to differences between quantum and classical vibrational distributions of H(2)(v(1)=high). Zero point energy violation has been found in the classical reactive-dissociative probabilities. Extension of these findings to full-dimensional treatments is examined.