Search

Context. Molecular hydrogen, H2, is the most abundant molecule in the Universe. Thanks to its widely spaced energy levels, it predominantly lights up in warm gas, T ≳ 102 K, such as shocked regions externally irradiated or not by interstellar UV photons, and it is one of the prime targets of James Webb Space Telescope (JWST) observations. These may include shocks from protostellar outflows, supernova remnants impinging on molecular clouds, all the way up to starburst galaxies and active galactic nuclei. Aims. Sophisticated shock models are able to simulate H2 emission from such shocked regions. We aim to explore H2 excitation using shock models, and to test over which parameter space distinct signatures are produced in H2 emission. Methods. We here present simulated H2 emission using the Paris-Durham shock code over an extensive grid of ~14 000 plane-parallel stationary shock models, a large subset of which are exposed to a semi-isotropic external UV radiation field. The grid samples six input parameters: the preshock density, shock velocity, transverse magnetic field strength, UV radiation field strength, the cosmic-ray-ionization rate, and the abundance of polycyclic aromatic hydrocarbons, PAHs. Physical quantities resulting from our self-consistent calculations, such as temperature, density, and width, have been extracted along with H2 integrated line intensities. These simulations and results are publicly available on the Interstellar Medium Services platform. Results. The strength of the transverse magnetic field, as quantified by the magnetic scaling factor, b, plays a key role in the excitation of H2. At low values of b (≲0.3, J-type shocks), H2 excitation is dominated by vibrationally excited lines; whereas, at higher values (b ≳ 1, C-type shocks), rotational lines dominate the spectrum for shocks with an external radiation field comparable to (or lower than) the solar neighborhood. Shocks with b ≫ 1 can potentially be spatially resolved with JWST for nea

(Abridged) H2 is the most abundant molecule in the Universe. Thanks to its widely spaced energy levels, it predominantly lights up in warm gas, T > 100 K, such as shocked regions, and it is one of the key targets of JWST observations. These include shocks from protostellar outflows, all the way up to starburst galaxies and AGN. Shock models are able to simulate H2 emission. We aim to explore H2 excitation using such models, and to test over which parameter space distinct signatures are produced in H2 emission. We present simulated H2 emission using the Paris-Durham shock code over an extensive grid of 14,000 plane-parallel stationary shock models, a large subset of which are exposed to an external UV radiation field. The grid samples 6 input parameters: preshock density, shock velocity, transverse magnetic field strength, UV radiation field strength, cosmic-ray-ionization rate, and PAH abundance. Physical quantities, such as temperature, density, and width, have been extracted along with H2 integrated line intensities. The strength of the transverse magnetic field, set by the scaling factor, b, plays a key role in the excitation of H2. At low values of b (<~ 0.3, J-type shocks), H2 excitation is dominated by vibrationally excited lines; at higher values (b >~ 1, C-type shocks), rotational lines dominate the spectrum for shocks with an external radiation field comparable to (or lower than) the solar neighborhood. Shocks with b >= 1 can be spatially resolved with JWST for nearby objects. When the input kinetic energy flux increases, the excitation and integrated intensity of H2 increases similarly. An external UV field mainly serves to increase the excitation, particularly for shocks where the input radiation energy is comparable to the input kinetic energy flux. These results provide an overview of the energetic reprocessing of input kinetic energy flux and the resulting H2 line emission., Comment: Published in A&A

Context. Molecular hydrogen, H2, is the most abundant molecule in the Universe. Thanks to its widely spaced energy levels, it predominantly lights up in warm gas, T ≳ 102 K, such as shocked regions externally irradiated or not by interstellar UV photons, and it is one of the prime targets of James Webb Space Telescope (JWST) observations. These may include shocks from protostellar outflows, supernova remnants impinging on molecular clouds, all the way up to starburst galaxies and active galactic nuclei. Aims. Sophisticated shock models are able to simulate H2 emission from such shocked regions. We aim to explore H2 excitation using shock models, and to test over which parameter space distinct signatures are produced in H2 emission. Methods. We here present simulated H2 emission using the Paris-Durham shock code over an extensive grid of ~14 000 plane-parallel stationary shock models, a large subset of which are exposed to a semi-isotropic external UV radiation field. The grid samples six input parameters: the preshock density, shock velocity, transverse magnetic field strength, UV radiation field strength, the cosmic-ray-ionization rate, and the abundance of polycyclic aromatic hydrocarbons, PAHs. Physical quantities resulting from our self-consistent calculations, such as temperature, density, and width, have been extracted along with H2 integrated line intensities. These simulations and results are publicly available on the Interstellar Medium Services platform. Results. The strength of the transverse magnetic field, as quantified by the magnetic scaling factor, b, plays a key role in the excitation of H2. At low values of b (≲0.3, J-type shocks), H2 excitation is dominated by vibrationally excited lines; whereas, at higher values (b ≳ 1, C-type shocks), rotational lines dominate the spectrum for shocks with an external radiation field comparable to (or lower than) the solar neighborhood. Shocks with b ≫ 1 can potentially be spatially resolved with JWST for nea

The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes of different phases of the ISM, as demonstrated by recent studies with the Herschel Space Observatory. The observations performed under the HyGAL program will allow us to address several questions related to the lifecycle of molecular material in the ISM and the physical processes that impact its phase transition, such as: (1) What is the distribution function of the H2 fraction in the ISM? (2) How does the ionization rate due to low-energy cosmic-rays vary within the Galaxy? (3) What is the nature of interstellar turbulence, and what mechanisms lead to its dissipation? This overview discusses the observing strategy, synergies with ancillary and archival observations, the data reduction and analysis schemes adopted; and presents the first results obtained toward three of the survey targets, W3(OH), W3IRS5 and NGC7538IRS1. Robust measurements of the column densities of these hydrides -- obtained through widespread observations of absorption lines-- help address the questions raised, and there is a timely synergy between these observations and the development of theoretical models, particularly pertaining to the formation of H2 within the turbulent ISM. The provision of enhanced HyGAL data products will therefore serve as a legacy for future ISM studies., Comment: 33 pages, 17 figures, 7 tables. Accepted for publication in the Astrophysical Journal

The origin of the cold gas in central galaxies in groups is still a matter of debate. We present Multi-Unit Spectroscopic Explorer (MUSE) observations of 18 optically selected local Brightest Group Galaxies (BGGs) to study the kinematics and distribution of the optical emission-line gas. MUSE observations reveal a distribution of gas morphologies including ten complex networks of filaments extending up to 10 kpc to two compact (<3 kpc) and five extended (>5 kpc) disk-dominated structures. Some rotating disks show rings and elongated structures arising from the central disk. The kinematics of the stellar component is mainly rotation-dominated, which is very different from the disturbed kinematics and distribution found in the filamentary sources. The ionized gas is kinematically decoupled from the stellar component for most systems, suggesting an external origin for the gas. We also find that the Halpha luminosity correlates with the cold molecular mass. By exploring the thermodynamical properties of the hot atmospheres, we find that the filamentary sources and compact disks are found in systems with small central entropy values and tcool/teddy ratios. This suggests that, like for Brightest Cluster Galaxies in cool core clusters, the ionized gas are likely formed from hot halo gas condensations, consistently with the Chaotic Cold Accretion simulations (as shown via the C-ratio, Tat, and k-plot). We note that gaseous rotating disks are more frequent than in BCGs. An explanation for the origin of the gas in those objects is a contribution to gas fueling by mergers or group satellites, as qualitatively hinted by some sources of the present sample. Nonetheless, we discuss the possibility that some extended disks could also be a transition stage in an evolutionary sequence including filaments, extended disks and compact disks, as described by hot gas condensation models of cooling flows., Comment: 38 pages, and 16 figures. Resubmitted to A&A after referee report

The HyGAL Stratospheric Observatory for Infrared Astronomy legacy program surveys six hydride molecules-ArH+, OH+, H2O+, SH, OH, and CH-and two atomic constituents-C+ and O-within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes of different phases of the ISM, as demonstrated by recent studies with the Herschel Space Observatory. The observations performed under the HyGAL program will allow us to address several questions related to the life cycle of molecular material in the ISM and the physical processes that impact the phase transition from atomic to molecular gas, such as: (1) What is the distribution function of the H-2 fraction in the ISM? (2) How does the ionization rate due to low-energy cosmic rays vary within the Galaxy? (3) What is the nature of interstellar turbulence (e.g., typical shear or shock velocities), and what mechanisms lead to its dissipation? In this overview, we discuss the observing strategy, the synergies with ancillary and archival observations of other small molecules, and the data reduction and analysis schemes we adopted; and we present the first results obtained toward three of the survey targets, W3(OH), W3 IRS5, and NGC 7538 IRS1. Robust measurements of the column densities of these hydrides-obtained through widespread observations of absorption lines-help address the questions raised, and there is a very timely synergy between these observations and the development of theoretical models, particularly pertaining to the formation of H-2 within the turbulent ISM. The provision of enhanced HyGAL data products will therefore serve as a legacy for future ISM studies.

As a complement to the HyGAL Stratospheric Observatory for Infrared Astronomy Legacy Program, we report the results of a ground-based absorption line survey of simple molecules in diffuse and translucent Galactic clouds. Using the Institut de Radioastronomie Millim\'etrique (IRAM) 30 m telescope, we surveyed molecular lines in the 2 mm and 3 mm wavelength ranges toward 15 millimeter continuum sources. These sources, which are all massive star-forming regions located mainly in the first and second quadrants of the Milky Way, form the subset of the HyGAL sample that can be observed by the IRAM 30 m telescope. We detected HCO$^+$ absorption lines toward 14 sightlines, toward which we identified 78 foreground cloud components, as well as lines from HCN, HNC, C$_2$H, and c-C$_3$H$_2$ toward most sightlines. In addition, CS and H$_2$S absorption lines are found toward at least half of the continuum sources. Static Meudon photodissociation region (PDR) isobaric models that consider ultraviolet-dominated chemistry were unable to reproduce the column densities of all seven molecular species by just a factor of a few, except for H$_2$S. The inclusion of other formation routes driven by turbulent dissipation could possibly explain the observed high column densities of these species in diffuse clouds. There is a tentative trend for H$_2$S and CS abundances relative to H$_2$ to be larger in diffuse clouds ($X$(H$_2$S) and $X$(CS) $\sim 10^{-8} - 10^{-7}$) than in translucent clouds ($X$(H$_2$S) and $X$(CS) $\sim 10^{-9} - 10^{-8}$) toward a small sample; however, a larger sample is required in order to confirm this trend. The derived H$_2$S column densities are higher than the values predicted from the isobaric PDR models, suggesting that chemical desorption of H$_2$S from sulfur-containing ice mantles may play a role in increasing the H$_2$S abundance., Comment: 41 pages, 25 figures, 10 tables, and 5 appendices. Accepted for publication in A&A

The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes of different phases of the ISM, as demonstrated by recent studies with the Herschel Space Observatory. The observations performed under the HyGAL program will allow us to address several questions related to the lifecycle of molecular material in the ISM and the physical processes that impact its phase transition, such as: (1) What is the distribution function of the H2 fraction in the ISM? (2) How does the ionization rate due to low-energy cosmic-rays vary within the Galaxy? (3) What is the nature of interstellar turbulence, and what mechanisms lead to its dissipation? This overview discusses the observing strategy, synergies with ancillary and archival observations, the data reduction and analysis schemes adopted; and presents the first results obtained toward three of the survey targets, W3(OH), W3IRS5 and NGC7538IRS1. Robust measurements of the column densities of these hydrides -- obtained through widespread observations of absorption lines-- help address the questions raised, and there is a timely synergy between these observations and the development of theoretical models, particularly pertaining to the formation of H2 within the turbulent ISM. The provision of enhanced HyGAL data products will therefore serve as a legacy for future ISM studies., Comment: 33 pages, 17 figures, 7 tables. Accepted for publication in the Astrophysical Journal

The origin of the cold gas in central galaxies in groups is still a matter of debate. We present Multi-Unit Spectroscopic Explorer (MUSE) observations of 18 optically selected local Brightest Group Galaxies (BGGs) to study the kinematics and distribution of the optical emission-line gas. MUSE observations reveal a distribution of gas morphologies including ten complex networks of filaments extending up to 10 kpc to two compact (<3 kpc) and five extended (>5 kpc) disk-dominated structures. Some rotating disks show rings and elongated structures arising from the central disk. The kinematics of the stellar component is mainly rotation-dominated, which is very different from the disturbed kinematics and distribution found in the filamentary sources. The ionized gas is kinematically decoupled from the stellar component for most systems, suggesting an external origin for the gas. We also find that the Halpha luminosity correlates with the cold molecular mass. By exploring the thermodynamical properties of the hot atmospheres, we find that the filamentary sources and compact disks are found in systems with small central entropy values and tcool/teddy ratios. This suggests that, like for Brightest Cluster Galaxies in cool core clusters, the ionized gas are likely formed from hot halo gas condensations, consistently with the Chaotic Cold Accretion simulations (as shown via the C-ratio, Tat, and k-plot). We note that gaseous rotating disks are more frequent than in BCGs. An explanation for the origin of the gas in those objects is a contribution to gas fueling by mergers or group satellites, as qualitatively hinted by some sources of the present sample. Nonetheless, we discuss the possibility that some extended disks could also be a transition stage in an evolutionary sequence including filaments, extended disks and compact disks, as described by hot gas condensation models of cooling flows., Comment: 38 pages, and 16 figures. Resubmitted to A&A after referee report

As material from an infalling protostellar envelope hits the forming disk, an accretion shock may develop which could (partially) alter the envelope material entering the disk. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) indicate that emission originating from warm SO and SO$_2$ might be good tracers of such accretion shocks. The goal of this work is to test under what shock conditions the abundances of gas-phase SO and SO$_2$ increase in an accretion shock at the disk-envelope interface. Detailed shock models including gas dynamics are computed using the Paris-Durham shock code for non-magnetized J-type accretion shocks in typical inner envelope conditions. The effect of pre-shock density, shock velocity, and strength of the ultraviolet (UV) radiation field on the abundance of warm SO and SO$_2$ is explored. Warm gas-phase chemistry is efficient in forming SO under most J-type shock conditions considered. In lower-velocity (~3 km/s) shocks, the abundance of SO is increased through subsequent reactions starting from thermally desorbed CH$_4$ toward H$_2$CO and finally SO. In higher velocity (>4 km/s) shocks, both SO and SO$_2$ are formed through reactions of OH and atomic S. The strength of the UV radiation field is crucial for SO and in particular SO$_2$ formation through the photodissociation of H$_2$O. Thermal desorption of SO and SO$_2$ ice is only relevant in high-velocity (>5 km/s) shocks at high densities ($>10^7$ cm$^{-3}$). Warm emission from SO and SO$_2$ is a possible tracer of accretion shocks at the disk-envelope interface as long as a local UV field is present. Additional observations with ALMA at high-angular resolution could provide further constraints. Moreover, the James Webb Space Telescope will give access to other possible slow, dense shock tracers such as H$_2$, H$_2$O, and [S I] 25$\mu$m., Comment: Accepted to A&A, 18 pages, 17 figures

We present a comparative analysis of the $\rm CH^+$(1-0) and $\rm Ly \alpha$ lines, observed with the Atacama Large Millimeter Array (ALMA) and Keck telescope respectively, in the field of the submillimetre-selected galaxy (SMG) SMM\,J02399$-$0136 at $z\sim2.8$, which comprises a heavily obscured starburst galaxy and a broad absorption line quasar, immersed in a large $\rm Ly \alpha$ nebula. This comparison highlights the critical role played by turbulence in channeling the energy across gas phases and scales, splitting the energy trail between hot/thermal and cool/turbulent phases in the circum-galactic medium (CGM). The unique chemical and spectroscopic properties of $\rm CH^+$ are used to infer the existence of a massive ($\sim 3.5 \times 10^{10}$ ${\rm M}_\odot$), highly turbulent reservoir of diffuse molecular gas of radius $\sim 20\,$kpc coinciding with the core of the $\rm Ly \alpha$ nebula. The whole cool and cold CGM is shown to be inflowing towards the galaxies at a velocity $\sim$ 400 km$\,s^{-1}$. Several kpc-scale shocks are detected tentatively in $\rm CH^+$ emission. Their specific location in space and velocity with respect to the high-velocity $\rm Ly \alpha$ emission suggests that they lie at the interface of the inflowing CGM and the high-velocity $\rm Ly \alpha$ emission, and signpost the feeding of CGM turbulence by AGN- and stellar-driven outflows. The mass and energy budgets of the CGM require net mass accretion at a rate commensurate with the star formation rate (SFR). From this similarity, we infer that the merger-driven burst of star formation and black-hole growth are ultimately fuelled by large-scale gas accretion., Comment: 24 pages, 17 figures. Accepted for publication in MNRAS

As material from an infalling protostellar envelope hits the forming disk, an accretion shock may develop which could (partially) alter the envelope material entering the disk. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) indicate that emission originating from warm SO and SO$_2$ might be good tracers of such accretion shocks. The goal of this work is to test under what shock conditions the abundances of gas-phase SO and SO$_2$ increase in an accretion shock at the disk-envelope interface. Detailed shock models including gas dynamics are computed using the Paris-Durham shock code for non-magnetized J-type accretion shocks in typical inner envelope conditions. The effect of pre-shock density, shock velocity, and strength of the ultraviolet (UV) radiation field on the abundance of warm SO and SO$_2$ is explored. Warm gas-phase chemistry is efficient in forming SO under most J-type shock conditions considered. In lower-velocity (~3 km/s) shocks, the abundance of SO is increased through subsequent reactions starting from thermally desorbed CH$_4$ toward H$_2$CO and finally SO. In higher velocity (>4 km/s) shocks, both SO and SO$_2$ are formed through reactions of OH and atomic S. The strength of the UV radiation field is crucial for SO and in particular SO$_2$ formation through the photodissociation of H$_2$O. Thermal desorption of SO and SO$_2$ ice is only relevant in high-velocity (>5 km/s) shocks at high densities ($>10^7$ cm$^{-3}$). Warm emission from SO and SO$_2$ is a possible tracer of accretion shocks at the disk-envelope interface as long as a local UV field is present. Additional observations with ALMA at high-angular resolution could provide further constraints. Moreover, the James Webb Space Telescope will give access to other possible slow, dense shock tracers such as H$_2$, H$_2$O, and [S I] 25$\mu$m., Comment: Accepted to A&A, 18 pages, 17 figures

We present a comparative analysis of the $\rm CH^+$(1-0) and $\rm Ly \alpha$ lines, observed with the Atacama Large Millimeter Array (ALMA) and Keck telescope respectively, in the field of the submillimetre-selected galaxy (SMG) SMM\,J02399$-$0136 at $z\sim2.8$, which comprises a heavily obscured starburst galaxy and a broad absorption line quasar, immersed in a large $\rm Ly \alpha$ nebula. This comparison highlights the critical role played by turbulence in channeling the energy across gas phases and scales, splitting the energy trail between hot/thermal and cool/turbulent phases in the circum-galactic medium (CGM). The unique chemical and spectroscopic properties of $\rm CH^+$ are used to infer the existence of a massive ($\sim 3.5 \times 10^{10}$ ${\rm M}_\odot$), highly turbulent reservoir of diffuse molecular gas of radius $\sim 20\,$kpc coinciding with the core of the $\rm Ly \alpha$ nebula. The whole cool and cold CGM is shown to be inflowing towards the galaxies at a velocity $\sim$ 400 km$\,s^{-1}$. Several kpc-scale shocks are detected tentatively in $\rm CH^+$ emission. Their specific location in space and velocity with respect to the high-velocity $\rm Ly \alpha$ emission suggests that they lie at the interface of the inflowing CGM and the high-velocity $\rm Ly \alpha$ emission, and signpost the feeding of CGM turbulence by AGN- and stellar-driven outflows. The mass and energy budgets of the CGM require net mass accretion at a rate commensurate with the star formation rate (SFR). From this similarity, we infer that the merger-driven burst of star formation and black-hole growth are ultimately fuelled by large-scale gas accretion., Comment: 24 pages, 17 figures. Accepted for publication in MNRAS

The excitation of the filamentary gas structures surrounding giant elliptical galaxies at the center of cool-core clusters, a.k.a BCGs (brightest cluster galaxies), is key to our understanding of active galactic nucleus feedback, and of the impact of environmental and local effects on star formation. We investigate the contribution of the thermal radiation from the cooling flow surrounding BCGs to the excitation of the filaments. We explore the effects of small levels of extra-heating (turbulence), and of metallicity, on the optical and infrared lines. Using the Cloudy code, we model the photoionization and photodissociation of a slab of gas of optical depth AV{\leq}30mag at constant pressure, in order to calculate self-consistently all of the gas phases, from ionized gas to molecular gas. The ionizing source is the EUV and soft X-ray radiation emitted by the cooling gas. We test these models comparing their predictions to the rich multi-wavelength observations, from optical to submillimeter. These models reproduce most of the multi-wavelength spectra observed in the nebulae surrounding the BCGs, not only the LINER-like optical diagnostics: [O iii]{\lambda} 5007 {\AA}/H\b{eta}, [N ii]{\lambda} 6583 {\AA}/H{\alpha} and ([S ii]{\lambda} 6716 {\AA}+[S ii]{\lambda} 6731 {\AA})/H{\alpha} but also the infrared emission lines from the atomic gas. The modeled ro-vib H2 lines also match observations, which indicates that near and mid-IR H2 lines are mostly excited by collisions between H2 molecules and secondary electrons produced naturally inside the cloud by the interaction between the X-rays and the cold gas in the filament. However, there is still some tension between ionized and molecular line tracers (i.e. CO), which requires to optimize the cloud structure and the density of the molecular zone., Comment: 16 pages, 9 figures, 3 tables. Accepted for publication in A&A

Recent ALMA observations suggest that the highest velocity part of molecular protostellar jets are launched from the dust-sublimation regions of the accretion disks (<0.3 au). However, formation and survival of molecules in inner protostellar disk winds, in the presence of a harsh FUV radiation field and the absence of dust, remain unexplored. We aim at determining if simple molecules can be synthesized and spared in fast and collimated dust-free disk winds or if a fraction of dust is necessary to explain the observed molecular abundances. This work is based on the Paris-Durham shock code designed to model irradiated environments. Fundamental properties of the dust-free chemistry are investigated from single point models. A laminar 1D disk wind model is then built using a parametric flow geometry. This model includes time-dependent chemistry and the attenuation of the radiation field by gas-phase photoprocesses. We show that a small fraction of H2 (< 1e-2), primarily formed through the H- route, can efficiently initiate molecule synthesis such as CO and SiO above TK ~ 800 K. The attenuation of the radiation field by atomic species (eg. C, Si, S) proceeds through continuum self-shielding. This process ensures efficient formation of CO, OH, SiO, H2O through neutral-neutral reactions, and the survival of these molecules. Class 0 dust-free winds with high mass-loss rates ($\dot{M}_w >$ 2e-6 Msun/yr) are predicted to be rich in molecules if warm (TK > 800 K). The molecular content of disk winds is very sensitive to the presence of dust and a mass-fraction of surviving dust as small as 1e-5 significantly increases the H2O and SiO abundances. Chemistry of high-velocity jets is a powerful tool to probe their content in dust and uncover their launching point. Models of internal shocks are required to fully exploit the current (sub-)millimeter observations and prepare future JWST observations., Comment: accepted for publication in A&A

Shocks are ubiquitous in the interstellar and intergalactic media, where their chemical and radiative signatures reveal the physical conditions in which they arise. Detailed astrochemical models of shocks at all velocities are necessary to understand the physics of many environments including protostellar outflows, supernova remnants, and galactic outflows. We present an accurate treatment of the self-generated UV radiation in intermediate velocity, stationary, weakly magnetised, J-type, molecular shocks. Shock solutions computed with the Paris-Durham shock code are post-processed using a multi-level accelerated $\Lambda$-iteration radiative transfer algorithm to compute Ly$\alpha$, Ly$\beta$, and 2-photon continuum emission. The subsequent impacts on the ionisation and dissociation of key atomic and molecular species as well as on the heating by the photoelectric effect take the wavelength dependent cross-sections and the fluid velocity profile into account. We analyse shock models with velocities $V=25-60$ km/s, propagating in dense ($n \geq 10^4$ ${\rm cm}^{-3}$), shielded gas. Self-absorption traps Ly$\alpha$ photons in a small region in the shock, though a large fraction escapes into the line wings. We find a critical velocity $V\sim 30$ km/s above which shocks produce a Ly$\alpha$ photon flux exceeding that of the standard ISRF. The escaping photons generate a warm slab (T~100 K) ahead of the shock as well as pre-ionise C and S. These shocks are traced by bright atomic fine structure (e.g. O and S) and metastable (e.g. O and C) lines, substantive molecular emission (e.g. H2, OH, and CO), enhanced column densities of several species (e.g. CH+ and HCO+), as well as a severe destruction of H2O. As much as 13-21% of the initial kinetic energy of the shock escapes in Ly$\alpha$ and Ly$\beta$ photons if the dust opacity in the radiative precursor allows it., Comment: Accepted for publication in A&A, 16 main body pages, 8 appendices pages

Recent ALMA observations suggest that the highest velocity part of molecular protostellar jets are launched from the dust-sublimation regions of the accretion disks (<0.3 au). However, formation and survival of molecules in inner protostellar disk winds, in the presence of a harsh FUV radiation field and the absence of dust, remain unexplored. We aim at determining if simple molecules can be synthesized and spared in fast and collimated dust-free disk winds or if a fraction of dust is necessary to explain the observed molecular abundances. This work is based on the Paris-Durham shock code designed to model irradiated environments. Fundamental properties of the dust-free chemistry are investigated from single point models. A laminar 1D disk wind model is then built using a parametric flow geometry. This model includes time-dependent chemistry and the attenuation of the radiation field by gas-phase photoprocesses. We show that a small fraction of H2 (< 1e-2), primarily formed through the H- route, can efficiently initiate molecule synthesis such as CO and SiO above TK ~ 800 K. The attenuation of the radiation field by atomic species (eg. C, Si, S) proceeds through continuum self-shielding. This process ensures efficient formation of CO, OH, SiO, H2O through neutral-neutral reactions, and the survival of these molecules. Class 0 dust-free winds with high mass-loss rates ($\dot{M}_w >$ 2e-6 Msun/yr) are predicted to be rich in molecules if warm (TK > 800 K). The molecular content of disk winds is very sensitive to the presence of dust and a mass-fraction of surviving dust as small as 1e-5 significantly increases the H2O and SiO abundances. Chemistry of high-velocity jets is a powerful tool to probe their content in dust and uncover their launching point. Models of internal shocks are required to fully exploit the current (sub-)millimeter observations and prepare future JWST observations., Comment: accepted for publication in A&A

(abridged) The far-UV wavelength range (912-2000A) provides access to atomic and molecular transitions of many species the interstellar medium (ISM), circumgalactic medium (CGM), and intergalactic medium, within phases spanning a wide range of ionization, density, temperature, and molecular gas fraction. Far-UV space telescopes have enabled detailed studies of the ISM in the Milky Way thanks to absorption features appearing in the UV spectra of hot stars and yielding fundamental insights into the composition and physical characteristics of all phases of the ISM along with the processes that influence them. However, we have yet to design a spectrometer able to observe the full UV domain at resolving power R>10^5 with a signal-to-noise ratio SNR>500. Such a resolution is necessary to resolve lines from both the cold molecular hydrogen and the warm metal ions with a turbulent velocity of about 1 km s-1, and to differentiate distinct velocity components. Future UV spectroscopic studies of the Milky Way ISM must revolutionize our understanding of the ISM as a dynamical, unstable, and magnetized medium, and rise to the challenge brought forward by current theories. Another interesting prospect is to transpose the same level of details that has been reached for the Milky Way to the ISM in external galaxies, in particular in metal-poor galaxies, where the ISM chemical composition, physical conditions, and topology change dramatically, with significant consequences on the star-formation properties. Finally, we need to be able to perform statistical analyses of background quasar lines of sight intersecting the CGM of galaxies at various redshifts and to comprehend the role of gas exchanges and flows for galaxy evolution., Comment: This ESA white paper draws from and expands on the following white papers previously submitted for the Astro 2020 Decadal Survey, Gry et al. (2019), Lebouteiller et al. (2019), and Yan et al. (2019). The original ESA white paper was reformatted as a 2-column paper for the present document. (arXiv admin note: substantial text overlap with arXiv:1903.07295, arXiv:1903.06250, arXiv:1903.11065)

Multi-phase filamentary structures around Brightest Cluster Galaxies are likely a key step of AGN-feedback. We observed molecular gas in 3 cool cluster cores: Centaurus, Abell S1101, and RXJ1539.5 and gathered ALMA and MUSE data for 12 other clusters. Those observations show clumpy, massive and long, 3--25 kpc, molecular filaments, preferentially located around the radio bubbles inflated by the AGN (Active Galactic Nucleus). Two objects show nuclear molecular disks. The optical nebula is certainly tracing the warm envelopes of cold molecular filaments. Surprisingly, the radial profile of the H$\alpha$/CO flux ratio is roughly constant for most of the objects, suggesting that (i) between 1.2 to 7 times more cold gas could be present and (ii) local processes must be responsible for the excitation. Projected velocities are between 100--400 km s$^{-1}$, with disturbed kinematics and sometimes coherent gradients. This is likely due to the mixing in projection of several thin unresolved filaments. The velocity fields may be stirred by turbulence induced by bubbles, jets or merger-induced sloshing. Velocity and dispersions are low, below the escape velocity. Cold clouds should eventually fall back and fuel the AGN. We compare the filament's radial extent, r$_{fil}$, with the region where the X-ray gas can become thermally unstable. The filaments are always inside the low-entropy and short cooling time region, where t$_{cool}$/t$_{ff}$<20 (9 of 13 sources). The range t$_{cool}$/t$_{ff}$, 8-23 at r$_{fil}$, is likely due to (i) a more complex gravitational potential affecting the free-fall time (e.g., sloshing, mergers); (ii) the presence of inhomogeneities or uplifted gas in the ICM, affecting the cooling time. For some of the sources, r$_{fil}$ lies where the ratio of the cooling time to the eddy-turnover time, t$_{cool}$/t$_{eddy}$, is approximately unity., Comment: 33 Pages, 17 Figures, submitted to A&A

Multi-phase filamentary structures around Brightest Cluster Galaxies are likely a key step of AGN-feedback. We observed molecular gas in 3 cool cluster cores: Centaurus, Abell S1101, and RXJ1539.5 and gathered ALMA and MUSE data for 12 other clusters. Those observations show clumpy, massive and long, 3--25 kpc, molecular filaments, preferentially located around the radio bubbles inflated by the AGN (Active Galactic Nucleus). Two objects show nuclear molecular disks. The optical nebula is certainly tracing the warm envelopes of cold molecular filaments. Surprisingly, the radial profile of the H$\alpha$/CO flux ratio is roughly constant for most of the objects, suggesting that (i) between 1.2 to 7 times more cold gas could be present and (ii) local processes must be responsible for the excitation. Projected velocities are between 100--400 km s$^{-1}$, with disturbed kinematics and sometimes coherent gradients. This is likely due to the mixing in projection of several thin unresolved filaments. The velocity fields may be stirred by turbulence induced by bubbles, jets or merger-induced sloshing. Velocity and dispersions are low, below the escape velocity. Cold clouds should eventually fall back and fuel the AGN. We compare the filament's radial extent, r$_{fil}$, with the region where the X-ray gas can become thermally unstable. The filaments are always inside the low-entropy and short cooling time region, where t$_{cool}$/t$_{ff}$<20 (9 of 13 sources). The range t$_{cool}$/t$_{ff}$, 8-23 at r$_{fil}$, is likely due to (i) a more complex gravitational potential affecting the free-fall time (e.g., sloshing, mergers); (ii) the presence of inhomogeneities or uplifted gas in the ICM, affecting the cooling time. For some of the sources, r$_{fil}$ lies where the ratio of the cooling time to the eddy-turnover time, t$_{cool}$/t$_{eddy}$, is approximately unity., Comment: 33 Pages, 17 Figures, submitted to A&A

(abridged) The far-UV wavelength range (912-2000A) provides access to atomic and molecular transitions of many species the interstellar medium (ISM), circumgalactic medium (CGM), and intergalactic medium, within phases spanning a wide range of ionization, density, temperature, and molecular gas fraction. Far-UV space telescopes have enabled detailed studies of the ISM in the Milky Way thanks to absorption features appearing in the UV spectra of hot stars and yielding fundamental insights into the composition and physical characteristics of all phases of the ISM along with the processes that influence them. However, we have yet to design a spectrometer able to observe the full UV domain at resolving power R>10^5 with a signal-to-noise ratio SNR>500. Such a resolution is necessary to resolve lines from both the cold molecular hydrogen and the warm metal ions with a turbulent velocity of about 1 km s-1, and to differentiate distinct velocity components. Future UV spectroscopic studies of the Milky Way ISM must revolutionize our understanding of the ISM as a dynamical, unstable, and magnetized medium, and rise to the challenge brought forward by current theories. Another interesting prospect is to transpose the same level of details that has been reached for the Milky Way to the ISM in external galaxies, in particular in metal-poor galaxies, where the ISM chemical composition, physical conditions, and topology change dramatically, with significant consequences on the star-formation properties. Finally, we need to be able to perform statistical analyses of background quasar lines of sight intersecting the CGM of galaxies at various redshifts and to comprehend the role of gas exchanges and flows for galaxy evolution., Comment: This ESA white paper draws from and expands on the following white papers previously submitted for the Astro 2020 Decadal Survey, Gry et al. (2019), Lebouteiller et al. (2019), and Yan et al. (2019). The original ESA white paper was reformatted as a 2-column paper for the present document. (arXiv admin note: substantial text overlap with arXiv:1903.07295, arXiv:1903.06250, arXiv:1903.11065)

The transition between atomic and molecular hydrogen is associated with important changes in the structure of interstellar clouds, and marks the beginning of interstellar chemistry. Because of the relatively simple networks controlling their abundances, molecular ions are usually good probes of the underlying physical conditions including for instance the fraction of gas in molecular form or the fractional ionization. In this paper we focus on three possible probes of the molecular hydrogen column density, HCO+, HOC+, and CF+. We presented high sensitivity ALMA absorption data toward a sample of compact HII regions and bright QSOs with prominent foreground absorption, in the ground state transitions of the molecular ions HCO+, HOC+, and CF+ and the neutral species HCN and HNC, and from the excited state transitions of C3H+(4-3) and 13CS(2-1). These data are compared with Herschel absorption spectra of the ground state transition of HF and p-H2O. We show that the HCO+, HOC+, and CF+ column densities are well correlated with each other. HCO+ and HOC+ are tightly correlated with p-H2O, while they exhibit a different correlation pattern with HF depending on whether the absorbing matter is located in the Galactic disk or in the central molecular zone. We report new detections of C3H+ confirming that this ion is ubiquitous in the diffuse matter, with an abundance relative to H2 of ~7E-11. We confirm that the CF+ abundance is lower than predicted by simple chemical models and propose that the rate of the main formation reaction is lower by a factor of about 3 than usually assumed. In the absence of CH or HF data, we recommend to use the ground state transitions of HCO+, CCH, and HOC+ to trace diffuse molecular hydrogen, with mean abundances relative to H2 of 3E-9, 4E-8 and 4E-11., Comment: 15 pages, 9 Figures, accepted by Astronomy and Astrophysics. slightly modified abstract

Cosmic rays are predominantly accelerated in shocks associated with star formation such as supernova remnants and stellar wind bubbles, so the cosmic-ray flux and thus cosmic-ray ionization rate, $\zeta_{\rm H}$, should correlate with the star-formation rate in a galaxy. Submillimeter bright galaxies (SMGs) are some of the most prolific star forming galaxies in the Universe, and gravitationally lensed SMGs provide bright continuum sources suitable for absorption line studies. Abundances of OH$^+$ and H$_2$O$^+$ are useful for inferring $\zeta_{\rm H}$ when combined with chemical models, and have been used for this purpose within the Milky Way. At redshifts $z\gtrsim2$ transitions out of the ground rotational states of OH$^+$ and H$_2$O$^+$ are observable with ALMA, and we present observations of both molecules in absorption toward the lensed SMGs SMM J2135$-$0102 and SDP 17b. These detections enable an exploration of $\zeta_{\rm H}$ in galaxies with extreme star formation and high supernova rates, both of which should significantly enhance cosmic-ray production. The observed OH$^+$ and H$_2$O$^+$ absorption is thought to arise in massive, extended halos of cool, diffuse gas that surround these galaxies. Using a chemical model designed to focus on the reaction network important to both species, we infer cosmic-ray ionization rates of $\zeta_{\rm H}\sim10^{-16}$-$10^{-14}$ s$^{-1}$ in these extended gaseous halos. As our estimates come from gas that is far away from the sites of cosmic-ray acceleration, they imply that cosmic-ray ionization rates in the compact regions where star formation occurs in these galaxies are orders of magnitude higher., Comment: 16 pages, 5 figures, 3 tables, accepted for publication in ApJ

Aims. We aim to investigate the chemistry of internal photon-dominated regions (PDRs) surrounding deeply embedded hypercompact (HC) and ultracompact (UC) HII regions. We search for specific tracers of this evolutionary stage of massive star formation that can be detected with current astronomical facilities. Methods. We modeled hot cores with embedded HC/UCHII regions (called HII region models in the article despite the fact that we do not model the HII region itself), by coupling the astrochemical code Saptarsy to a radiative transfer framework obtaining the spatio-temporal evolution of abundances as well as time-dependent synthetic spectra. In these models where we focused on the internal PDR surrounding the HII region, the gas temperature is set to the dust temperature and we do not include dynamics thus the density structure is fixed. We compared this to hot molecular core (HMC) models and studied the effect on the chemistry of the radiation field which is included in the HII region models only during the computation of abundances. In addition, we investigated the chemical evolution of the gas surrounding HII regions with models of different densities at the ionization front, different sizes of the ionized cavity and different initial abundances. Results. We obtain the time evolution of synthetic spectra for a dozen of selected species as well as ratios of their integrated intensities. We find that some molecules such as C, N2H+, CN, and HCO do not trace the inner core and so are not good tracers to distinguish the HII/PDR regions to the HMCs phase. On the contrary, C+ and O trace the internal PDRs, in the two models starting with different initial abundances, but are unfortunately currently unobservable with the current achievable spatial resolution because of the very thin internal PDR (Delta r(PDR) < 100 AU). The emission of these two tracers is very dependent on the size of the HII region and on the density in the PDR. In addition, we find that the abundance

Aims. We aim to investigate the chemistry of internal photon-dominated regions (PDRs) surrounding deeply embedded hypercompact (HC) and ultracompact (UC) HII regions. We search for specific tracers of this evolutionary stage of massive star formation that can be detected with current astronomical facilities. Methods. We modeled hot cores with embedded HC/UCHII regions (called HII region models in the article despite the fact that we do not model the HII region itself), by coupling the astrochemical code Saptarsy to a radiative transfer framework obtaining the spatio-temporal evolution of abundances as well as time-dependent synthetic spectra. In these models where we focused on the internal PDR surrounding the HII region, the gas temperature is set to the dust temperature and we do not include dynamics thus the density structure is fixed. We compared this to hot molecular core (HMC) models and studied the effect on the chemistry of the radiation field which is included in the HII region models only during the computation of abundances. In addition, we investigated the chemical evolution of the gas surrounding HII regions with models of different densities at the ionization front, different sizes of the ionized cavity and different initial abundances. Results. We obtain the time evolution of synthetic spectra for a dozen of selected species as well as ratios of their integrated intensities. We find that some molecules such as C, N2H+, CN, and HCO do not trace the inner core and so are not good tracers to distinguish the HII/PDR regions to the HMCs phase. On the contrary, C+ and O trace the internal PDRs, in the two models starting with different initial abundances, but are unfortunately currently unobservable with the current achievable spatial resolution because of the very thin internal PDR (Delta r(PDR) < 100 AU). The emission of these two tracers is very dependent on the size of the HII region and on the density in the PDR. In addition, we find that the abundance

The transition between atomic and molecular hydrogen is associated with important changes in the structure of interstellar clouds, and marks the beginning of interstellar chemistry. Because of the relatively simple networks controlling their abundances, molecular ions are usually good probes of the underlying physical conditions including for instance the fraction of gas in molecular form or the fractional ionization. In this paper we focus on three possible probes of the molecular hydrogen column density, HCO+, HOC+, and CF+. We presented high sensitivity ALMA absorption data toward a sample of compact HII regions and bright QSOs with prominent foreground absorption, in the ground state transitions of the molecular ions HCO+, HOC+, and CF+ and the neutral species HCN and HNC, and from the excited state transitions of C3H+(4-3) and 13CS(2-1). These data are compared with Herschel absorption spectra of the ground state transition of HF and p-H2O. We show that the HCO+, HOC+, and CF+ column densities are well correlated with each other. HCO+ and HOC+ are tightly correlated with p-H2O, while they exhibit a different correlation pattern with HF depending on whether the absorbing matter is located in the Galactic disk or in the central molecular zone. We report new detections of C3H+ confirming that this ion is ubiquitous in the diffuse matter, with an abundance relative to H2 of ~7E-11. We confirm that the CF+ abundance is lower than predicted by simple chemical models and propose that the rate of the main formation reaction is lower by a factor of about 3 than usually assumed. In the absence of CH or HF data, we recommend to use the ground state transitions of HCO+, CCH, and HOC+ to trace diffuse molecular hydrogen, with mean abundances relative to H2 of 3E-9, 4E-8 and 4E-11., Comment: 15 pages, 9 Figures, accepted by Astronomy and Astrophysics. slightly modified abstract

Cosmic rays are predominantly accelerated in shocks associated with star formation such as supernova remnants and stellar wind bubbles, so the cosmic-ray flux and thus cosmic-ray ionization rate, $\zeta_{\rm H}$, should correlate with the star-formation rate in a galaxy. Submillimeter bright galaxies (SMGs) are some of the most prolific star forming galaxies in the Universe, and gravitationally lensed SMGs provide bright continuum sources suitable for absorption line studies. Abundances of OH$^+$ and H$_2$O$^+$ are useful for inferring $\zeta_{\rm H}$ when combined with chemical models, and have been used for this purpose within the Milky Way. At redshifts $z\gtrsim2$ transitions out of the ground rotational states of OH$^+$ and H$_2$O$^+$ are observable with ALMA, and we present observations of both molecules in absorption toward the lensed SMGs SMM J2135$-$0102 and SDP 17b. These detections enable an exploration of $\zeta_{\rm H}$ in galaxies with extreme star formation and high supernova rates, both of which should significantly enhance cosmic-ray production. The observed OH$^+$ and H$_2$O$^+$ absorption is thought to arise in massive, extended halos of cool, diffuse gas that surround these galaxies. Using a chemical model designed to focus on the reaction network important to both species, we infer cosmic-ray ionization rates of $\zeta_{\rm H}\sim10^{-16}$-$10^{-14}$ s$^{-1}$ in these extended gaseous halos. As our estimates come from gas that is far away from the sites of cosmic-ray acceleration, they imply that cosmic-ray ionization rates in the compact regions where star formation occurs in these galaxies are orders of magnitude higher., Comment: 16 pages, 5 figures, 3 tables, accepted for publication in ApJ

Starburst galaxies at the peak of cosmic star formation are among the most extreme starforming engines in the universe, producing stars over ~100 Myr. The star formation rates of these galaxies, which exceed 100 $M_\odot$ per year, require large reservoirs of cold molecular gas to be delivered to their cores, despite strong feedback from stars or active galactic nuclei. Starburst galaxies are therefore ideal targets to unravel the critical interplay between this feedback and the growth of a galaxy. The methylidyne cation, CH$^+$, is a most useful molecule for such studies because it cannot form in cold gas without supra-thermal energy input, so its presence highlights dissipation of mechanical energy or strong UV irradiation. Here, we report the detection of CH$^+$(J=1-0) emission and absorption lines in the spectra of six lensed starburst galaxies at redshifts z~2.5. This line has such a high critical density for excitation that it is emitted only in very dense ($>10^5$ cm$^{-3}$) gas, and is absorbed in low-density gas. We find that the CH$^+$ emission lines, which are broader than 1000 km s$^{-1}$, originate in dense shock waves powered by hot galactic winds. The CH$^+$ absorption lines reveal highly turbulent reservoirs of cool ($T\sim 100$K), low-density gas, extending far outside (>10 kpc) the starburst cores (radii <1 kpc). We show that the galactic winds sustain turbulence in the 10 kpc-scale environments of the starburst cores, processing these environments into multi-phase, gravitationally bound reservoirs. However, the mass outflow rates are found to be insufficient to balance the star formation rates. Another mass input is therefore required for these reservoirs, which could be provided by on-going mergers or cold stream accretion. Our results suggest that galactic feedback, coupled jointly to turbulence and gravity, extends the starburst phase instead of quenching it., Comment: Published as a Nature Letter

Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]$_{\rm 63 \mu m}$ emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. We present observations of the OI $^3$P$_1 \rightarrow$ $^3$P$_2$, OH between $^2\Pi_{1/2}$ $J = 3/2$ and $J = 1/2$ at 1837.8 GHz, and CO (16-15) lines with SOFIA-GREAT at three positions in the Cep E outflow: mm (the driving protostar), BI (in the southern lobe), and BII (the terminal position in the southern lobe). The CO line is detected at all three positions. The OI line is detected in BI and BII, whereas the OH line is not detected. In BII, we identify three kinematical components in OI and CO, already detected in CO: the jet, the HH377 terminal bow-shock, and the outflow cavity. The OI column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is where the abundance ratio of OI to CO is the lowest (about 0.2), whereas the jet component is atomic (ratio $\sim$2.7). In the jet, we compare the OI observations with shock models that successfully fit the integrated intensity of 10 CO lines: these models do not fit the OI data. The high intensity of OI emission points towards the propagation of additional dissociative or alternative FUV-irradiated shocks, where the illumination comes from the shock itself. From the sample of low-to-high mass protostellar outflows where similar observations have been performed, the effects of illumination seem to increase with the mass of the protostar., Comment: 9 pages, 4 figures, 6 tables, Astronomy and astrophysics, in press

Context. Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jet and outflow systems are launched, and to quantify their chemical and energetic impacts on the surrounding medium. Aims. We performed a high-spectral resolution study of the [OI](63 mu m) emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. The goal is to determine the structure of the outflow, to constrain the chemical conditions in the various components, and to understand the nature of the underlying shocks, thus probing the origin of the mass-loss phenomenon. Methods. We present observations of the O I P-3(1) -> P-3(2), OH between (2)Pi(1/2) J = 3/2 and J = 1/2 at 1837.8 GHz, and CO (16-15) lines with the GREAT receiver onboard SOFIA towards three positions in the Cep E protostellar outflow: Cep E-mm (the driving protostar), Cep E-BI (in the southern lobe), and Cep E-BII (the terminal position in the southern lobe). Results. The CO (16-15) line is detected at all three positions. The [OI](63 mu m) line is detected in Cep E-BI and BII, whereas the OH line is not detected. In Cep E-BII, we identify three kinematical components in O I and CO. These were already detected in CO transitions and relate to spatial components: the jet, the HH377 terminal bow-shock, and the outflow cavity. We measure line temperature and line integrated intensity ratios for all components. The O I column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is the region where the abundance ratio of O I to CO is the lowest (about 0.2), whereas the jet component is atomic (N(O I) /N(CO) similar to 2.7). In the jet, we compare the [OI](63 mu m) observations with shock models that successfully fit the integrated intensity of 10 CO lines. We find that these models m

Context. The methylidyne cation (CH+) and hydroxyl (OH) are key molecules in the warm interstellar chemistry, but their formation and excitation mechanisms are not well understood. Their abundance and excitation are predicted to be enhanced by the presence of vibrationally excited H-2 or hot gas (similar to 500-1000 K) in photodissociation regions (PDRs) with high incident far-ultraviolet (FUV) radiation field. The excitation may also originate in dense gas (>10(5) cm(-3)) followed by nonreactive collisions with H-2, H, and electrons. Previous observations of the Orion Bar suggest that the rotationally excited CH+ and OH correlate with the excited CO, which is a tracer of dense and warm gas, and that formation pumping contributes to CH+ excitation. Aims. Our goal is to examine the spatial distribution of the rotationally excited CH+ and OH emission lines in the Orion Bar to establish their physical origin and main formation and excitation mechanisms. Methods. We present spatially sampled maps of the CH+ J = 3-2 transition at 119.8 mu m and the OH Lambda doublet at 84 mu m in the Orion Bar over an area of 110 '' Chi 110 '' with Herschel/PACS. We compare the spatial distribution of these molecules with those of their chemical precursors, C-broken vertical bar, O and H-2, and tracers of warm and dense gas (high-J CO). We assess the spatial variation of the CH+ J = 2-1 velocity-resolved line profile at 1669 GHz with Herschel/HIFI spectrometer observations. Results. The OH and especially CH+ lines correlate well with the high-J CO emission and delineate the warm and dense molecular region at the edge of the Bar. While notably similar, the differences in the CH+ and OH morphologies indicate that CH+ formation and excitation are strongly related to the observed vibrationally excited H-2. This, together with the observed broad CH+ line widths, indicates that formation pumping contributes to the excitation of this reactive molecular ion. Interestingly, the peak of the rotati

Context. The methylidyne cation (CH+) and hydroxyl (OH) are key molecules in the warm interstellar chemistry, but their formation and excitation mechanisms are not well understood. Their abundance and excitation are predicted to be enhanced by the presence of vibrationally excited H-2 or hot gas (similar to 500-1000 K) in photodissociation regions (PDRs) with high incident far-ultraviolet (FUV) radiation field. The excitation may also originate in dense gas (>10(5) cm(-3)) followed by nonreactive collisions with H-2, H, and electrons. Previous observations of the Orion Bar suggest that the rotationally excited CH+ and OH correlate with the excited CO, which is a tracer of dense and warm gas, and that formation pumping contributes to CH+ excitation. Aims. Our goal is to examine the spatial distribution of the rotationally excited CH+ and OH emission lines in the Orion Bar to establish their physical origin and main formation and excitation mechanisms. Methods. We present spatially sampled maps of the CH+ J = 3-2 transition at 119.8 mu m and the OH Lambda doublet at 84 mu m in the Orion Bar over an area of 110 '' Chi 110 '' with Herschel/PACS. We compare the spatial distribution of these molecules with those of their chemical precursors, C-broken vertical bar, O and H-2, and tracers of warm and dense gas (high-J CO). We assess the spatial variation of the CH+ J = 2-1 velocity-resolved line profile at 1669 GHz with Herschel/HIFI spectrometer observations. Results. The OH and especially CH+ lines correlate well with the high-J CO emission and delineate the warm and dense molecular region at the edge of the Bar. While notably similar, the differences in the CH+ and OH morphologies indicate that CH+ formation and excitation are strongly related to the observed vibrationally excited H-2. This, together with the observed broad CH+ line widths, indicates that formation pumping contributes to the excitation of this reactive molecular ion. Interestingly, the peak of the rotati

Context. Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jet and outflow systems are launched, and to quantify their chemical and energetic impacts on the surrounding medium. Aims. We performed a high-spectral resolution study of the [OI](63 mu m) emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. The goal is to determine the structure of the outflow, to constrain the chemical conditions in the various components, and to understand the nature of the underlying shocks, thus probing the origin of the mass-loss phenomenon. Methods. We present observations of the O I P-3(1) -> P-3(2), OH between (2)Pi(1/2) J = 3/2 and J = 1/2 at 1837.8 GHz, and CO (16-15) lines with the GREAT receiver onboard SOFIA towards three positions in the Cep E protostellar outflow: Cep E-mm (the driving protostar), Cep E-BI (in the southern lobe), and Cep E-BII (the terminal position in the southern lobe). Results. The CO (16-15) line is detected at all three positions. The [OI](63 mu m) line is detected in Cep E-BI and BII, whereas the OH line is not detected. In Cep E-BII, we identify three kinematical components in O I and CO. These were already detected in CO transitions and relate to spatial components: the jet, the HH377 terminal bow-shock, and the outflow cavity. We measure line temperature and line integrated intensity ratios for all components. The O I column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is the region where the abundance ratio of O I to CO is the lowest (about 0.2), whereas the jet component is atomic (N(O I) /N(CO) similar to 2.7). In the jet, we compare the [OI](63 mu m) observations with shock models that successfully fit the integrated intensity of 10 CO lines. We find that these models m

Starburst galaxies at the peak of cosmic star formation are among the most extreme starforming engines in the universe, producing stars over ~100 Myr. The star formation rates of these galaxies, which exceed 100 $M_\odot$ per year, require large reservoirs of cold molecular gas to be delivered to their cores, despite strong feedback from stars or active galactic nuclei. Starburst galaxies are therefore ideal targets to unravel the critical interplay between this feedback and the growth of a galaxy. The methylidyne cation, CH$^+$, is a most useful molecule for such studies because it cannot form in cold gas without supra-thermal energy input, so its presence highlights dissipation of mechanical energy or strong UV irradiation. Here, we report the detection of CH$^+$(J=1-0) emission and absorption lines in the spectra of six lensed starburst galaxies at redshifts z~2.5. This line has such a high critical density for excitation that it is emitted only in very dense ($>10^5$ cm$^{-3}$) gas, and is absorbed in low-density gas. We find that the CH$^+$ emission lines, which are broader than 1000 km s$^{-1}$, originate in dense shock waves powered by hot galactic winds. The CH$^+$ absorption lines reveal highly turbulent reservoirs of cool ($T\sim 100$K), low-density gas, extending far outside (>10 kpc) the starburst cores (radii <1 kpc). We show that the galactic winds sustain turbulence in the 10 kpc-scale environments of the starburst cores, processing these environments into multi-phase, gravitationally bound reservoirs. However, the mass outflow rates are found to be insufficient to balance the star formation rates. Another mass input is therefore required for these reservoirs, which could be provided by on-going mergers or cold stream accretion. Our results suggest that galactic feedback, coupled jointly to turbulence and gravity, extends the starburst phase instead of quenching it., Comment: Published as a Nature Letter

Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]$_{\rm 63 \mu m}$ emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. We present observations of the OI $^3$P$_1 \rightarrow$ $^3$P$_2$, OH between $^2\Pi_{1/2}$ $J = 3/2$ and $J = 1/2$ at 1837.8 GHz, and CO (16-15) lines with SOFIA-GREAT at three positions in the Cep E outflow: mm (the driving protostar), BI (in the southern lobe), and BII (the terminal position in the southern lobe). The CO line is detected at all three positions. The OI line is detected in BI and BII, whereas the OH line is not detected. In BII, we identify three kinematical components in OI and CO, already detected in CO: the jet, the HH377 terminal bow-shock, and the outflow cavity. The OI column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is where the abundance ratio of OI to CO is the lowest (about 0.2), whereas the jet component is atomic (ratio $\sim$2.7). In the jet, we compare the OI observations with shock models that successfully fit the integrated intensity of 10 CO lines: these models do not fit the OI data. The high intensity of OI emission points towards the propagation of additional dissociative or alternative FUV-irradiated shocks, where the illumination comes from the shock itself. From the sample of low-to-high mass protostellar outflows where similar observations have been performed, the effects of illumination seem to increase with the mass of the protostar., Comment: 9 pages, 4 figures, 6 tables, Astronomy and astrophysics, in press

The abundance of CH+ and OH and excitation are predicted to be enhanced by the presence of vibrationally excited H2 or hot gas (~500-1000 K) in PDRs with high incident FUV radiation field. The excitation may also originate in dense gas (>10^5 cm-3) followed by nonreactive collisions. Previous observations suggest that the CH+ and OH correlate with dense and warm gas, and formation pumping contributes to CH+ excitation. We examine the spatial distribution of the CH+ and OH emission in the Orion Bar to establish their physical origin and main formation and excitation mechanisms. We present spatially sampled maps of the CH+ J=3-2 transition at 119.8 {\mu}m and the OH {\Lambda}-doublet at 84 {\mu}m in the Orion Bar over an area of 110"x110" with Herschel (PACS). We compare the spatial distribution of these molecules with those of their chemical precursors, C+, O and H2, and tracers of warm and dense gas. We assess the spatial variation of CH+ J=2-1 velocity-resolved line profile observed with Herschel (HIFI). The OH and CH+ lines correlate well with the high-J CO emission and delineate the warm and dense molecular region. While similar, the differences in the CH+ and OH morphologies indicate that CH+ formation and excitation are related to the observed vibrationally excited H2. This indicates that formation pumping contributes to the excitation of CH+. Interestingly, the peak of the rotationally excited OH 84 {\mu}m emission coincides with a bright young object, proplyd 244-440, which shows that OH can be an excellent tracer of UV-irradiated dense gas. The spatial distribution of CH+ and OH revealed in our maps is consistent with previous modeling studies. Both formation pumping and nonreactive collisions in a UV-irradiated dense gas are important CH+ J=3-2 excitation processes. The excitation of the OH {\Lambda}-doublet at 84 {\mu}m is mainly sensitive to the temperature and density., Comment: Accepted to A&A

More complete knowledge of galaxy evolution requires understanding the process of star formation and interaction between the interstellar radiation field and the interstellar medium in galactic environments traversing a wide range of physical parameter space. Here we focus on the impact of massive star formation on the surrounding low metallicity ISM in 30 Doradus in the Large Magellanic Cloud. A low metal abundance, as is the case of some galaxies of the early universe, results in less ultra-violet shielding for the formation of the molecular gas necessary for star formation to proceed. The half-solar metallicity gas in this region is strongly irradiated by the super star cluster R136, making it an ideal laboratory to study the structure of the ISM in an extreme environment. Our spatially resolved study investigates the gas heating and cooling mechanisms, particularly in the photo-dissociation regions where the chemistry and thermal balance are regulated by far-ultraviolet photons (6 eV< h\nu <13.6 eV). We present Herschel observations of far-infrared fine-structure lines obtained with PACS and SPIRE/FTS. We have combined atomic fine-structure lines from Herschel and Spitzer observations with ground-based CO data to provide diagnostics on the properties and the structure of the gas by modeling it with the Meudon PDR code. We derive the spatial distribution of the radiation field, the pressure, the size, and the filling factor of the photodissociated gas and molecular clouds. We find a range of pressure of ~ 10^5 - 1.7x10^6 cm^{-3} K and a range of incident radiation field G_UV ~ 10^2 - 2.5x10^4 through PDR modeling. Assuming a plane-parallel geometry and a uniform medium, we find a total extinction of 1-3 mag , which correspond to a PDR cloud size of 0.2 to 3pc, with small CO depth scale of 0.06 to 0.5pc. We also determine the three dimensional structure of the gas. (Abridged), Comment: 20 pages, 23 figures, accepted in A&A

The abundance of CH+ and OH and excitation are predicted to be enhanced by the presence of vibrationally excited H2 or hot gas (~500-1000 K) in PDRs with high incident FUV radiation field. The excitation may also originate in dense gas (>10^5 cm-3) followed by nonreactive collisions. Previous observations suggest that the CH+ and OH correlate with dense and warm gas, and formation pumping contributes to CH+ excitation. We examine the spatial distribution of the CH+ and OH emission in the Orion Bar to establish their physical origin and main formation and excitation mechanisms. We present spatially sampled maps of the CH+ J=3-2 transition at 119.8 {\mu}m and the OH {\Lambda}-doublet at 84 {\mu}m in the Orion Bar over an area of 110"x110" with Herschel (PACS). We compare the spatial distribution of these molecules with those of their chemical precursors, C+, O and H2, and tracers of warm and dense gas. We assess the spatial variation of CH+ J=2-1 velocity-resolved line profile observed with Herschel (HIFI). The OH and CH+ lines correlate well with the high-J CO emission and delineate the warm and dense molecular region. While similar, the differences in the CH+ and OH morphologies indicate that CH+ formation and excitation are related to the observed vibrationally excited H2. This indicates that formation pumping contributes to the excitation of CH+. Interestingly, the peak of the rotationally excited OH 84 {\mu}m emission coincides with a bright young object, proplyd 244-440, which shows that OH can be an excellent tracer of UV-irradiated dense gas. The spatial distribution of CH+ and OH revealed in our maps is consistent with previous modeling studies. Both formation pumping and nonreactive collisions in a UV-irradiated dense gas are important CH+ J=3-2 excitation processes. The excitation of the OH {\Lambda}-doublet at 84 {\mu}m is mainly sensitive to the temperature and density., Comment: Accepted to A&A

More complete knowledge of galaxy evolution requires understanding the process of star formation and interaction between the interstellar radiation field and the interstellar medium in galactic environments traversing a wide range of physical parameter space. Here we focus on the impact of massive star formation on the surrounding low metallicity ISM in 30 Doradus in the Large Magellanic Cloud. A low metal abundance, as is the case of some galaxies of the early universe, results in less ultra-violet shielding for the formation of the molecular gas necessary for star formation to proceed. The half-solar metallicity gas in this region is strongly irradiated by the super star cluster R136, making it an ideal laboratory to study the structure of the ISM in an extreme environment. Our spatially resolved study investigates the gas heating and cooling mechanisms, particularly in the photo-dissociation regions where the chemistry and thermal balance are regulated by far-ultraviolet photons (6 eV< h\nu <13.6 eV). We present Herschel observations of far-infrared fine-structure lines obtained with PACS and SPIRE/FTS. We have combined atomic fine-structure lines from Herschel and Spitzer observations with ground-based CO data to provide diagnostics on the properties and the structure of the gas by modeling it with the Meudon PDR code. We derive the spatial distribution of the radiation field, the pressure, the size, and the filling factor of the photodissociated gas and molecular clouds. We find a range of pressure of ~ 10^5 - 1.7x10^6 cm^{-3} K and a range of incident radiation field G_UV ~ 10^2 - 2.5x10^4 through PDR modeling. Assuming a plane-parallel geometry and a uniform medium, we find a total extinction of 1-3 mag , which correspond to a PDR cloud size of 0.2 to 3pc, with small CO depth scale of 0.06 to 0.5pc. We also determine the three dimensional structure of the gas. (Abridged), Comment: 20 pages, 23 figures, accepted in A&A

We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright continuum sources. We have used the GREAT instrument on SOFIA to observe the 1383 GHz $^2\Pi_{3/2} J=5/2-3/2$ transitions of SH towards the star-forming regions W31C, G29.96-0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards all five sources; and the EMIR receivers on the IRAM 30m telescope at Pico Veleta to detect the H$_2$S 1(10)-1(01), CS J=2-1 and SO 3(2)-2(1) transitions. In nine foreground absorption components detected towards these sources, the inferred column densities of the four detected molecules showed relatively constant ratios, with N(SH)/N(H$_2$S) in the range 1.1 - 3.0, N(CS)/N(H$_2$S) in the range 0.32 - 0.61, and N(SO)/N(H$_2$S) in the range 0.08 - 0.30. The observed SH/H$_2$ ratios - in the range (0.5-2.6) $\times 10^{-8}$ - indicate that SH (and other sulphur-bearing molecules) account for << 1% of the gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules, however, greatly exceed those predicted by standard models of cold diffuse molecular clouds, providing further evidence for the enhancement of endothermic reaction rates by elevated temperatures or ion-neutral drift. We have considered the observed abundance ratios in the context of shock and turbulent dissipation region (TDR) models. Using the TDR model, we find that the turbulent energy available at large scale in the diffuse ISM is sufficient to explain the observed column densities of SH and CS. Standard shock and TDR models, however, fail to reproduce the column densities of H$_2$S and SO by a factor of about 10; more elaborate shock models - in which account is taken of the velocity drift, relative to H$_2$, of SH molecules produced by the dissociative recombination of H$_3$S$^+$ - reduce this discrepancy to a factor ~ 3., Comment: 30 pages, accepted for publication in A&A

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen-bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H-2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta H) and molecular hydrogen fraction (f(H2)). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042 +/- 0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H-2 fractions. We also infer zeta H throughout our sample, and find a lognormal distribution with mean log(zeta(H)) = -15.75 (zeta(H) = 1.78 x 10(-16) s(-1)) and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H-3(+) observations. Ionization rates in the Galactic center tend to be 10-100 times larger than found in the Galactic disk, also in accord with prior studies.

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen-bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H-2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta H) and molecular hydrogen fraction (f(H2)). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042 +/- 0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H-2 fractions. We also infer zeta H throughout our sample, and find a lognormal distribution with mean log(zeta(H)) = -15.75 (zeta(H) = 1.78 x 10(-16) s(-1)) and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H-3(+) observations. Ionization rates in the Galactic center tend to be 10-100 times larger than found in the Galactic disk, also in accord with prior studies.

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen-bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H-2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta H) and molecular hydrogen fraction (f(H2)). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042 +/- 0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H-2 fractions. We also infer zeta H throughout our sample, and find a lognormal distribution with mean log(zeta(H)) = -15.75 (zeta(H) = 1.78 x 10(-16) s(-1)) and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H-3(+) observations. Ionization rates in the Galactic center tend to be 10-100 times larger than found in the Galactic disk, also in accord with prior studies.

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen-bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H-2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta H) and molecular hydrogen fraction (f(H2)). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042 +/- 0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H-2 fractions. We also infer zeta H throughout our sample, and find a lognormal distribution with mean log(zeta(H)) = -15.75 (zeta(H) = 1.78 x 10(-16) s(-1)) and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H-3(+) observations. Ionization rates in the Galactic center tend to be 10-100 times larger than found in the Galactic disk, also in accord with prior studies.

We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright continuum sources. We have used the GREAT instrument on SOFIA to observe the 1383 GHz $^2\Pi_{3/2} J=5/2-3/2$ transitions of SH towards the star-forming regions W31C, G29.96-0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards all five sources; and the EMIR receivers on the IRAM 30m telescope at Pico Veleta to detect the H$_2$S 1(10)-1(01), CS J=2-1 and SO 3(2)-2(1) transitions. In nine foreground absorption components detected towards these sources, the inferred column densities of the four detected molecules showed relatively constant ratios, with N(SH)/N(H$_2$S) in the range 1.1 - 3.0, N(CS)/N(H$_2$S) in the range 0.32 - 0.61, and N(SO)/N(H$_2$S) in the range 0.08 - 0.30. The observed SH/H$_2$ ratios - in the range (0.5-2.6) $\times 10^{-8}$ - indicate that SH (and other sulphur-bearing molecules) account for << 1% of the gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules, however, greatly exceed those predicted by standard models of cold diffuse molecular clouds, providing further evidence for the enhancement of endothermic reaction rates by elevated temperatures or ion-neutral drift. We have considered the observed abundance ratios in the context of shock and turbulent dissipation region (TDR) models. Using the TDR model, we find that the turbulent energy available at large scale in the diffuse ISM is sufficient to explain the observed column densities of SH and CS. Standard shock and TDR models, however, fail to reproduce the column densities of H$_2$S and SO by a factor of about 10; more elaborate shock models - in which account is taken of the velocity drift, relative to H$_2$, of SH molecules produced by the dissociative recombination of H$_3$S$^+$ - reduce this discrepancy to a factor ~ 3., Comment: 30 pages, accepted for publication in A&A

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta_H) and molecular hydrogen fraction, f(H2). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042+-0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H2 fractions. We also infer zeta_H throughout our sample, and find a log-normal distribution with mean log(zeta_H)=-15.75, (zeta_H=1.78x10^-16 s^-1), and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H3+ observations. Ionization rates in the Galactic center tend to be 10--100 times larger than found in the Galactic disk, also in accord with prior studies., Comment: 76 pages, 25 figures, 6 tables; accepted for publication in ApJ

Ionized carbon is the main gas-phase reservoir of carbon in the neutral diffuse interstellar medium and its 158 micron fine structure transition [CII] is the most important cooling line of the diffuse interstellar medium (ISM). We combine [CII] absorption and emission spectroscopy to gain an improved understanding of physical conditions in the different phases of the ISM. We present high resolution [CII] spectra obtained with the Herschel/HIFI instrument towards bright dust continuum sources regions in the Galactic plane, probing simultaneously the diffuse gas along the line of sight and the background high-mass star forming regions. These data are complemented by observations of the 492 and 809 GHz fine structure lines of atomic carbon and by medium spectral resolution spectral maps of the fine structure lines of atomic oxygen at 63 and 145 microns with Herschel/PACS. We show that the presence of foreground absorption may completely cancel the emission from the background source in medium spectral resolution data and that high spectral resolution spectra are needed to interpret the [CII] and [OI] emission and the [CII]/FIR ratio. This phenomenon may explain part of the [CII]/FIR deficit seen in external luminous infrared galaxies. The C+ and C excitation in the diffuse gas is consistent with a median pressure of 5900 Kcm-3 for a mean TK ~100 K. The knowledge of the gas density allows us to determine the filling factor of the absorbing gas along the selected lines of sight: the median value is 2.4 %, in good agreement with the CNM properties. The mean excitation temperature is used to derive the average cooling due to C+ in the Galactic plane : 9.5 x 10^{-26} erg/s/H. Along the observed lines of sight, the gas phase carbon abundance does not exhibit a strong gradient as a function of Galacto-centric radius and has a weighted average of C/H = 1.5 +/- 0.4 x 10^{-4}., Comment: 25 pages

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

We present new observations of the (6,6) and (9,9) inversion transitions of the hydronium ion toward Sagittarius B2(N) and W31C. Sensitive observations toward Sagittarius B2(N) show that the high, ~ 500 K, rotational temperatures characterizing the population of the highly-excited metastable H3O+ rotational levels are present over a wide range of velocities corresponding to the Sagittarius B2 envelope, as well as the foreground gas clouds between the Sun and the source. Observations of the same lines toward W31C, a line of sight that does not intersect the Central Molecular Zone, but instead traces quiescent gas in the Galactic disk, also imply a high rotational temperature of ~ 380 K, well in excess of the kinetic temperature of the diffuse Galactic interstellar medium. While it is plausible that some fraction of the molecular gas may be heated to such high temperatures in the active environment of the Galactic center, characterized by high X-ray and cosmic ray fluxes, shocks and high degree of turbulence, this is unlikely in the largely quiescent environment of the Galactic disk clouds. We suggest instead that the highly-excited states of the hydronium ion are populated mainly by exoergic chemical formation processes and temperature describing the rotational level population does not represent the physical temperature of the medium. The same arguments may be applicable to other symmetric top rotors, such as ammonia. This offers a simple explanation to the long-standing puzzle of the presence of a pervasive, hot molecular gas component in the central region of the Milky Way. Moreover, our observations suggest that this is a universal process, not limited to the active environments associated with galactic nuclei., Comment: Accepted for publication in The Astrophysical Journal, 10 pages, 4 figures, 3 tables

Aims. We describe the assignment of a previously unidentified interstellar absorption line to ArH+ and discuss its relevance in the context of hydride absorption in diffuse gas with a low H-2 fraction. The confidence of the assignment to ArH+ is discussed, and the column densities are determined toward several lines of sight. The results are then discussed in the framework of chemical models, with the aim of explaining the observed column densities. Methods. We fitted the spectral lines with multiple velocity components, and determined column densities from the line-to-continuum ratio. The column densities of ArH+ were compared to those of other species, tracing interstellar medium (ISM) components with different H-2 abundances. We constructed chemical models that take UV radiation and cosmic ray ionization into account. Results. Thanks to the detection of two isotopologues, (ArH+)-Ar-36 and (ArH+)-Ar-38, we are confident about the carrier assignment to NeH+. is not detected with a limit of [NeH+]/[ArH+] <= 0.1. The derived column densities agree well with the predictions of chemical models. ArH+ is a unique tracer of gas with a fractional H-2 abundance of 10(-4)-10(-3) and shows little correlation to H2O+, which traces gas with a fractional H2 abundance of Conclusions. A careful analysis of variations in the ArH+, OH, H2O+, and HF column densities promises to be a faithful tracer of the distribution of the H-2 fractional abundance by providing unique information on a poorly known phase in the cycle of interstellar matter and on its transition from atomic diffuse gas to dense molecular gas traced by CO emission. Abundances of these species put strong observational constraints upon magnetohydrodynamical (MHD) simulations of the interstellar medium, and potentially could evolve into a tool characterizing the ISM. Paradoxically, the ArH+ molecule is a better tracer of almost purely atomic hydrogen gas than Hi itself. since Hi can also be present in gas with a significa