SH$^+$ is a surprisingly widespread molecular ion in diffuse interstellar clouds. There, it plays an important role triggering the sulfur chemistry. In addition, SH$^+$ emission lines have been detected at the UV-illuminated edges of dense molecular clouds, \mbox{so-called} photo-dissociation regions (PDRs), and toward high-mass protostars. An accurate determination of the SH$^+$ abundance and of the physical conditions prevailing in these energetic environments relies on knowing the rate coefficients of inelastic collisions between SH$^+$ molecules and hydrogen atoms, hydrogen molecules, and electrons. In this paper, we derive SH$^+$--H fine and hyperfine-resolved rate coefficients from the recent quantum calculations for the SH$^+$--H collisions, including inelastic, exchange and reactive processes. The method used is based on the infinite order sudden approach. State-to-state rate coefficients between the first 31 fine levels and 61 hyperfine levels of SH$^+$ were obtained for temperatures ranging from 10 to 1000 K. Fine structure-resolved rate coefficients present a strong propensity rule in favour of $\Delta j = \Delta N$ transitions. The $\Delta j = \Delta F$ propensity rule is observed for the hyperfine transitions. {The new rate coefficients will help significantly in the interpretation of SH$^+$ spectra from PDRs and UV-irradiated shocks where the abundance of hydrogen atoms with respect to hydrogen molecules can be significant.

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

The rate constants for the formation, destruction, and collisional excitation of SH$^+$ are calculated from quantum mechanical approaches using two new SH$_2^+$ potential energy surfaces (PESs) of $^4A''$ and $^2A''$ electronic symmetry. The PESs were developed to describe all adiabatic states correlating to the SH$^+$ ($^3\Sigma^-$) + H($^2S$) channel. The formation of SH$^+$ through the S$^+$ + H$_2$ reaction is endothermic by $\approx$ 9860 K, and requires at least two vibrational quanta on the H$_2$ molecule to yield significant reactivity. Quasi-classical calculations of the total formation rate constant for H$_2$($v=2$) are in very good agreement with the quantum results above 100K. Further quasi-classical calculations are then performed for $v=3$, 4, and 5 to cover all vibrationally excited H$_2$ levels significantly populated in dense photodissociation regions (PDR). The new calculated formation and destruction rate constants are two to six times larger than the previous ones and have been introduced in the Meudon PDR code to simulate the physical and illuminating conditions in the Orion bar prototypical PDR. New astrochemical models based on the new molecular data produce four times larger SH$^+$ column densities, in agreement with those inferred from recent ALMA observations of the Orion bar.

GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) to determine the fractional abundances of CO, HCO+, HCN, CS, SO, HCS+, and N2H+ in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from Av~3 to ~20mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of (1-5)x10$^3$ cm$^{-3}$; and ii) the dense phase, located at Av>10mag, with molecular hydrogen densities >10$^4$ cm$^{-3}$. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone where C/H~8x10$^{-5}$ and C/O~0.8-1, until the beginning of the dense phase at Av~10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H~(0.4 - 2.2) x10$^{-6}$, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H~8x10$^{-8}$). Based on our chemical modeling, we constrain the value of $\zeta_{\rm H_2}$ to ~(0.5 - 1.8) x10$^{-16}$ s$^{-1}$ in the translucent cloud.

The presence of SiS in space seems to be restricted to a few selected types of astronomical environments. It is long known to be present in circumstellar envelopes around evolved stars and it has also been detected in a handful of star-forming regions with evidence of outflows, like Sgr\,B2, Orion\,KL and more recently L1157-B1. The kinetics of reactions involving SiS is very poorly known and here we revisit the chemistry of SiS in space by studying some potentially important reactions of formation and destruction of this molecule. We calculated {\it ab initio} potential energy surfaces of the SiOS system and computed rate coefficients in the temperature range 50-2500 K for the reaction of destruction of SiS, in collisions with atomic O, and of its formation, through the reaction between Si and SO. We find that both reactions are rapid, with rate coefficients of a few times 10$^{-10}$ cm$^3$ s$^{-1}$, almost independent of temperature. In the reaction between Si and SO, SiO production is 5-7 times more efficient than SiS formation. The reaction of SiS with O atoms can play an important role in destroying SiS in envelopes around evolved stars. We built a simple chemical model of a postshock gas to study the chemistry of SiS in protostellar outflows and we found that SiS forms with a lower abundance and later than SiO, that SiS is efficiently destroyed through reaction with O, and that the main SiS-forming reactions are Si + SO and Si + SO$_2$.

We present a detailed theoretical study of the rotational excitation of CH$^+$ due to reactive and nonreactive collisions involving C$^+(^2P)$, H$_2$, CH$^+$, H and free electrons. Specifically, the formation of CH$^+$ proceeds through the reaction between C$^+(^2P)$ and H$_2(\nu_{\rm H_2}=1, 2)$, while the collisional (de)excitation and destruction of CH$^+$ is due to collisions with hydrogen atoms and free electrons. State-to-state and initial-state-specific rate coefficients are computed in the kinetic temperature range 10-3000~K for the inelastic, exchange, abstraction and dissociative recombination processes using accurate potential energy surfaces and the best scattering methods. Good agreement, within a factor of 2, is found between the experimental and theoretical thermal rate coefficients, except for the reaction of CH$^+$ with H atoms at kinetic temperatures below 50~K. The full set of collisional and chemical data are then implemented in a radiative transfer model. Our Non-LTE calculations confirm that the formation pumping due to vibrationally excited H$_2$ has a substantial effect on the excitation of CH$^+$ in photon-dominated regions. In addition, we are able to reproduce, within error bars, the far-infrared observations of CH$^+$ toward the Orion Bar and the planetary nebula NGC~7027. Our results further suggest that the population of $\nu_{\rm H_2}=2$ might be significant in the photon-dominated region of NGC~7027.

As many organic molecules, formic acid (HCOOH) has two conformers (trans and cis). The energy barrier to internal conversion from trans to cis is much higher than the thermal energy available in molecular clouds. Thus, only the most stable conformer (trans) is expected to exist in detectable amounts. We report the first interstellar detection of cis-HCOOH. Its presence in ultraviolet (UV) irradiated gas exclusively (the Orion Bar photodissociation region), with a low trans-to-cis abundance ratio of 2.8+-1.0, supports a photoswitching mechanism: a given conformer absorbs a stellar photon that radiatively excites the molecule to electronic states above the interconversion barrier. Subsequent fluorescent decay leaves the molecule in a different conformer form. This mechanism, which we specifically study with ab initio quantum calculations, was not considered in Space before but likely induces structural changes of a variety of interstellar molecules submitted to UV radiation.

Because of the very peculiar conditions of chemistry in many astrophysical gases (low densities, mostly low temperatures, kinetics-dominated chemical evolution), great efforts have been devoted to study molecular signatures and chemical evolution. While experiments are being performed in many laboratories, it appears that the efforts directed towards theoretical works are not as strong. This report deals with the present status of chemical physics/physical chemistry theory, for the qualitative and quantitative understanding of kinetics of molecular scattering, being it reactive or inelastic. By gathering several types of expertise, from applied mathematics to physical chemistry, dialog is made possible, as a step towards new and more adapted theoretical frameworks, capable of meeting the theoretical, methodological and numerical challenges of kinetics-dominated gas phase chemistry in astrophysical environments. A state of the art panorama is presented, alongside present-day strengths and shortcomings. However, coverage is not complete, being limited in this report to actual attendance of the workshop. Some paths towards relevant progress are proposed.

Barnard B1b has revealed as one of the most interesting globules from the chemical and dynamical point of view. It presents a rich molecular chemistry characterized by large abundances of deuterated and complex molecules. Furthermore, it hosts an extremely young Class 0 object and one candidate to First Hydrostatic Core (FHSC). Our aim was to determine the cosmic ray ionization rate and the depletion factors in this extremely young star forming region. We carried out a spectral survey towards Barnard 1b as part of the IRAM Large program ASAI using the IRAM 30-m telescope at Pico Veleta (Spain). This provided a very complete inventory of neutral and ionic C-, N- and S- bearing species with, up to our knowledge, the first secure detections of the deuterated ions DCS+ and DOCO+. We used a state-of-the-art pseudo-time-dependent gas-phase chemical model to determine the value of the cosmic ray ionization rate and the depletion factors. The observational data were well fitted with $\zeta_{H_2}$ between 3E-17 s$^{-1}$ and 1E-16 s$^{-1}$. Elemental depletions were estimated to be ~10 for C and O, ~1 for N and ~25 for S. Barnard B1b presents similar depletions of C and O than those measured in pre-stellar cores. The depletion of sulfur is higher than that of C and O but not as extreme as in cold cores. In fact, it is similar to the values found in some bipolar outflows, hot cores and photon-dominated regions. Several scenarios are discussed to account for these peculiar abundances. We propose that it is the consequence of the initial conditions (important outflows and enhanced UV fields in the surroundings) and a rapid collapse (~0.1 Myr) that permits to maintain most S- and N-bearing species in gas phase to great optical depths. The interaction of the compact outflow associated with B1b-S with the surrounding material could enhance the abundances of S-bearing molecules, as well.

The rate constants required to model the OH$^+$ observations in different regions of the interstellar medium have been determined using state of the art quantum methods. First, state-to-state rate constants for the H$_2(v=0,J=0,1)$+ O$^+$($^4S$) $\rightarrow$ H + OH$^+(X ^3\Sigma^-, v', N)$ reaction have been obtained using a quantum wave packet method. The calculations have been compared with time-independent results to asses the accuracy of reaction probabilities at collision energies of about 1 meV. The good agreement between the simulations and the existing experimental cross sections in the $0.01-$1 eV energy range shows the quality of the results. The calculated state-to-state rate constants have been fitted to an analytical form. Second, the Einstein coefficients of OH$^+$ have been obtained for all astronomically significant ro-vibrational bands involving the $X^3\Sigma^-$ and/or $A^3\Pi$ electronic states. For this purpose the potential energy curves and electric dipole transition moments for seven electronic states of OH$^+$ are calculated with {\it ab initio} methods at the highest level and including spin-orbit terms, and the rovibrational levels have been calculated including the empirical spin-rotation and spin-spin terms. Third, the state-to-state rate constants for inelastic collisions between He and OH$^+(X ^3\Sigma^-)$ have been calculated using a time-independent close coupling method on a new potential energy surface. All these rates have been implemented in detailed chemical and radiative transfer models. Applications of these models to various astronomical sources show that inelastic collisions dominate the excitation of the rotational levels of OH$^+$. In the models considered the excitation resulting from the chemical formation of OH$^+$ increases the line fluxes by about 10 % or less depending on the density of the gas.

Depolarizing collisions are elastic or quasielastic collisions that equalize the populations and destroy the coherence between the magnetic sublevels of atomic levels. In astrophysical plasmas, the main depolarizing collider is neutral hydrogen. We consider depolarizing rates on the lowest levels of neutral and singly ionized alkaly-earths Mg I, Sr I, Ba I, Mg II, Ca II, and Ba II, due to collisions with H. We compute ab initio potential curves of the atom-H system and solve the quantum mechanical dynamics. From the scattering amplitudes we calculate the depolarizing rates for Maxwellian distributions of colliders at temperatures T <10000 K. A comparative analysis of our results and previous calculations in the literature is done. We discuss the effect of these rates on the formation of scattering polarization patterns of resonant lines of alkali-earths in the solar atmosphere, and their effect on Hanle effect diagnostics of solar magnetic fields.