We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

The Leiden Atomic and Molecular Database (LAMDA) collects spectroscopic information and collisional rate coefficients for molecules, atoms, and ions of astrophysical and astrochemical interest. We describe the developments of the database since its inception in 2005, and outline our plans for the near future. Such a database is constrained both by the nature of its uses and by the availability of accurate data: we suggest ways to improve the synergies among users and suppliers of data. We summarize some recent developments in computation of collisional cross sections and rate coefficients. We consider atomic and molecular data that are needed to support astrophysics and astrochemistry with upcoming instruments that operate in the mid- and far-infrared parts of the spectrum.

The measured nitrogen-to-carbon ratio in comets is lower than for the Sun, a discrepancy which could be alleviated if there is an unknown reservoir of nitrogen in comets. The nucleus of comet 67P/Churyumov-Gerasimenko exhibits an unidentified broad spectral reflectance feature around 3.2 micrometers, which is ubiquitous across its surface. On the basis of laboratory experiments, we attribute this absorption band to ammonium salts mixed with dust on the surface. The depth of the band indicates that semivolatile ammonium salts are a substantial reservoir of nitrogen in the comet, potentially dominating over refractory organic matter and more volatile species. Similar absorption features appear in the spectra of some asteroids, implying a compositional link between asteroids, comets, and the parent interstellar cloud.

We present the first accurate rate coefficients for the rotational excitation of CO by H2O in the kinetic temperature range 5-100 K. The statistical adiabatic channel method (SACM) is combined with a high-level rigid-rotor CO-H2O intermolecular potential surface. Transitions among the first 11 rotational levels of CO and the first 8 rotational levels of both para-H2O and ortho-H2O are considered. Our rate coefficients are compared to previous data from the literature and they are also incorporated in a simple non-LTE model of cometary coma including collision-induced transitions, solar radiative pumping and radiative decay. We find that the uncertainties in the collision data have significant influence on the CO population distribution for H2O densities in the range 10^3-10^8 cm^-3. We also show that the rotational distribution of H2O plays an important role in CO excitation (owing to correlated energy transfer in both CO and H2O), while the impact of the ortho-to-para ratio of H2O is found to be negligible.

Using APEX-1 and APEX-2 observations, we have detected and studied the rotational lines of the HC$_3$N molecule (cyanoacetylene) in the powerful outflow/hot molecular core G331.512-0.103. We identified thirty-one rotational lines at $J$ levels between 24 and 39; seventeen of them in the ground vibrational state $v$=0 (9 lines corresponding to the main C isotopologue and 8 lines corresponding to the $^{13}$C isotopologues), and fourteen in the lowest vibrationally excited state $v_7$=1. Using LTE-based population diagrams for the beam-diluted $v$=0 transitions, we determined $T_{\rm exc}$=85$\pm$4 K and $N$(HC$_3$N)=(6.9$\pm$0.8)$\times$10$^{14}$ cm$^{-2}$, while for the beam-diluted $v_7$=1 transitions we obtained $T_{\rm exc}$=89$\pm$10 K and $N$(HC$_3$N)=2$\pm$1$\times$10$^{15}$ cm$^{-2}$. Non-LTE calculations using H$_2$ collision rates indicate that the HC$_3$N emission is in good agreement with LTE-based results. From the non-LTE method we estimated $T_{\rm kin}$ $\simeq$90~K, $n$(H$_2$)$\simeq$2$\times$10$^7$~cm$^{-3}$ for a central core of 6 arcsec in size. A vibrational temperature in the range from 130~K to 145~K was also determined, values which are very likely lower limits. Our results suggest that rotational transitions are thermalized, while IR radiative pumping processes are probably more efficient than collisions in exciting the molecule to the vibrationally excited state $v_7$=1. Abundance ratios derived under LTE conditions for the $^{13}$C isotopologues suggest that the main formation pathway of HC$_3$N is ${\rm C}_2{\rm H}_2 + {\rm CN} \rightarrow {\rm HC}_3{\rm N} + {\rm H}$.

We investigate the deuteration of methanol towards the high-mass star forming region NGC 7538-IRS1. We have carried out a multi-transition study of CH$_3$OH, $^{13}$CH$_3$OH and of the deuterated fllavors, CH$_2$DOH and CH$_3$OD, between 1.0--1.4 mm with the IRAM-30~m antenna. In total, 34 $^{13}$CH$_3$OH, 13 CH$_2$DOH lines and 20 CH$_3$OD lines spanning a wide range of upper-state energies (E$_{up}$) were detected. From the detected transitions, we estimate that the measured D/H does not exceed 1$\%$, with a measured CH$_2$DOH/CH$_3$OH and CH$_3$OD/CH$_3$OH of about (32$\pm$8)$\times$10$^{-4}$ and (10$\pm$4)$\times$10$^{-4}$, respectively. This finding is consistent with the hypothesis of a short-time scale formation during the pre-stellar phase. We find a relative abundance ratio CH$_2$DOH/CH$_3$OD of 3.2 $\pm$ 1.5. This result is consistent with a statistical deuteration. We cannot exclude H/D exchanges between water and methanol if water deuteration is of the order 0.1$\%$, as suggested by recent Herschel observations.

We report extensive theoretical calculations on the rotation-inversion excitation of interstellar ammonia (NH3) due to collisions with atomic and molecular hydrogen (both para- and ortho-H2). Close-coupling calculations are performed for total energies in the range 1-2000 cm-1 and rotational cross sections are obtained for all transitions among the lowest 17 and 34 rotation-inversion levels of ortho- and para-NH3, respectively. Rate coefficients are deduced for kinetic temperatures up to 200 K. Propensity rules for the three colliding partners are discussed and we also compare the new results to previous calculations for the spherically symmetrical He and para-H2 projectiles. Significant differences are found between the different sets of calculations. Finally, we test the impact of the new rate coefficients on the calibration of the ammonia thermometer. We find that the calibration curve is only weakly sensitive to the colliding partner and we confirm that the ammonia thermometer is robust.

We present an analysis of gas densities in the central R=300 parsecs of the Milky Way, focusing on three clouds: GCM-0.02-0.07 (the 50 km/s cloud), GCM-0.13-0.08 (the 20 km/s cloud), and GCM0.25+0.01 (the "Brick"). Densities are determined using observations of the J=(3-2), (4-3), (5-4), (10-9), (18-17), (19-18), (21-20), and (24-23) transitions of the molecule HC3N. We find evidence of at least two excitation regimes for HC3N and constrain the low-excitation component to have a density less than 10^4 cm^-3 and the high-excitation component to have a density between 10^5 and 10^6 cm^-3. This is much less than densities of 10^7 cm^-3 that are found in Sgr B2, the most actively star-forming cloud in the Galactic center. This is consistent with the requirement of a higher density threshold for star formation in the Galactic center than is typical in the Galactic disk. We are also able to constrain the column density of each component in order to determine the mass fraction of 'dense' (n>10^5 cm^-3) gas for these clouds. We find that this is ~15% for all three clouds. Applying the results of our models to ratios of the (10-9) and (3-2) line across the entire central R=300 pc, we find that the fraction of gas with n>10^4 cm^-3 increases inward of a radius of ~140 pc, consistent with the predictions of recent models for the gas dynamics in this region. Our observations show that HC3N is an excellent molecule for probing the density structure of clouds in the Galactic center.

The origin of the heavily fractionated reservoir of nitrogen in comets remains an issue in the theory of their formation and hence of the solar system. Whether the fractionated reservoir traced by comets is inherited from the interstellar cloud or is the product of processes taking place in the protostar, or in the protoplanetary disk, remains unclear. So far, observations of nitrogen isotopic ratios in protostars or prestellar cores have not securely identified such a fractionated reservoir owing to the intrinsic difficulty of direct isotopic ratios measurements. In this article, we report the detection of 5 rotational lines of \ce{HC3N}, {including the weaker components of the hyperfine multiplets}, and two rotational lines of its $^{15}$N isotopologue, towards the L1544 prestellar core. Based on a MCMC/non-LTE multi-line analysis at the hyperfine level, we derive the column densities of \ce{HC3N} ($8.0\pm0.4\tdix{13}$\cc) and \ce{HC3^{15}N} ($2.0\pm0.4\tdix{11}$\cc) and derive an isotopic ratio of 400$\pm$20(1$\sigma$). This value suggests that \ce{HC3N} is slightly depleted in $^{15}$N in L1544 with respect to the elemental $^{14}$N/$^{15}$N ratio {of $\approx$330} in the present-day local interstellar medium. Our study also stresses the need for radiative calculations at the hyperfine level. Finally, the comparison of the derived ratio with those obtained in CN and HCN in the same core seems to favor CN+C$_2$H$_2$ as the dominant formation route to HC$_3$N. However, uncertainties in the isotopic ratios preclude definitive conclusions.

The inelastic cross sections and rate coefficients for the rotational excitation of methanimine (CH$_2$NH) by cold H$_2$ have been determined quantum mechanically based on a new highly correlated five-dimensional potential energy surface. This surface was fitted to more than 60 000 ab initio points with a root mean square error of $\sim$ 1-2 cm$^{-1}$ in the region of the potential well whose depth is -374.0 cm$^{-1}$. The rotationally inelastic rate coefficients have been combined with spectroscopic data to generate the non-equilibrium spectrum of CH$_2$NH towards the interstellar molecular cloud Sgr~B2(N). The transition $1_{10}\to 1_{11}$ at 5.29 GHz is found to be inverted and the predicted weak maser emission spectrum is in excellent agreement with new observations performed with the 100-meter Green Bank Telescope.

We study the gravitational collapse of prestellar sources and the associated evolution of their chemical composition. We use the University of Grenoble Alpes Astrochemical Network (UGAN), which includes reactions involving the different nuclear--spin states of H2, H3+, and of the hydrides of carbon, nitrogen, oxygen, and sulfur, for reactions involving up to seven protons. In addition, species-to-species rate coefficients are provided for the ortho/para interconversion of the H3+ + H2 system and isotopic variants. The composition of the medium is followed from an initial steady state through the early phase of isothermal gravitational collapse. Both the freeze--out of the molecules on to grains and the coagulation of the grains were incorporated in the model. The predicted abundances and column densities of the spin isomers of ammonia and its deuterated forms are compared with those measured recently towards the prestellar cores H-MM1, L16293E, and Barnard B1. We find that gas--phase processes alone account satisfactorily for the observations, without recourse to grain-surface reactions. In particular, our model reproduces both the isotopologue abundance ratios and the ortho:para ratios of NH2D and NHD2 within observational uncertainties. More accurate observations are necessary to distinguish between full scrambling processes---as assumed in our gas-phase network---and direct nucleus- or atom-exchange reactions.

The isotopic ratio of nitrogen in nearby protoplanetary disks, recently measured in CN and HCN, indicates that a fractionated reservoir of nitrogen is available at the stage of comet formation. This reservoir presents a 3:1 enrichment in 15N relative to the elemental ratio of 330, identical to that between comets and the sun, suggesting that similar processes are responsible for the fractionation in the protosolar nebula (PSN) and in PSN analogs. However, how the fractionation of nitrogen takes place is an open question. Previously obtained HCN/HC15N abundance ratios suggest that HCN may be enriched in 15N in prestellar cores, although doubts remain on these measurements, which rely on the double-isotopologue method. Here we present direct measurements of the HCN/H13CN and HCN/HC15N abundance ratios in the L1498 prestellar core based on spatially resolved spectra of HCN(1-0), (3-2), H13CN(1-0), and HC15N(1-0). We use state-of-the-art radiative transfer calculations using ALICO, a 1D radiative transfer code capable of treating hyperfine (hf) overlaps. From an analysis of dust emission maps of L1498, we derive a new physical structure. We also use new, high-accuracy HCN-H2 hf collisional rates, which enable us to quantitatively reproduce all the features of the HCN(1-0) and HCN(3-2) spectra, especially the anomalous hf line ratios. Special attention is devoted to derive meaningful uncertainties on the abundance ratios. The obtained values, HCN/H13CN=45(3) and HCN/HC15N=338(28), indicate that carbon is fractionated in HCN, but nitrogen is not. For the HCN/HC15N abundance ratio, our study validates to some extent analyses based on the single excitation temperature assumption. Comparisons with other measurements suggest significant core-to-core variability. Furthermore, the 13C enrichment we found in HCN could explain the superfractionation of nitrogen measured in solar system chondrites.

R-matrix calculations combined with the adiabatic-nuclei-rotation and Coulomb-Born approximations are used to compute electron-impact rotational rate coefficients for two open-shell diatomic cations of astrophysical interest: the hydoxyl and sulphanyl ions, OH$^+$ and SH$^+$. Hyperfine resolved rate coefficients are deduced using the infinite-order-sudden approximation. The propensity rule $\Delta F=\Delta j=\Delta N=\pm 1$ is observed, as is expected for cations with a large dipole moment. A model for OH$^+$ excitation in the Orion Bar photon-dominated region (PDR) is presented which nicely reproduces Herschel observations for an electron fraction $x_e=10^{-4}$ and an OH$^+$ column density of $3\times 10^{13}$~cm$^{-2}$. Electron impact electronic excitation cross sections and rate coefficients for the ions are also presented.

Our knowledge about the "cold" Universe often relies on molecular spectra. A general property of such spectra is that the energy level populations are rarely at local thermodynamic equilibrium. Solving the radiative transfer thus requires the availability of collisional rate coefficients with the main colliding partners over the temperature range 10-1000 K. These rate coefficients are notoriously difficult to measure and expensive to compute. In particular, very few reliable collisional data exist for collisions involving reactive radicals or ions. Here we explore the use of a fast quantum statistical method to determine molecular collisional excitation rate coefficients. The method is benchmarked against accurate (but costly) close-coupling calculations. For collisions proceeding through the formation of a strongly-bound complex, the method is found to be highly satisfactory up to room temperature. Its accuracy decreases with the potential well depth and with increasing temperature, as expected. This new method opens the way to the determination of accurate inelastic collisional data involving key reactive species such as H3+, H2O+, and H3O+ for which exact quantum calculations are currently not feasible.

Isotopic ratios are keys to understanding the origin and early evolution of the solar system in the context of Galactic nucleosynthesis. The large range of measured $^{14}$N/$^{15}$N isotopic ratios in the solar system reflects distinct reservoirs of nitrogen whose origins remain to be determined. We have directly measured a C$^{14}$N/C$^{15}$N abundance ratio of 323$\pm$30 in the disk orbiting the nearby young star TW Hya. This value, which is in good agreement with nitrogen isotopic ratios measured for prestellar cores, likely reflects the primary present-day reservoir of nitrogen in the solar neighbourhood. These results support models invoking novae as primary $^{15}$N sources as well as outward migration of the Sun over its lifetime, and suggest that comets sampled a secondary, $^{15}$N-rich reservoir during solar system formation.

Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity and spectral resolution to further investigate physical and chemical processes in pre-stellar cores. Here we present the velocity-resolved Herschel-HIFI observations of the ortho-NH_3(1_0-0_0) line at 572 GHz and study the abundance profile of ammonia across the pre-stellar core L1544 to test current theories of its physical and chemical structure. Recently calculated collisional coefficients have been included in our non-LTE radiative transfer code to reproduce Herschel observations. A gas-grain chemical model, including spin-state chemistry and applied to the (static) physical structure of L1544 is also used to infer the abundance profile of ortho-NH_3 . The hyperfine structure of ortho-NH_3(1_0-0_0) is resolved for the first time in space. All the hyperfine components are strongly self-absorbed. The profile can be reproduced if the core is contracting in quasi-equilibrium, consistent with previous work, and if the NH_3 abundance is slightly rising toward the core centre, as deduced from previous interferometric observations of para-NH_3(1,1). The chemical model overestimates the NH_3 abundance at radii between ~ 4000 and 15000 AU by about two orders of magnitude and underestimates the abundance toward the core centre by more than one order of magnitude. Our observations show that chemical models applied to static clouds have problems in reproducing NH_3 observations.

Cold gas-phase water has recently been detected in a cold dark cloud, Barnard 5 located in the Perseus complex, by targeting methanol peaks as signposts for ice mantle evaporation. Observed morphology and abundances of methanol and water are consistent with a transient non-thermal evaporation process only affecting the outermost ice mantle layers, possibly triggering a more complex chemistry. We present the detection of the Complex Organic Molecules (COMs) acetaldehyde and methyl formate as well as formic acid and ketene, and the tentative detection of di-methyl ether towards the methanol hotspot of Barnard 5 located between two dense cores using the single dish OSO 20m, IRAM 30m, and NRO 45m telescopes. The high energy cis- conformer of formic acid is detected, suggesting that formic acid is mostly formed at the surface of interstellar grains and then evaporated. The detection of multiple transitions for each species allows us to constrain their abundances through LTE and non-LTE methods. All the considered COMs show similar abundances between $\sim 1$ and $\sim 10$ % relative to methanol depending on the assumed excitation temperature. The non-detection of glycolaldehyde, an isomer of methyl formate, with a [glycolaldehyde]/[methyl formate] abundance ratio lower than 6 %, favours gas phase formation pathways triggered by methanol evaporation. According to their excitation temperatures derived in massive hot cores, formic acid, ketene, and acetaldehyde have been designated as "lukewarm" COMs whereas methyl formate and di-methyl ether were defined as "warm" species. Comparison with previous observations of other types of sources confirms that "lukewarm" and "warm" COMs show similar abundances in low-density cold gas whereas the "warm" COMs tend to be more abundant than the "lukewarm" species in warm protostellar cores.

We report the discovery of water maser emission at frequencies above 1 THz. Using the GREAT instrument on SOFIA, we have detected emission in the 1.296411 THz 8(27)-7(34) transition of water toward three oxygen-rich evolved stars: W Hya, U Her, and VY CMa. An upper limit on the 1.296 THz line flux was obtained toward R Aql. Near-simultaneous observations of the 22.23508 GHz 6(16)-5(23) water maser transition were carried out towards all four sources using the Effelsberg 100m telescope. The measured line fluxes imply 22 GHz / 1.296 THz photon luminosity ratios of 0.012, 0.12, and 0.83 respectively for W Hya, U Her, and VY CMa, values that confirm the 22 GHz maser transition to be unsaturated in W Hya and U Her. We also detected the 1.884888 THz 8(45)-7(53) transition toward W Hya and VY CMa, and the 1.278266 THz 7(43)-6(52) transition toward VY CMa. Like the 22 GHz maser transition, all three of the THz emission lines detected here originate from the ortho-H2O spin isomer. Based upon a model for the circumstellar envelope of W Hya, we estimate that stimulated emission is responsible for ~ 85% of the observed 1.296 THz line emission, and thus that this transition may be properly described as a terahertz-frequency maser. In the case of the 1.885 THz transition, by contrast, our W Hya model indicates that the observed emission is dominated by spontaneous radiative decay, even though a population inversion exists.

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.

Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud.

Cyanopolyynes are chains of carbon atoms with an atom of hydrogen and a CN group on either side. They are detected almost everywhere in the ISM, as well as in comets. In the past, they have been used to constrain the age of some molecular clouds, since their abundance is predicted to be a strong function of time. We present an extensive study of the cyanopolyynes distribution in the solar-type protostar IRAS16293-2422 based on TIMASSS IRAM-30m observations. The goals are (i) to obtain a census of the cyanopolyynes in this source and of their isotopologues; (ii) to derive how their abundance varies across the protostar envelope; and (iii) to obtain constraints on the history of IRAS16293-2422. We detect several lines from HC3N and HC5N, and report the first detection of DC3N, in a solar-type protostar. We found that the HC3N abundance is roughly constant (~1.3x10^(-11)) in the outer cold envelope of IRAS16293-2422, and it increases by about a factor 100 in the inner region where Tdust>80K. The HC5N has an abundance similar to HC3N in the outer envelope and about a factor of ten lower in the inner region. The HC3N abundance derived in the inner region, and where the increase occurs, also provide strong constraints on the time taken for the dust to warm up to 80K, which has to be shorter than ~10^3-10^4yr. Finally, the cyanoacetylene deuteration is about 50\% in the outer envelope and <5$\% in the warm inner region. The relatively low deuteration in the warm region suggests that we are witnessing a fossil of the HC3N abundantly formed in the tenuous phase of the pre-collapse and then frozen into the grain mantles at a later phase. The accurate analysis of the cyanopolyynes in IRAS16293-2422 unveils an important part of its past story. It tells us that IRAS16293-2422 underwent a relatively fast (<10^5yr) collapse and a very fast (<10^3-10^4yr) warming up of the cold material to 80K.

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.

New calculations for rotational excitation of cyanoacetylene by collisions with hydrogen molecules are performed to include the lowest 38 rotational levels of HC3N and kinetic temperatures to 300 K. Calculations are based on the interaction potential of Wernli et al. A&A, 464, 1147 (2007) whose accuracy is checked against spectroscopic measurements of the HC3N-H2 complex. The quantum coupled-channel approach is employed and complemented by quasi-classical trajectory calculations. Rate coefficients for ortho-H2 are provided for the first time. Hyperfine resolved rate coefficients are also deduced. Collisional propensity rules are discussed and comparisons between quantum and classical rate coefficients are presented. This collisional data should prove useful in interpreting HC3N observations in the cold and warm ISM, as well as in protoplanetary disks.

Understanding the processes that could lead to enrichment of molecules in $^{15}$N atoms is of particular interest in order to shed light on the relatively large variations observed in the $^{14}$N/$^{15}$N ratio in various solar system environments. Currently, the sample of molecular clouds where $^{14}$N/$^{15}$N ratios have been measured is small and has to be enlarged in order to allow statistically significant studies. In particular, the N$_2$H$^+$ molecule currently shows the largest spread of $^{14}$N/$^{15}$N ratios in high-mass star forming regions. However, the $^{14}$N/$^{15}$N ratio in N$_2$H$^+$ was obtained in only two low-mass star forming regions (L1544 and B1b). The current work extends this sample to a third dark cloud. We targeted the 16293E prestellar core, where the N$^{15}$NH$^+$ $J$=1-0 line was detected. Using a model previously developed for the physical structure of the source, we solved the molecular excitation with a non-local radiative transfer code. For that purpose, we computed specific collisional rate coefficients for the N$^{15}$NH$^+$-H$_2$ collisional system. As a first step of the analysis, the N$_2$H$^+$ abundance profile was constrained by reproducing the N$_2$H$^+$ $J$=1-0 and 3-2 maps. A scaling factor was then applied to this profile to match the N$^{15}$NH$^+$ $J$=1-0 spectrum. We derive a column density ratio N$_2$H$^+$ / N$^{15}$NH$^+$ = $330^{+170}_{-100}$. The current estimate $\sim$330 agrees with the value typical of the elemental isotopic ratio in the local ISM. It is however lower than in some other cores, where values as high as 1300 have been reported.

Complex organic molecules (COMs) are detected in many regions of the interstellar medium, including prestellar cores. However, their formation mechanisms in cold (~10 K) cores remain to this date poorly understood. The formyl radical HCO is an important candidate precursor for several O-bearing terrestrial COMs in cores, as an abundant building block of many of these molecules. Several chemical routes have been proposed to account for its formation, both on grain surfaces, as an incompletely hydrogenated product of H addition to frozen-out CO molecules, or in the gas phase, either the product of the reaction between H2CO and a radical, or as a product of dissociative recombination of protonated formaldehyde H2COH+. The detection and abundance determination of H2COH+, if present, could provide clues as to whether this latter scenario might apply. We searched for protonated formaldehyde H2COH+ in the prestellar core L1689B using the IRAM 30m telescope. The H2COH+ ion is unambiguously detected, for the first time in a cold (~10 K) source. The derived abundance agrees with a scenario in which the formation of H2COH+ results from the protonation of formaldehyde. We use this abundance value to constrain the branching ratio of the dissociative recombination of H2COH+ towards the HCO channel to ~10-30%. This value could however be smaller if HCO can be efficiently formed from gas-phase neutral-neutral reactions, and we stress the need for laboratory measurements of the rate constants of these reactions at 10 K. Given the experimental difficulties in measuring branching ratios experimentally, observations can bring valuable constraints on these values, and provide a useful input for chemical networks.

The availability of collisional rate coefficients is a prerequisite for an accurate interpretation of astrophysical observations, since the observed media often harbour densities where molecules are populated under non--LTE conditions. In the current study, we present calculations of rate coefficients suitable to describe the various spin isomers of multiply deuterated ammonia, namely the ND$_2$H and ND$_3$ isotopologues. These calculations are based on the most accurate NH$_3$--H$_2$ potential energy surface available, which has been modified to describe the geometrical changes induced by the nuclear substitutions. The dynamical calculations are performed within the close--coupling formalism and are carried out in order to provide rate coefficients up to a temperature of $T$ = 50K. For the various isotopologues/symmetries, we provide rate coefficients for the energy levels below $\sim$ 100 cm$^{-1}$. Subsequently, these new rate coefficients are used in astrophysical models aimed at reproducing the NH$_2$D, ND$_2$H and ND$_3$ observations previously reported towards the prestellar cores B1b and 16293E. We thus update the estimates of the corresponding column densities and find a reasonable agreement with the previous models. In particular, the ortho--to--para ratios of NH$_2$D and NHD$_2$ are found to be consistent with the statistical ratios.

The recent unexpected detection of terrestrial complex organic molecules in the cold (~ 10 K) gas has cast doubts on the commonly accepted formation mechanisms of these species. Standard gas-phase mechanisms are inefficient and tend to underproduce these molecules, and many of the key reactions involved are unconstrained. Grain-surface mechanisms, which were presented as a viable alternative, suffer from the fact that they rely on grain surface diffusion of heavy radicals, which is not possible thermally at very low temperatures. One of the simplest terrestrial complex organic molecules, methanol is believed to form on cold grain surfaces following from successive H atom additions on CO. Unlike heavier species, H atoms are very mobile on grain surfaces even at 10 K. Intermediate species involved in grain surface methanol formation by CO hydrogenation are the radicals HCO and CH3O, as well as the stable species formaldehyde H2CO. These radicals are thought to be precursors of complex organic molecules on grain surfaces. We present new observations of the HCO and CH3O radicals in a sample of prestellar cores and carry out an analysis of the abundances of the species HCO, H2CO, CH3O, and CH3OH, which represent the various stages of grain-surface hydrogenation of CO to CH3OH. The abundance ratios between the various intermediate species in the hydrogenation reaction of CO on grains are similar in all sources of our sample, HCO:H2CO:CH3O:CH3OH ~ 10:100:1:100. We argue that these ratios may not be representative of the primordial abundances on the grains but, rather, suggest that the radicals HCO and CH3O are gas-phase products of the precursors H2CO and CH3OH, respectively. Gas-phase pathways are considered and simple estimates of HCO and CH3O abundances are compared to the observations. Critical reaction rate constants, branching ratios, and intermediate species are finally identified.

The 1$_{11}$-1$_{01}$ lines of ortho and para--NH$_2$D (o/p-NH$_2$D), respectively at 86 and 110 GHz, are commonly observed to provide constraints on the deuterium fractionation in the interstellar medium. In cold regions, the hyperfine structure due to the nitrogen ($^{14}$N) nucleus is resolved. To date, this splitting is the only one which is taken into account in the NH$_2$D column density estimates. We investigate how the inclusion of the hyperfine splitting caused by the deuterium (D) nucleus affects the analysis of the rotational lines of NH$_2$D. We present 30m IRAM observations of the above mentioned lines, as well as APEX o/p-NH$_2$D observations of the 1$_{01}$-0$_{00}$ lines at 333 GHz. The hyperfine spectra are first analyzed with a line list that only includes the hyperfine splitting due to the $^{14}$N nucleus. We find inconsistencies between the line widths of the 1$_{01}$-0$_{00}$ and 1$_{11}$-1$_{01}$ lines, the latter being larger by a factor of $\sim$1.6$\pm0.3$. Such a large difference is unexpected given the two sets of lines are likely to originate from the same region. We next employ a newly computed line list for the o/p-NH$_2$D transitions, where the hyperfine structure induced by both nitrogen and deuterium nuclei is included. With this new line list, the analysis of the previous spectra leads to linewidths which are compatible. Neglecting the hyperfine structure owing to D leads to overestimate the linewidths of the o/p-NH$_2$D lines at 3 mm. The error for a cold molecular core is about 50\%. This error propagates directly to the column density estimate. It is therefore recommended to take into account the hyperfine splittings caused by both the $^{14}$N and D nuclei in any analysis relying on these lines.

We have used the Herschel-HIFI instrument to observe both nuclear spin symmetries of amidogen (NH2) towards the high-mass star-forming regions W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4) and G34.3+0.1. The aim is to investigate the ratio of nuclear spin types, the ortho-to-para ratio (OPR), of NH2. The excited NH2 transitions are used to construct radiative transfer models of the hot cores and surrounding envelopes in order to investigate the excitation and possible emission of the ground state rotational transitions of ortho-NH2 N_(K_a,K_c} J=1_(1,1) 3/2 - 0_(0,0) 1/2 and para-NH2 2_(1,2) 5/2 - 1_(0,1) 3/2$ used in the OPR calculations. Our best estimate of the average OPR in the envelopes lie above the high temperature limit of three for W49N, specifically 3.5 with formal errors of \pm0.1, but for W31C, W51, and G34.3+0.1 we find lower values of 2.5\pm0.1, 2.7\pm0.1, and 2.3\pm0.1, respectively. Such low values are strictly forbidden in thermodynamical equilibrium since the OPR is expected to increase above three at low temperatures. In the translucent interstellar gas towards W31C, where the excitation effects are low, we find similar values between 2.2\pm0.2 and 2.9\pm0.2. In contrast, we find an OPR of 3.4\pm0.1 in the dense and cold filament connected to W51, and also two lower limits of >4.2 and >5.0 in two other translucent gas components towards W31C and W49N. At low temperatures (T \lesssim 50 K) the OPR of H2 is <10^-1, far lower than the terrestrial laboratory normal value of three. In such a "para-enriched H2" gas, our astrochemical models can reproduce the variations of the observed OPR, both below and above the thermodynamical equilibrium value, by considering nuclear-spin gas-phase chemistry. The models suggest that values below three arise in regions with temperatures >20-25 K, depending on time, and values above three at lower temperatures.

The deuterium fractionation of gas-phase molecules in hot cores is believed to reflect the composition of interstellar ices. The deuteration of methanol is a major puzzle, however, because the isotopologue ratio [CH2DOH]/[CH3OD], which is predicted to be equal to 3 by standard grain chemistry models, is much larger (~20) in low-mass hot corinos and significantly lower (~1) in high-mass hot cores. This dichotomy in methanol deuteration between low-mass and massive protostars is currently not understood. In this study, we report a simplified rate equation model of the deuterium chemistry occurring in the icy mantles of interstellar grains. We apply this model to the chemistry of hot corinos and hot cores, with IRAS 16293-2422 and the Orion~KL Compact Ridge as prototypes, respectively. The chemistry is based on a statistical initial deuteration at low temperature followed by a warm-up phase during which thermal hydrogen/deuterium (H/D) exchanges occur between water and methanol. The exchange kinetics is incorporated using laboratory data. The [CH2DOH]/[CH3OD] ratio is found to scale inversely with the D/H ratio of water, owing to the H/D exchange equilibrium between the hydroxyl (-OH) functional groups of methanol and water. Our model is able to reproduce the observed [CH2DOH]/[CH3OD] ratios provided that the primitive fractionation of water ice [HDO]/[H2O] is ~ 2% in IRAS 16293-2422 and ~0.6% in Orion~KL. We conclude that the molecular D/H ratios measured in hot cores may not be representative of the original mantles because molecules with exchangeable deuterium atoms can equilibrate with water ice during the warm-up phase.

{\bf R}-matrix calculations combined with the adiabatic nuclei approximation are used to compute electron-impact rotiational excitation rates for three closed-shell diatomic cations, HeH$^+$, CH$^+$, ArH$^+$. Comparisons with previous studies show that an improved treatment of threshold effects leads to significant changes in the low temperature rates, furthermore the new calculations suggest that excitation of CH$^+$ is dominated by $\Delta J =1$ transitions as is expected for cations with a large dipole moment. A model for ArH$^+$ excitation in the Crab Nebula is presented which gives results consistent with the observations for electron densities in the range $2-3\times 10^3$~cm$^{-3}$.

Context: The solution of the nonlocal thermodynamical equilibrium (non-LTE) radiative transfer equation usually relies on stationary iterative methods, which may falsely converge in some cases. Furthermore, these methods are often unable to handle large-scale systems, such as molecular spectra emerging from, for example, cool stellar atmospheres. Aims: Our objective is to develop a new method, which aims to circumvent these problems, using nonstationary numerical techniques and taking advantage of parallel computers. Methods: The technique we develop may be seen as a generalization of the coupled escape probability method. It solves the statistical equilibrium equations in all layers of a discretized model simultaneously. The numerical scheme adopted is based on the generalized minimum residual method. Result:. The code has already been applied to the special case of the water spectrum in a red supergiant stellar atmosphere. This demonstrates the fast convergence of this method, and opens the way to a wide variety of astrophysical problems.

C$_2$H is one of the first radicals to be detected in the interstellar medium. Its higher rotational transitions have recently become available with the Herschel Space Observatory. We aim to constrain the physical parameters of the C$_2$H emitting gas toward the Orion Bar. We analyse the C$_2$H line intensities measured toward the Orion Bar CO$^+$ Peak and Herschel/HIFI maps of C$_2$H, CH, and HCO$^+$, and a NANTEN map of [CI]. We interpret the observed C$_2$H emission using radiative transfer and PDR models. Five rotational transitions of C$_2$H have been detected in the HIFI frequency range toward the CO$^+$ peak. A single component rotational diagram gives a rotation temperature of ~64 K and a beam-averaged C$_2$H column density of 4$\times$10$^{13}$ cm$^{-2}$. The measured transitions cannot be explained by any single parameter model. According to a non-LTE model, most of the C$_2$H column density produces the lower-$N$ C$_2$H transitions and traces a warm ($T_{\rm{kin}}$ ~ 100-150 K) and dense ($n$(H$_2$)~10$^5$-10$^6$ cm$^{-3}$) gas. A small fraction of the C$_2$H column density is required to reproduce the intensity of the highest-$N$ transitions ($N$=9-8 and N=10-9) originating from a high density ($n$(H$_2$)~5$\times$10$^6$ cm$^{-3}$) hot ($T_{\rm{kin}}$ ~ 400 K) gas. The total beam-averaged C$_2$H column density in the model is 10$^{14}$ cm$^{-2}$. Both the non-LTE radiative transfer model and a simple PDR model representing the Orion Bar with a plane-parallel slab of gas and dust suggest, that C$_2$H cannot be described by a single pressure component, unlike the reactive ion CH$^+$, which was previously analysed toward the Orion Bar CO$^+$ peak. The physical parameters traced by the higher rotational transitions ($N$=6-5,...,10-9) of C$_2$H may be consistent with the edges of dense clumps exposed to UV radiation near the ionization front of the Orion Bar.

Chemical models used to study the chemical composition of the gas and the ices in the interstellar medium are based on a network of chemical reactions and associated rate coefficients. These reactions and rate coefficients are partially compiled from data in the literature, when available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public gas-phase network first released in 2012. In addition to a description of the many specific updates, we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as compared with the results of the previous version of our network, kida.uva.2011.

We present quantum dynamical calculations that describe the rotational excitation of H$_2$O due to collisions with H atoms. We used a recent, high accuracy potential energy surface, and solved the collisional dynamics with the close-coupling formalism, for total energies up to 12 000 cm$^{-1}$. From these calculations, we obtained collisional rate coefficients for the first 45 energy levels of both ortho- and para-H$_2$O and for temperatures in the range T = 5-1500 K. These rate coefficients are subsequently compared to the values previously published for the H$_2$O / He and H$_2$O / H$_2$ collisional systems. It is shown that no simple relation exists between the three systems and that specific calculations are thus mandatory.

The availability of collisional rate coefficients with H$_2$ is a pre-requisite for interpretation of observations of molecules whose energy levels are populated under non local thermodynamical equilibrium conditions. In the current study, we present collisional rate coefficients for the NH$_2$D / para--H$_2$($J_2 = 0,2$) collisional system, for energy levels up to $J_\tau = 7_7$ ($E_u$$\sim$735 K) and for gas temperatures in the range $T = 5-300$K. The cross sections are obtained using the essentially exact close--coupling (CC) formalism at low energy and at the highest energies, we used the coupled--states (CS) approximation. For the energy levels up to $J_\tau = 4_2$ ($E_u$$\sim$215 K), the cross sections obtained through the CS formalism are scaled according to a few CC reference points. These reference points are subsequently used to estimate the accuracy of the rate coefficients for higher levels, which is mainly limited by the use of the CS formalism. Considering the current potential energy surface, the rate coefficients are thus expected to be accurate to within 5\% for the levels below $J_\tau = 4_2$, while we estimate an accuracy of 30\% for higher levels.

A non-LTE radiative transfer treatment of cis-methyl formate (HCOOCH3) rotational lines is presented for the first time using a set of theoretical collisional rate coefficients. These coefficients have been computed in the temperature range 5-30 K by combining coupled-channel scattering calculations with a high accuracy potential energy surface for HCOOCH3-He. The results are compared to observations toward the Sagittarius B2(N) molecular cloud using the publicly available PRIMOS survey from the Green Bank Telescope. A total of 49 low-lying transitions of methyl formate, with upper levels below 25 K, are identified. These lines are found to probe a presumably cold (~30 K), moderately dense (~1e4 cm-3) and extended region surrounding Sgr B2(N). The derived column density of ~4e14 cm-2 is only a factor of ~10 larger than the column density of the trans conformer in the same source. Provided that the two conformers have the same spatial distribution, this result suggests that strongly non-equilibrium processes must be involved in their synthesis. Finally, our calculations show that all detected emission lines with a frequency below 30 GHz are (collisionally pumped) weak masers amplifying the continuum of Sgr B2(N). This result demonstrates the importance and generality of non-LTE effects in the rotational spectra of complex organic molecules at centimetre wavelengths.

The aim of the present work is to perform a comprehensive analysis of the interstellar chemistry of nitrogen, focussing on the gas-phase formation of the smallest polyatomic species and in particular nitrogen hydrides. We present a new chemical network in which the kinetic rates of critical reactions have been updated based on recent experimental and theoretical studies, including nuclear spin branching ratios. Our network thus treats the different spin symmetries of the nitrogen hydrides self-consistently together with the ortho and para forms of molecular hydrogen. This new network is used to model the time evolution of the chemical abundances in dark cloud conditions. The steady-state results are analysed, with special emphasis on the influence of the overall amounts of carbon, oxygen, and sulphur. Our calculations are also compared with Herschel/HIFI observations of NH, NH$_2$, and NH$_3$ detected towards the external envelope of the protostar IRAS 16293-2422. The observed abundances and abundance ratios are reproduced for a C/O gas-phase elemental abundance ratio of $\sim0.8$, provided that the sulphur abundance is depleted by a factor larger than 2. The ortho-to-para ratio of H$_2$ in these models is $\sim10^{-3}$. Our models also provide predictions for the ortho-to-para ratios of NH$_2$ and NH$_3$ of $\sim2.3$ and $\sim0.7$ respectively. We conclude that the abundances of nitrogen hydrides in dark cloud conditions are consistent with the gas-phase synthesis predicted with our new chemical network.

We report the first detection of a Galactic source of OH+ line emission: the Orion Bar, a bright nearby photon-dominated region. Line emission is detected over ~1' (0.12 pc), tracing the Bar itself as well as the Southern tip of the Orion Ridge. The line width of ~4 km/s suggests an origin of the OH+ emission close to the PDR surface, at a depth of A_V ~0.3-0.5 into the cloud where most hydrogen is in atomic form. Steady-state collisional and radiative excitation models require unrealistically high OH+ column densities to match the observed line intensity, indicating that the formation of OH+ in the Bar is rapid enough to influence its excitation. Our best-fit OH+ column density of ~1x10^14 cm^-2 is similar to that in previous absorption line studies, while our limits on the ratios of OH+/H2O+ (>~40) and OH+/H3O+ (>~15) are higher than seen before. The column density of OH+ is consistent with estimates from a thermo-chemical model for parameters applicable to the Orion Bar, given the current uncertainties in the local gas pressure and the spectral shape of the ionizing radiation field. The unusually high OH+/H2O+ and OH+/H3O+ ratios are probably due to the high UV radiation field and electron density in this object. In the Bar, photodissociation and electron recombination are more effective destroyers of OH+ than the reaction with H2, which limits the production of H2O+. The appearance of the OH+ lines in emission is the result of the high density of electrons and H atoms in the Orion Bar, since for these species, inelastic collisions with OH+ are faster than reactive ones. In addition, chemical pumping, far-infrared pumping by local dust, and near-UV pumping by Trapezium starlight contribute to the OH+ excitation. Similar conditions may apply to extragalactic nuclei where OH+ lines are seen in emission.

Aims. We investigate the dominant formation mechanism of H2CO and CH3OH in the Horsehead PDR and its associated dense core. Methods. We performed deep integrations of several H2CO and CH3OH lines at two positions in the Horsehead, namely the PDR and dense core, with the IRAM-30m telescope. In addition, we observed one H2CO higher frequency line with the CSO telescope at both positions. We determine the H2CO and CH3OH column densities and abundances from the single-dish observations complemented with IRAM-PdBI high-angular resolution maps (6") of both species. We compare the observed abundances with PDR models including either pure gas-phase chemistry or both gas-phase and grain surface chemistry. Results. We derive CH3OH abundances relative to total number of hydrogen atoms of ~1.2e-10 and ~2.3e-10 in the PDR and dense core positions, respectively. These abundances are similar to the inferred H2CO abundance in both positions (~2e-10). We find an abundance ratio H2CO/CH3OH of ~2 in the PDR and ~1 in the dense core. Pure gas-phase models cannot reproduce the observed abundances of either H2CO or CH3OH at the PDR position. Both species are therefore formed on the surface of dust grains and are subsequently photodesorbed into the gas-phase at this position. At the dense core, on the other hand, photodesorption of ices is needed to explain the observed abundance of CH3OH, while a pure gas-phase model can reproduce the observed H2CO abundance. The high-resolution observations show that CH3OH is depleted onto grains at the dense core. CH3OH is thus present in an envelope around this position, while H2CO is present in both the envelope and the dense core itself. Conclusions. Photodesorption is an efficient mechanism to release complex molecules in low FUV-illuminated PDRs, where thermal desorption of ice mantles is ineffective.

Understanding the origin of the composition of solar system cosmomaterials is a central question, not only in the cosmochemistry and astrochemistry fields, and requires various approaches to be combined. Measurements of isotopic ratios in cometary materials provide strong constraints on the content of the protosolar nebula. Their relation with the composition of the parental dark clouds is, however, still very elusive. In this paper, we bring new constraints based on the isotopic composition of nitrogen in dark clouds, with the aim of understanding the chemical processes that are responsible for the observed isotopic ratios. We have observed and detected the fundamental rotational transition of C$^{15}$N towards two starless dark clouds, L1544 and L1498. We were able to derive the column density ratio of C$^{15}$N over $^{13}$CN towards the same clouds, and obtain the CN/C$^{15}$N isotopic ratios, which were found to be $500\pm75$ for both L1544 and L1498. These values are therefore marginally consistent with the protosolar value of 441. Moreover, this ratio is larger than the isotopic ratio of nitrogen measured in HCN. In addition, we present model calculations of the chemical fractionation of nitrogen in dark clouds, which make it possible to understand how CN can be deprived of $^{15}$N and HCN can simultaneously be enriched in heavy nitrogen. The non-fractionation of N2H+, however, remains an open issue and we propose some chemical way of alleviating the discrepancy between model predictions and the observed ratios.

Utilising previous work by the authors on the spin-coupled rotational cross-sections for electron-CN collisions, data for the associated rate coefficients is presented. Data on rotational, fine-structure and hyperfine-structure transition involving rotational levels up to $N$=20 are computed for temperatures in the range 10 -- 1000~K. Rates are calculated by combining Born-corrected R-matrix calculations with the infinite-order-sudden (IOS) approximation. The dominant hyperfine transitions are those with $\Delta N=\Delta j= \Delta F=1$. For dipole-allowed transitions, electron-impact rates are shown to exceed those for excitation of CN by para-H$_2$($j=0$) by five orders of magnitude. The role of electron collisions in the excitation of CN in diffuse clouds, where local excitation competes with the cosmic microwave background (CMB) photons, is considered. Radiative transfer calculations are performed and the results compared to observations. These comparisons suggest that electron density lies in the range $n(e)\sim 0.01-0.06$~cm$^{-3}$ for typical physical conditions present in diffuse clouds.

Complex (iso-)nitrile molecules, such as CH3CN and HC3N, are relatively easily detected in our Galaxy and in other galaxies. We constrain their chemistry through observations of two positions in the Horsehead edge: the photo-dissociation region (PDR) and the dense, cold, and UV-shielded core just behind it. We systematically searched for lines of CH3CN, HC3N, C3N, and some of their isomers in our sensitive unbiased line survey at 3, 2, and 1mm. We derived column densities and abundances through Bayesian analysis using a large velocity gradient radiative transfer model. We report the first clear detection of CH3NC at millimeter wavelength. We detected 17 lines of CH3CN at the PDR and 6 at the dense core position, and we resolved its hyperfine structure for 3 lines. We detected 4 lines of HC3N, and C3N is clearly detected at the PDR position. We computed new electron collisional rate coefficients for CH3CN, and we found that including electron excitation reduces the derived column density by 40% at the PDR position. While CH3CN is 30 times more abundant in the PDR than in the dense core, HC3N has similar abundance at both positions. The isomeric ratio CH3NC/CH3CN is 0.15+-0.02. In the case of CH3CN, pure gas phase chemistry cannot reproduce the amount of CH3CN observed in the UV-illuminated gas. We propose that CH3CN gas phase abundance is enhanced when ice mantles of grains are destroyed through photo-desorption or thermal-evaporation in PDRs, and through sputtering in shocks. (abridged)

We compute the rotational quenching rates of the first 81 rotational levels of ortho- and para-H2CO in collision with ortho- and para-H2, for a temperature range of 10-300 K. We make use of the quantum close-coupling and coupled-states scattering methods combined with the high accuracy potential energy surface of Troscompt et al. (2009a). Rates are significantly different from the scaled rates of H2CO in collision with He; consequently, critical densities are noticeably lower. We compare a full close- coupling computation of pressure broadening cross sections with experimental data and show that our results are compatible with the low temperature measurements of Mengel & De Lucia (2000), for a spin temperature of H2 around 50 K.

The interpretation of water lines in red supergiant stellar atmospheres has been much debated over the past decade. The introduction of the so-called MOLspheres to account for near-infrared "extra" absorption has been controversial. We propose that non-LTE effects should be taken into account before considering any extra-photospheric contribution. After a brief introduction on the radiative transfer treatment and the inadequacy of classical treatments in the case of large-scale systems such as molecules, we present a new code, based on preconditioned Krylov subspace methods. Preliminary results suggest that NLTE effects lead to deeper water bands, as well as extra cooling.

The line intensities of the fundamental rotational transitions of H13CN and HC15N were observed towards two prestellar cores, L183 and L1544, and lead to molecular isotopic ratios 140 6 14N/15N 6 250 and 140 6 14N/15N 6 360, respectively. The range of values reflect genuine spatial variations within the cores. A comprehensive analysis of the available measurements of the nitrogen isotopic ratio in prestellar cores show that molecules carrying the nitrile functional group appear to be systematically 15N-enriched com- pared to those carrying the amine functional group. A chemical origin for the differential 15N-enhance- ment between nitrile- and amine-bearing interstellar molecules is proposed. This sheds new light on several observations of Solar System objects: (i) the similar N isotopic fractionation in Jupiter's NH3 and solar wind N+; (ii) the 15N-enrichments in cometary HCN and CN (that might represent a direct inter- stellar inheritance); and (iii) 15N-enrichments observed in organics in primitive cosmomaterials. The large variations in the isotopic composition of N-bearing molecules in Solar System objects might then simply reflect the different interstellar N reservoirs from which they are originating.

Water is a key molecule in many astrophysical studies. Its high dipole moment makes this molecule to be subthermally populated under the typical conditions of most astrophysical objects. This motivated the calculation of various sets of collisional rate coefficients (CRC) for H$_2$O (with He or H$_2$) which are necessary to model its rotational excitation and line emission. We performed accurate non--local non--LTE radiative transfer calculations using different sets of CRC in order to predict the line intensities from transitions that involve the lowest energy levels of H$_2$O (E $<$ 900 K). The results obtained from the different CRC sets are then compared using line intensity ratio statistics. For the whole range of physical conditions considered in this work, we obtain that the intensities based on the quantum and QCT CRC are in good agreement. However, at relatively low H$_2$ volume density ($n$(H$_2$) $<$ 10$^7$ cm$^{-3}$) and low water abundance ($\chi$(H$_2$O) $<$ 10$^{-6}$), these physical conditions being relevant to describe most molecular clouds, we find differences in the predicted line intensities of up to a factor of $\sim$ 3 for the bulk of the lines. Most of the recent studies interpreting early Herschel Space Observatory spectra used the QCT CRC. Our results show that although the global conclusions from those studies will not be drastically changed, each case has to be considered individually, since depending on the physical conditions, the use of the QCT CRC may lead to a mis--estimate of the water vapour abundance of up to a factor of $\sim$ 3.

The HDO/H2O ratio is a powerful diagnostic to understand the evolution of water from the first stages of star formation to the formation of planets and comets. Our aim is to determine precisely the abundance distribution of HDO towards the low-mass protostar IRAS16293-2422 and learn more about the water formation mechanisms by determining the HDO/H2O abundance ratio. A spectral survey of the source IRAS16293-2422 was carried out in the framework of the CHESS Herschel Key program with the HIFI instrument, allowing the detection of numerous HDO lines. Other transitions have been observed previously with ground-based telescopes. The spherical Monte Carlo radiative transfer code RATRAN was used to reproduce the observed line profiles of HDO by assuming an abundance jump. To determine the H2O abundance throughout the envelope, a similar study was made of the H2-18O observed lines, as the H2O main isotope lines are contaminated by the outflows. We derive an inner HDO abundance of 1.7e-7 and an outer HDO abundance of 8e-11. To reproduce the HDO absorption lines, it is necessary to add an absorbing layer in front of the envelope. It may correspond to a water-rich layer created by the photodesorption of the ices at the edges of the molecular cloud. The HDO/H2O ratio is ~1.4-5.8% in the hot corino whereas it is ~0.2-2.2% in the outer envelope. It is estimated at ~4.8% in the added absorbing layer. Although it is clearly higher than the cosmic D/H abundance, the HDO/H2O ratio remains lower than the D/H ratio derived for other deuterated molecules observed in the same source. The similarity of the ratios derived in the hot corino and in the added absorbing layer suggests that water formed before the gravitational collapse of the protostar, contrary to formaldehyde and methanol, which formed later once the CO molecules had depleted on the grains.

We present a novel chemical database for gas-phase astrochemistry. Named the KInetic Database for Astrochemistry (KIDA), this database consists of gas-phase reactions with rate coefficients and uncertainties that will be vetted to the greatest extent possible. Submissions of measured and calculated rate coefficients are welcome, and will be studied by experts before inclusion into the database. Besides providing kinetic information for the interstellar medium, KIDA is planned to contain such data for planetary atmospheres and for circumstellar envelopes. Each year, a subset of the reactions in the database (kida.uva) will be provided as a network for the simulation of the chemistry of dense interstellar clouds with temperatures between 10 K and 300 K. We also provide a code, named Nahoon, to study the time-dependent gas-phase chemistry of 0D and 1D interstellar sources.

The clumpy density structure of photon-dominated regions is well established, but the physical properties of the clumps and of the surrounding interclump medium are only approximately known. The aim of this paper is to constrain the physical and chemical conditions in the Orion Bar, a prototypical nearby photon-dominated region. We present observations of the HF J=1-0 line, which appears in emission toward the Orion Bar, and compare the brightness of the line to non-LTE radiative transfer calculations. The large width of the HF line suggests an origin of the emission in the interclump gas, but collisional excitation by H2 in the interclump gas underpredicts the observed line intensity by factors of 3-5. In contrast, an origin of the line in the dense clumps requires a density of ~10^9 cm^-3, 10-100 times higher than previous estimates, which is unlikely. However, electron impact excitation reproduces our observations for T = 100 K and n(e) = 10 cm^-3, as expected for the interclump gas. We conclude that HF emission is a signpost of molecular gas with a high electron density. Similar conditions may apply to active galactic nuclei where HF also appears in emission.

Nitrogen bearing species are common tracers of the physical conditions in a wide variety of objects, and most remarkably in dark clouds. The reservoir of gaseous nitrogen is expected to be atomic or molecular, but none of the two species are observable in the dark gas. Their abundances therefore derive indirectly from those of N-bearing species through chemical modelling. The recent years have accumulated data which stress our incomplete understanding of the nitrogen chemistry in dark cloud conditions. To tackle this problem of the nitrogen chemistry in cold gas, we have revised the formation of nitrogen hydrides, which is initiated by the key reaction \ce{N+ + H2 -> NH+ + H}. We propose a new rate for this reaction which depends on the ortho-to-para ratio of H2. This new rate allows to reproduce the abundance ratios of the three nitrogen hydrides, NH, \ce{NH2}, and \ce{NH3}, observed towards IRAS16293-2422, provided that the channel leading to NH from the dissociative recombination of \ce{N2H+} is not closed at low temperature. The ortho-to-para ratio of H2 is constrained to O/P=$10^{-3}$ by the abundance ratio NH:NH2, which provides a new method to measure O/P. This work stresses the need for reaction rates at the low temperatures of dark clouds, and for branching ratios of critical dissociative recombination reactions.

We present the first detection of hydrogen chlorine in a protostellar shock, by observing the fundamental transition at 626 GHz with the Herschel HIFI spectrometer. We detected two of the three hyperfine lines, from which we derived a line opacity < 1. Using a non-LTE LVG code, we constrained the HCl column density, temperature and density of the emitting gas. The hypothesis that the emission originates in the molecular cloud is ruled out, as it would imply a too dense gas. Conversely, assuming that the emission originates in the 10"-15" size shocked gas previously observed at the IRAM PdB interferometer, we obtain: N(HCl)=0.7-2 x 10(13) cm-2, temperature > 15 K and density > 3 x 10(5) cm-3}. Combining with the Herschel HIFI CO(5-4) observations allows to further constrain the gas density and temperature, 10(5)-10(6) cm-3 and 120-250 K, as well as the HCl column density, 2 x 10(13) cm-2, and, finally, abundance: 3-6 x 10(-9). The estimated HCl abundance is consistent with that previously observed in low- and high- mass protostars. This puzzling result in the L1157-B1 shock, where species from volatile and refractory grains components are enhanced, suggests either that HCl is not the main reservoir of chlorine in the gas phase, against previous chemical models predictions, or that the elemental chlorine abundance is low in L1157-B1. Astrochemical modelling suggests that HCl is in fact formed in the gas phase, at low temperatures, prior to the occurance of the shock, and that the latter does not enhance its abundance.

Collisional de-excitation rates of partially deuterated molecules are different from the fully hydrogenated species because of lowering of symmetry. We compute the collisional (de)excitation rates of ND2H by ground state para-H2, extending the previous results for He- lium. We describe the changes in the potential energy surface of NH3- H2 involved by the pres- ence of two deuterium nuclei. Cross sections are calculated within the full close-coupling ap- proach and augmented with coupled-state calculations. Collisional rate coefficients are given between 5 and 35 K, a range of temperatures which is relevant to cold interstellar conditions. We find that the collisional rates of ND2H by H2 are about one order of magnitude higher than those obtained with Helium as perturber. These results are essential to radiative transfer modelling and will allow to interpret the millimeter and submillimeter detections of ND2H with better constraints than previously.

Nitrogen is the fifth most abundant element in the Universe, yet the gas-phase chemistry of N-bearing species remains poorly understood. Nitrogen hydrides are key molecules of nitrogen chemistry. Their abundance ratios place strong constraints on the production pathways and reaction rates of nitrogen-bearing molecules. We observed the class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI, covering most of the frequency range from 0.48 to 1.78~THz at high spectral resolution. The hyperfine structure of the amidogen radical o-NH2 is resolved and seen in absorption against the continuum of the protostar. Several transitions of ammonia from 1.2 to 1.8~THz are also seen in absorption. These lines trace the low-density envelope of the protostar. Column densities and abundances are estimated for each hydride. We find that NH:NH2:NH3=5:1:300. {Dark clouds chemical models predict steady-state abundances of NH2 and NH3 in reasonable agreement with the present observations, whilst that of NH is underpredicted by more than one order of magnitude, even using updated kinetic rates. Additional modelling of the nitrogen gas-phase chemistry in dark-cloud conditions is necessary before having recourse to heterogen processes.

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

Shock modelling predicts an electron density enhancement within the magnetic precursor of C-shocks. Previous observations of SiO, H13CO+, HN13C and H13CN toward the young L1448-mm outflow showed an over-excitation of the ion fluid that was attributed to an electron density enhancement in the precursor. We re-visit this interpretation and test if it still holds when we consider different source morphologies and kinetic temperatures for the observed molecules, and also give some insight on the spatial extent of the electron density enhancement around L1448-mm. We estimate the opacities of H13CO+ and HN13C by observing the J=3\to2 lines of rarer isotopologues to confirm that the emission is optically thin. To model the excitation of the molecules, we use the large velocity gradient (LVG) approximation with updated collisional coefficients to i) re- analyse the observations toward the positions where the over-excitation of H13CO+ has previously been observed [i.e. toward L1448- mm at offsets (0,0) and (0,-10)], and ii) to investigate if the electron density enhancement is still required for the cases of extended and compact emission, and for kinetic temperatures of up to 400 K. We also report several lines of SiO, HN13C and H13CO+ toward new positions around this outflow, to investigate the spatial extent of the over-excitation of the ions in L1448-mm. From the isotopologue observations, we find that the emission of H13CO+ and HN13C from the precursor is optically thin if this emission is extended. Using the new collisional coefficients, an electron density enhancement is still needed to explain the excitation of H13CO+ for extended emission and for gas temperatures of\le 400 K toward L1448-mm (0,-10), and possibly also toward L1448-mm (0,0). For compact emission the data cannot be fitted. We do not find any evidence for the over-excitation of the ion fluid toward the newly observed positions around L1448-mm. The observed line emission of SiO, H13CO+ and HN13C toward L1448-mm (0,0) and (0,-10) is consistent with an electron density enhancement in the precursor component, if this emission is spatially extended. This is also true for the case of high gas temperatures (\le400 K) toward the (0,-10) offset. The electron density enhancement seems to be restricted to the southern, redshifted lobe of the L1448-mm outflow. Interferometric images of the line emission of these molecules are needed to confirm the spatial extent of the over-excitation of the ions and thus, of the electron density enhancement in the magnetic precursor of L1448-mm.

The internal energy available in vibrationally excited H2 molecules can be used to overcome or diminish the activation barrier of various chemical reactions of interest for molecular astrophysics. In this article we investigate in detail the impact on the chemical composition of interstellar clouds of the reactions of vibrationally excited H2 with C+, He+, O, OH, and CN, based on the available chemical kinetics data. It is found that the reaction of H2 (v>0) and C+ has a profound impact on the abundances of some molecules, especially CH+, which is a direct product and is readily formed in astronomical regions with fractional abundances of vibrationally excited H2, relative to ground state H2, in excess of 10^(-6), independently of whether the gas is hot or not. The effects of these reactions on the chemical composition of the diffuse clouds zeta Oph and HD 34078, the dense PDR Orion Bar, the planetary nebula NGC 7027, and the circumstellar disk around the B9 star HD 176386 are investigated through PDR models. We find that formation of CH+ is especially favored in dense and highly FUV illuminated regions such as the Orion Bar and the planetary nebula NGC 7027, where column densities in excess of 10^(13) cm^(-2) are predicted. In diffuse clouds, however, this mechanism is found to be not efficient enough to form CH+ with a column density close to the values derived from astronomical observations.

Modeling of molecular emission from interstellar clouds requires the calculation of rates for excitation by collisions with the most abundant species. The present paper focuses on the calculation of rate coefficients for rotational excitation of the HCN and HNC molecules in their ground vibrational state in collision with He. The calculations are based on new two-dimensional potential energy surfaces obtained from highly correlated \textit{ab initio} calculations. Calculations of pure rotational (de)excitation cross sections of HCN and HNC by He were performed using the essentially exact close-coupling method. Cross sections for transitions among the 8 first rotational levels of HCN and HNC were calculated for kinetic energies up to 1000 cm$^{-1}$. These cross sections were used to determine collisional rate constants for temperatures ranging from 5 K to 100 K. A propensity for even $\Delta j$ transitions is observed in the case of HCN--He collisions whereas a propensity for odd $\Delta j$ transitions is observed in the case of HNC--He collisions. The consequences for astrophysical models are evaluated and it is shown that the use of HCN rate coefficients to interpret HNC observations can lead to significant inaccuracies in the determination of the HNC abundance, in particular in cold dark clouds for which the new HNC rates show that the $j=1-0$ line of this species will be more easily excited by collisions than HCN. An important result of the new HNC-He rates is that the HNC/HCN abundance ratio derived from observations in cold clouds has to be revised from $>$1 to $\simeq$1, in good agreement with detailed chemical models available in the literature.

Aims. Intermediate Mass (IM) stars are an important component of our Galaxy, as they significantly contribute to the interstellar FUV field and, consequently, play an important role in the energy balance of the ISM. Despite their importance, very little is known about their formation process and only a few studies have been devoted to characterize the first phases in the evolution of intermediate mass protostars. Here we consider in great detail the case of the brightest and closest known young IM protostar: FIR4 in the OMC2 component of the Orion molecular cloud complex. Methods. We analyzed the available continuum emission (maps and SED) through one-dimensional dust radiative transfer calculations. We ran large grids of models to find the envelope model that best fits the data. The derived dust density and temperature profiles have been then used to compute the gas temperature profile, equating gas cooling and heating terms across the envelope. Last, we computed the water line spectrum for various possible values of water abundance. Results. The luminosity of FIR4 has been reevaluated to 1000 Lo, making FIR4 definitively an Intermediate Mass protostar. The envelope surrounding FIR4 has a relatively shallow density power law index, ~ 0.6. The most surprising result is that the gas and dust are thermally decoupled in the interior of the envelope, where the dust ices sublimate at 100 K. This has important consequences in the interpretation of the line data. We provide the predictions for the water spectrum, and discuss in detail the lines which will be observed by the Herschel Space Observatory .

Ammonia inversion lines are often used as probes of the physical conditions in the dense ISM. The excitation temperature between the first two para metastable (rotational) levels is an excellent probe of the gas kinetic temperature. However, the calibration of this ammonia thermometer depends on the accuracy of the collisional rates with H2. Here we present new collisional rates for ortho-NH3 and para-NH3 colliding with para-H2 (J=0) and we investigate the effects of these new rates on the excitation of ammonia. Scattering calculations employ a new, high accuracy, potential energy surface computed at the coupled-cluster CCSD(T) level with a basis set extrapolation procedure. Rates are obtained for all transitions involving ammonia levels with J <= 3 and for kinetic temperatures in the range 5-100 K. We find that the calibration curve of the ammonia thermometer -- which relates the observed excitation temperature between the first two para metastable levels to the gas kinetic temperature -- does not change significantly when these new rates are used. Thus, the calibration of ammonia thermometer appears to be robust. Effects of the new rates on the excitation temperature of inversion and rotation-inversion transitions are also found to be small.

The molecular gas composition in the inner 1 kpc disk of the starburst galaxy M82 resembles that of Galactic Photon Dominated Regions (PDRs). In particular, large abundances of the reactive ions HOC+ and CO+ have been measured in the nucleus of this galaxy. To investigate the origin of the large abundances of reactive ions in M82, we have completed our previous 30m HOC+ J=1-0 observations with the higher excitation HCO+ and HOC+ J=4-3 and 3-2 rotational lines. In addition, we have obtained with the IRAM Plateau de Bure Interferometer (PdBI) a 4" resolution map of the HOC+ 1-0 emission, the first ever obtained in a Galactic or extragalactic source. Our HOC+ interferometric image shows that the emission of the HOC+ 1-0 line is mainly restricted to the nuclear disk, with the maxima towards the E. and W. molecular peaks. In addition, line excitation calculations imply that the HOC+ emission arises in dense gas. Therefore, the HOC+ emission is arising in the dense PDRs embedded in the M82 nuclear disk, rather than in the intercloud phase and/or wind. We have improved our previous chemical model of M82 by (i) using the new version of the Meudon PDR code, (ii) updating the chemical network, and (iii) considering two different types of clouds (with different thickness). Most molecular observations (HCO+, HOC+, CO+, CN, HCN, H3O+) are well explained assuming that ~ 87 % of the mass of the molecular gas is forming small clouds (Av=5 mag) while only ~ 13 % of the mass is in large molecular clouds (Av=50 mag). Such small number of large molecular clouds suggests that M82 is an old starburst, where star formation has almost exhausted the molecular gas reservoir.

Rotational excitation of isotopologues of HCN and HNC by thermal electron-impact is studied using the molecular {\bf R}-matrix method combined with the adiabatic-nuclei-rotation (ANR) approximation. Rate coefficients are obtained for electron temperatures in the range 5$-$6000 K and for transitions among all levels up to J=8. Hyperfine rates are also derived using the infinite-order-sudden (IOS) scaling method. It is shown that the dominant rotational transitions are dipole allowed, that is those for which $\Delta J=1$. The hyperfine propensity rule $\Delta J=\Delta F$ is found to be stronger than in the case of He$-$HCN collisions. For dipole allowed transitions, electron-impact rates are shown to exceed those for excitation of HCN by He atoms by 6 orders of magnitude. As a result, the present rates should be included in any detailed population model of isotopologues of HCN and HNC in sources where the electron fraction is larger than 10$^{-6}$, for example in interstellar shocks and comets.

The interpretation of water line emission from existing observations and future HIFI/Herschel data requires a detailed knowledge of collisional rate coefficients. Among all relevant collisional mechanisms, the rotational (de)excitation of H2O by H2 molecules is the process of most interest in interstellar space. To determine rate coefficients for rotational de-excitation among the lowest 45 para and 45 ortho rotational levels of H2O colliding with both para and ortho-H2 in the temperature range 20-2000 K. Rate coefficients are calculated on a recent high-accuracy H2O-H2 potential energy surface using quasi-classical trajectory calculations. Trajectories are sampled by a canonical Monte-Carlo procedure. H2 molecules are assumed to be rotationally thermalized at the kinetic temperature. By comparison with quantum calculations available for low lying levels, classical rates are found to be accurate within a factor of 1-3 for the dominant transitions, that is those with rates larger than a few 10^{-12}cm^{3}s^{-1}. Large velocity gradient modelling shows that the new rates have a significant impact on emission line fluxes and that they should be adopted in any detailed population model of water in warm and hot environments.

Rates for rotational excitation of HC3N by collisions with He atoms and H2 molecules are computed for kinetic temperatures in the range 5-20K and 5-100K, respectively. These rates are obtained from extensive quantum and quasi-classical calculations using new accurate potential energy surfaces (PES).

Magnetic precursors of C-shocks accelerate, compress and heat molecular ions, modifying the kinematics and the physical conditions of the ion fluid with respect to the neutral one. Electron densities are also expected to be significantly enhanced in shock precursors. In this Letter, we present observations of strongly polar ion and neutral molecules such as SiO, H13CO+, HN13C and H13CN, which reveal the electron density enhancements associated with the precursor of the young L1448-mm outflow. While in the ambient gas the excitation of the ions and neutrals is explained by collisional excitation by H2 with a single density of 10E5cmE-3, H13CO+ shows an over excitation in the shock precursor component that requires H2 densities of a factor of >10 larger than those derived from the neutral species. This over excitation in H13CO+ can be explained if we consider an additional excitation by collisions with electrons and an electron density enhancement in the precursor stage by a factor of 500, i.e. a fractional ionization of 5x10E-5. These results show that multiline observations can be used to study the evolution of the ion and electron fluids at the first stages of the C-shock interaction.

Cross sections for the rotational (de)excitation of CO by ground state para- and ortho-H_2 are obtained using quantum scattering calculations for collision energies between 1 and 520 cm^{-1}. A new CO-H_2 potential energy surface is employed and its quality is assessed by comparison with explicitly correlated CCSD(T)-R12 calculations. Rate constants for rotational levels of CO up to 5 and temperatures in the range 5-70 K are deduced. The new potential is found to have a strong influence on the resonance structure of the cross sections at very low collision energies. As a result, the present rates at 10 K differ by up to 50% with those obtained by \citet{flower01} on a previous, less accurate, potential energy surface.

We present observations of methanol lines in a sample of Class 0 low mass protostars. Using a 1-D radiative transfer model, we derive the abundances in the envelopes. In two sources of the sample, the observations can only be reproduced by the model if the methanol abundance is enhanced by about two order of magnitude in the inner hot region of the envelope. Two other sources show similar jumps, although at a lower confidence level. The observations for the other three sources are well reproduced with a constant abundance, but the presence of a jump cannot be ruled out. The observed methanol abundances in the warm gas around low mass protostars are orders of magnitude higher than gas phase chemistry models predict. Hence, in agreement with other evidences, this suggest that the high methanol abundance reflects recent evaporation of ices due to the heating by the newly formed star. The observed abundance ratios of CH3 OH, H2 CO, and CO are in good agreement with grain surface chemistry models. However, the absolute abundances are more difficult to reproduce and may point towards the presence of multiple ice components in these regions.

We present SO and SO2 observations in the region of the low mass protostar IRAS2/NGC1333. The East-West outflow originating from this source has been mapped in four transitions of both SO and SO2. In addition, CS observations published in the literature have been used. We compute the SO, SO2 and CS column densities and the physical conditions at several positions of the outflow using LTE and a non-LTE LVG approximations. The SO2/SO and CS/SO abundance ratios are compared with the theoretical predictions of a chemical model adapted to the physical conditions in the IRAS2 outflow. The SO2/SO abundance ratios are constant in the two lobes of the outflow whereas CS/SO is up to 6 times lower in the shocked gas of the East lobe than in the West one. The comparison with the chemical model allows us to constrain the age of the outflow produced by IRAS2 to >5e3yr. We find low densities and temperatures for the outflow of IRAS2 (< 1e6 cm-3 and >70 K) from SO and SO2 emission probably because the two molecules trace the cooled entrained material. The East lobe of the outflow shows denser gas compared to the West lobe. We also discuss some constraints on the depleted form of sulphur.