Rotational excitation of NS+ by H2 revisited: A new global potential energy surface and rate coefficients.

Due to the lack of specific collisional data, the abundance of NS+ in cold dense interstellar clouds was determined using collisional rate coefficients of CS as a substitute. To better understand the chemistry of sulfur in the interstellar medium, further abundance modeling using the actual NS+ collisional rate coefficients is needed. For this purpose, we have computed the first full 4D potential energy surface of the NS+-H2 van der Waals complex using the explicitly correlated coupled cluster approach with single, double, and non-iterative triple excitation in conjunction with the augmented-correlation consistent-polarized valence triple zeta basis set. The potential energy surface exhibits a global minimum of 848.24 cm-1 for a planar configuration of the complex. The long-range interaction energy, described using multipolar moments, is sensitive to the orientation of H2 up to radial distances of ∼50 a0. From this new interaction potential, we derived excitation cross sections, induced by collision with ortho- and para-H2, for the 15 low-lying rotational levels of NS+ using the quantum mechanical close-coupling approach. By thermally averaging these data, we determined downward rate coefficients for temperatures up to 50 K. By comparing them with the previous NS+-H2 data, we demonstrated that reduced dimensional approaches are not suited for this system. In addition, we found that the CS collisional data underestimate our results by up to an order of magnitude. The differences clearly indicate that the abundance of NS+, in cold dense clouds retrieved from observational spectra, must be reassessed using these new collisional rate coefficients.

Rotranslational dynamics of confined water. II. Spectroscopic evidence of confinement effects on the far-infrared spectra of water isotopologues in argon and krypton matrices.

Water molecules trapped in rare gas matrices exhibit conspicuous shifts in their far-infrared (FIR), rotranslational spectral features compared with the corresponding transitions observed in the gas phase. These confinement-induced perturbations have been related not only to the quantization of translational motion but also to the coupling between the orientational and positional degrees of freedom: the rotation-translation coupling (RTC). As the propensity displayed by the nuclear spin isomers (NSI) of water to undergo interconversion in confinement is intimately related to how its nuclear spin degrees of freedom are coupled with those for intra- and intermolecular motions, confinement-induced RTC should also strongly impact the NSI interconversion mechanisms and rates. Insight into the rotranslational dynamics for H2 16O, H2 17O, and H2 18O, confined in argon and krypton matrices, is provided here based on the evolution of rotranslational spectra induced by NSI interconversion while a definitive assignment is provided from the transition energies and intensities calculated using the confined rotor model [Paper I, Wespiser et al., J. Chem. Phys. 156, 074304 (2021)]. In order to build a complete rotranslational energy diagram of confined water, which is fundamental to understand the NSI interconversion rates, the energy difference between the ground ortho and para rotranslational states is derived from the temperature dependence of the intensity ratio of mid-infrared lines emerging from these states. These investigations should provide deeper insight of the factors that control NSI interconversion of water isotopologues under extreme confinement.

Interaction of H2O with CO: potential energy surface, bound states and scattering calculations.

Collisions between H2O and CO play a crucial role in the gaseous component of comets and protoplanetary disks. We present here a five-dimensional potential energy surface (PES) for the H2O-CO collisional complex. Ab initio calculations were carried out using the explicitly-correlated closed-shell single- and double-excitation coupled cluster approach with the non-iterative perturbative treatment of triple excitations [CCSD(T)-F12a] method with the augmented correlation-consistent aug-cc-pVTZ basis sets. The most stable configuration of the complex, where the carbon atom of CO is pointing towards the OH bond of water, has a binding energy De = 646.1 cm-1. The end-over-end rotational constant of the H2O-CO complex was extracted from bound state calculations and it was found to be B0 = 0.0916 cm-1, in excellent agreement with experimental measurements. Finally, cross sections for the rotational excitation of CO by H2O are computed for s-wave (J = 0) scattering at the full close-coupling level of theory. These results will serve as a benchmark for future studies.

Ab initio potential and rotational spectra of the CO-N2 complex.

Ab initio calculations of the intermolecular potential energy surface (PES) of CO-N2 have been carried out using the closed-shell single- and double-excitation coupled cluster approach with a non-iterative perturbative treatment of triple excitations method and the augmented correlation-consistent quadruple-zeta (aug-cc-pVQZ) basis set supplemented with midbond functions. The global minimum (De = 117.35 cm-1) of the four-dimensional PES corresponds to an approximately T-shaped structure with the N2 subunit forming the leg and CO the top. The bound rovibrational levels of the CO-N2 complex were calculated for total angular momenta J = 0-8 on this intermolecular potential surface. The calculated dissociation energies D0 are 75.60 and 76.79 cm-1 for the ortho-N2 (A-symmetry) and para-N2 (B-symmetry) nuclear spin modifications of CO-N2, respectively. Guided by these bound state calculations, a new millimeter-wave survey for the CO-N2 complex in the frequency range of 110-145 GHz was performed using the intracavity OROTRON jet spectrometer. Transitions not previously observed were detected and assigned to the subbands connecting the K = 0 and 1, (jCO, jN2 ) = (1, 0) states with a new K = 1, (jCO, jN2 ) = (2, 0) state. Finally, the measured rotational energy levels of the CO-N2 complex were compared to the theoretical bound state results, thus providing a critical test of the quality of the PES presented. The computed rovibrational wave functions were analyzed to characterize the nature of the different bound states observed for the two nuclear spin species of CO-N2.

Rotational study of the NH3-CO complex: millimeter-wave measurements and ab initio calculations.

The rotational spectrum of the van der Waals complex NH3-CO has been measured with the intracavity OROTRON jet spectrometer in the frequency range of 112-139 GHz. Newly observed and assigned transitions belong to the K = 0-0, K = 1-1, K = 1-0, and K = 2-1 subbands correlating with the rotationless (jk)NH3 = 00 ground state of free ortho-NH3 and the K = 0-1 and K = 2-1 subbands correlating with the (jk)NH3 = 11 ground state of free para-NH3. The (approximate) quantum number K is the projection of the total angular momentum J on the intermolecular axis. Some of these transitions are continuations to higher J values of transition series observed previously [C. Xia et al., Mol. Phys. 99, 643 (2001)], the other transitions constitute newly detected subbands. The new data were analyzed together with the known millimeter-wave and microwave transitions in order to determine the molecular parameters of the ortho-NH3-CO and para-NH3-CO complexes. Accompanying ab initio calculations of the intermolecular potential energy surface (PES) of NH3-CO has been carried out at the explicitly correlated coupled cluster level of theory with single, double, and perturbative triple excitations and an augmented correlation-consistent triple zeta basis set. The global minimum of the five-dimensional PES corresponds to an approximately T-shaped structure with the N atom closest to the CO subunit and binding energy De = 359.21 cm(-1). The bound rovibrational levels of the NH3-CO complex were calculated for total angular momentum J = 0-6 on this intermolecular potential surface and compared with the experimental results. The calculated dissociation energies D0 are 210.43 and 218.66 cm(-1) for ortho-NH3-CO and para-NH3-CO, respectively.

Rotational study of the CH4-CO complex: Millimeter-wave measurements and ab initio calculations.

The rotational spectrum of the van der Waals complex CH4-CO has been measured with the intracavity OROTRON jet spectrometer in the frequency range of 110-145 GHz. Newly observed and assigned transitions belong to the K = 2-1 subband correlating with the rotationless jCH4 = 0 ground state and the K = 2-1 and K = 0-1 subbands correlating with the jCH4 = 2 excited state of free methane. The (approximate) quantum number K is the projection of the total angular momentum J on the intermolecular axis. The new data were analyzed together with the known millimeter-wave and microwave transitions in order to determine the molecular parameters of the CH4-CO complex. Accompanying ab initio calculations of the intermolecular potential energy surface (PES) of CH4-CO have been carried out at the explicitly correlated coupled cluster level of theory with single, double, and perturbative triple excitations [CCSD(T)-F12a] and an augmented correlation-consistent triple zeta (aVTZ) basis set. The global minimum of the five-dimensional PES corresponds to an approximately T-shaped structure with the CH4 face closest to the CO subunit and binding energy De = 177.82 cm(-1). The bound rovibrational levels of the CH4-CO complex were calculated for total angular momentum J = 0-6 on this intermolecular potential surface and compared with the experimental results. The calculated dissociation energies D0 are 91.32, 94.46, and 104.21 cm(-1) for A (jCH4 = 0), F (jCH4 = 1), and E (jCH4 = 2) nuclear spin modifications of CH4-CO, respectively.

Hints of a rotating spiral structure in the innermost regions around IRC +10216.

The Atacama Large Millimeter/submillimeter Array (ALMA) is allowing us to study the innermost regions of the circumstellar envelopes of evolved stars with un-precedented precision and sensitivity. Key processes in the ejection of matter and dust from these objects occur in their inner zones. In this work, we present sub-arcsecond interferometric maps of transitions of metal-bearing molecules towards the prototypical C-rich evolved star IRC +10216. While Al-bearing molecules seem to be present as a roughly spherical shell, the molecular emission from the salts NaCl and KCl presents an elongation in the inner regions, with a central minimum. In order to accurately analyze the emission from the NaCl rotational lines, we present new calculations of the collisional rates for this molecule based on new spectroscopic constants. The most plausible interpretation for the spatial distribution of the salts is a spiral with a NaCl mass of 0.08M☉. Alternatively, a torus of gas and dust would result in similar structures as those observed. From the torus scenario we derive a mass of ~ 1.1 × 10-4M☉. In both cases, the spiral and the torus, the NaCl structure presents an inner minimum of 27 AU. In the case of the torus, the outer radius is 73 AU. The kinematics of both the spiral and the torus suggests that they are slowly expanding and rotating. Alternative explanations for the presence of the elongation are explored. The presence of these features only in KCl and NaCl might be a result of their comparatively high dipole moment with respect to the Al-bearing species.