An unbiased NOEMA 2.6 to 4 mm survey of the GG Tau ring: First detection of CCS in a protoplanetary disk
Context. Molecular line surveys are among the main tools to probe the structure and physical conditions in protoplanetary disks (PPDs), the birthplace of planets. The large radial and vertical temperature as well as density gradients in these PPDs lead to a complex chemical composition, making chemistry an important step to understand the variety of planetary systems. Aims. We aimed to study the chemical content of the protoplanetary disk surrounding GG Tau A, a well-known triple T Tauri system. Methods. We used NOEMA with the new correlator PolyFix to observe rotational lines at ∼2.6 to 4 mm from a few dozen molecules. We analysed the data with a radiative transfer code to derive molecular densities and the abundance relative to 13CO, which we compare to those of the TMC1 cloud and LkCa 15 disk. Results. We report the first detection of CCS in PPDs. We also marginally detect OCS and find 16 other molecules in the GG Tauri outer disk. Ten of them had been found previously, while seven others (13CN, N2H+, HNC, DNC, HC3N, CCS, and C34S) are new detections in this disk. Conclusions. The analysis confirms that sulphur chemistry is not yet properly understood. The D/H ratio, derived from DCO+/HCO+, DCN/HCN, and DNC/HNC ratios, points towards a low temperature chemistry. The detection of the rare species CCS confirms that GG Tau is a good laboratory to study the protoplanetary disk chemistry, thanks to its large disk size and mass.
Theoretical study of the HCS+–H2 van der Waals complex: potential energy surface, rovibrational bound states, and rotationally inelastic collisional cross sections
An accurate ground-state intermolecular potential energy surface (PES) was calculated for the HCS (Formula presented.) –H (Formula presented.) complex. The surface was constructed from 3023 ab initio energies, computed with explicitly correlated coupled-cluster theory, CCSD(T)-F12b, with extrapolation to the complete basis set limit (VTZ-F12/VQZ-F12). The new 4D PES was used to compute rovibrational energies, rotational constants, and close-coupling quantum scattering calculations at low collision energies. A symmetry-adapted Lanczos algorithm was used to variationally compute the low-lying rovibrational bound states. Rotational constants for states of two isomers were determined from the energy levels and reported in hopes of motivating experiments. The rotationally inelastic state-to-state cross sections of HCS (Formula presented.) by collision with para- and ortho-H (Formula presented.) were computed and compared, and trends and propensities are discussed.
Rotational relaxation of HCO+and DCO+by collision with H2
The HCO+ and DCO+ molecules are commonly used as tracers in the interstellar medium. Therefore, accurate rotational rate coefficients of these systems with He and H2 are crucial in non-local thermal equilibrium models. We determine in this work the rotational de-excitation rate coefficients of HCO+ in collision with both para- and ortho-H2, and also analyse the isotopic effects by studying the case of DCO+. A new four-dimensional potential energy surface from ab initio calculations was developed for the HCO+-H2 system, and adapted to the DCO+-H2 case. These surfaces are then employed in close-coupling calculations to determine the rotational de-excitation cross-sections and rate coefficients for the lower rotational states of HCO+ and DCO+. The new rate coefficients for HCO+ + para-H2 were compared with the available data, and a set of rate coefficients for HCO+ + ortho-H2 is also reported. The difference between the collision rates with ortho- and para-H2 is found to be small. These calculations confirm that the use of the rate coefficients for HCO+ + para-H2 for estimating those for HCO+ + ortho-H2 as well as for DCO+ + para-H2 is a good approximation.
An explicitly correlated six-dimensional potential energy surface for the SiCSi + H2 complex
The first six-dimensional potential energy surface (PES) for the SiCSi + H2 complex is presented in this work. This surface is developed from a large number of ab initio energies computed at the explicitly correlated coupled-cluster level of theory together with the augmented correlation-consistent polarized valence triple zeta basis set (CCSD(T)-F12/aug-cc-pVTZ). These energies are fitted to an analytical function through a procedure that combines spline, least-squares, and kernel-based methods. Two minimums of similar depths were found at the equilibrium geometry of the SiCSi molecule. The dependence of the PES on the bending angle is analyzed. Furthermore, a reduced four-dimensional PES averaged over the H2 orientation is presented. Finally, the six-dimensional PES is used for computing the second virial coefficient of the SiCSi + H2 pair using classical and semi-classical methods.
Quantum study of the bending relaxation of H2O by collision with H
Vibrationally excited levels of the H2O molecule are currently detected in various environments of the interstellar medium (ISM), and collisional data for H2O, including vibration with the main colliders of the ISM, are needed. The present study focuses on the bending relaxation of H2O by collision with H when taking bending-rotation coupling explicitly into account with the rigid-bender close-coupling (RB-CC) method. With this aim, a new four-dimensional potential energy surface including the H2O bending mode is developed from a large grid of ab initio energies computed using a high level of theory. For purely rotational transitions, our RB-CC rates show very good agreement with rigid-rotor calculations performed using our new potential energy surface (PES) and with those available in the literature. Calculations for pure rotational transitions inside the excited bending level ν2 = 1 of H2O are performed and compared with their equivalents inside ν2 = 0. Vibrational quenching of H2O is also calculated and found to be much more efficient through collision with H rather than with He. © 2022 The Author(s)
A Rigid Bender Study of the Bending Relaxation of H2O and D2O by Collisions with Ar
The bending relaxation of H2O and D2O by collisions with Ar is studied at the Close Coupling level. Two new 4D PES are developed for these two systems. They are tested by performing rigid rotor calculations as well as by computing the D2O−Ar bound states. The results are compared with available theoretical and experimental data. Propensity rules for the dynamics are discussed and compared to those of H2O colliding with Ne or He. The bending relaxation cross sections and rates are then calculated for these two systems. The results are analysed and compared with available experimental data. © 2024 The Authors. ChemPhysChem published by Wiley-VCH GmbH.
Rotational Relaxation of AlNC and AlCN by para-H2(j = 0) at Low Temperatures
An accurate determination of the abundances of metal-containing molecules in the interstellar medium or circumstellar gas requires knowledge of molecular data, including the collisional rate coefficients. This work is focused on the study of the collision of the aluminum isocyanide (AlNC) molecule, as well as its isomer AlCN, with para-H2 (j = 0). For the AlNC + H2 and AlCN + H2 complexes, averaged potential energy surfaces are developed from ab initio energies computed at the coupled cluster with the single, double, and perturbative triple excitation level of theory. Such surfaces are used in close-coupling calculations. The rate coefficients at low temperature are compared with those for the collisions with He. The use of the mass scaling procedure is a good approximation in the case of AlCN. However, for the collision with AlNC, a different propensity rule is found between the rates with He and para-H2 (j = 0). Finally, rotational rate coefficients for the lowest 26 rotational levels of both molecules, AlCN and AlNC, by collision with para-H2 (j = 0) are reported.
A benchmark for the size of the QM system required for accurate hybrid QM/MM calculations on the metal site of the protein copper, zinc superoxide dismutase
Study of the CN(X2ς+) + N(4S) Reaction at High Temperatures: Potential Energy Surface and Thermal Rate Coefficients
Reactions involving C and N play an essential role in the chemistry around the surface of a hypersonic spacecraft during its atmospheric re-entry. The collision of CN with other molecules and atoms has particular interest in aerothermodynamic modeling. This work focuses on the study of the CN + N → N2 + C reaction in the triplet manifold 3A″ of CN2. A high-level full-dimensional potential energy surface for this system is developed from ab initio calculations at the MRCI-F12 + Q level of theory. This surface is employed in quasiclassical trajectory calculations, and thermal rate coefficients from 100 to 20,000 K are computed. The rates for the formation of N2 are compared with the available experimental data, and good agreement is found. At low and intermediate temperatures, the N2 formation is more efficient than the N-exchange process, while at high temperatures, the rates for both processes are comparable. Finally, analytically modified Arrhenius expressions for the reaction rates of N2 formation and N-exchange are reported.
Deexcitation rate coefficients of C3by collision with H2at low temperatures
Context. An accurate analysis of the physical-chemical conditions in the regions of the interstellar medium in which C3 is observed requires knowing the collisional rate coefficients of this molecule with He, H2, electrons, and H. Aims. The main goals of this study are to present the first potential energy surface for the C3 +H2 complex, to study the dynamics of the system, and to report a set of rate coefficients at low temperature for the lower rotational states of C3 with para- and ortho-H2. Methods. A large grid of ab initio energies was computed at the explicitly correlated coupled-cluster with single-, double-, and perturbative triple-excitation level of theory, together with the augmented correlation-consistent quadruple zeta basis set (CCSD(T)-F12a/aug-cc-pVQZ). This grid of energies was fit to an analytical function. The potential energy surface was employed in close-coupling calculations at low collisional energies. Results. We present a high-level four-dimensional potential energy surface (PES) for studying the collision of C3 with H2. The global minimum of the surface is found in the linear HH-CCC configuration. Rotational deexcitation state-to-state cross sections of C3 by collision with para- and ortho-H2 are computed. Furthermore, a reduced two-dimensional surface is developed by averaging the surface over the orientation of H2. The cross sections for the collision with para-H2 using this approximation and those from the four-dimensional PES agree excellently. Finally, a set of rotational rate coefficients for the collision of C3 with para- and ortho-H2 at low temperatures are reported.