Issue 
A&A
Volume 638, June 2020



Article Number  A31  
Number of page(s)  5  
Section  Atomic, molecular, and nuclear data  
DOI  https://doi.org/10.1051/00046361/202037821  
Published online  05 June 2020 
Rotational relaxation of H_{2}S by collision with He
^{1}
Núcleo de Astroquímica y Astrofísica, Instituto de Ciencias Químicas Aplicadas, Facultad de Ingeniería, Universidad Autónoma de Chile, Av. Pedro de Valdivia 425, Providencia, Santiago, Chile
email: otoniel.denis@uautonoma.cl
^{2}
Institut des Sciences Moléculaires, Université de Bordeaux, CNRS UMR 5255, 33405 Talence Cedex, France
email: thierry.stoecklin@ubordeaux.fr
Received:
25
February
2020
Accepted:
17
April
2020
Context. The H_{2}S molecule has been detected in several regions of the interstellar medium (ISM). The use of nonLTE models requires knowledge of accurate collisional rate coefficients of the molecules detected with the most common collider in the ISM.
Aims. The main goal of this work is to study the collision of H_{2}S with He.
Methods. A grid of ab initio energies was computed at the coupled cluster level of theory including single, double, and perturbative triple excitations (CCSD(T)) and using the augmented correlation consistent polarized quadruple zeta (augccpVQZ) basis set supplemented by a set of midbond functions. These energies were fitted to an analytical function, which was employed to study the dynamics of the system. Close coupling calculations were performed to study the collision of H_{2}S with He.
Results. The rate coefficients determined from the close coupling calculation were compared with those of the collision with H_{2}O+He, and large differences were found. Finally, the rate coefficients for the lower rotational deexcitation of H_{2}S by collision with He are reported.
Key words: astrochemistry / molecular data / molecular processes / scattering / ISM: molecules
© ESO 2020
1. Introduction
Sulfur molecules play an important role in biological systems on Earth, which increases the interest in following the chemical history of these molecules in the Universe (Fuente et al. 2017). Approximately 20 molecules containing sulfur have been detected in the interstellar medium (ISM), as can be seen in the Cologne Database for Molecular Spectroscopy (CDMS; Müller et al. 2005). Owing to their properties and abundance, these systems are useful tracers of interstellar material (Duley et al. 1980; Nagy et al. 2013; Fuente et al. 2016; Goicoechea et al. 2017; Gorai 2018; Cernicharo et al. 2018; Drozdovskaya et al. 2018).
Among all sulfurbearing molecules, we focus on hydrogen sulfide (H_{2}S). This system has been observed in a variety of astrophysical environments (Thaddeus et al. 1972; Ukita & Morris 1983; Minh et al. 1989, 1990; Martin et al. 2005; Li et al. 2015; Holdship et al. 2016; Minh 2016; Phuong et al. 2018). H_{2}S has been used as a probe of sulfur chemistry in the ISM and also as a probe of the physical state of the gas (Crockett et al. 2014).
The local thermal equilibrium (LTE) can be employed to determine the physicochemical conditions of the region where the molecules are observed, for example, density, excitation temperature, and molecular abundance. In such a model, the rotational population of the molecule follows a Boltzmann distribution. However, in typical molecular clouds the density is low and collisions are rare. Also, the LTE does not apply, and nonLTE models should be considered. Such models require the input of both the Einstein coefficients, which are usually known, and accurate statetostate collisional rate coefficients with the most common colliders of the ISM (e.g., He, H_{2}, H, and e^{−}).
Recently, Crockett et al. (2014) studied the H_{2}S molecule in the Orion KL nebula and performed nonLTE analysis. In the absence of collisional rates for this system, theses authors used two approximations. First, they employed the experimental depopulation cross section of Ball et al. (1999) for the 1_{1, 0} → 1_{0, 1} transition of H_{2}S+He to derive a base rate coefficient of 1.35 × 10^{−11} cm^{3} s^{−1} at 10 K. The collision rates were then scaled in proportion to radiative line strengths so that the sum of all downward rates from each upper state is equal to the base rate (Crockett et al. 2014). Second, Crockett et al. (2014) used the rates for the collisions of H_{2}O with H_{2} to scale those for the collisions with H_{2}S. To conclude, it is important to point out that the rate coefficients available for H_{2}S in the LAMDA database (Schöier et al. 2005) are also those scaled from the collision with H_{2}O, while the spectroscopic characteristics of the H_{2}O and H_{2}S molecules are very different, and differences in the interaction potential energy surfaces (PESs) of He with H_{2}S and H_{2}O are most probable. A full set of collisional rates for H_{2}S is then required for evaluating the validity of these two approximations.
A very recent work dedicated to the collisions of H_{2} with H_{2}S (Dagdigian 2020) was published online at the same time as the present work was submitted to publication. A brief comparison of the rates for collision of H_{2}S with He and H_{2} are therefore used to discuss the validity of the mass scaling procedure for estimating the rates for collisions with H_{2} from those with He.
The main goal of this paper is to provide a new set of reliable rotational rate coefficients for the collision of H_{2}S with He calculated at the close coupling level on a new ab initio PES. This article is organized as follows: In the next section, the methods used are presented. In Sect. 3, the results are discussed, while the conclusions can be found in Sect. 4.
2. Methods
2.1. Ab initio calculations
A large grid of interaction ab initio energies were computed to study the H_{2}S+He complex. The calculations were performed at the coupled cluster level of theory including single, double, and perturbative triple excitations (CCSD(T); Raghavachari et al. 1989) and using the augmented correlation consistent polarized quadruple zeta (augccpVQZ; Dunning 1989) basis set supplemented by a set of midbond functions (Cybulski & Toczyłowski 1999). The basis set superposition error was corrected using the counterpoise procedure (Boys & Bernardi 1970). The calculations were performed with the MOLPRO package (Werner et al. 2012). The center of the coordinates of the system was located at the center of mass of H_{2}S molecule, H_{2}S is in the plane xz, and the scattering coordinates R, θ, and φ are defined in Fig. 1.
Fig. 1.
Internal coordinates used to describe the H_{2}S+He system. 
The rigid rotor approximation is used to describe the H_{2}S molecule and the HS bond length r_{S − H} = 1.35 Å and bond angle ∢HSH = 92.13° are fixed to their vibrationally averaged value in the rovibrational ground state Cook et al. (1975). The radial grid includes 25 values of R ranging from 1.5 Å to 9 Å, while θ and ϕ vary from 0° to 180° in steps of 10° and 15°, respectively. A total of 6175 ab initio energies were computed.
2.2. Analytical representation of the PES
The ab initio energies were fitted to the following analytical form:
where the angular part is represented by a product of normalized associated Legendre and cosine functions, and m is even. For each R value, the ab initio energies were fitted to this analytical form using a leastsquares procedure. Then, the f_{l, m} radial coefficients were obtained using a reproducing kernel Hilbert space (RKHS) representation,
with
where N is the number of radial points of the grid and q^{5, 2}(R, R_{k}) are the onedimensional kernels defined by Ho and Rabitz (Ho & Rabitz 1996). The coefficients are found by solving the linear system f_{lm}(R_{i}) = α^{l, m}q(R_{i}, R_{j}), where i and j label different radial grid points, while the use of the q^{5, 2}(R, R_{k}) kernel ensures a R^{−6} longrange behavior (Soldán & Hutson 2000).
2.3. Dynamics
The dynamics of the system was studied by separately solving the close coupling equation for the two nuclear spin modifications of H_{2}S in the spacefixed frame. The calculations were performed with the most recent version of the Newmat code Stoecklin et al. (2019) and the logderivative propagator.
The rotational constants of H_{2}S were fixed to their experimental values A_{0} = 10.3601 cm^{−1}, B_{0}= 9.0156 cm^{−1}, and C_{0} = 4.7318 cm^{−1} (Cook et al. 1975), while the rotational basis set was made of 10 values of the rotational quantum number j. This included 50 ortho rotational states up to 908.085 cm^{−1} and 50 para states up to 908.065 cm^{−1}. In what follows, we use the notation j_{KAKC} to designate the asymmetric top rotational energies of H_{2}S, where K_{A} and K_{C} are, respectively, the absolute values of the Z molecular axis projections of the rotational angular momentum j in the limit of a prolate or an oblate top. This molecule is a nearoblate asymmetric top, and the oblatelimit projection quantum number K_{C} is then an approximately good quantum number. This is in contrast with water, which is an intermediate case closer to a prolate limit for which neither K_{A} nor K_{C} are nearly conserved. Scattering calculations were carried out for energies up to 4000 cm^{−1} for each H_{2}S–He pair of nuclear spin modifications. This range of energies could involve exited bending and stretching levels of H_{2}S, while the H_{2}S vibration is not taken into account in the present work. We know however from previous work (Stoecklin et al. 2019) that the statetostate rotational excitation uncertainty resulting from the use of this approximation above the excited vibrational level thresholds of H_{2}S is very small as rovibrational cross sections are two orders of magnitude smaller than purely rotational cross section. Partial waves up to total angular momentum J = 90 were included in the calculations to achieve a 10^{−3} relative criterion for the convergence of the stateselected quenching cross section as a function of the maximal value of the total angular momentum quantum number J. The maximum propagation distance was 50 bohr, and convergence was checked as a function of the propagator step size. The computed statetostate integral cross sections were eventually employed to calculate the corresponding stateselected rate coefficients between 1 and 1000 K by Boltzmann averaging over the collision energy E_{c}, that is,
3. Results
3.1. PES
The grid of ab initio energies was fitted to the analytical form of Eq. (1). The two panels of Fig. 2 shows contour plots of the resulting PES for φ = 0° and φ = 90°, respectively. The global minimum of the PES (−36.49 cm^{−1}) appearing on the left panel of this figure is associated with the bent configuration (R = 3.46 Å, θ = 72°, φ = 0°). A secondary minimum, at R = 3.84 Å and θ = 180° (φ is undefined for θ = 0° and 180°), can be seen on this figure. The global minimum energy of H_{2}S+He PES is then very close to that of the H_{2}O+He system (−34.9 cm^{−1}) (Patkowski et al. 2002), which is also associated with a similar geometry of the complex (R = 3.13 Å, θ = 75°, φ = 0°). As a result, the magnitudes of the cross sections should approximately follow the scaling law based on the ratio of the reduced masses of the two systems used in the Lambda database (Schöier et al. 2005). Conversely, the angular anisotropies of the two PES are very different because the H_{2}O+He PES does not exhibit any secondary minimum, as shown in Fig. 3 of Patkowski et al. (2002). As a result and because of the consequences of the differences between the values of the rotational constants of H_{2}O and H_{2}S discussed in Sect. 2.3, the propensity rules followed by the rotational transitions of the two collisional systems are expected to differ.
Fig. 2.
Contour plots of the H_{2}S+He system at φ = 0° (panel A) and φ = 90° (panel B). Red contour lines correspond to positive energies ranging from 0 to 100 cm^{−1} in steps of 10 cm^{−1}, and in steps of 100 cm^{−1} in the [100, 2000] cm^{−1} interval. Blue contour lines show the negative energies varying in steps of 2 cm^{−1}. 
3.2. Dynamics
Figures 3 and 4 show deexcitation cross sections from examples of para (K_{a} + K_{c} even) and ortho (K_{a} + K_{c} odd) states of the j = 2 and j = 4 multiplets. In each figure, the lower and upper panels are dedicated to examples of initial para and ortho states, respectively. These are on one side the (j = 2, K_{a} = 2, K_{c} = 0) and (j = 4, K_{a} = 4, K_{c} = 0) para states and on the other the (j = 2, K_{a} = 2, K_{c} = 1) and (j = 4, K_{a} = 4, K_{c} = 1) ortho state of H_{2}S. Before analyzing these figures we note that all these transitions tend to follow the general propensity rule demonstrated by Alexander (1982). This early paper states that transitions between states with the same sign of the parity index ϵ = (−)^{j + KA + KC} are favored, and the collisional selectivity becomes increasingly rigorous as both j and the collisional energy increase. In order to simplify the notation, in the following we only give the three integers mnp to designate the j = m, K_{a} = n, K_{c} = p level.
Fig. 3.
Rotational deexcitation cross sections of H_{2}S by collision with He. The upper and lower panels show transitions issued from the 2 2 1 and 2 2 0 initial levels of H_{2}S, respectively. 
Fig. 4.
Rotational deexcitation cross sections of H_{2}S by collision with He. The upper and lower panels show transitions issued from the 4 4 1 and 4 4 0 initial levels of H_{2}S, respectively. 
If we first look at the transitions issued from the 2 2 1 ortho level of H_{2}S represented in the upper panel of Fig. 3, we see that the 2 1 2 final level gives the largest cross section in the [1–100] cm^{−1} interval. This is expected as this transition is elastic in j (Δj = 0), while these two levels have the same parity index and are the closest in energy. They however differ in K_{c} (ΔK_{c} = 1), whereas H_{2}S is a nearoblate symmetric top and K_{c} should then nearly be a good quantum number. At higher collisional energy, the (ΔK_{c} = 0) propensity rule is conversely respected as the 1 0 1 final level gives the largest cross section. Similar behavior can be observed for the transitions issued from the initial 2 2 0 para state represented in the lower panel of the same Fig. 3). The 2 1 1 level followed by the 2 0 2 level are closest in energy, have the same parity index as 2 2 0, and are j elastic. Again these two transitions give the largest cross section in the [1–100] cm^{−1} interval. At higher energy the (ΔK_{c} = 0) propensity rule makes the 0 0 0 transition progressively increase, whereas it is inelastic in j (Δj = 2) and is the most distant in energy.
If we now look at the transitions issued from the 4 4 1 ortho level represented in the upper panel of Fig. 4, we see that again the 4 3 2 and 4 2 3 levels, which are the closest in energy with the same parity index, give the highest cross section up to 300 cm^{−1}. However, the 3 2 1 level, whose parity index is opposite, is inelastic in j (Δj = 1) and is quite distant in energy. But this level follows the (ΔK_{c} = 0) propensity rule giving cross sections that are very close in magnitude and become largest at higher energy. Similar behavior can be observed in the lower panel of Fig. 4 for the transitions issued from the 4 4 0 para level. The 2 2 0 level, which is quite distant in energy and is inelastic in j (Δj = 2), gives the largest cross section above 100 cm^{−1}. We conclude that the (ΔK_{c} = 0) propensity rule becomes increasingly rigorous as j and the collision energy both increase, while for collision energy below the well depth, it is rather the conservation of the parity index and the lowest energy difference that lead to the highest transition cross sections.
We now turn to the comparison between the H_{2}S+He and H_{2}O+He collisional systems as we have seen that the well depths of the two associated PES are very close, while the ratio of the reduced masses of the two systems nearly equals one. These are two reasons that lead several authors to assume that the collisional rate coefficients would also be equal. The rotational deexcitation rate coefficients of H_{2}S by collision with He from the initial states 2 2 1 and 2 2 0 of H_{2}S are shown in the upper and lower panels of Fig. 5, respectively, where they are compared with those reported for the H_{2}O+He system by Yang et al. (2013), which can be found in the Basecol database (Dubernet et al. 2013). All three temperature dependence, magnitude, and propensity rules differ strongly for the two systems. These strong differences result from both the differences in the angular anisotropies of the two PESs and the rotational constants; H_{2}S has a nearly oblate symmetric top, while H_{2}O is a good example of an asymmetric top intermediate case. The data in this work computed will be sent to the Basecol database soon.
Fig. 5.
Rotational deexcitation rate coefficients of H_{2}S by collision with He. The upper and lower panels show transitions issued from the 2 2 1 and 2 2 0 initial levels of water, respectively. The rates are compared with the corresponding rates for the H_{2}S + paraH_{2} collisions (designated by pH_{2}) and with those for H_{2}O+He collisions (designated by H_{2}O). 
We then compare our results with the only experimental data available of Ball et al. (1999) for the 1_{10} → 1_{01} transition of H_{2}S colliding with He. As a matter of fact, PES calculations have reached a very high level of accuracy, which is illustrated by the very good agreement between theory and experiment obtained in several recent studies dedicated to the H_{2}CO (Faure et al. 2016), D_{2}CO (Stoecklin et al. 2017), ArCO (Mertens et al. 2017), NO+O_{2} (Gao et al. 2018), and ND_{3}D_{2} (Gao et al. 2019) systems. Ball et al. (1999) did not obtain statetostate rate coefficients but gave instead depopulation Boltzmann averaged cross sections for the rotational transition 1_{01} → 1_{10} as a function of temperature. They neglected the contribution of higher notational levels of H_{2}S and used an approximate twolevel model, which is expected to be valid up to approximately 10 K. They reported a value of 6.84 Å^{2} for the depopulation cross section of the 110 level of H_{2}S at 11.1 K or equivalently a value of 1.75 10^{−11} cm^{3} molecule^{−1} s^{−1} for the corresponding rate at this temperature. Our calculated value is about three times larger (5.7 10^{−11} cm^{3} molecule^{−1} s^{−1}). Such an agreement between theory and experiment, if not satisfactory, is understandable knowing that energy resonances in this domain of energy may drastically change the magnitude of the cross sections in a very narrow energy interval, while possible sources of experimental error are mentioned by these authors. In the same publication Ball et al. (1999) also performed close coupling calculations and obtained a good agreement with experiment at this temperature. This agreement however appears to be fortuitous as they used an approximate model of PES whose well depth is only a half of ours. Other experimental data are then needed to decide if there is a real disagreement between theory and experiment.
To conclude, we also compared in the upper and lower panels of Fig. 5 the rotational deexcitation rate coefficients of H_{2}S by collision with He from the initial states 2 2 1 and 2 2 0 of H_{2}S with those recently published by Dagdigian (Dagdigian 2020) for the collisions with H_{2}. The well depth of the H_{2}H_{2}S complex (146.8 cm^{−1}) calculated by Dagdigian (2020) is five times deeper than that of HeH_{2}S, while the relative masses of the former is about a half of the later. In Fig. 5, the rate coefficients for the two systems are seen to be about of the same magnitude, while differing both in the law variation as a function of temperature and in the propensity rules. The scaling law based on the square root of the inverse ratio of the relative masses therefore does not appear to be a good option to predict the rate coefficients for the collisions of H_{2}S with paraH_{2} from those with He, while the excitation rate coefficients with He need to be included in the modeling of the rotational excitation of H_{2}S.
4. Summary
A threedimensional PES for the collision of H_{2}S+He was developed. The global minimum (−36.49 cm^{−1}) is obtained for a bent geometry and is comparable in both aspects to that of the HeH_{2}O system, while a secondary minimum is identified that does not exist for HeH_{2}O. Close coupling calculations were performed for the 50 lowest rotational states of the two nuclear spin modifications of H_{2}S. At low and intermediate energy, the transition toward levels that are the closest in energy and have the same parity index gives the largest cross section, while for higher energies the nearoblate (ΔK_{C} = 0) propensity rule progressively takes the lead. Our calculation gives a value that is the triple of the only experimental data available for the 1_{01} → 1_{10} depopulation rate coefficient at 10 K. However, this comparison is not conclusive and more experimental data will be needed.
While the existing databases suggested up to now to use the HeH_{2}O statetostate rate coefficients for the HeH_{2}S collisions, we demonstrated that the three temperature dependence, magnitude and propensity rules differ strongly for the two systems. These strong differences result from both the differences in the angular anisotropies of the two PESs and in the rotational constants, where H_{2}S has a nearly oblate symmetric top while H_{2}O is a good example of asymmetric top intermediate case. We also compared our results with those by Dagdigian (2020) for the collisions with H_{2}. The rate coefficients for the two systems were found to be about of the same magnitude, but these rates differ in the law variation as a function of temperature and in the propensity rules. We then recommend using both the new statetostate rate coefficients with He and H_{2} in any nonLTE analysis of astronomical observation of H_{2}S.
Acknowledgments
We acknowledge the support from the ECOSSUDCONICYT project number C17E06 (Programa de Cooperación Científica ECOSCONICYT ECOS170039). We also thank Paul Dagdigian for the valuable discussion and for providing the data for the collision H_{2}S+H_{2}. Computer time for this study was provided by the Mésocentre de Calcul Intensif Aquitain, which is the computing facility of Université de Bordeaux et Université de Pau et des Pays de l’Adour.
References
 Alexander, M. H. 1982, J. Chem. Phys., 77, 1855 [NASA ADS] [CrossRef] [Google Scholar]
 Ball, C. D., Mengel, M., De Lucia, F. C., & Woon, D. E. 1999, J. Chem. Phys., 111, 8893 [NASA ADS] [CrossRef] [Google Scholar]
 Boys, S. F., & Bernardi, F. 1970, Mol. Phys., 19, 553 [NASA ADS] [CrossRef] [Google Scholar]
 Cernicharo, J., Lefloch, B., Agúndez, M., et al. 2018, ApJ, 853, L22 [NASA ADS] [CrossRef] [Google Scholar]
 Cook, R. L., De Lucia, F. C., & Helminger, P. 1975, J. Mol. Struct., 28, 237 [NASA ADS] [CrossRef] [Google Scholar]
 Crockett, N., Bergin, E., Neill, J., et al. 2014, ApJ, 781, 114 [NASA ADS] [CrossRef] [Google Scholar]
 Cybulski, S. M., & Toczyłowski, R. 1999, J. Chem. Phys., 111, 10520 [NASA ADS] [CrossRef] [Google Scholar]
 Dagdigian, P. J. 2020, J. Chem. Phys., 152, 074307 [NASA ADS] [CrossRef] [Google Scholar]
 Drozdovskaya, M. N., van Dishoeck, E. F., Jørgensen, J. K., et al. 2018, MNRAS, 476, 4949 [NASA ADS] [CrossRef] [Google Scholar]
 Dubernet, M., Alexander, M., Ba, Y., et al. 2013, A&A, 553, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Duley, W., Millar, T., & Williams, D. 1980, MNRAS, 192, 945 [NASA ADS] [Google Scholar]
 Dunning, T. H. 1989, J. Chem. Phys., 90, 1007 [NASA ADS] [CrossRef] [Google Scholar]
 Faure, A., Jankowski, P., Stoecklin, T., & Szalewicz, K. 2016, Sci. Rep., 6, 28449 [NASA ADS] [CrossRef] [Google Scholar]
 Fuente, A., Cernicharo, J., Roueff, E., et al. 2016, A&A, 593, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fuente, A., Goicoechea, J. R., Pety, J., et al. 2017, ApJ, 851, L49 [NASA ADS] [CrossRef] [Google Scholar]
 Gao, Z., Karman, T., Tang, G., et al. 2018, Phys. Chem. Chem. Phys., 20, 12444 [CrossRef] [Google Scholar]
 Gao, Z., Loreau, J., van der Avoird, A., & van de Meerakker, S. Y. 2019, Phys. Chem. Chem. Phys., 21, 14033 [CrossRef] [Google Scholar]
 Goicoechea, J. R., Cuadrado, S., Pety, J., et al. 2017, A&A, 601, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gorai, P. 2018, Exploring the Universe: From Near Space to ExtraGalactic (Springer), p467–475 [NASA ADS] [CrossRef] [Google Scholar]
 Ho, T. S., & Rabitz, H. 1996, J. Phys. Chem., 104, 2584 [NASA ADS] [CrossRef] [Google Scholar]
 Holdship, J., Viti, S., JimenezSerra, I., et al. 2016, MNRAS, 463, 802 [NASA ADS] [CrossRef] [Google Scholar]
 Li, J., Wang, J., Zhu, Q., Zhang, J., & Li, D. 2015, ApJ, 802, 40 [NASA ADS] [CrossRef] [Google Scholar]
 Martin, S., MartínPintado, J., Mauersberger, R., Henkel, C., & GarcíaBurillo, S. 2005, ApJ, 620, 210 [NASA ADS] [CrossRef] [Google Scholar]
 Mertens, L. A., Labiad, H., DenisAlpizar, O., et al. 2017, Chem. Phys. Lett., 683, 521 Ahmed Zewail (1946–2016) Commemoration Issue of Chemical Physics Letters [NASA ADS] [CrossRef] [Google Scholar]
 Minh, Y. 2016, J. Phys. Conf. Ser., 728, 052007 [CrossRef] [Google Scholar]
 Minh, Y., Irvine, W. M., & Ziurys, L. M. 1989, ApJ, 345, L63 [NASA ADS] [CrossRef] [Google Scholar]
 Minh, Y., Ziurys, L. M., Irvine, W. M., & McGonagle, D. 1990, ApJ, 360, 136 [NASA ADS] [CrossRef] [Google Scholar]
 Müller, H. S., Schlöder, F., Stutzki, J., & Winnewisser, G. 2005, J. Mol. Struct., 742, 215 [NASA ADS] [CrossRef] [Google Scholar]
 Nagy, Z., Van der Tak, F., Ossenkopf, V., et al. 2013, A&A, 550, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Patkowski, K., Korona, T., Moszynski, R., Jeziorski, B., & Szalewicz, K. 2002, J. Mol. Struct. (Theochem.), 591, 231 [CrossRef] [Google Scholar]
 Phuong, N., Chapillon, E., Majumdar, L., et al. 2018, A&A, 616, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Raghavachari, K., Trucks, G. W., Pople, J. A., & HeadGordon, M. 1989, Chem. Phys. Lett., 157, 479 [NASA ADS] [CrossRef] [Google Scholar]
 Schöier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & Black, J. H. 2005, A&A, 432, 369 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Soldán, P., & Hutson, J. M. 2000, J. Chem. Phys., 112, 4415 [NASA ADS] [CrossRef] [Google Scholar]
 Stoecklin, T., Faure, A., Jankowski, P., et al. 2017, Phys. Chem. Chem. Phys., 19, 189 [NASA ADS] [CrossRef] [Google Scholar]
 Stoecklin, T., DenisAlpizar, O., Clergerie, A., et al. 2019, J. Phys. Chem. A, 123, 5704 [CrossRef] [Google Scholar]
 Thaddeus, P., Kutner, M., Penzias, A., Wilson, R., & Jefferts, K. 1972, ApJ, 176, L73 [NASA ADS] [CrossRef] [Google Scholar]
 Ukita, N., & Morris, M. 1983, A&A, 121, 15 [NASA ADS] [Google Scholar]
 Werner, H.J., Knowles, P. J., Knizia, G., Manby, F. R., & Schütz, M. 2012, WIREs Comput. Mol. Sci., 2, 242 [CrossRef] [Google Scholar]
 Yang, B., Nagao, M., Satomi, W., Kimura, M., & Stancil, P. 2013, ApJ, 765, 77 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1.
Internal coordinates used to describe the H_{2}S+He system. 

In the text 
Fig. 2.
Contour plots of the H_{2}S+He system at φ = 0° (panel A) and φ = 90° (panel B). Red contour lines correspond to positive energies ranging from 0 to 100 cm^{−1} in steps of 10 cm^{−1}, and in steps of 100 cm^{−1} in the [100, 2000] cm^{−1} interval. Blue contour lines show the negative energies varying in steps of 2 cm^{−1}. 

In the text 
Fig. 3.
Rotational deexcitation cross sections of H_{2}S by collision with He. The upper and lower panels show transitions issued from the 2 2 1 and 2 2 0 initial levels of H_{2}S, respectively. 

In the text 
Fig. 4.
Rotational deexcitation cross sections of H_{2}S by collision with He. The upper and lower panels show transitions issued from the 4 4 1 and 4 4 0 initial levels of H_{2}S, respectively. 

In the text 
Fig. 5.
Rotational deexcitation rate coefficients of H_{2}S by collision with He. The upper and lower panels show transitions issued from the 2 2 1 and 2 2 0 initial levels of water, respectively. The rates are compared with the corresponding rates for the H_{2}S + paraH_{2} collisions (designated by pH_{2}) and with those for H_{2}O+He collisions (designated by H_{2}O). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.