Rotational relaxation of CS by collision with ortho- and para-H2 molecules

Abundances of sulphur molecules in the Horsehead nebula First NS+ detection in a photodissociation region

CONTEXT

Sulphur is one of the most abundant elements in the Universe (S/H∼1.3×10 -5 ) and plays a crucial role in biological systems on Earth. The understanding of its chemistry is therefore of major importance.

AIMS

Our goal is to complete the inventory of S-bearing molecules and their abundances in the prototypical photodissociation region (PDR) the Horsehead nebula to gain insight into sulphur chemistry in UV irradiated regions. Based on the WHISPER (Wide-band High-resolution Iram-30m Surveys at two positions with Emir Receivers) millimeter (mm) line survey, our goal is to provide an improved and more accurate description of sulphur species and their abundances towards the core and PDR positions in the Horsehead.

METHODS

The Monte Carlo Markov Chain (MCMC) methodology and the molecular excitation and radiative transfer code RADEX were used to explore the parameter space and determine physical conditions and beam-averaged molecular abundances.

RESULTS

A total of 13 S-bearing species (CS, SO, SO2, OCS, H2CS - both ortho and para - HDCS, C2S, HCS+, SO+, H2S, S2H, NS and NS+) have been detected in the two targeted positions. This is the first detection of SO+ in the Horsehead and the first detection of NS+ in any PDR. We find a differentiated chemical behaviour between C-S and O-S bearing species within the nebula. The C-S bearing species C2S and o-H2CS present fractional abundances a factor of > two higher in the core than in the PDR. In contrast, the O-S bearing molecules SO, SO2, and OCS present similar abundances towards both positions. A few molecules, SO+, NS, and NS+, are more abundant towards the PDR than towards the core, and could be considered as PDR tracers.

CONCLUSIONS

This is the first complete study of S-bearing species towards a PDR. Our study shows that CS, SO, and H2S are the most abundant S-bearing molecules in the PDR with abundances of ∼ a few 10-9. We recall that SH, SH+, S, and S+ are not observable at the wavelengths covered by the WHISPER survey. At the spatial scale of our observations, the total abundance of S atoms locked in the detected species is < 10-8, only ∼0.1% of the cosmic sulphur abundance.

Inelastic vibrational dynamics of CS in collision with H2 using a full-dimensional potential energy surface

We report a six-dimensional (6D) potential energy surface (PES) for the CS-H2 system computed using high-level electronic structure theory and fitted using a hybrid invariant polynomial method. Full-dimensional quantum close-coupling scattering calculations have been carried out using this potential for rotational and, for the first time, vibrational quenching transitions of CS induced by H2. State-to-state cross sections and rate coefficients for rotational transitions in CS from rotational levels j1 = 0-5 in the ground vibrational state are compared with previous theoretical results obtained using a rigid-rotor approximation. For vibrational quenching, state-to-state and total cross sections and rate coefficients were calculated for the vibrational transitions in CS(v1 = 1,j1) + H2(v2 = 0,j2) → CS(v1' = 0,j1') + H2(v2' = 0,j2') collisions, for j1 = 0-5. Cross sections for collision energies in the range 1 to 3000 cm-1 and rate coefficients in the temperature range of 5 to 600 K are obtained for both para-H2 (j2 = 0) and ortho-H2 (j2 = 1) collision partners. Application of the computed results in astrophysics is also discussed.

Comparative experimental and theoretical study of the rotational excitation of CO by collision with ortho- and para-D2 molecules

A joint crossed beam and quantum mechanical investigation of the rotationally inelastic collisions of CO with ortho- and para-D₂ molecules is reported. A new 4D potential energy surface (PES) averaged over the ground vibrational states of D₂ and CO is used to calculate the rovibrational bound states of the ortho-D₂-CO complexes. Close coupling calculations are then performed in the rigid rotor approximation for ortho- and para-D₂ colliding with CO for the experimentally investigated transition of CO (j = 0 → 1) and for collision energies ranging from 0.1 to 25 cm^-1. The agreement between theory and experiment is found to be very good for both the bound state energies of the ortho-D₂-CO complexes and for the inelastic scattering cross-sections showing that the free rotation of two rigid rotors is a very good model of the D₂-CO system in this low collision energy domain.

On the importance of full-dimensionality in low-energy molecular scattering calculations

Scattering of H₂ on CO is of great importance in astrophysics and also is a benchmark system for comparing theory to experiment. We present here a new 6-dimensional potential energy surface for the ground electronic state of H₂-CO with an estimated uncertainty of about 0.6 cm⁻¹ in the global minimum region, several times smaller than achieved earlier. This potential has been used in nearly exact 6-dimensional quantum scattering calculations to compute state-to-state cross-sections measured in low-energy crossed-beam experiments. Excellent agreement between theory and experiment has been achieved in all cases. We also show that the fully 6-dimensional approach is not needed with the current accuracy of experimental data since an equally good agreement with experiment was obtained using only a 4-dimensional treatment, which validates the rigid-rotor approach widely used in scattering calculations. This finding, which disagrees with some literature statements, is important since for larger systems full-dimensional scattering calculations are currently not possible.