Showing 1–50 of 142 results for author: Neufeld, D A

JWST/MIRI detection of suprathermal OH rotational emissions: probing the dissociation of the water by Lyman alpha photons near the protostar HOPS 370

Authors: David A. Neufeld, P. Manoj, Himanshu Tyagi, Mayank Narang, Dan M. Watson, S. Thomas Megeath, Ewine F. Van Dishoeck, Robert A. Gutermuth, Thomas Stanke, Yao-Lun Yang, Adam E. Rubinstein, Guillem Anglada, Henrik Beuther, Alessio Caratti o Garatti, Neal J. Evans II, Samuel Federman, William J. Fischer, Joel Green, Pamela Klaassen, Leslie W. Looney, Mayra Osorio, Pooneh Nazari, John J. Tobin, Lukasz Tychoniec, Scott Wolk

Abstract: Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative str… ▽ More Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative strengths of the observed OH lines provide a powerful diagnostic of the ultraviolet radiation field in a heavily-extinguished region ($A_V \sim 10 - 20$) where direct UV observations are impossible. To high precision, the OH line strengths are consistent with a picture in which the suprathermal OH states are populated following the photodissociation of water in its $\tilde B - X$ band by ultraviolet radiation produced by fast ($\sim 80\,\rm km\,s^{-1}$) shocks along the jet. The observed dominance of emission from symmetric ($A^\prime$) OH states over that from antisymmetric ($A^{\prime\prime}$) states provides a distinctive signature of this particular population mechanism. Moreover, the variation of intensity with rotational quantum number suggests specifically that Ly$α$ radiation is responsible for the photodissociation of water, an alternative model with photodissociation by a 10$^4$ K blackbody being disfavored at a high level of significance. Using measurements of the Br$α$ flux to estimate the Ly$α$ production rate, we find that $\sim 4\%$ of the Ly$α$ photons are absorbed by water. Combined with direct measurements of water emissions in the $ν_2 = 1 -0$ band, the OH observations promise to provide key constraints on future models for the diffusion of Ly$α$ photons in the vicinity of a shock front. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 30 pages, 7 figures. Accepted for publication in ApJ Letters

JWST detection of extremely excited outflowing CO and H2O in VV 114 E SW: a possible rapidly accreting IMBH

Authors: Eduardo González-Alfonso, Ismael García-Bernete, Miguel Pereira-Santaella, David A. Neufeld, Jacqueline Fischer, Fergus R. Donnan

Abstract: Mid-infrared (mid-IR) gas-phase molecular bands are powerful diagnostics of the warm interstellar medium. We report the James Webb Space Telescope detection of the CO v=1-0 (4.4-5.0 um) and H2O nu2=1-0 (5.0-7.8um) ro-vibrational bands, both in absorption, toward the ``s2'' core in the southwest nucleus of the merging galaxy VV 114 E. All ro-vibrational CO lines up to J_low=33 (E_low~3000 K) are de… ▽ More Mid-infrared (mid-IR) gas-phase molecular bands are powerful diagnostics of the warm interstellar medium. We report the James Webb Space Telescope detection of the CO v=1-0 (4.4-5.0 um) and H2O nu2=1-0 (5.0-7.8um) ro-vibrational bands, both in absorption, toward the ``s2'' core in the southwest nucleus of the merging galaxy VV 114 E. All ro-vibrational CO lines up to J_low=33 (E_low~3000 K) are detected, as well as a forest of H2O lines up to 13_{0,13} (E_low~2600 K). The highest-excitation lines are blueshifted by ~180 km s^{-1} relative to the extended molecular cloud, which is traced by the rotational CO J=3-2 346 GHz line observed with the Atacama Large Millimeter/submillimeter Array. The bands also show absorption in a low-velocity component (blueshifted by ~30 km s^{-1}) with lower excitation. The analysis shows that the bands are observed against a continuum with effective temperature of T_bck~550 K extinguished with tau_6um^ext~ 2.5-3 (A_k~6.9-8.3 mag). The high-excitation CO and H2O lines are consistent with v=0 thermalization with T_rot~450 K and column densities of N_CO~(1.7-3.5)x10^{19} cm^{-2} and N_H2O~(1.5-3.0)x10^{19} cm$^{-2}$. Thermalization of the v=0 levels of H2O requires either an extreme density of n_H2>~10^9 cm^{-3}, or radiative excitation by the mid-IR field in a very compact (<1 pc) optically thick source emitting ~10^{10} L_sun. The latter alternative is favored, implying that the observed absorption probes the very early stages of a fully enshrouded active black hole (BH). On the basis of a simple model for BH growth and applying a lifetime constraint to the s2 core, an intermediate-mass BH (IMBH, M_BH~4.5x10^4 M_sun) accreting at super-Eddington rates is suggested, where the observed feedback has not yet been able to break through the natal cocoon. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Dissociative Recombination of Rotationally Cold OH$^+$ and Its Implications for the Cosmic Ray Ionization Rate in Diffuse Clouds

Authors: Ábel Kálosi, Lisa Gamer, Manfred Grieser, Robert von Hahn, Leonard W. Isberner, Julia I. Jäger, Holger Kreckel, David A. Neufeld, Daniel Paul, Daniel W. Savin, Stefan Schippers, Viviane C. Schmidt, Andreas Wolf, Mark G. Wolfire, Oldřich Novotný

Abstract: Observations of OH$^+$ are used to infer the interstellar cosmic ray ionization rate in diffuse atomic clouds, thereby constraining the propagation of cosmic rays through and the shielding by interstellar clouds, as well as the low energy cosmic ray spectrum. In regions where the H$_2$ to H number density ratio is low, dissociative recombination (DR) is the dominant destruction process for OH$^+$… ▽ More Observations of OH$^+$ are used to infer the interstellar cosmic ray ionization rate in diffuse atomic clouds, thereby constraining the propagation of cosmic rays through and the shielding by interstellar clouds, as well as the low energy cosmic ray spectrum. In regions where the H$_2$ to H number density ratio is low, dissociative recombination (DR) is the dominant destruction process for OH$^+$ and the DR rate coefficient is important for predicting the OH$^+$ abundance and inferring the cosmic ray ionization rate. We have experimentally studied DR of electronically and vibrationally relaxed OH$^+$ in its lowest rotational levels, using an electron--ion merged-beams setup at the Cryogenic Storage Ring. From these measurements, we have derived a kinetic temperature rate coefficient applicable to diffuse cloud chemical models, i.e., for OH$^+$ in its electronic, vibrational, and rotational ground level. At typical diffuse cloud temperatures, our kinetic temperature rate coefficient is a factor of $\sim 5$ times larger than the previous experimentally derived value and a factor of $\sim 33$ times larger than the value calculated by theory. Our combined experimental and modelling results point to a significant increase for the cosmic ray ionization rate inferred from observations of OH$^+$ and H$_2$O$^+$, corresponding to a geometric mean of $(6.6 \pm 1.0) \times 10^{-16}\,\mathrm{s}^{-1}$, which is more than a factor of two larger than the previously inferred values of the cosmic ray ionization rate in diffuse atomic clouds. Combined with observations of diffuse and dense molecular clouds, these findings indicate a greater degree of cosmic ray shielding in interstellar clouds than has been previously inferred. △ Less

Submitted 29 September, 2023; originally announced September 2023.

Comments: Main paper: PDFLaTeX with 10 pages, 4 figures. Appendix starting on page 11: PDFLaTeX with 10 pages, 3 figures, 2 tables. This is the Accepted Manuscript version of an article accepted for publication in The Astrophysical Journal Letters. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules II. The absorption line survey with the IRAM 30 m telescope

Authors: W. -J. Kim, P. Schilke, D. A. Neufeld, A. M. Jacob, Á. Sánchez-Monge, D. Seifried, B. Godard, K. M. Menten, S. Walch, E. Falgarone, V. S. Veena, S. Bialy, T. Möller, F. Wyrowski

Abstract: 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étrique (IRAM) 30 m telescope, we surveyed molecular lines in the 2 mm and 3 mm wavelength ranges toward 15 millimeter continuum source… ▽ More 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étrique (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. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 41 pages, 25 figures, 10 tables, and 5 appendices. Accepted for publication in A&A

Journal ref: A&A 670, A111 (2023)

Experimental determination of the dissociative recombination rate coefficient for rotationally-cold CH$^{+}$ and its implications for the diffuse cloud chemistry

Authors: Daniel Paul, Manfred Grieser, Florian Grussie, Robert von Hahn, Leonard W. Isberner, Ábel Kálosi, Claude Krantz, Holger Kreckel, Damian Müll, David A. Neufeld, Daniel W. Savin, Stefan Schippers, Patrick Wilhelm, Andreas Wolf, Mark G. Wolfire, Oldřich Novotný

Abstract: Observations of CH$^+$ are used to trace the physical properties of diffuse clouds, but this requires an accurate understanding of the underlying CH$^+$ chemistry. Until this work, the most uncertain reaction in that chemistry was dissociative recombination (DR) of CH$^+$. Using an electron-ion merged-beams experiment at the Cryogenic Storage Ring, we have determined the DR rate coefficient of the… ▽ More Observations of CH$^+$ are used to trace the physical properties of diffuse clouds, but this requires an accurate understanding of the underlying CH$^+$ chemistry. Until this work, the most uncertain reaction in that chemistry was dissociative recombination (DR) of CH$^+$. Using an electron-ion merged-beams experiment at the Cryogenic Storage Ring, we have determined the DR rate coefficient of the CH$^+$ electronic, vibrational, and rotational ground state applicable for different diffuse cloud conditions. Our results reduce the previously unrecognized order-of-magnitude uncertainty in the CH$^+$ DR rate coefficient to $\sim \pm 20\%$ and are applicable at all temperatures relevant to diffuse clouds, ranging from quiescent gas to gas locally heated by processes such as shocks and turbulence. Based on a simple chemical network, we find that DR can be an important destruction mechanism at temperatures relevant to quiescent gas. As the temperature increases locally, DR can continue to be important up to temperatures of $ \sim 600\,\mathrm{K} $ if there is also a corresponding increase in the electron fraction of the gas. Our new CH$^+$ DR rate coefficient data will increase the reliability of future studies of diffuse cloud physical properties via CH$^+$ abundance observations. △ Less

Submitted 31 August, 2022; originally announced August 2022.

Comments: Main paper: PDFLaTeX with 7 pages, 3 figures. Appendix starting on page 7: PDFLaTeX with 11 pages, 2 figures, 4 tables. This article has been accepted by The Astrophysical Journal

Infrared line emissions from atoms and atomic ions in NGC 7027: improved wavelength determinations for infrared metal lines and a probable detection of Zn$^{5+}$

Authors: David A. Neufeld

Abstract: An infrared L- and M-band spectral survey, performed toward the young planetary nebula NGC 7027 with the iSHELL instrument on NASA's Infrared Telescope Facility (IRTF), has revealed more than 20 vibrational lines of the molecules HeH$^+$, H$_2$, and CH$^+$ and more than 50 spectral lines of atoms and atomic ions. The present paper focuses on the atomic line emissions, the molecular lines having be… ▽ More An infrared L- and M-band spectral survey, performed toward the young planetary nebula NGC 7027 with the iSHELL instrument on NASA's Infrared Telescope Facility (IRTF), has revealed more than 20 vibrational lines of the molecules HeH$^+$, H$_2$, and CH$^+$ and more than 50 spectral lines of atoms and atomic ions. The present paper focuses on the atomic line emissions, the molecular lines having been discussed in two previous publications (Neufeld et al. 2020, 2021). The atomic lines detected with high confidence in the 2.951 - 5.24 micron region covered (incompletely) by this survey comprise (1) six collisionally-excited lines of metal ions that had previously been identified in astrophysical nebulae but for which the present observations provide the most accurate wavelength determinations obtained to date; (2) a spectral line at 4.6895 micron, not previously reported, for which the probable identification is the $^4F_{7/2} - $ $^4F_{9/2}$ fine structure transition of [Zn VI]; (3) 39 recombination lines of H and He$^+$, with upper states of principal quantum number up to 38 (H) or 24 (He$^+$); (4) 10 recombination lines of the multielectron species He, C$^{2+}$, and C$^{3+}$. △ Less

Submitted 11 May, 2022; originally announced May 2022.

Comments: 34 pages, 13 figures, accepted for publication in the Astrophysical Journal

HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules -- I. Survey description and a first look toward W3(OH), W3 IRS5 and NGC 7538 IRS1

Authors: A. M. Jacob, D. A. Neufeld, P. Schilke, H. Wiesemeyer, W. Kim, S. Bialy, M. Busch, D. Elia, E. Falgarone, M. Gerin, B. Godard, R. Higgins, P. Hennebelle, N. Indriolo, D. C. Lis, K. M. Menten, A. Sanchez-Monge, V. Ossenkopf-Okada, M. R. Rugel, D. Seifried, P. Sonnentrucker, S. Walch, M. Wolfire, F. Wyrowski, V. Valdivia

Abstract: 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… ▽ More 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. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 33 pages, 17 figures, 7 tables. Accepted for publication in the Astrophysical Journal

ArH+ and H2O+ absorption towards luminous galaxies

Authors: Arshia M. Jacob, Karl M. Menten, Friedrich Wyrowski, Benjamin Winkel, David A. Neufeld, Baerbel S. Koribalski

Abstract: Along several sight lines within the Milky Way ArH+ has been ubiquitously detected with only one detection in extragalactic environments, namely along two sight lines in the red shift z=0.89 absorber towards the lensed blazar PKS 1830-211. Being formed in predominantly atomic gas by reactions between Ar+, which were initially ionised by cosmic rays and molecular hydrogen, ArH+ has been shown to be… ▽ More Along several sight lines within the Milky Way ArH+ has been ubiquitously detected with only one detection in extragalactic environments, namely along two sight lines in the red shift z=0.89 absorber towards the lensed blazar PKS 1830-211. Being formed in predominantly atomic gas by reactions between Ar+, which were initially ionised by cosmic rays and molecular hydrogen, ArH+ has been shown to be an excellent tracer of atomic gas as well as the impinging cosmic-ray ionisation rates. In this work, we attempt to extend the observations of ArH+in extragalactic sources to examine its use as a tracer of the atomic interstellar medium (ISM) in these galaxies. We report the detection of ArH+ towards two luminous nearby galaxies, NGC 253 and NGC 4945, and the non-detection towards Arp 220 observed using the SEPIA660 receiver on the APEX 12 m telescope. In addition, the two sidebands of this receiver allowed us to observe the NKaKc=1_1,0-1_0,1 transitions of another atomic gas tracer p-H2O+ at 607.227 GHz with the ArH+ line, simultaneously. We modelled the optically thin spectra of both species and compared their observed line profiles with that of other well-known atomic gas tracers such as OH+ and o-H2O+ and diffuse and dense molecular gas tracers HF and CO, respectively. By further assuming that the observed absorption from the ArH+, OH+, and H2O+ molecules are affected by the same flux of cosmic rays, we investigate the properties of the different cloud layers based on a steady-state analysis of the chemistry of these three species. △ Less

Submitted 31 December, 2021; originally announced December 2021.

Comments: Accepted for publication in Astronomy and Astrophysics 16 pages, 11 figures, 3 tables

Journal ref: A&A 659, A152 (2022)

The chemistry of chlorine-bearing species in the diffuse interstellar medium, and new SOFIA/GREAT observations of HCl$^+$

Authors: David A. Neufeld, Helmut Wiesemeyer, Mark J. Wolfire, Arshia Jacob, Christof Buchbender, Maryvonne Gerin, Harshal Gupta, Rolf Güsten, Peter Schilke

Abstract: We have revisited the chemistry of chlorine-bearing species in the diffuse interstellar medium with new observations of the HCl$^+$ molecular ion and new astrochemical models. Using the GREAT instrument on board SOFIA, we observed the $^2Π_{3/2}\, J = 5/2 - 3/2$ transition of HCl$^+$ near 1444 GHz toward the bright THz continuum source W49N. We detected absorption by diffuse foreground gas unassoc… ▽ More We have revisited the chemistry of chlorine-bearing species in the diffuse interstellar medium with new observations of the HCl$^+$ molecular ion and new astrochemical models. Using the GREAT instrument on board SOFIA, we observed the $^2Π_{3/2}\, J = 5/2 - 3/2$ transition of HCl$^+$ near 1444 GHz toward the bright THz continuum source W49N. We detected absorption by diffuse foreground gas unassociated with the background source, and were able to thereby measure the distribution of HCl$^+$ along the sight-line. We interpreted the observational data using an updated version of an astrochemical model used previously in a theoretical study of Cl-bearing interstellar molecules. The abundance of HCl$^+$ was found to be almost constant relative to the related H$_2$Cl$^+$ ion, but the observed $n({\rm H_2Cl^+})/n({\rm HCl^+})$ abundance ratio exceeds the predictions of our astrochemical model by an order-of-magnitude. This discrepancy suggests that the rate of the primary destruction process for ${\rm H_2Cl^+}$, dissociative recombination, has been significantly overestimated. For HCl$^+$, the model predictions can provide a satisfactory fit to the observed column densities along the W49N sight-line while simultaneously accounting for the ${\rm OH^+}$ and ${\rm H_2O^+}$ column densities. △ Less

Submitted 1 June, 2021; originally announced June 2021.

Comments: 34 pages. 8 figures. Accepted for publication in the Astrophysical Journal

Anomalous intensities in the infrared emission of CH$^+$ explained by quantum nuclear motion and electric dipole calculations

Authors: P. Bryan Changala, David A. Neufeld, Benjamin Godard

Abstract: The unusual infrared emission patterns of CH$^+$, recently detected in the planetary nebula NGC 7027, are examined theoretically with high-accuracy rovibrational wavefunctions and $ab$ $initio$ dipole moment curves. The calculated transition dipole moments quantitatively reproduce the observed $J$-dependent intensity variation, which is ascribed to underlying centrifugal distortion-induced interfe… ▽ More The unusual infrared emission patterns of CH$^+$, recently detected in the planetary nebula NGC 7027, are examined theoretically with high-accuracy rovibrational wavefunctions and $ab$ $initio$ dipole moment curves. The calculated transition dipole moments quantitatively reproduce the observed $J$-dependent intensity variation, which is ascribed to underlying centrifugal distortion-induced interference effects. We discuss the implications of this anomalous behavior for astrochemical modeling of CH$^+$ production and excitation, and provide a simple expression to estimate the magnitude of this effect for other light diatomic molecules with small dipole derivatives. △ Less

Submitted 28 May, 2021; originally announced May 2021.

Comments: 21 pages, 6 figures, accepted for publication in the Astrophysical Journal

Observations and analysis of CH$^+$ vibrational emissions from the young, carbon-rich planetary nebula NGC 7027: a textbook example of chemical pumping

Authors: David A. Neufeld, Benjamin Godard, P. Bryan Changala, Alexandre Faure, T. R. Geballe, Rolf Güsten, Karl M. Menten, Helmut Wiesemeyer

Abstract: We discuss the detection of 14 rovibrational lines of CH$^+$, obtained with the iSHELL spectrograph on NASA's Infrared Telescope Facility (IRTF) on Maunakea. Our observations in the 3.49 - 4.13 $μ$m spectral region, obtained with a 0.375" slit width that provided a spectral resolving power $λ/Δλ\sim 80,000$, have resulted in the unequivocal detection of the $R(0) - R(3)$ and $P(1)-P(10)$ transitio… ▽ More We discuss the detection of 14 rovibrational lines of CH$^+$, obtained with the iSHELL spectrograph on NASA's Infrared Telescope Facility (IRTF) on Maunakea. Our observations in the 3.49 - 4.13 $μ$m spectral region, obtained with a 0.375" slit width that provided a spectral resolving power $λ/Δλ\sim 80,000$, have resulted in the unequivocal detection of the $R(0) - R(3)$ and $P(1)-P(10)$ transitions within the $v=1-0$ band of CH$^+$. The $R$-branch transitions are anomalously weak relative to the $P$-branch transitions, a behavior that is explained accurately by rovibronic calculations of the transition dipole moment reported in a companion paper (Changala et al. 2021). Nine infrared transitions of H$_2$ were also detected in these observations, comprising the $S(8)$, $S(9)$, $S(13)$ and $S(15)$ pure rotational lines; the $v=1-0$ $O(4) - O(7)$ lines, and the $v=2-1$ $O(5)$ line. We present a photodissociation region model, constrained by the CH$^+$ and H$_2$ line fluxes that we measured, that includes a detailed treatment of the excitation of CH$^+$ by inelastic collisions, optical pumping, and chemical ("formation") pumping. The latter process is found to dominate the excitation of the observed rovibrational lines of CH$^+$, and the model is remarkably successful in explaining both the absolute and relative strengths of the CH$^+$ and H$_2$ lines. △ Less

Submitted 28 May, 2021; originally announced May 2021.

Comments: 38 pages, 15 figures, accepted for publication in the Astrophysical Journal

Terahertz Water Masers: II. Further SOFIA/GREAT Detections toward Circumstellar Outflows, and a Multitransition Analysis

Authors: David A. Neufeld, Karl M. Menten, Carlos Durán, Rolf Güsten, Michael J. Kaufman, Alex Kraus, Parichay Mazumdar, Gary J. Melnick, Gisela Ortiz-León, Helmut Wiesemeyer, Friedrich Wyrowski

Abstract: Following up on our discovery of terahertz water masers, reported in 2017, we report two further detections of water maser emission at frequencies above 1 THz. Using the GREAT instrument on SOFIA, we have detected emission in the 1.296411 THz $8_{27}-7_{34}$ transition of water toward two additional oxygen-rich evolved stars, omicron Ceti (Mira) and R Crateris, and obtained an upper limit on the 1… ▽ More Following up on our discovery of terahertz water masers, reported in 2017, we report two further detections of water maser emission at frequencies above 1 THz. Using the GREAT instrument on SOFIA, we have detected emission in the 1.296411 THz $8_{27}-7_{34}$ transition of water toward two additional oxygen-rich evolved stars, omicron Ceti (Mira) and R Crateris, and obtained an upper limit on the 1.296 THz line emission from U Orionis. Toward these three sources, and toward the red supergiant star VY Canis Majorae from which 1.296 THz line emission was reported previously, we have also observed several lower-frequency (sub)millimeter water maser transitions using the APEX 12-m telescope along with the 22 GHz transition using the Effelsberg 100-m telescope. We have used a simple model to analyse the multi-transition data thereby obtained. Adopting, as a prior, independent literature estimates of the mass-loss-rates in these four sources and in W Hydrae, we infer water abundances in a remarkably narrow range: $n({\rm H_2O})/n({\rm H_2}) = 1.4 - 2.5 \times 10^{-4}$. For o Cet, VY CMa, and W Hya, the model is successful in predicting the maser line fluxes to within a typical factor $\sim 1.6 - 3$. For R Crt and U Ori, the model is less successful, with typical line flux predictions lying an order of magnitude above or below the observations; such discrepancies are perhaps unsurprising given the exponential nature of maser amplification. △ Less

Submitted 19 November, 2020; v1 submitted 3 November, 2020; originally announced November 2020.

Comments: 37 pages, including 10 figures. Accepted for publication in the Astrophysical Journal. Two typos corrected: in two instances the ratio n(H2O)/n(H2) had erroneously been inverted in version 1

Extending the view of ArH+ chemistry in diffuse clouds

Authors: Arshia M. Jacob, Karl M. Menten, Friedrich Wyrowski, Benjamin Winkel, David A. Neufeld

Abstract: One of the surprises of the Herschel mission was the detection of ArH+ towards the Crab Nebula in emission and in absorption towards strong Galactic background sources. Although these detections were limited to the first quadrant of the Galaxy, the existing data suggest that ArH+ ubiquitously and exclusively probes the diffuse atomic regions of the ISM. In this study, we extend the coverage of ArH… ▽ More One of the surprises of the Herschel mission was the detection of ArH+ towards the Crab Nebula in emission and in absorption towards strong Galactic background sources. Although these detections were limited to the first quadrant of the Galaxy, the existing data suggest that ArH+ ubiquitously and exclusively probes the diffuse atomic regions of the ISM. In this study, we extend the coverage of ArH+ to other parts of the Galaxy with new observations of its J = 1-0 transition along seven Galactic sight lines towards bright sub-mm continuum sources. We aim to benchmark its efficiency as a tracer of purely atomic gas by evaluating its correlation (or lack there of) with other well-known atomic and molecular gas tracers. The observations of ArH+ near 617.5 GHz were made feasible with the new, sensitive SEPIA660 receiver on the APEX 12 m telescope. The two sidebands of this receiver allowed us to observe p-H2O+ transitions of at 607.227 GHz simultaneously with the ArH+ line. By analysing the steady state chemistry of OH+ and o-H2O+, we derive on average a cosmic-ray ionisation rate (CRIR), of 2.3e-16 s^-1 towards the sight lines studied in this work. Using the derived values of the CRIR and the observed ArH+ abundances we constrain the molecular fraction of the gas traced by ArH+ to lie below 2e-2 with a median value of 8.8e-4. Combined, our observations of ArH+, OH+, H2O+, and CH probe different regimes of the ISM, from diffuse atomic to diffuse and translucent molecular clouds. Over Galactic scales, we see that the distribution of N(ArH+) is associated with that of N(H), particularly in the inner Galaxy with potentially even contributions from the warm neutral medium phase of atomic gas at larger galactocentric distances. We derive an average o/p-ratio for H2O+ of 2.1, which corresponds to a nuclear spin temperature of 41 K, consistent with the typical gas temperatures of diffuse clouds. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: Accepted for publication in Astronomy & Astrophysics (24 pages, 30 figures)

Journal ref: A&A 643, A91 (2020)

High Resolution Infrared Spectroscopy of Hot Molecular Gas in AFGL 2591 and AFGL 2136: Accretion in the Inner Regions of Disks Around Massive Young Stellar Objects

Authors: Andrew G. Barr, Adwin Boogert, Curtis N. DeWitt, Edward Montiel, Matthew J. Richter, John H. Lacy, David A. Neufeld, Nick Indriolo, Yvonne Pendleton, Jean Chiar, Alexander G. G. M. Tielens

Abstract: We have performed a high resolution 4-13 $μm$ spectral survey of the hot molecular gas associated with the massive protostars AFGL 2591 and AFGL 2136, utilising the Echelon-Cross-Echelle-Spectrograph (EXES) on-board the Stratospheric Observatory for Infrared Astronomy (SOFIA), and the iSHELL instrument and Texas Echelon Cross Echelle Spectrograph (TEXES) on the NASA Infrared Telescope Facility (IR… ▽ More We have performed a high resolution 4-13 $μm$ spectral survey of the hot molecular gas associated with the massive protostars AFGL 2591 and AFGL 2136, utilising the Echelon-Cross-Echelle-Spectrograph (EXES) on-board the Stratospheric Observatory for Infrared Astronomy (SOFIA), and the iSHELL instrument and Texas Echelon Cross Echelle Spectrograph (TEXES) on the NASA Infrared Telescope Facility (IRTF). Here we present results of this survey with analysis of CO, HCN, C$_2$H$_2$, NH$_3$ and CS, deriving the physical conditions for each species. Also from the IRTF, iSHELL data at 3 $μm$ for AFGL 2591 are presented that show HCN and C$_2$H$_2$ in emission. In the EXES and TEXES data, all species are detected in absorption, and temperatures and abundances are found to be high (600 K and 10$^{-6}$, respectively). Differences of up to an order of magnitude in the abundances of transitions that trace the same ground state level are measured for HCN and C$_2$H$_2$. The mid-infrared continuum is known to originate in a disk, hence we attribute the infrared absorption to arise in the photosphere of the disk. As absorption lines require an outwardly decreasing temperature gradient, we conclude that the disk is heated in the mid-plane by viscous heating due to accretion. We attribute the near-IR emission lines to scattering by molecules in the upper layers of the disk photosphere. The absorption lines trace the disk properties at 50 AU where a high temperature gas-phase chemistry is taking place. Abundances are consistent with chemical models of the inner disk of Herbig disks. △ Less

Submitted 22 July, 2020; originally announced July 2020.

The H$_2$O Spectrum of the Massive Protostar AFGL 2136 IRS 1 from 2 to 13 $μ$m at High Resolution: Probing the Circumstellar Disk

Authors: Nick Indriolo, D. A. Neufeld, A. G. Barr, A. C. A. Boogert, C. N. DeWitt, A. Karska, E. J. Montiel, M. J. Richter, A. G. G. M. Tielens

Abstract: We have observed the massive protostar AFGL 2136 IRS 1 in multiple wavelength windows in the near-to-mid-infrared at high ($\sim3$ km s$^{-1}$) spectral resolution using VLT+CRIRES, SOFIA+EXES, and Gemini North+TEXES. There is an abundance of H$_2$O absorption lines from the $ν_1$ and $ν_3$ vibrational bands at 2.7 $μ$m, from the $ν_2$ vibrational band at 6.1 $μ$m, and from pure rotational transit… ▽ More We have observed the massive protostar AFGL 2136 IRS 1 in multiple wavelength windows in the near-to-mid-infrared at high ($\sim3$ km s$^{-1}$) spectral resolution using VLT+CRIRES, SOFIA+EXES, and Gemini North+TEXES. There is an abundance of H$_2$O absorption lines from the $ν_1$ and $ν_3$ vibrational bands at 2.7 $μ$m, from the $ν_2$ vibrational band at 6.1 $μ$m, and from pure rotational transitions near 10-13 $μ$m. Analysis of state-specific column densities derived from the resolved absorption features reveals that an isothermal absorbing slab model is incapable of explaining the relative depths of different absorption features. In particular, the strongest absorption features are much weaker than expected, indicating optical depth effects resulting from the absorbing gas being well-mixed with the warm dust that serves as the "background" continuum source at all observed wavelengths. The velocity at which the strongest H$_2$O absorption occurs coincides with the velocity centroid along the minor axis of the compact disk in Keplerian rotation recently observed in H$_2$O emission with ALMA. We postulate that the warm regions of this dust disk dominate the continuum emission at near-to-mid infrared wavelengths, and that H$_2$O and several other molecules observed in absorption are probing this disk. Absorption line profiles are not symmetric, possibly indicating that the warm dust in the disk that produces the infrared continuum has a non-uniform distribution similar to the substructure observed in 1.3 mm continuum emission. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Comments: 30 pages, 12 figures, 3 tables, ApJ accepted

Journal ref: ApJ 894 (2020) 107

Distribution of Water Vapor in Molecular Clouds. II

Authors: Gary J. Melnick, Volker Tolls, Ronald L. Snell, Michael J. Kaufman, Edwin A. Bergin, Javier R. Goicoechea, Paul F. Goldsmith, Eduardo González-Alfonso, David J. Hollenbach, Dariusz C. Lis, David A. Neufeld

Abstract: The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temp… ▽ More The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temperatures greater than 20K. Where water is abundant as ice on grain surfaces, and unavailable to cool the gas, significant amounts of oxygen are removed from the gas phase, suppressing the gas-phase chemical reactions that lead to a number of oxygen-bearing species, including O2. Models of FUV-illuminated clouds predict that the gas-phase water abundance peaks in the range Av ~3 and 8mag of the cloud surface, depending on the gas density and FUV field strength. Deeper within such clouds, water is predicted to exist mainly as ice on grain surfaces. More broadly, these models are used to analyze a variety of other regions, including outflow cavities associated with young stellar objects and the surface layers of protoplanetary disks. In this paper, we report the results of observational tests of FUV-illuminated cloud models toward the Orion Molecular Ridge and Cepheus B using data obtained from the Herschel Space Observatory and the Five College Radio Astronomy Observatory. Toward Orion, 2220 spatial positions were observed along the face-on Orion Ridge in the H2O 110-101 557GHz and NH3 J,K=1,0-0,0 572GHz lines. Toward Cepheus B, two strip scans were made in the same lines across the edge-on ionization front. These new observations demonstrate that gas-phase water exists primarily within a few magnitudes of dense cloud surfaces, strengthening the conclusions of an earlier study based on a much smaller data set, and indirectly supports the prediction that water ice is quite abundant in dense clouds. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: 59 pages, including 23 figures, Accepted for publication in ApJ

Detection of vibrational emissions from the helium hydride ion (HeH$^+$) in the planetary nebula NGC 7027

Authors: David A. Neufeld, Miwa Goto, T. R. Geballe, Rolf Güsten, Karl M. Menten, Helmut Wiesemeyer

Abstract: We report the detection of emission in the v=1-0 P(1) (3.51629 micron) and P(2) (3.60776 micron) rovibrational lines of the helium hydride cation (HeH+) from the planetary nebula NGC 7027. These detections were obtained with the iSHELL spectrograph on NASA's Infrared Telescope Facility (IRTF) on Maunakea. The confirm the discovery of HeH+ reported recently by Guesten et al. (2019), who used the GR… ▽ More We report the detection of emission in the v=1-0 P(1) (3.51629 micron) and P(2) (3.60776 micron) rovibrational lines of the helium hydride cation (HeH+) from the planetary nebula NGC 7027. These detections were obtained with the iSHELL spectrograph on NASA's Infrared Telescope Facility (IRTF) on Maunakea. The confirm the discovery of HeH+ reported recently by Guesten et al. (2019), who used the GREAT instrument on the SOFIA airborne observatory to observe its pure rotational J=1-0 transition at 149.137 micron. The flux measured for the HeH+ v=1-0 P(1) line is in good agreement with our model for the formation, destruction and excitation of HeH+ in NGC 7027. The measured strength of the J=1-0 pure rotational line, however, exceeds the model prediction significantly, as does that of the v=1-0 P(2) line, by factors of 2.9 and 2.3 respectively. Possible causes of these discrepancies are discussed. Our observations of NGC 7027, covering the 3.26 - 3.93 micron spectral region, have led to the detection of more than sixty spectral lines including nine rovibrational emissions from CH+. The latter are detected for the first time in an astronomical source. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: 49 pages, including 17 figures. Accepted for publication in ApJ

The Structure of Dark Molecular Gas in the Galaxy -- II. Physical State of "CO-Dark" Gas in the Perseus Arm

Authors: Michael P. Busch, Ronald J. Allen, Philip D. Engelke, David E. Hogg, David A. Neufeld, Mark G. Wolfire

Abstract: We report the results from a new, highly sensitive ($ΔT_{mb} \sim 3 $mK) survey for thermal OH emission at 1665 and 1667 MHz over a dense, 9 x 9-pixel grid covering a $1°$ x $1°$ patch of sky in the direction of $l = 105°, b = +2.50°$ towards the Perseus spiral arm of our Galaxy. We compare our Green Bank Telescope (GBT) 1667 MHz OH results with archival CO J=1-0 observations from the Five College… ▽ More We report the results from a new, highly sensitive ($ΔT_{mb} \sim 3 $mK) survey for thermal OH emission at 1665 and 1667 MHz over a dense, 9 x 9-pixel grid covering a $1°$ x $1°$ patch of sky in the direction of $l = 105°, b = +2.50°$ towards the Perseus spiral arm of our Galaxy. We compare our Green Bank Telescope (GBT) 1667 MHz OH results with archival CO J=1-0 observations from the Five College Radio Astronomy Observatory (FCRAO) Outer Galaxy Survey within the velocity range of the Perseus Arm at these galactic coordinates. Out of the 81 statistically-independent pointings in our survey area, 86% show detectable OH emission at 1667 MHz, and 19% of them show detectable CO emission. We explore the possible physical conditions of the observed features using a set of diffuse molecular cloud models. In the context of these models, both OH and CO disappear at current sensitivity limits below an A$_{\rm v}$ of 0.2, but the CO emission does not appear until the volume density exceeds 100-200 cm$^{-3}$. These results demonstrate that a combination of low column density A$_{\rm v}$ and low volume density $n_{H}$ can explain the lack of CO emission along sight lines exhibiting OH emission. The 18-cm OH main lines, with their low critical density of $n^{*}$ $ \sim 1 $ cm$^{-3}$, are collisionally excited over a large fraction of the quiescent galactic environment and, for observations of sufficient sensitivity, provide an optically-thin radio tracer for diffuse H$_2$. △ Less

Submitted 13 August, 2019; originally announced August 2019.

Comments: Accepted to The Astrophysical Journal, 12 pages, 7 figures, 2 tables

SOFIA/EXES observations of warm H$_2$ at high spectral resolution: witnessing para-to-ortho conversion behind a molecular shock wave in HH7

Authors: David A. Neufeld, Curtis DeWitt, Pierre Lesaffre, Sylvie Cabrit, Antoine Gusdorf, Le Ngoc Tram, Matthew Richter

Abstract: Spectrally-resolved observations of three pure rotational lines of H$_2$, conducted with the EXES instrument on SOFIA toward the classic bow shock HH7, reveal systematic velocity shifts between the S(5) line of ortho-H$_2$ and the two para-H$_2$ lines [S(4) and S(6)] lying immediately above and below it on the rotational ladder. These shifts, reported here for the first time, imply that we are wit… ▽ More Spectrally-resolved observations of three pure rotational lines of H$_2$, conducted with the EXES instrument on SOFIA toward the classic bow shock HH7, reveal systematic velocity shifts between the S(5) line of ortho-H$_2$ and the two para-H$_2$ lines [S(4) and S(6)] lying immediately above and below it on the rotational ladder. These shifts, reported here for the first time, imply that we are witnessing the conversion of para-H$_2$ to ortho-H$_2$ within a shock wave driven by an outflow from a young stellar object. The observations are in good agreement with the predictions of models for non-dissociative, C-type molecular shocks. They provide a clear demonstration of the chemical changes wrought by interstellar shock waves, in this case the conversion of para-H$_2$ to ortho-H$_2$ in reactive collisions with atomic hydrogen, and provide among the most compelling evidence yet obtained for C-type shocks in which the flow velocity changes continuously. △ Less

Submitted 24 May, 2019; originally announced May 2019.

Comments: 19 pages, 4 figures. Accepted for publication in the Astrophysical Journal Letters

Dark Matter that Interacts with Baryons: Experimental Limits on the Interaction Cross-section for 27 Atomic Nuclei, and Resultant Constraints on the Particle Properties

Authors: David A. Neufeld, Daniel J. Brach-Neufeld

Abstract: To constrain the properties of dark matter (DM) that interacts with nucleons, we have conducted an experimental search for any anomalous heating of ordinary baryonic matter at 77 K. Our tabletop experiment is motivated by the possibility (discussed in a previous paper) that DM particles with masses in the $\sim 1 - 2 m_{\rm p}$ range could be captured by and concentrated within the Earth. For suit… ▽ More To constrain the properties of dark matter (DM) that interacts with nucleons, we have conducted an experimental search for any anomalous heating of ordinary baryonic matter at 77 K. Our tabletop experiment is motivated by the possibility (discussed in a previous paper) that DM particles with masses in the $\sim 1 - 2 m_{\rm p}$ range could be captured by and concentrated within the Earth. For suitable parameters, this phenomenon could lead to a substantial density $\sim 10^{14}\,\rm cm^{-3}$ of thermalized (300 K) DM particles at Earth's surface that would heat cooler baryonic matter. Our experiment involves precise differential measurements of the evaporation rate of liquid nitrogen in a storage dewar within which various materials are immersed. The results revealed no statistically-significant detections of heating in the 27 elements with molar fractions $> 10^{-5}$ in Earth's crust. For material with the average composition of Earth's crust, our measurements imply a $3 σ$ upper limit of $1.32 \times 10^{-27}\, n_{\rm 14}^{-1} (m_{\rm DM}/2\,m_{\rm p})^{-1/2}\,\rm cm^2$ on the mean cross-section for scattering with thermal HIDM at 300 K, where $10^{14}\,n_{14} \, \rm cm^{-3}$ is the particle density at Earth's surface. In combination with a lower limit on the scattering cross-section, obtained from a consideration of the heat flow through the Earth's crust (Neufeld, Farrar & McKee 2018), our experiment places an upper limit of $1.6 \times 10^{13}\,\rm cm^{-3}$ on the density of DM at the Earth's surface. This in turn, significantly constrains the properties for any DM candidate that interacts with baryons. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Comments: 23 pages, including 4 figures and 2 tables. Accepted for publication in the Astrophysical Journal

High-Resolution SOFIA/EXES Spectroscopy of $\mathrm{SO}_2$ Gas in the Massive Young Stellar Object MonR2 IRS3: Implications for the Sulfur Budget

Authors: Ryan Dungee, Adwin Boogert, Curtis N. DeWitt, Edward Montiel, Matthew J. Richter, Andrew G. Barr, Geoffrey A. Blake, Steven B. Charnley, Nick Indriolo, Agata Karska, David A. Neufeld, Rachel L. Smith, Alexander G. G. M. Tielens

Abstract: Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the $ν_3$ rovibrational band of SO$_2$ obtained with the Echelon-Cross-Echelle Spectrograp… ▽ More Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the $ν_3$ rovibrational band of SO$_2$ obtained with the Echelon-Cross-Echelle Spectrograph on the Stratospheric Observatory for Infrared Astronomy. Using local thermodynamic equilibrium models we derive physical parameters for the SO$_2$ gas in the massive YSO MonR2 IRS3. This yields a SO$_2$/$\mathrm{H}$ abundance lower limit of $5.6\pm0.5\times10^{-7}$, or $>\!4\%$ of the cosmic sulfur budget, and an intrinsic line width (Doppler parameter) of $b<3.20\;\mathrm{km\;s}^{-1}$. The small line widths and high temperature ($T_\mathrm{ex}=234\pm15\;\mathrm{K}$) locate the gas in a relatively quiescent region near the YSO, presumably in the hot core where ices have evaporated. This sublimation unlocks a volatile sulfur reservoir (e.g., sulfur allotropes as detected abundantly in comet 67P/Churyumov--Gerasimenko), which is followed by SO$_2$ formation by warm, dense gas-phase chemistry. The narrowness of the lines makes formation of SO$_2$ from sulfur sputtered off grains in shocks less likely toward MonR2 IRS3. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Comments: 9 pages, 2 figures, accepted for publication in ApJL

Infrared Detection of Abundant CS in the Hot Core AFGL 2591 at High Spectral Resolution with SOFIA/EXES

Authors: Andrew G. Barr, Adwin Boogert, Curtis N. DeWitt, Edward Montiel, Matthew J. Richter, Nick Indriolo, David A. Neufeld, Yvonne Pendleton, Jean Chiar, Ryan Dungee, Alexander G. G. M. Tielens

Abstract: We have performed a 5-8 $μ$m spectral line survey of the hot molecular core associated with the massive protostar AFGL 2591, using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We have supplemented these data with a ground based study in the atmospheric M band around 4.5 $μ$m using the iSHELL instrument on the Infrared Telescope Faci… ▽ More We have performed a 5-8 $μ$m spectral line survey of the hot molecular core associated with the massive protostar AFGL 2591, using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We have supplemented these data with a ground based study in the atmospheric M band around 4.5 $μ$m using the iSHELL instrument on the Infrared Telescope Facility (IRTF), and the full N band window from 8-13 $μ$m using the Texas Echelon Cross Echelle Spectrograph (TEXES) on the IRTF. Here we present the first detection of ro-vibrational transitions of CS in this source. The absorption lines are centred on average around -10 kms$^{-1}$ and the line widths of CS compare well with the hot component of $^{13}$CO (around 10 kms$^{-1}$). Temperatures for CS, hot $^{13}$CO and $^{12}$CO v=1-2 agree well and are around 700 K. We derive a CS abundance of 8$\times$10$^{-3}$ and 2$\times$10$^{-6}$ with respect to CO and H$_2$ respectively. This enhanced CS abundance with respect to the surrounding cloud (1$\times$10$^{-8}$) may reflect sublimation of H$_2$S ice followed by gas-phase reactions to form CS. Transitions are in LTE and we derive a density of $>$10$^7$ cm$^{-3}$, which corresponds to an absorbing region of $<$0.04$''$. EXES observations of CS are likely to probe deeply into the hot core, to the base of the outflow. Submillimeter and infrared observations trace different components of the hot core as revealed by the difference in systemic velocities, line widths and temperatures, as well as the CS abundance. △ Less

Submitted 25 October, 2018; originally announced October 2018.

High spectral resolution observations toward Orion BN at 6 $μ$m: no evidence for hot water

Authors: Nick Indriolo, Jonathan C. Tan, A. C. A. Boogert, C. N. DeWitt, E. J. Montiel, D. A. Neufeld, M. J. Richter

Abstract: Orion BN has a large proper motion and radial velocity with respect to the gas and other stars in the region where it is presumed to have formed. Multiple dynamical interaction scenarios have been proposed to explain this motion. In one case BN is thought to have interacted with stars in the Trapezium cluster, while in another it is thought to have interacted with source I while deeply embedded in… ▽ More Orion BN has a large proper motion and radial velocity with respect to the gas and other stars in the region where it is presumed to have formed. Multiple dynamical interaction scenarios have been proposed to explain this motion. In one case BN is thought to have interacted with stars in the Trapezium cluster, while in another it is thought to have interacted with source I while deeply embedded in molecular gas. If there is dense gas that has been retained in close proximity to BN, it may be evidence that the latter scenario is favored. We observed BN at high spectral resolution in three windows near 6 $μ$m using SOFIA/EXES targeting the $ν_2$ vibrational band of H$_2$O. Absorption from only three transitions of H$_2$O is detected, and through kinematic analysis is associated with cool, dense foreground gas, not BN itself. We find no evidence for H$_2$O absorption or emission at the systemic velocity of BN. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Comments: 8 pages, 3 figures, 1 table, accepted for publication in ApJL SOFIA special issue

Journal ref: ApJL 865 (2018) L18

Constraints on the Cosmic-Ray Ionization Rate in the $z\sim2.3$ Lensed Galaxies SMM J2135$-$0102 and SDP 17b from Observations of OH$^+$ and H$_2$O$^+$

Authors: Nick Indriolo, E. A. Bergin, E. Falgarone, B. Godard, M. A. Zwaan, D. A. Neufeld, M. G. Wolfire

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, $ζ_{\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… ▽ More 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, $ζ_{\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 $ζ_{\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 $ζ_{\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 $ζ_{\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. △ Less

Submitted 14 August, 2018; originally announced August 2018.

Comments: 16 pages, 5 figures, 3 tables, accepted for publication in ApJ

Journal ref: ApJ 865 (2018) 127

Dark Matter that Interacts with Baryons: Density Distribution within the Earth and New Constraints on the Interaction Cross-section

Authors: David A. Neufeld, Glennys R. Farrar, Christopher F. McKee

Abstract: For dark matter (DM) particles with masses in the 0.6 - 6 m_p range, we set stringent constraints on the interaction cross-sections for scattering with ordinary baryonic matter. These constraints follow from the recognition that such particles can be captured by - and thermalized within - the Earth, leading to a substantial accumulation and concentration of DM that interact with baryons. Here, we… ▽ More For dark matter (DM) particles with masses in the 0.6 - 6 m_p range, we set stringent constraints on the interaction cross-sections for scattering with ordinary baryonic matter. These constraints follow from the recognition that such particles can be captured by - and thermalized within - the Earth, leading to a substantial accumulation and concentration of DM that interact with baryons. Here, we discuss the probability that DM intercepted by the Earth will be captured, the number of DM particles thereby accumulated over Earth's lifetime, the fraction of such particles retained in the face of evaporation, and the density distribution of such particles within the Earth. In the latter context, we note that a previous treatment of the density distribution of DM, presented by Gould and Raffelt and applied subsequently to DM in the Sun, is inconsistent with considerations of hydrostatic equilibrium. Our analysis provides an estimate of the DM particle density at Earth's surface, which may exceed 1.E+14 cm-3 for the mass range under consideration. Based upon our determination of the DM density at Earth's surface, we derive constraints on the scattering cross-sections. These constraints are placed by four considerations: (1) the lifetime of the relativistic proton beam at the Large Hadron collider (LHC); (2) the orbital decay of spacecraft in low Earth orbit (LEO); (3) the vaporization rate of cryogenic liquids in well-insulated storage dewars; and (4) the thermal conductivity of Earth's crust. As an example application of our results, we show that for the scattering cross-sections that were invoked recently in Barkana's original explanation for the anomalously deep 21 cm absorption reported by EDGES, DM particle masses in the 0.6 - 4 m_p range are ruled out. △ Less

Submitted 30 July, 2018; v1 submitted 22 May, 2018; originally announced May 2018.

Comments: Accepted for publication in ApJ. This version has several minor changes to improve clarity. A new appendix had been added to justify our treatment of the dark matter density distribution within the Earth, where we note that a previous analysis of the density distribution of DM in the Sun by Gould and Raffelt is inconsistent with considerations of hydrostatic equilibrium

SOFIA/GREAT Discovery of Terahertz Water Masers

Authors: David A. Neufeld, Gary J. Melnick, Michael J. Kaufman, Helmut Wiesemeyer, Rolf Güsten, Alex Kraus, Karl M. Menten, Oliver Ricken, Alexandre Faure

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

Submitted 26 May, 2017; originally announced May 2017.

Comments: Accepted for oublication in ApJ

The cosmic ray ionization rate in the Galactic disk, as determined from observations of molecular ions

Authors: David A. Neufeld, Mark G. Wolfire

Abstract: We have obtained estimates for the cosmic-ray ionization rate (CRIR) in the Galactic disk, using a detailed model for the physics and chemistry of diffuse interstellar gas clouds to interpret previously-published measurements of the abundance of four molecular ions: ArH$^+$, OH$^+$, H$_2$O$^+$ and H$_3^+$. For diffuse $atomic$ clouds at Galactocentric distances in the range $R_g \sim 4 - 9$ kpc, o… ▽ More We have obtained estimates for the cosmic-ray ionization rate (CRIR) in the Galactic disk, using a detailed model for the physics and chemistry of diffuse interstellar gas clouds to interpret previously-published measurements of the abundance of four molecular ions: ArH$^+$, OH$^+$, H$_2$O$^+$ and H$_3^+$. For diffuse $atomic$ clouds at Galactocentric distances in the range $R_g \sim 4 - 9$ kpc, observations of ArH$^+$, OH$^+$, and H$_2$O$^+$ imply a mean primary CRIR of $(2.2 \pm 0.3) \exp [(R_0-R_g)/4.7\,\rm{kpc}] \times 10^{-16} \rm \, s^{-1}$ per hydrogen atom, where $R_0=8.5$ kpc. Within diffuse $molecular$ clouds observed toward stars in the solar neighborhood, measurements of H$_3^+$ and H$_2$ imply a primary CRIR of $(2.3 \pm 0.6) \times 10^{-16}\,\,\rm s^{-1}$ per H atom, corresponding to a total ionization rate per H$_2$ molecule of $(5.3 \pm 1.1) \times 10^{-16}\,\,\rm s^{-1},$ in good accord with previous estimates. These estimates are also in good agreement with a rederivation, presented here, of the CRIR implied by recent observations of carbon and hydrogen radio recombination lines along the sight-line to Cas A. Here, our best-fit estimate for the primary CRIR is $2.9 \times 10^{-16}\,\,\rm s^{-1}$ per H atom. Our results show marginal evidence that the CRIR in diffuse molecular clouds decreases with cloud extinction, $A_{\rm V}({\rm tot})$, with a best-fit dependence $\propto A_{\rm V}({\rm tot})^{-1}$ for $A_{\rm V}({\rm tot}) \ge 0.5$. △ Less

Submitted 12 April, 2017; originally announced April 2017.

Comments: 48 pages, including 3 tables (at end) and 13 figures. Accepted for publication in the Astrophysical Journal

The chemistry of interstellar argonium and other probes of the molecular fraction in diffuse clouds

Authors: David A. Neufeld, Mark G. Wolfire

Abstract: We present a general parameter study, in which the abundance of interstellar argonium (ArH$^+$) is predicted using a model for the physics and chemistry of diffuse interstellar gas clouds. Results have been obtained as a function of UV radiation field, cosmic-ray ionization rate, and cloud extinction. No single set of cloud parameters provides an acceptable fit to the typical ArH$^+$, OH$^+$ and… ▽ More We present a general parameter study, in which the abundance of interstellar argonium (ArH$^+$) is predicted using a model for the physics and chemistry of diffuse interstellar gas clouds. Results have been obtained as a function of UV radiation field, cosmic-ray ionization rate, and cloud extinction. No single set of cloud parameters provides an acceptable fit to the typical ArH$^+$, OH$^+$ and $\rm H_2O^+$ abundances observed in diffuse clouds within the Galactic disk. Instead, the observed abundances suggest that ArH$^+$ resides primarily in a separate population of small clouds of total visual extinction of at most 0.02 mag per cloud, within which the column-averaged molecular fraction is in the range $10^{-5} - 10^{-2}$, while OH$^+$ and $\rm H_2O^+$ reside primarily in somewhat larger clouds with a column-averaged molecular fraction $\sim 0.2$. This analysis confirms our previous suggestion that the argonium molecular ion is a unique tracer of almost purely atomic gas. △ Less

Submitted 1 July, 2016; originally announced July 2016.

Comments: Accepted for publication in ApJ. 26 pages, 8 figures

Herschel/HIFI observations of the circumstellar ammonia lines in IRC+10216

Authors: M. R. Schmidt, J. H. He, R. Szczerba, V. Bujarrabal, J. Alcolea, J. Cernicharo, L. Decin, K. Justtanont, D. Teyssier, K. M. Menten, D. A. Neufeld, H. Olofsson, P. Planesas, A. P. Marston, A. M. Sobolev, A. de Koter, F. L. Schöier

Abstract: New high-resolution far-infrared (FIR) observations of both ortho- and para-NH3 transitions toward IRC+10216 were obtained with Herschel, with the goal of determining the ammonia abundance and constraining the distribution of NH3 in the envelope of IRC+10216. We used the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel to observe all rotational transitions up to the J=3 level (t… ▽ More New high-resolution far-infrared (FIR) observations of both ortho- and para-NH3 transitions toward IRC+10216 were obtained with Herschel, with the goal of determining the ammonia abundance and constraining the distribution of NH3 in the envelope of IRC+10216. We used the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel to observe all rotational transitions up to the J=3 level (three ortho- and six para-NH3 lines). We conducted non-LTE multilevel radiative transfer modelling, including the effects of near-infrared (NIR) radiative pumping through vibrational transitions. We found that NIR pumping is of key importance for understanding the excitation of rotational levels of NH3. The derived NH3 abundances relative to molecular hydrogen were (2.8+-0.5)x10^{-8} for ortho-NH3 and (3.2^{+0.7}_{-0.6})x10^{-8} for para-NH3, consistent with an ortho/para ratio of 1. These values are in a rough agreement with abundances derived from the inversion transitions, as well as with the total abundance of NH3 inferred from the MIR absorption lines. To explain the observed rotational transitions, ammonia must be formed near to the central star at a radius close to the end of the wind acceleration region, but no larger than about 20 stellar radii (1 sigma confidence level). △ Less

Submitted 6 June, 2016; originally announced June 2016.

Comments: 13 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 592, A131 (2016)

Herschel/HIFI Spectral Mapping of C$^+$, CH$^+$, and CH in Orion BN/KL: The Prevailing Role of Ultraviolet Irradiation in CH$^+$ Formation

Authors: Patrick W. Morris, Harshal Gupta, Zsofia Nagy, John C. Pearson, Volker Ossenkopf-Okada, Edith Falgarone, Dariusz C. Lis, Maryvonne Gerin, Gary Melnick, David A. Neufeld, Edwin A. Bergin

Abstract: The CH$^+$ ion is a key species in the initial steps of interstellar carbon chemistry. Its formation in diverse environments where it is observed is not well understood, however, because the main production pathway is so endothermic (4280 K) that it is unlikely to proceed at the typical temperatures of molecular clouds. We investigation CH$^+$ formation with the first velocity-resolved spectral ma… ▽ More The CH$^+$ ion is a key species in the initial steps of interstellar carbon chemistry. Its formation in diverse environments where it is observed is not well understood, however, because the main production pathway is so endothermic (4280 K) that it is unlikely to proceed at the typical temperatures of molecular clouds. We investigation CH$^+$ formation with the first velocity-resolved spectral mapping of the CH$^+$ $J=1-0, 2-1$ rotational transitions, three sets of CH $Λ$-doubled triplet lines, $^{12}$C$^+$ and $^{13}$C$^+$, and CH$_3$OH 835~GHz E-symmetry Q branch transitions, obtained with Herschel/HIFI over $\approx$12 arcmin$^2$ centered on the Orion BN/KL source. We present the spatial morphologies and kinematics, cloud boundary conditions, excitation temperatures, column densities, and $^{12}$C$^+$ optical depths. Emission from C$^+$, CH$^+$, and CH is indicated to arise in the diluted gas, outside of the explosive, dense BN/KL outflow. Our models show that UV-irradiation provides favorable conditions for steady-state production of CH$^+$ in this environment. Surprisingly, no spatial or kinematic correspondences of these species are found with H$_2$ S(1) emission tracing shocked gas in the outflow. We propose that C$^+$ is being consumed by rapid production of CO to explain the lack of C$^+$ and CH$^+$ in the outflow, and that fluorescence provides the reservoir of H$_2$ excited to higher ro-vibrational and rotational levels. Hence, in star-forming environments containing sources of shocks and strong UV radiation, a description of CH$^+$ formation and excitation conditions is incomplete without including the important --- possibly dominant --- role of UV irradiation. △ Less

Submitted 30 April, 2016; v1 submitted 19 April, 2016; originally announced April 2016.

Comments: Accepted for publication in The Astrophysical Journal

Far-infrared study of tracers of oxygen chemistry in diffuse clouds

Authors: H. Wiesemeyer, R. Güsten, S. Heyminck, H. W. Hübers, K. M. Menten, D. A. Neufeld, H. Richter, R. Simon, J. Stutzki, B. Winkel, F. Wyrowski

Abstract: Context. The chemistry of the diffuse interstellar medium rests upon three pillars: exothermic ion-neutral reactions (" cold chemistry "), endothermic neutral-neutral reactions with significant activation barriers (" warm chemistry "), and reactions on the surfaces of dust grains. While warm chemistry becomes important in the shocks associated with turbulent dissipation regions, the main path for… ▽ More Context. The chemistry of the diffuse interstellar medium rests upon three pillars: exothermic ion-neutral reactions (" cold chemistry "), endothermic neutral-neutral reactions with significant activation barriers (" warm chemistry "), and reactions on the surfaces of dust grains. While warm chemistry becomes important in the shocks associated with turbulent dissipation regions, the main path for the formation of interstellar OH and H2O is that of cold chemistry. Aims. The aim of this study is to observationally confirm the association of atomic oxygen with both atomic and molecular gas phases, and to understand the measured abundances of OH and OH + as a function of the available reservoir of H2. Methods. We obtained absorption spectra of the ground states of OH, OH+ and OI with high-velocity resolution, with GREAT on-board SOFIA, and with the THz receiver at the APEX. We analyzed them along with ancillary spectra of HF and CH from HIFI. To deconvolve them from the hyperfine structure and to separate the blend that is due to various velocity components on the sightline, we fit model spectra consisting of an appropriate number of Gaussian profiles using a method combining simulated annealing with downhill simplex minimization. Together with HF and/or CH as a surrogate for H2, and HI $λ$21 cm data, the molecular hydrogen fraction f^N\_H2 = N(H 2)/(N(H) + 2N(H 2)) can be determined. We then investigated abundance ratios as a function of f^N\_H2. Results. The column density of OI is correlated at a high significance with the amount of available molecular and atomic hydrogen, with an atomic oxygen abundance of $3 \times 10 ^{-4}$ relative to H nuclei. While the velocities of the absorption features of OH and OH+ are loosely correlated and reflect the spiral arm crossings on the sightline, upon closer inspection they display an anticorrespondence. The arm-to-interarm density contrast is found to be higher in OH than in OH+. While both species can coexist, with a higher abundance in OH than in OH+, the latter is found less frequently in absence of OH than the other way around, which is a direct consequence of the rapid destruction of OH+ by dissociative recombination when not enough H2 is available. This conjecture has been substantiated by a comparison between the OH/OH+ ratio with f^N\_H2, showing a clear correlation. The hydrogen abstraction reaction chain OH+ (H2,H) H2O+ (H2,H)H3O+ is confirmed as the pathway for the production of OH and H 2 O. Our estimate of the branching ratio of the dissociative recombination of H3O+ to OH and H2O is confined within the interval of 84 to 91%, which matches laboratory measurements (74 to 83%). -- A correlation between the linewidths and column densities of OH+ features is found to be significant with a false-alarm probability below 5%. Such a correlation is predicted by models of interstellar MHD turbulence. For OH the same correlation is found to be insignificant because there are more narrow absorption features. Conclusions. While it is difficult to assess the contributions of warm neutral-neutral chemistry to the observed abundances, it seems fair to conclude that the predictions of cold ion-neutral chemistry match the abundance patterns we observed. △ Less

Submitted 10 February, 2016; originally announced February 2016.

Journal ref: Astronomy and Astrophysics, EDP Sciences, 2016, 585, pp.A76

Interstellar Hydrides

Authors: Maryvonne Gerin, David A. Neufeld, Javier R. Goicoechea

Abstract: Interstellar hydrides -- that is, molecules containing a single heavy element atom with one or more hydrogen atoms -- were among the first molecules detected outside the solar system. They lie at the root of interstellar chemistry, being among the first species to form in initially-atomic gas, along with molecular hydrogen and its associated ions. Because the chemical pathways leading to the forma… ▽ More Interstellar hydrides -- that is, molecules containing a single heavy element atom with one or more hydrogen atoms -- were among the first molecules detected outside the solar system. They lie at the root of interstellar chemistry, being among the first species to form in initially-atomic gas, along with molecular hydrogen and its associated ions. Because the chemical pathways leading to the formation of interstellar hydrides are relatively simple, the analysis of the observed abundances is relatively straightforward and provides key information about the environments where hydrides are found. Recent years have seen rapid progress in our understanding of interstellar hydrides, thanks largely to far-IR and submillimeter observations performed with the Herschel Space Observatory. In this review, we will discuss observations of interstellar hydrides, along with the advanced modeling approaches that have been used to interpret them, and the unique information that has thereby been obtained. △ Less

Submitted 12 January, 2016; originally announced January 2016.

Comments: Accepted for publication in Annual Review of Astronomy and Astrophysics 2016, Vol. 54

Evolution of Mass Outflow in Protostars

Authors: Dan M. Watson, Nuria P. Calvet, William J. Fischer, W. J. Forrest, P. Manoj, S. Thomas Megeath, Gary J. Melnick, Joan Najita, David A. Neufeld, Patrick D. Sheehan, Amelia M. Stutz, John J. Tobin

Abstract: We have surveyed 84 Class 0, Class I, and flat-spectrum protostars in mid-infrared [Si II], [Fe II] and [S I] line emission, and 11 of these in far-infrared [O I] emission. We use the results to derive their mass outflow rates. Thereby we observe a strong correlation of mass outflow rates with bolometric luminosity, and with the inferred mass accretion rates of the central objects, which continues… ▽ More We have surveyed 84 Class 0, Class I, and flat-spectrum protostars in mid-infrared [Si II], [Fe II] and [S I] line emission, and 11 of these in far-infrared [O I] emission. We use the results to derive their mass outflow rates. Thereby we observe a strong correlation of mass outflow rates with bolometric luminosity, and with the inferred mass accretion rates of the central objects, which continues through the Class 0 range the trend observed in Class II young stellar objects. Along this trend from large to small mass-flow rates, the different classes of young stellar objects lie in the sequence Class 0 -- Class I/flat-spectrum -- Class II, indicating that the trend is an evolutionary sequence in which mass outflow and accretion rates decrease together with increasing age, while maintaining rough proportionality. The survey results include two which are key tests of magnetocentrifugal outflow-acceleration mechanisms: the distribution of the outflow/accretion branching ratio b, and limits on the distribution of outflow speeds. Neither rule out any of the three leading outflow-acceleration, angular-momentum-ejection mechanisms, but they provide some evidence that disk winds and accretion-powered stellar winds (APSWs) operate in many protostars. An upper edge observed in the branching-ratio distribution is consistent with the upper bound of b = 0.6 found in models of APSWs, and a large fraction (0.31) of the sample have branching ratio sufficiently small that only disk winds, launched on scales as large as several AU, have been demonstrated to account for them. △ Less

Submitted 15 November, 2015; originally announced November 2015.

Comments: Version submitted to ApJ: 36 pages, 3 tables, 8 figures

Detection of extragalactic argonium, ArH$^+$, toward PKS~1830$-$211

Authors: Holger S. P. Müller, Sébastien Muller, Peter Schilke, Edwin A. Bergin, John H. Black, Maryvonne Gerin, Dariusz C. Lis, David A. Neufeld, Sümeyye Suri

Abstract: Argonium has recently been detected as a ubiquitous molecule in our Galaxy. Model calculations indicate that its abundance peaks at molecular fractions in the range of 1E-4 to 1E-3 and that the observed column densities require high values of the cosmic ray ionization rate. Therefore, this molecular cation may serve as an excellent tracer of the very diffuse interstellar medium (ISM), as well as a… ▽ More Argonium has recently been detected as a ubiquitous molecule in our Galaxy. Model calculations indicate that its abundance peaks at molecular fractions in the range of 1E-4 to 1E-3 and that the observed column densities require high values of the cosmic ray ionization rate. Therefore, this molecular cation may serve as an excellent tracer of the very diffuse interstellar medium (ISM), as well as an indicator of the cosmic ray ionization rate. We attempted to detect ArH+ in extragalactic sources to evaluate its diagnostic power as a tracer of the almost purely atomic ISM in distant galaxies. We obtained ALMA observations of a foreground galaxy at z = 0.89 in the direction of the lensed blazar PKS 1830-211. Two isotopologs of argonium, 36ArH+ and 38ArH+, were detected in absorption along two different lines of sight toward PKS 1830-211, known as the SW and NE images of the background blazar. The argonium absorption is clearly enhanced on the more diffuse line of sight (NE) compared to other molecular species. The isotopic ratio 36Ar/38Ar is 3.46 +- 0.16 toward the SW image, i.e., significantly lower than the solar value of 5.5. Our results demonstrate the suitability of argonium as a tracer of the almost purely atomic, diffuse ISM in high-redshift sources. The evolution of the isotopic ratio with redshift may help to constrain nucleosynthetic scenarios in the early Universe. △ Less

Submitted 23 September, 2015; originally announced September 2015.

Comments: Astronomy and Astrophysics, accepted; 5 pages, 3 figures, 1 table

Journal ref: A&A 582, L4 (2015)

Widespread Galactic CF+ absorption: detection toward W49 with the Plateau de Bure Interferometer

Authors: H. S. Liszt, V. V. Guzman, J. Pety, M. Gerin, D. A. Neufeld, P. Gratier

Abstract: To study the usefulness of \CFP\ as a tracer of the regions where C\p\ and \HH\ coexist in the interstellar medium. We used the Plateau de Bure Interferometer to synthesize \CFP\ J=1-0 absorption at 102.6 GHz toward the core of the distant HII region W49N at l = 43.2\degr, b=0.0\degr, and we modeled the fluorine chemistry in diffuse/translucent molecular gas. We detected \CFP\ absorption over a br… ▽ More To study the usefulness of \CFP\ as a tracer of the regions where C\p\ and \HH\ coexist in the interstellar medium. We used the Plateau de Bure Interferometer to synthesize \CFP\ J=1-0 absorption at 102.6 GHz toward the core of the distant HII region W49N at l = 43.2\degr, b=0.0\degr, and we modeled the fluorine chemistry in diffuse/translucent molecular gas. We detected \CFP\ absorption over a broad range of velocity showing that \CFP\ is widespread in the \HH-bearing Galactic disk gas. Originally detected in dense gas in the Orion Bar and Horsehead PDR, \CFP\ was subsequently detected in absorption from diffuse and translucent clouds seen toward \bll\ and 3C111. Here we showed that \CFP\ is distributed throughout the diffuse and translucent molecular disk gas with N(\CFP)/N(\HH) $= 1.5-2.0\times10^{-10}$, increasing to N(\CFP)/N(\HH) $= 3.5\times10^{-10}$ in one cloud at 39 \kms\ having higher N(\HH) $\approx 3\times10^{21}\pcc$. Models of the fluorine chemistry reproduce the observed column densities and relative abundance of HF, from which \CFP\ forms, but generally overpredict the the column density of \CFP\ by factors of 1.4-4. We show that a free space photodissociation rate $Γ\ga 10^{-9}\ps$, comparable to that of CH, might account for much of the discrepancy but a recent calculation finds a value about ten times smaller. In the heavily blended and kinematically complex spectra seen toward W49, \CFP\ absorption primarily traces the peaks of the \HH\ distribution. △ Less

Submitted 29 May, 2015; originally announced May 2015.

Comments: A&A in press

Journal ref: A&A 579, A12 (2015)

A Herschel/HIFI Legacy Survey of HF and H2O in the Galaxy: Probing Diffuse Molecular Cloud Chemistry

Authors: P. Sonnentrucker, M. Wolfire, D. A. Neufeld, N. Flagey, M. Gerin, P. Goldsmith, D. Lis, R. Monje

Abstract: We combine Herschel observations of a total of 12 sources to construct the most uniform survey of HF and H2O in our Galactic disk. Both molecules are detected in absorption along all sight lines. The high spectral resolution of the Heterodyne Instrument for the Far-Infrared (HIFI) allows us to compare the HF and H2O distributions in 47 diffuse cloud components sampling the disk. We find that the H… ▽ More We combine Herschel observations of a total of 12 sources to construct the most uniform survey of HF and H2O in our Galactic disk. Both molecules are detected in absorption along all sight lines. The high spectral resolution of the Heterodyne Instrument for the Far-Infrared (HIFI) allows us to compare the HF and H2O distributions in 47 diffuse cloud components sampling the disk. We find that the HF and H2O velocity distributions follow each other almost perfectly and establish that HF and H2O probe the same gas-phase volume. Our observations corroborate theoretical predictions that HF is a sensitive tracer of H2 in diffuse clouds, down to molecular fractions of only a few percent. Using HF to trace H2 in our sample, we find that the N(H2O)-to-N(HF) ratio shows a narrow distribution with a median value of 1.51. Our results further suggest that H2O might be used as a tracer of H2 -within a factor 2.5- in the diffuse interstellar medium. We show that the measured factor of ~2.5 variation around the median is driven by true local variations in the H2O abundance relative to H2 throughout the disk. The latter variability allows us to test our theoretical understanding of the chemistry of oxygen-bearing molecules in the diffuse gas. We show that both gas-phase and grain-surface chemistry are required to reproduce our H2O observations. This survey thus confirms that grain surface reactions can play a significant role in the chemistry occurring in the diffuse interstellar medium n_H < 1000 cm^-3. △ Less

Submitted 27 May, 2015; originally announced May 2015.

Comments: 53 pages; 12 figures, accepted for publication in ApJ main journal

Herschel observations of interstellar chloronium. II - Detections toward G29.96-0.02, W49N, W51, and W3(OH), and determinations of the ortho-to-para and $^{35}$Cl/$^{37}$Cl isotopic ratios

Authors: David A. Neufeld, John H. Black, Maryvonne Gerin, Javier R. Goicoechea, Paul F. Goldsmith, Cecile Gry, Harshal Gupta, Eric Herbst, Nick Indriolo, Dariusz Lis, Karl M. Menten, Raquel Monje, Bhaswati Mookerjea, Carina Persson, Paule Sonnentrucker, Mark G. Wolfire

Abstract: We report additional detections of the chloronium molecular ion, H$_2$Cl$^+$, toward four bright submillimeter continuum sources: G29.96, W49N, W51, and W3(OH). With the use of the HIFI instrument on the Herschel Space Observatory, we observed the $2_{12}-1_{01}$ transition of ortho-H$_2^{35}$Cl$^+$ at 781.627 GHz in absorption toward all four sources. Much of the detected absorption arises in dif… ▽ More We report additional detections of the chloronium molecular ion, H$_2$Cl$^+$, toward four bright submillimeter continuum sources: G29.96, W49N, W51, and W3(OH). With the use of the HIFI instrument on the Herschel Space Observatory, we observed the $2_{12}-1_{01}$ transition of ortho-H$_2^{35}$Cl$^+$ at 781.627 GHz in absorption toward all four sources. Much of the detected absorption arises in diffuse foreground clouds that are unassociated with the background continuum sources and in which our best estimates of the $N({\rm H_2Cl^+})/N({\rm H})$ ratio lie in the range $(0.9 - 4.8) \times 10^{-9}$. These chloronium abundances relative to atomic hydrogen can exceed the predictions of current astrochemical models by up to a factor of 5. Toward W49N, we have also detected the $2_{12}-1_{01}$ transition of ortho-H$_2^{37}$Cl$^+$ at 780.053 GHz and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ at 485.418 GHz. These observations imply $\rm H_2^{35}Cl^+/H_2^{37}Cl^+$ column density ratios that are consistent with the solar system $^{35}$Cl/$^{37}$Cl isotopic ratio of 3.1, and chloronium ortho-to-para ratios consistent with 3, the ratio of spin statistical weights. △ Less

Submitted 4 May, 2015; originally announced May 2015.

Comments: 31 pages, including 7 figures. Accepted for publication in the ApJ

SOFIA/EXES Observations of Water Absorption in the Protostar AFGL 2591 at High Spectral Resolution

Authors: Nick Indriolo, D. A. Neufeld, C. N. DeWitt, M. J. Richter, A. C. A. Boogert, G. M. Harper, D. T. Jaffe, K. R. Kulas, M. E. McKelvey, N. Ryde, W. Vacca

Abstract: We present high spectral resolution (~3 km/s) observations of the nu_2 ro-vibrational band of H2O in the 6.086--6.135 micron range toward the massive protostar AFGL 2591 using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). Ten absorption features are detected in total, with seven caused by transitions in the nu_2 band of H2O, two by t… ▽ More We present high spectral resolution (~3 km/s) observations of the nu_2 ro-vibrational band of H2O in the 6.086--6.135 micron range toward the massive protostar AFGL 2591 using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). Ten absorption features are detected in total, with seven caused by transitions in the nu_2 band of H2O, two by transitions in the first vibrationally excited nu_2 band of H2O, and one by a transition in the nu_2 band of H2{18}O. Among the detected transitions is the nu_2 1(1,1)--0(0,0) line which probes the lowest lying rotational level of para-H2O. The stronger transitions appear to be optically thick, but reach maximum absorption at a depth of about 25%, suggesting that the background source is only partially covered by the absorbing gas, or that the absorption arises within the 6 micron emitting photosphere. Assuming a covering fraction of 25%, the H2O column density and rotational temperature that best fit the observed absorption lines are N(H2O)=(1.3+-0.3)*10^{19} cm^{-2} and T=640+-80 K. △ Less

Submitted 23 February, 2015; originally announced February 2015.

Comments: 6 pages, 3 figures, 1 table, accepted for publication in ApJL

Journal ref: ApJ 802 (2015) L14

Sulphur-bearing molecules in diffuse molecular clouds: new results from SOFIA/GREAT and the IRAM 30 m telescope

Authors: D. A. Neufeld, B. Godard, M. Gerin, G. Pineau des Forêts, C. Bernier, E. Falgarone, U. U. Graf, R. Güsten, E. Herbst, P. Lesaffre, P. Schilke, P. Sonnentrucker, H. Wiesemeyer

Abstract: 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Π_{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 E… ▽ More 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Π_{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. △ Less

Submitted 19 February, 2015; originally announced February 2015.

Comments: 30 pages, accepted for publication in A&A

Journal ref: A&A 577, A49 (2015)

Herschel Survey of Galactic OH+, H2O+, and H3O+: Probing the Molecular Hydrogen Fraction and Cosmic-Ray Ionization Rate

Authors: Nick Indriolo, D. A. Neufeld, M. Gerin, P. Schilke, A. O. Benz, B. Winkel, K. M. Menten, E. T. Chambers, John H. Black, S. Bruderer, E. Falgarone, B. Godard, J. R. Goicoechea, H. Gupta, D. C. Lis, V. Ossenkopf, C. M. Persson, P. Sonnentrucker, F. F. S. van der Tak, E. F. van Dishoeck, Mark G. Wolfire, F. Wyrowski

Abstract: 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 ra… ▽ More 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. △ Less

Submitted 2 December, 2014; originally announced December 2014.

Comments: 76 pages, 25 figures, 6 tables; accepted for publication in ApJ

Journal ref: ApJ 800 (2015) 40

Herschel Far-Infrared Spectral-mapping of Orion BN/KL Outflows: Spatial distribution of excited CO, H2O, OH, O and C+ in shocked gas

Authors: Javier R. Goicoechea, Luis Chavarria, Jose Cernicharo, David A. Neufeld, Roland Vavrek, Edwin A. Bergin, Sara Cuadrado, Pierre Encrenaz, Mireya Etxaluze, Gary J. Melnick, Edward Polehampton

Abstract: We present ~2'x2' spectral-maps of Orion BN/KL outflows taken with Herschel at ~12'' resolution. For the first time in the far-IR domain, we spatially resolve the emission associated with the bright H2 shocked regions "Peak 1" and "Peak 2" from that of the Hot Core and ambient cloud. We analyze the ~54-310um spectra taken with the PACS and SPIRE spectrometers. More than 100 lines are detected, mos… ▽ More We present ~2'x2' spectral-maps of Orion BN/KL outflows taken with Herschel at ~12'' resolution. For the first time in the far-IR domain, we spatially resolve the emission associated with the bright H2 shocked regions "Peak 1" and "Peak 2" from that of the Hot Core and ambient cloud. We analyze the ~54-310um spectra taken with the PACS and SPIRE spectrometers. More than 100 lines are detected, most of them rotationally excited lines of 12CO (up to J=48-47), H2O, OH, 13CO, and HCN. Peaks 1/2 are characterized by a very high L(CO)/L(FIR)~5x10^{-3} ratio and a plethora of far-IR H2O emission lines. The high-J CO and OH lines are a factor ~2 brighter toward Peak 1 whereas several excited H2O lines are ~50% brighter toward Peak 2. A simplified non-LTE model allowed us to constrain the dominant gas temperature components. Most of the CO column density arises from Tk~200-500 K gas that we associate with low-velocity shocks that fail to sputter grain ice mantles and show a maximum gas-phase H2O/CO~10^{-2} abundance ratio. In addition, the very excited CO (J>35) and H2O lines reveal a hotter gas component (Tk~2500 K) from faster (v_S>25 km/s) shocks that are able to sputter the frozen-out H2O and lead to high H2O/CO>~1 abundance ratios. The H2O and OH luminosities cannot be reproduced by shock models that assume high (undepleted) abundances of atomic oxygen in the preshock gas and/or neglect the presence of UV radiation in the postshock gas. Although massive outflows are a common feature in other massive star-forming cores, Orion BN/KL seems more peculiar because of its higher molecular luminosities and strong outflows caused by a recent explosive event. △ Less

Submitted 11 November, 2014; originally announced November 2014.

Comments: 21 pages, 11 figures. Accepted for publication in ApJ (english not edited)

Shockingly low water abundances in Herschel / PACS observations of low-mass protostars in Perseus

Authors: A. Karska, L. E. Kristensen, E. F. van Dishoeck, M. N. Drozdovskaya, J. C. Mottram, G. J. Herczeg, S. Bruderer, S. Cabrit, N. J. Evans II, D. Fedele, A. Gusdorf, J. K. Jorgensen, M. J. Kaufman, G. J. Melnick, D. A. Neufeld, B. Nisini, G. Santangelo, M. Tafalla, S. F. Wampfler

Abstract: Protostars interact with their surroundings through jets and winds impacting on the envelope and creating shocks, but the nature of these shocks is still poorly understood. Our aim is to survey far-infrared molecular line emission from a uniform and significant sample of deeply-embedded low-mass young stellar objects in order to characterize shocks and the possible role of ultraviolet radiation in… ▽ More Protostars interact with their surroundings through jets and winds impacting on the envelope and creating shocks, but the nature of these shocks is still poorly understood. Our aim is to survey far-infrared molecular line emission from a uniform and significant sample of deeply-embedded low-mass young stellar objects in order to characterize shocks and the possible role of ultraviolet radiation in the immediate protostellar environment. Herschel/PACS spectral maps of 22 objects in the Perseus molecular cloud were obtained as part of the `William Herschel Line Legacy' survey. Line emission from H$_\mathrm{2}$O, CO, and OH is tested against shock models from the literature. Observed line ratios are remarkably similar and do not show variations with source physical parameters. Observations show good agreement with the shock models when line ratios of the same species are compared. Ratios of various H$_\mathrm{2}$O lines provide a particularly good diagnostic of pre-shock gas densities, $n_\mathrm{H}\sim10^{5}$ cm$^{-3}$, in agreement with typical densities obtained from observations of the post-shock gas. The corresponding shock velocities, obtained from comparison with CO line ratios, are above 20 km\,s$^{-1}$. However, the observations consistently show one-to-two orders of magnitude lower H$_\mathrm{2}$O-to-CO and H$_\mathrm{2}$O-to-OH line ratios than predicted by the existing shock models. The overestimated model H$_\mathrm{2}$O fluxes are most likely caused by an overabundance of H$_\mathrm{2}$O in the models since the excitation is well-reproduced. Illumination of the shocked material by ultraviolet photons produced either in the star-disk system or, more locally, in the shock, would decrease the H$_\mathrm{2}$O abundances and reconcile the models with observations. Detections of hot H$_\mathrm{2}$O and strong OH lines support this scenario. △ Less

Submitted 16 September, 2014; originally announced September 2014.

Comments: 28 pages, 12 figures, accepted to Astronomy & Astrophysics

Journal ref: A&A 572, A9 (2014)

Ubiquitous argonium (ArH$^+$) in the diffuse interstellar medium -- a molecular tracer of almost purely atomic gas

Authors: Peter Schilke, David A. Neufeld, Holger S. P. Müller, Claudia Comito, Edwin A. Bergin, Dariusz C. Lis, Maryvonne Gerin, John H. Black, Mark Wolfire, Nick Indriolo, John. C. Pearson, Karl M. Menten, Benjamin Winkel, Alvaro Sanchez-Monge, Thomas Moeller, Benjamin Godard, Edith Falgarone

Abstract: 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 column densities along several lines of sight are determined and discussd in the framework of chemical models. The column densities of ArH$^+$ are compared to those of other species, tracing inte… ▽ More 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 column densities along several lines of sight are determined and discussd in the framework of chemical models. The column densities of ArH$^+$ are compared to those of other species, tracing interstellar medium (ISM) components with different H$_2$ abundances. Chemical models are constructed, taking UV radiation and cosmic ray ionization into account. Due to the detection of two isotopologues, $^{36}$ArH$^+$ and $^{38}$ArH$^+$, we are confident about the carrier assignment to ArH$^+$. NeH$^+$ is not detected with a limit of [NeH$^+$]/[ArH$^+$] $\le$ 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 with H$_2$O$^+$, which traces gas with a fractional H$_2$ abundance of $\approx $0.1. A careful analysis of variations in the ArH$^+$, OH$^+$, H$_2$O$^+$ and HF column densities promises to be a faithful tracer of the distribution of the H$_2$ fractional abundance, providing unique information on a poorly known phase in the cycle of interstellar matter, 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 to characterize the ISM. Paradoxically, the ArH$^+$ molecule is a better tracer of \new{almost} purely atomic hydrogen gas than H{\sc i} itself, since H{\sc i} can also be present in gas with a significant molecular content, but ArH$^+$ singles out gas that is $>99.9$\% atomic. △ Less

Submitted 14 April, 2014; v1 submitted 31 March, 2014; originally announced March 2014.

Comments: 12 pages, 11 figures, Astronomy and Astrophysics, accepted revised version: previous version had column density of SgrB2 sources wrong

Journal ref: Astron. Astrophys. 566 (2014) Art. No. A29

Hydrogen Fluoride toward Luminous Nearby Galaxies: NGC 253 and NGC 4945

Authors: R. R. Monje, S. Lord, E. Falgarone, D. C. Lis, D. A. Neufeld, T. G. Phillips, R. Güsten

Abstract: We present the detection of hydrogen fluoride, HF, in two luminous nearby galaxies NGC 253 and NGC 4945 using the Heterodyne Instrument for the Far-Infrared (HIFI) on board the Herschel Space Observatory. The HF line toward NGC 253 has a P-Cygni profile, while an asymmetric absorption profile is seen toward NGC 4945. The P-Cygni profile in NGC 253 suggests an outflow of molecular gas with a mass o… ▽ More We present the detection of hydrogen fluoride, HF, in two luminous nearby galaxies NGC 253 and NGC 4945 using the Heterodyne Instrument for the Far-Infrared (HIFI) on board the Herschel Space Observatory. The HF line toward NGC 253 has a P-Cygni profile, while an asymmetric absorption profile is seen toward NGC 4945. The P-Cygni profile in NGC 253 suggests an outflow of molecular gas with a mass of M(H$_2$)$_{out}$ $\sim$ 1 $\times$ 10$^7$ M$_\odot$ and an outflow rate as large as M $\sim$ 6.4 M$_\odot$ yr$^{-1}$. In the case of NGC 4945, the axisymmetric velocity components in the HF line profile is compatible with the interpretation of a fast-rotating nuclear ring surrounding the nucleus and the presence of inflowing gas. The gas falls into the nucleus with an inflow rate of $\le$ 1.2 M$_\odot$ yr$^{-1}$, inside a inner radius of $\le$ 200 pc. The gas accretion rate to the central AGN is much smaller, suggesting that the inflow can be triggering a nuclear starburst. From these results, the HF $J = 1-0$ line is seen to provide an important probe of the kinematics of absorbing material along the sight-line to nearby galaxies with bright dust continuum and a promising new tracer of molecular gas in high-redshift galaxies. △ Less

Submitted 26 February, 2014; originally announced February 2014.

Journal ref: 2014 ApJ 785 22

The water abundance behind interstellar shocks: results from $Herschel$/PACS and $Spitzer$/IRS observations of H$_2$O, CO, and H$_2$

Authors: David A. Neufeld, Antoine Gusdorf, Rolf Güsten, Greg J. Herczeg, Lars Kristensen, Gary J. Melnick, Brunella Nisini, Volker Ossenkopf, Mario Tafalla, Ewine F. van Dishoeck

Abstract: We have investigated the water abundance in shock-heated molecular gas, making use of $Herschel$ measurements of far-infrared CO and H$_2$O line emissions in combination with $Spitzer$ measurements of mid-IR H$_2$ rotational emissions. We present far-infrared line spectra obtained with $Herschel$'s PACS instrument in range spectroscopy mode towards two positions in the protostellar outflow NGC 207… ▽ More We have investigated the water abundance in shock-heated molecular gas, making use of $Herschel$ measurements of far-infrared CO and H$_2$O line emissions in combination with $Spitzer$ measurements of mid-IR H$_2$ rotational emissions. We present far-infrared line spectra obtained with $Herschel$'s PACS instrument in range spectroscopy mode towards two positions in the protostellar outflow NGC 2071 and one position each in the supernova remnants W28 and 3C391. These spectra provide unequivocal detections, at one or more positions, of 12 rotational lines of water, 14 rotational lines of CO, 8 rotational lines of OH (4 lambda doublets), and 7 fine-structure transitions of atoms or atomic ions. We first used a simultaneous fit to the CO line fluxes, along with H$_2$ rotational line fluxes measured previously by $Spitzer$, to constrain the temperature and density distribution within the emitting gas; and we then investigated the water abundances implied by the observed H$_2$O line fluxes. The water line fluxes are in acceptable agreement with standard theoretical models for nondissociative shocks that predict the complete vaporization of grain mantles in shocks of velocity $v \sim 25$ km/s, behind which the characteristic gas temperature is $\sim 1300$ K and the H$_2$O/CO ratio is 1.2 △ Less

Submitted 23 December, 2013; originally announced December 2013.

Comments: 42 pages, 15 figures, accepted for publication in the Astrophysical Journal

Interstellar water chemistry: from laboratory to observations

Authors: Ewine F. van Dishoeck, Eric Herbst, David A. Neufeld

Abstract: Water is observed throughout the universe, from diffuse interstellar clouds to protoplanetary disks around young stars, and from comets in our own solar system and exoplanetary atmospheres to galaxies at high redshifts. This review summarizes the spectroscopy and excitation of water in interstellar space as well as the basic chemical processes that form and destroy water under interstellar conditi… ▽ More Water is observed throughout the universe, from diffuse interstellar clouds to protoplanetary disks around young stars, and from comets in our own solar system and exoplanetary atmospheres to galaxies at high redshifts. This review summarizes the spectroscopy and excitation of water in interstellar space as well as the basic chemical processes that form and destroy water under interstellar conditions. Three major routes to water formation are identified: low temperature ion-molecule chemistry, high-temperature neutral-neutral chemistry and gas-ice chemistry. The rate coefficients of several important processes entering the networks are discussed in detail; several of them have been determined only in the last decade through laboratory experiments and theoretical calculations. Astronomical examples of each of the different chemical routes are presented using data from powerful new telescopes, in particular the Herschel Space Observatory. Basic chemical physics studies remain critically important to analyze astronomical data. △ Less

Submitted 17 December, 2013; originally announced December 2013.

Comments: Authors' manuscript 138 pages, 34 figures, 4 tables, published in a Thematic Issue "Astrochemistry" in Chemical Reviews (December 2013), volume 113, 9043-9085 following peer review by the American Chemical Society. The published paper is available as open access at http://pubs.acs.org/doi/abs/10.1021/cr4003177

Journal ref: Chemical Reviews, 113, 9043-9085 (2013)

Fundamental Vibrational Transitions of HCl Detected in CRL 2136

Authors: Miwa Goto, T. Usuda, T. R. Geballe, K. M. Menten, N. Indriolo, D. A. Neufeld

Abstract: We would like to understand the chemistry of dense clouds and their hot cores more quantitatively by obtaining more complete knowledge of the chemical species present in them. We have obtained high-resolution infrared absorption spectroscopy at 3-4 um toward the bright infrared source CRL 2136. The fundamental vibration-rotation band of HCl has been detected within a dense cloud for the first time… ▽ More We would like to understand the chemistry of dense clouds and their hot cores more quantitatively by obtaining more complete knowledge of the chemical species present in them. We have obtained high-resolution infrared absorption spectroscopy at 3-4 um toward the bright infrared source CRL 2136. The fundamental vibration-rotation band of HCl has been detected within a dense cloud for the first time. The HCl is probably located in the warm compact circumstellar envelope or disk of CRL 2136. The fractional abundance of HCl is (4.9-8.7)e-8, indicating that approximately 20 % of the elemental chlorine is in gaseous HCl. The kinetic temperature of the absorbing gas is 250 K, half the value determined from infrared spectroscopy of 13CO and water. The percentage of chlorine in HCl is approximately that expected for gas at this temperature. The reason for the difference in temperatures between the various molecular species is unknown. △ Less

Submitted 15 October, 2013; originally announced October 2013.

Comments: 6 pages, 3 figures, A&A in press

HIFISTARS Herschel/HIFI observations of VY Canis Majoris. Molecular-line inventory of the envelope around the largest known star

Authors: J. Alcolea, V. Bujarrabal, P. Planesas, D. Teyssier, J. Cernicharo, E. De Beck, L. Decin, C. Dominik, K. Justtanont, A. de Koter, A. P. Marston, G. Melnick, K. M. Menten, D. A. Neufeld, H. Olofsson, M. Schmidt, F. L. Schöier, R. Szczerba, L. B. F. M. Waters

Abstract: To gain insight into the physical conditions and kinematics of the warm (100-1000 K) gas around the red hyper-giant VY CMa, we performed sensitive high spectral resolution observations of molecular lines in the sub-mm/FIR using the HIFI instrument of the Herschel Space Observatory. We observed CO, H2O, and other molecular species, sampling excitation energies from a few tens to a few thousand K. T… ▽ More To gain insight into the physical conditions and kinematics of the warm (100-1000 K) gas around the red hyper-giant VY CMa, we performed sensitive high spectral resolution observations of molecular lines in the sub-mm/FIR using the HIFI instrument of the Herschel Space Observatory. We observed CO, H2O, and other molecular species, sampling excitation energies from a few tens to a few thousand K. These observations are part of the Herschel Guaranteed Time Key Program HIFISTARS. We detected the J=6-5, J=10-9, and J=16-15 lines of 12CO and 13CO at about 100, 300, and 750K above the ground state (and the 13CO J=9-8 line). These lines are crucial for improving the modelling of the internal layers of the envelope around VY CMa. We also detected 27 lines of H2O and its isotopomers, and 96 lines of species such as NH3, SiO, SO, SO2 HCN, OH and others, some of them originating from vibrationally excited levels. Three lines were not unambiguously assigned. Our observations confirm that VY CMa's envelope must consist of two or more detached components. The molecular excitation in the outer layers is significantly lower than in the inner ones, resulting in strong self-absorbed profiles in molecular lines that are optically thick in this outer envelope, for instance, low-lying lines of H2O. Except for the most abundant species, CO and H2O, most of the molecular emission detected at these sub-mm/FIR wavelengths arise from the central parts of the envelope. The spectrum of VY CMa is very prominent in vibrationally excited lines, which are caused by the strong IR pumping present in the central regions. Compared with envelopes of other massive evolved stars, VY CMa's emission is particularly strong in these vibrationally excited lines, as well as in the emission from less abundant species such as H13CN, SO, and NH3. △ Less

Submitted 9 October, 2013; originally announced October 2013.

Comments: 25 pags. 18 figs

Water Absorption from Gas Very Near the Massive Protostar AFGL 2136 IRS 1

Authors: Nick Indriolo, D. A. Neufeld, A. Seifahrt, M. J. Richter

Abstract: We present ground-based observations of the nu_1 and nu_3 fundamental bands of H2O toward the massive protostar AFGL 2136 IRS 1, identifying absorption features due to 47 different ro-vibrational transitions between 2.468 micron and 2.561 micron. Analysis of these features indicates the absorption arises in warm (T=506+-25 K), very dense (n(H2)>5*10^9 cm^-3) gas, suggesting an origin close to the… ▽ More We present ground-based observations of the nu_1 and nu_3 fundamental bands of H2O toward the massive protostar AFGL 2136 IRS 1, identifying absorption features due to 47 different ro-vibrational transitions between 2.468 micron and 2.561 micron. Analysis of these features indicates the absorption arises in warm (T=506+-25 K), very dense (n(H2)>5*10^9 cm^-3) gas, suggesting an origin close to the central protostar. The total column density of warm water is estimated to be N(H2O)=(1.02+-0.02)*10^19 cm^-2, giving a relative abundance of N(H2O)/N(H2)~10^-4. Our study represents the first extensive use of water vapor absorption lines in the near infrared, and demonstrates the utility of such observations in deriving physical parameters. △ Less

Submitted 13 August, 2013; originally announced August 2013.

Comments: 17 pages, 4 figures, 1 table, accepted for publication in The Astrophysical Journal

Journal ref: ApJ 776 (2013) 8

Far-Infrared Spectroscopy of the Galactic Center. Hot Molecular Gas: Shocks versus Radiation near SgrA*

Authors: Javier R. Goicoechea, M. Etxaluze, J. Cernicharo, M. Gerin, D. A. Neufeld, A. Contursi, T. A. Bell, M. De Luca, P. Encrenaz, N. Indriolo, D. C. Lis, E. T. Polehampton, P. Sonnentrucker

Abstract: We present a 52-671um spectral scan toward SgrA* taken with the PACS and SPIRE spectrometers onboard Herschel. The achieved angular resolution allows us to separate, for the first time at far-IR wavelengths, the emission toward the central cavity (gas in the inner central parsec of the galaxy) from that of the surrounding circum-nuclear disk. The spectrum toward SgrA* is dominated by strong [OIII]… ▽ More We present a 52-671um spectral scan toward SgrA* taken with the PACS and SPIRE spectrometers onboard Herschel. The achieved angular resolution allows us to separate, for the first time at far-IR wavelengths, the emission toward the central cavity (gas in the inner central parsec of the galaxy) from that of the surrounding circum-nuclear disk. The spectrum toward SgrA* is dominated by strong [OIII], [OI], [CII], [NIII], [NII], and [CI] fine structure lines (in decreasing order of luminosity) arising in gas irradiated by UV-photons from the central stellar cluster. In addition, rotationally excited lines of 12CO (from J=4-3 to 24-23), 13CO, H2O, OH, H3O+, HCO+ and HCN, as well as ground-state absorption lines of OH+, H2O+, H3O+, CH+, H2O, OH, HF, CH and NH are detected. The excitation of the 12CO ladder is consistent with a hot isothermal component at Tk ~ 10^{3.1} K and n(H2)< 10^4 cm^{-3}. It is also consistent with a distribution of temperature components at higher density with most CO at Tk<300 K. The detected molecular features suggest that, at present, neither very enhanced X-ray, nor cosmic-ray fluxes play a dominant role in the heating of the hot molecular gas. The hot CO component (either the bulk of the CO column or just a small fraction depending on the above scenario) results from a combination of UV- and shock-driven heating. If irradiated dense clumps/clouds do not exist, shocks likely dominate the heating of the hot molecular gas. This is consistent with the high-velocity gas detected toward SgrA*. △ Less

Submitted 6 May, 2013; originally announced May 2013.

Comments: 20 pages, 3 figures, 1 table. Accepted for publication in Astrophysical Journal Letters