Showing 1–42 of 42 results for author: Kaufman, M J

Sagittarius B1 -- A Patchwork of H II Regions and PhotoDissociation Regions

Authors: Janet P. Simpson, Sean W. J. Colgan, Angela S. Cotera, Michael J. Kaufman, Susan R. Stolovy

Abstract: Sgr B1 is a luminous H II region in the Galactic Center immediately next to the massive star-forming giant molecular cloud Sgr B2 and apparently connected to it from their similar radial velocities. In 2018 we showed from SOFIA FIFI-LS observations of the [O III] 52 and 88 micron lines that there is no central exciting star cluster and that the ionizing stars must be widely spread throughout the r… ▽ More Sgr B1 is a luminous H II region in the Galactic Center immediately next to the massive star-forming giant molecular cloud Sgr B2 and apparently connected to it from their similar radial velocities. In 2018 we showed from SOFIA FIFI-LS observations of the [O III] 52 and 88 micron lines that there is no central exciting star cluster and that the ionizing stars must be widely spread throughout the region. Here we present SOFIA FIFI-LS observations of the [O I] 146 and [C II] 158 micron lines formed in the surrounding photodissociation regions (PDRs). We find that these lines correlate neither with each other nor with the [O III] lines although together they correlate better with the 70 micron Herschel PACS images from Hi-GAL. We infer from this that Sgr B1 consists of a number of smaller H II regions plus their associated PDRs, some seen face-on and the others seen more or less edge-on. We used the PDR Toolbox to estimate densities and the far-ultraviolet intensities exciting the PDRs. Using models computed with Cloudy, we demonstrate possible appearances of edge-on PDRs and show that the density difference between the PDR densities and the electron densities estimated from the [O III] line ratios is incompatible with pressure equilibrium unless there is a substantial pressure contribution from either turbulence or magnetic field or both. We likewise conclude that the hot stars exciting Sgr B1 are widely spaced throughout the region at substantial distances from the gas with no evidence of current massive star formation. △ Less

Submitted 13 February, 2021; originally announced February 2021.

Comments: 18 pages, 8 figures. Accepted by 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

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

SOFIA FIFI-LS Observations of Sgr B1: Ionization Structure and Sources of Excitation

Authors: Janet P. Simpson, Sean W. J. Colgan, Angela S. Cotera, Michael J. Kaufman, Susan R. Stolovy

Abstract: The current paradigm of Galactic Center (GC) gas motions and star formation envisions sequential star formation in streams of gas as they pass near the supermassive black hole, Sgr A*. This is based on the relative positions of dense molecular clouds, the very young star-forming region Sgr B2, the much older region Sgr C, and the several Myr old Arches and Quintuplet Clusters. Because Sgr B1 is fo… ▽ More The current paradigm of Galactic Center (GC) gas motions and star formation envisions sequential star formation in streams of gas as they pass near the supermassive black hole, Sgr A*. This is based on the relative positions of dense molecular clouds, the very young star-forming region Sgr B2, the much older region Sgr C, and the several Myr old Arches and Quintuplet Clusters. Because Sgr B1 is found with Sgr B2 in a common envelope of molecular gas and far-infrared emission, the two sources are thought to be physically related, even though there are indicators of a significantly greater age for Sgr B1. To clarify the status of Sgr B1, we have mapped it with the FIFI-LS spectrometer on SOFIA in the far-infrared lines of [O III] 52 and 88 micron. From the ratios of these lines and lines measured with the Spitzer Infrared Spectrograph, we find that there are at least eight separate sub-regions that must contain the stars that excite the gas. We infer spectral energy distributions (SEDs) of the ionizing sources from models and find they are in agreement only with SEDs of late O stars augmented at the highest frequencies with interstellar X-rays from fast shocks. We suggest that although the gas, from its velocity structure, must be part of the very young Sgr B2 complex, the stars that are ionizing the gas were not formed there but are the remnants of a previous generation of star formation in the GC. △ Less

Submitted 25 October, 2018; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: 9 pages, 3 figures. Accepted by ApJL. New version fixes typos and updates references

The Herschel-PACS legacy of low-mass protostars: Properties of warm and hot gas and its origin in far-UV illuminated shocks

Authors: Agata Karska, Michael J. Kaufman, Lars E. Kristensen, Ewine F. van Dishoeck, Gregory J. Herczeg, Joseph C. Mottram, Łukasz Tychoniec, Johan E. Lindberg, Neal J. Evans II, Joel D. Green, Yao-Lun Yang, Antoine Gusdorf, Dominika Itrich, Natasza Siódmiak

Abstract: Recent observations from Herschel allow the identification of important mechanisms responsible for the heating of gas surrounding low-mass protostars and its subsequent cooling in the far-infrared (FIR). Shocks are routinely invoked to reproduce some properties of the far-IR spectra, but standard models fail to reproduce the emission from key molecules, e.g. H$_2$O. Here, we present the Herschel-P… ▽ More Recent observations from Herschel allow the identification of important mechanisms responsible for the heating of gas surrounding low-mass protostars and its subsequent cooling in the far-infrared (FIR). Shocks are routinely invoked to reproduce some properties of the far-IR spectra, but standard models fail to reproduce the emission from key molecules, e.g. H$_2$O. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars (Class 0/I). The Herschel-PACS spectral maps covering $\sim55-210$ $μ$m with a field-of-view of $\sim$50'' are used to quantify the gas excitation conditions and spatial extent using rotational transitions of H$_{2}$O, high-$J$ CO, and OH, as well as [O I] and [C II]. We confirm that a warm ($\sim$300 K) CO reservoir is ubiquitous and that a hotter component ($760\pm170$ K) is frequently detected around protostars. The line emission is extended beyond $\sim$1000 AU spatial scales in 40/90 objects, typically in molecular tracers in Class 0 and atomic tracers in Class I objects. High-velocity emission ($\gtrsim90$ km s$^{-1}$) is detected in only 10 sources in the [O I] line, suggesting that the bulk of [O I] arises from gas that is moving slower than typical jets. Line flux ratios show an excellent agreement with models of $C$-shocks illuminated by UV photons for pre-shock densities of $\sim$$10^5$ cm$^{-3}$ and UV fields 0.1-10 times the interstellar value. The far-IR molecular and atomic lines are a unique diagnostic of feedback from UV emission and shocks in envelopes of deeply embedded protostars. △ Less

Submitted 9 February, 2018; originally announced February 2018.

Comments: Accepted by ApJS, 55 pages, 170 figures

Shocks and UV radiation around low-mass protostars: the Herschel-PACS legacy

Authors: Agata Karska, Michael J. Kaufman, Lars E. Kristensen, Ewine F. van Dishoeck

Abstract: Far-infrared spectroscopy reveals gas cooling and its underlying heating due to physical processes taking place in the surroundings of protostars. These processes are reflected in both the chemistry and excitation of abundant molecular species. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars from the WISH (van Dishoeck et al. 2011), DIGIT (Green et al. 201… ▽ More Far-infrared spectroscopy reveals gas cooling and its underlying heating due to physical processes taking place in the surroundings of protostars. These processes are reflected in both the chemistry and excitation of abundant molecular species. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars from the WISH (van Dishoeck et al. 2011), DIGIT (Green et al. 2013), and WILL surveys (Mottram et al. 2017). The $5\times5$ spectra covering the $\sim50''\times50''$ field-of-view include rotational transitions of CO, H$_2$O, and OH lines, as well as fine-structure [O I] and [C II] in the $\sim$50-200 $μ$m range. The CO rotational temperatures (for $J_\mathrm{u}\geq14)$ are typically $\sim$300 K, with some sources showing additional components with temperatures as high as $\sim$1000 K. The H$_2$O / CO and H$_2$O / OH flux ratios are low compared to stationary shock models, suggesting that UV photons may dissociate some H$_2$O and decrease its abundance. Comparison to C shock models illuminated by UV photons shows good agreement between the line emission and the models for pre-shock densities of $10^5$ cm$^{-3}$ and UV fields 0.1-10 times the interstellar value. The far-infrared molecular and atomic lines are the unique diagnostic of shocks and UV fields in deeply-embedded sources. △ Less

Submitted 8 November, 2017; originally announced November 2017.

Comments: To appear in "Astrochemistry VII -- Through the Cosmos from Galaxies to Planets", proceedings of the IAU Symposium No. 332, 2017, Puerto Varas, Chile; eds. M. Cunningham, T. Millar and Y. Aikawa (4 pages, 4 figures)

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

Mid-J CO Shock Tracing Observations of Infrared Dark Clouds III: SLED fitting

Authors: A. Pon, M. J. Kaufman, D. Johnstone, P. Caselli, F. Fontani, M. J. Butler, I. Jiménez-Serra, A. Palau, J. C. Tan

Abstract: Giant molecular clouds contain supersonic turbulence that can locally heat small fractions of gas to over 100 K. We run shock models for low-velocity, C-type shocks propagating into gas with densities between 10^3 and 10^5 cm^(-3) and find that CO lines are the most important cooling lines. Comparison to photodissociation region (PDR) models indicates that mid-J CO lines (J = 8-7 and higher) shoul… ▽ More Giant molecular clouds contain supersonic turbulence that can locally heat small fractions of gas to over 100 K. We run shock models for low-velocity, C-type shocks propagating into gas with densities between 10^3 and 10^5 cm^(-3) and find that CO lines are the most important cooling lines. Comparison to photodissociation region (PDR) models indicates that mid-J CO lines (J = 8-7 and higher) should be dominated by emission from shocked gas. In Papers I and II we presented CO J = 3-2, 8-7, and 9-8 observations toward four primarily quiescent clumps within infrared dark clouds. Here we fit PDR models to the combined spectral line energy distributions and show that the PDR models that best fit the low-J CO emission underpredict the mid-J CO emission by orders of magnitude, strongly hinting at a hot gas component within these clumps. The low-J CO data clearly show that the integrated intensities of both the CO J = 8-7 and 9-8 lines are anomalously high, such that the line ratio can be used to characterize the hot gas component. Shock models are reasonably consistent with the observed mid-J CO emission, with models with densities near 10^(4.5) cm^(-3) providing the best agreement. Where this mid-J CO is detected, the mean volume filling factor of the hot gas is 0.1%. Much of the observed mid-J CO emission, however, is also associated with known protostars and may be due to protostellar feedback. △ Less

Submitted 25 July, 2016; v1 submitted 9 June, 2016; originally announced June 2016.

Comments: 28 pages, 18 figures, accepted by ApJ, grammatical and typesetting corrections

Journal ref: ApJ 827, 107 (2016)

O$_2$ Emission Toward Orion H$_2$ Peak 1 and the Role of FUV-Illuminated C-Shocks

Authors: Gary J. Melnick, Michael J. Kaufman

Abstract: Molecular oxygen, O_2, has been the target of ground-based and space-borne searches for decades. Of the thousands of lines of sight surveyed, only those toward Rho Oph and Orion H_2 Peak 1 have yielded detections of any statistical significance. The detection of the O_2 N_J =3_3 -1_2 and 5_4 - 3_4 lines at 487.249 GHz and 773.840 GHz, respectively, toward Rho Ophiuchus has been attributed to a sho… ▽ More Molecular oxygen, O_2, has been the target of ground-based and space-borne searches for decades. Of the thousands of lines of sight surveyed, only those toward Rho Oph and Orion H_2 Peak 1 have yielded detections of any statistical significance. The detection of the O_2 N_J =3_3 -1_2 and 5_4 - 3_4 lines at 487.249 GHz and 773.840 GHz, respectively, toward Rho Ophiuchus has been attributed to a short-lived peak in the time-dependent, cold-cloud O_2 abundance, while the detection of the O_2 N_J =3_3 - 1_2, 5_4 - 3_4 lines, plus the 7_6 - 5_6 line at 1120.715 GHz, toward Orion has been ascribed to time-dependent preshock physical and chemical evolution and low-velocity (12 km/s) non-dissociative C-type shocks, both of which are fully shielded from far-ultraviolet (FUV) radiation, plus a postshock region that is exposed to a FUV field. We report a re-interpretation of the Orion O_2 detection based on new C-type shock models that fully incorporate the significant effects the presence of even a weak FUV field can have on the preshock gas, shock structure and postshock chemistry. In particular, we show that a family of solutions exists, depending on the FUV intensity, that reproduces both the observed O_2 intensities and O_2 line ratios. The solution in closest agreement with the shock parameters inferred for H_2 Peak 1 from other gas tracers assumes a 23 km/s shock impacting gas with a preshock density of 8x10^4 cm^-3 and G_0 =1, substantially different from that inferred for the fully-shielded shock case. As pointed out previously, the similarity between the LSR velocity of all three O_2 lines (~11 km/s) and recently measured H_2O 5_32 - 4_41 maser emission at 620.701 GHz toward H_2 Peak 1 suggests that the O_2 emission arises behind the same shocks responsible for the maser emission, though the O_2 emission is almost certainly more extended than the localized high density maser spots △ Less

Submitted 14 May, 2015; originally announced May 2015.

Comments: 29 pages, 11 figures

A New Determination of the Binding Energy of Atomic Oxygen on Dust Grain Surfaces: Experimental Results and Simulations

Authors: Jiao He, Jianming Shi, Tyler Hopkins, Gianfranco Vidali, Michael J. Kaufman

Abstract: The energy to desorb atomic oxygen from an interstellar dust grain surface, $E_{\rm des}$, is an important controlling parameter in gas-grain models; its value impacts the temperature range over which oxygen resides on a dust grain. However, no prior measurement has been done of the desorption energy. We report the first direct measurement of $E_{\rm des}$ for atomic oxygen from dust grain analogs… ▽ More The energy to desorb atomic oxygen from an interstellar dust grain surface, $E_{\rm des}$, is an important controlling parameter in gas-grain models; its value impacts the temperature range over which oxygen resides on a dust grain. However, no prior measurement has been done of the desorption energy. We report the first direct measurement of $E_{\rm des}$ for atomic oxygen from dust grain analogs. The values of $E_{\rm des}$ are $1660\pm 60$~K and $1850\pm 90$~K for porous amorphous water ice and for a bare amorphous silicate film, respectively, or about twice the value previously adopted in simulations of the chemical evolution of a cloud. We use the new values to study oxygen chemistry as a function of depth in a molecular cloud. For $n=10^4$ cm$^{-3}$ and $G_0$=10$^2$ ($G_0$=1 is the average local interstellar radiation field), the main result of the adoption of the higher oxygen binding energy is that H$_2$O can form on grains at lower visual extinction $A_{\rm V}$, closer to the cloud surface. A higher binding energy of O results in more formation of OH and H$_2$O on grains, which are subsequently desorbed by FUV radiation, with consequences for gas-phase chemistry. For higher values of $n$ and $G_0$, the higher binding energy can lead to a large increase in the column of H$_2$O but a decrease in the column of O$_2$. △ Less

Submitted 29 January, 2015; v1 submitted 15 January, 2015; originally announced January 2015.

Comments: 7 pages, 8 figures. Accepted by ApJ

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)

Mid-J CO observations of Perseus B1-East 5: evidence for turbulent dissipation via low-velocity shocks

Authors: Andy Pon, Doug Johnstone, Michael J. Kaufman, Paola Caselli, René Plume

Abstract: Giant molecular clouds contain supersonic turbulence and magnetohydrodynamic simulations predict that this turbulence should decay rapidly. Such turbulent dissipation has the potential to create a warm (T ~100 K) gas component within a molecular cloud. We present observations of the CO J = 5-4 and 6-5 transitions, taken with the Herschel Space Observatory, towards the Perseus B1-East 5 region. We… ▽ More Giant molecular clouds contain supersonic turbulence and magnetohydrodynamic simulations predict that this turbulence should decay rapidly. Such turbulent dissipation has the potential to create a warm (T ~100 K) gas component within a molecular cloud. We present observations of the CO J = 5-4 and 6-5 transitions, taken with the Herschel Space Observatory, towards the Perseus B1-East 5 region. We combine these new observations with archival measurements of lower rotational transitions and fit photodissociation region models to the data. We show that Perseus B1-E5 has an anomalously large CO J = 6-5 integrated intensity, consistent with a warm gas component existing within the region. This excess emission is consistent with predictions for shock heating due to the dissipation of turbulence in low velocity shocks with the shocks having a volume filling factor of 0.15 per cent. We find that B1-E has a turbulent energy dissipation rate of 3.5 x 10$^{32}$ erg / s and a dissipation time-scale that is only a factor of 3 larger than the flow crossing time-scale. △ Less

Submitted 6 March, 2015; v1 submitted 8 September, 2014; originally announced September 2014.

Comments: 18 pages, 15 figures, 4 tables, accepted by MNRAS, fixed errors described in erratum

Journal ref: MNRAS 445 (2014) 1508-1520

The Chemistry of Interstellar OH+, H2O+, and H3O+: Inferring the Cosmic Ray Ionization Rates from Observations of Molecular Ions

Authors: David Hollenbach, M. J. Kaufman, D. Neufeld, M. Wolfire, J. R. Goicoechea

Abstract: We model the production of OH+, H2O+, and H3O+ in interstellar clouds, using a steady state photodissociation region code that treats the freeze-out of gas species, grain surface chemistry, and desorption of ices from grains. The code includes PAHs, which have important effects on the chemistry. All three ions generally have two peaks in abundance as a function of depth into the cloud, one at A_V<… ▽ More We model the production of OH+, H2O+, and H3O+ in interstellar clouds, using a steady state photodissociation region code that treats the freeze-out of gas species, grain surface chemistry, and desorption of ices from grains. The code includes PAHs, which have important effects on the chemistry. All three ions generally have two peaks in abundance as a function of depth into the cloud, one at A_V<~1 and one at A_V~3-8, the exact values depending on the ratio of incident ultraviolet flux to gas density. For relatively low values of the incident far ultraviolet flux on the cloud (χ<~ 1000; χ= 1= local interstellar value), the columns of OH+ and H2O+ scale roughly as the cosmic ray primary ionization rate ζ(crp) divided by the hydrogen nucleus density n. The H3O+ column is dominated by the second peak, and we show that if PAHs are present, N(H3O+) ~ 4x10^{13} cm^{-2} independent of ζ(crp) or n. If there are no PAHs or very small grains at the second peak, N(H3O+) can attain such columns only if low ionization potential metals are heavily depleted. We also model diffuse and translucent clouds in the interstellar medium, and show how observations of N(OH+)/N(H) and N(OH+)/N(H2O+) can be used to estimate ζ(crp)/n, χ/n and A_V in them. We compare our models to Herschel observations of these two ions, and estimate ζ(crp) ~ 4-6 x 10^-16 (n/100 cm^-3) s^-1 and χ/n = 0.03 cm^3 for diffuse foreground clouds towards W49N. △ Less

Submitted 29 May, 2012; originally announced May 2012.

Herschel Search for O2 Toward the Orion Bar

Authors: Gary J. Melnick, Volker Tolls, Paul F. Goldsmith, Michael J. Kaufman, David J. Hollenbach, John H. Black, Pierre Encrenaz, Edith Falgarone, Maryvonne Gerin, Åke Hjalmarson, Di Li, Dariusz C. Lis, René Liseau, David A. Neufeld, Laurent Pagani, Ronald L. Snell, Floris van der Tak, Ewine F. van Dishoeck

Abstract: We report the results of a search for molecular oxygen (O2) toward the Orion Bar, a prominent photodissociation region at the southern edge of the HII region created by the luminous Trapezium stars. We observed the spectral region around the frequency of the O2 N_J = 3_3 - 1_2 transition at 487 GHz and the 5_4 - 3_4 transition at 774 GHz using the Heterodyne Instrument for the Far Infrared on the… ▽ More We report the results of a search for molecular oxygen (O2) toward the Orion Bar, a prominent photodissociation region at the southern edge of the HII region created by the luminous Trapezium stars. We observed the spectral region around the frequency of the O2 N_J = 3_3 - 1_2 transition at 487 GHz and the 5_4 - 3_4 transition at 774 GHz using the Heterodyne Instrument for the Far Infrared on the Herschel Space Observatory. Neither line was detected, but the 3sigma upper limits established here translate to a total line-of-sight O2 column density < 1.5 10^16 cm^-2 for an emitting region whose temperature is between 30K and 250 K, or < 1 10^16 cm^-2 if the O2 emitting region is primarily at a temperature of ~< 100 K. Because the Orion Bar is oriented nearly edge-on relative to our line of sight, the observed column density is enhanced by a factor estimated to be between 4 and 20 relative to the face-on value. Our upper limits imply that the face-on O2 column density is less than 4 10^15 cm^-2, a value that is below, and possibly well below, model predictions for gas with a density of 10^4 - 10^5 cm^-3 exposed to a far ultraviolet flux 10^4 times the local value, conditions inferred from previous observations of the Orion Bar. The discrepancy might be resolved if: (1) the adsorption energy of O atoms to ice is greater than 800 K; (2) the total face-on Av of the Bar is less than required for O2 to reach peak abundance; (3) the O2 emission arises within dense clumps with a small beam filling factor; or, (4) the face-on depth into the Bar where O2 reaches its peak abundance, which is density dependent, corresponds to a sky position different from that sampled by our Herschel beams. △ Less

Submitted 25 April, 2012; originally announced April 2012.

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

Multi-line detection of O2 toward rho Oph A

Authors: R. Liseau, P. F. Goldsmith, B. Larsson, L. Pagani, P. Bergman, J. Le Bourlot, T. A. Bell, A. O. Benz, E. A. Bergin, P. Bjerkeli, J. H. Black, S. Bruderer, P. Caselli, E. Caux, J. -H. Chen, M. de Luca, P. Encrenaz, E. Falgarone, M. Gerin, J. R. Goicoechea, Å. Hjalmarson, D. J. Hollenbach, K. Justtanont, M. J. Kaufman, F. Le Petit , et al. (14 additional authors not shown)

Abstract: Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a bea… ▽ More Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a beam size of 10 arcmin. A larger telescope aperture such as that of the Herschel Space Observatory is required to resolve the O2 emission and to pinpoint its origin. We use the Heterodyne Instrument for the Far Infrared aboard Herschel to obtain high resolution O2 spectra toward selected positions in rho Oph A. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. The 487.2GHz line was clearly detected toward all three observed positions in rho Oph A. In addition, an oversampled map of the 773.8GHz transition revealed the detection of the line in only half of the observed area. Based on their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (> 50 K) adjacent to a much colder region, where temperatures are below 30 K. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on the derived gas temperatures, but fall into the range N(O2) = (3 to >6)e15/cm^2. Beam averaged O2 abundances are about 5e-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz of about 4 - 5 arcmin, encompassing the entire rho Oph A core. △ Less

Submitted 25 February, 2012; originally announced February 2012.

Comments: 10 pages, 8 figures intended for publication in A&A

Molecular Tracers of Turbulent Shocks in Giant Molecular Clouds

Authors: A. Pon, D. Johnstone, M. J. Kaufman

Abstract: Giant molecular clouds contain supersonic turbulence and simulations of magnetohydrodynamic turbulence show that these supersonic motions decay in roughly a crossing time, which is less than the estimated lifetimes of molecular clouds. Such a situation requires a significant release of energy. We run models of C-type shocks propagating into gas with densities around 10^3 cm^(-3) at velocities of a… ▽ More Giant molecular clouds contain supersonic turbulence and simulations of magnetohydrodynamic turbulence show that these supersonic motions decay in roughly a crossing time, which is less than the estimated lifetimes of molecular clouds. Such a situation requires a significant release of energy. We run models of C-type shocks propagating into gas with densities around 10^3 cm^(-3) at velocities of a few km / s, appropriate for the ambient conditions inside of a molecular cloud, to determine which species and transitions dominate the cooling and radiative energy release associated with shock cooling of turbulent molecular clouds. We find that these shocks dissipate their energy primarily through CO rotational transitions and by compressing pre-existing magnetic fields. We present model spectra for these shocks and by combining these models with estimates for the rate of turbulent energy dissipation, we show that shock emission should dominate over emission from unshocked gas for mid to high rotational transitions (J >5) of CO. We also find that the turbulent energy dissipation rate is roughly equivalent to the cosmic-ray heating rate and that the ambipolar diffusion heating rate may be significant, especially in shocked gas. △ Less

Submitted 25 January, 2012; v1 submitted 9 January, 2012; originally announced January 2012.

Comments: 14 pages, 5 figures, accepted for publication in ApJ; minor grammatical errors corrected

Journal ref: ApJ 748 (2012) 25-39

PDR Model Mapping of Physical Conditions via Spitzer/IRS Spectroscopy of H2: Theoretical Success toward NGC 2023-South

Authors: Y. Sheffer, M. G. Wolfire, D. J. Hollenbach, M. J. Kaufman, M. Cordier

Abstract: We use the IRS on Spitzer to observe the southern part of the reflection nebula NGC 2023, including the Southern Ridge, which is a photodissociation region par excellence excited by HD 37903. Five pure-rotational H2 emission lines are detected and mapped over and around the Southern Ridge in order to compare with predicted level column densities from theoretical PDR models. We find very good agree… ▽ More We use the IRS on Spitzer to observe the southern part of the reflection nebula NGC 2023, including the Southern Ridge, which is a photodissociation region par excellence excited by HD 37903. Five pure-rotational H2 emission lines are detected and mapped over and around the Southern Ridge in order to compare with predicted level column densities from theoretical PDR models. We find very good agreement between PDR model predictions and emission line intensities and ratios measured with Spitzer, leading us to conclude that grain photoelectric heating sufficiently warms the gas to produce the observed H2 line emission via collisional excitation. On the Southern Ridge, we infer a hydrogen nucleus density n_H approx 2x10^5 cm^-3 and radiation field strength χ approx 10^4 relative to the local Galactic interstellar radiation field. This high value for χ independently predicts a distance toward HD 37903 of 300 pc, and is consistent with the most recent HIPPARCOS results. Over the map we find that both n_H and χ vary by a factor of sim 3. Such 2-D variations provide clues about the underlying 3-D structure of the Southern Ridge field, which appears to be the tip of a molecular cloud. We also map variations in excitation temperature and the ortho-to-para ratio, the latter attaining values of sim 1.5---2.0 on the Southern Ridge, and find that PDR modeling can readily reproduce observed ortho-to-para ratios that are < 3 for rotational excitation dominated by collisional processes. Last, the stars Sellgren C and G are discovered to be resolved on archival HST images into two point sources each, with separations of lesssim 0.5". △ Less

Submitted 20 October, 2011; originally announced October 2011.

Comments: 16 pages, 13 figures, to be published in 2011 Nov 1 issue of ApJ

Journal ref: ApJ 741, 45 (2011)

Herschel Measurements of Molecular Oxygen in Orion

Authors: Paul F. Goldsmith, Rene Liseau, Tom A. Bell, John H. Black, Jo-Hsin Chen, David Hollenbach, Michael J. Kaufman, Di Li, Dariusz C. Lis, Gary Melnick, David Neufeld, Laurent Pagani, Ronald Snell, Arnold O. Benz, Edwin Bergin, Simon Bruderer, Paola Caselli, Emmanuel Caux, Pierre Encrenaz, Edith Falgarone, Maryvonne Gerin, Javier R. Goicoechea, Ake Hjalmarson, Bengt Larsson, Jacques Le Bourlot , et al. (9 additional authors not shown)

Abstract: We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 c… ▽ More We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2. △ Less

Submitted 1 August, 2011; originally announced August 2011.

Distribution of Water Vapor in Molecular Clouds

Authors: Gary J. Melnick, Volker Tolls, Ronald L. Snell, Edwin A. Bergin, David J. Hollenbach, Michael J. Kaufman, Di Li, David A. Neufeld

Abstract: We report the results of a large-area study of water vapor along the Orion Molecular Cloud ridge, the purpose of which was to determine the depth-dependent distribution of gas-phase water in dense molecular clouds. We find that the water vapor measured toward 77 spatial positions along the face-on Orion ridge, excluding positions surrounding the outflow associated with BN/KL and IRc2, display inte… ▽ More We report the results of a large-area study of water vapor along the Orion Molecular Cloud ridge, the purpose of which was to determine the depth-dependent distribution of gas-phase water in dense molecular clouds. We find that the water vapor measured toward 77 spatial positions along the face-on Orion ridge, excluding positions surrounding the outflow associated with BN/KL and IRc2, display integrated intensities that correlate strongly with known cloud surface tracers such as CN, C2H, 13CO J =5-4, and HCN, and less well with the volume tracer N2H+. Moreover, at total column densities corresponding to Av < 15 mag., the ratio of H2O to C18O integrated intensities shows a clear rise approaching the cloud surface. We show that this behavior cannot be accounted for by either optical depth or excitation effects, but suggests that gas-phase water abundances fall at large Av. These results are important as they affect measures of the true water-vapor abundance in molecular clouds by highlighting the limitations of comparing measured water vapor column densities with such traditional cloud tracers as 13CO or C18O. These results also support cloud models that incorporate freeze-out of molecules as a critical component in determining the depth-dependent abundance of water vapor. △ Less

Submitted 14 November, 2010; originally announced November 2010.

Water in low-mass star-forming regions with Herschel: HIFI spectroscopy of NGC1333

Authors: L. E. Kristensen, R. Visser, E. F. van Dishoeck, U. A. Yıldız, S. D. Doty, G. J. Herczeg, F. -C. Liu, B. Parise, J. K. Jørgensen, T. A. van Kempen, C. Brinch, S. F. Wampfler, S. Bruderer, A. O. Benz, M. R. Hogerheijde, E. Deul, R. Bachiller, A. Baudry, M. Benedettini, E. A. Bergin, P. Bjerkeli, G. A. Blake, S. Bontemps, J. Braine, P. Caselli , et al. (42 additional authors not shown)

Abstract: 'Water In Star-forming regions with Herschel' (WISH) is a key programme dedicated to studying the role of water and related species during the star-formation process and constraining the physical and chemical properties of young stellar objects. The Heterodyne Instrument for the Far-Infrared (HIFI) on the Herschel Space Observatory observed three deeply embedded protostars in the low-mass star-for… ▽ More 'Water In Star-forming regions with Herschel' (WISH) is a key programme dedicated to studying the role of water and related species during the star-formation process and constraining the physical and chemical properties of young stellar objects. The Heterodyne Instrument for the Far-Infrared (HIFI) on the Herschel Space Observatory observed three deeply embedded protostars in the low-mass star-forming region NGC1333 in several H2-16O, H2-18O, and CO transitions. Line profiles are resolved for five H16O transitions in each source, revealing them to be surprisingly complex. The line profiles are decomposed into broad (>20 km/s), medium-broad (~5-10 km/s), and narrow (<5 km/s) components. The H2-18O emission is only detected in broad 1_10-1_01 lines (>20 km/s), indicating that its physical origin is the same as for the broad H2-16O component. In one of the sources, IRAS4A, an inverse P Cygni profile is observed, a clear sign of infall in the envelope. From the line profiles alone, it is clear that the bulk of emission arises from shocks, both on small (<1000 AU) and large scales along the outflow cavity walls (~10 000 AU). The H2O line profiles are compared to CO line profiles to constrain the H2O abundance as a function of velocity within these shocked regions. The H2O/CO abundance ratios are measured to be in the range of ~0.1-1, corresponding to H2O abundances of ~10-5-10-4 with respect to H2. Approximately 5-10% of the gas is hot enough for all oxygen to be driven into water in warm post-shock gas, mostly at high velocities. △ Less

Submitted 18 July, 2010; originally announced July 2010.

Comments: Accepted for publication in the A&A HIFI special issue

Origin of the hot gas in low-mass protostars: Herschel-PACS spectroscopy of HH 46

Authors: T. A. van Kempen, L. E. Kristensen, G. J. Herczeg, R. Visser, E. F. van Dishoeck, S. F. Wampfler, S. Bruderer, A. O. Benz, S. D. Doty, C. Brinch, M. R. Hogerheijde, J. K. Jørgensen, M. Tafalla, D. Neufeld, R. Bachiller, A. Baudry, M. Benedettini, E. A. Bergin, P. Bjerkeli, G. A. Blake, S. Bontemps, J. Braine, P. Caselli, J. Cernicharo, C. Codella , et al. (40 additional authors not shown)

Abstract: 'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CI… ▽ More 'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CII] lines located between 63 and 186 um. The excitation and spatial distribution of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, [OI] emission is seen in low-velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km/s in both the red- and blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [OI] emission is attributed to J-type shocks in dense gas close to the protostar. In the scenario described here, the combined cooling by far-IR lines within the central spatial pixel is estimated to be 2 \times 10-2 L_sun, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope. △ Less

Submitted 12 May, 2010; originally announced May 2010.

Comments: Accepted for publication in Astronomy and Astrophysics (Herschel special issue)

Spitzer spectral line mapping of protostellar outflows: I. Basic data and outflow energetics

Authors: David A. Neufeld, Brunella Nisini, Teresa Giannini, Gary J. Melnick, Edwin A. Bergin, Yuan Yuan, Sebastien Maret, Volker Tolls, Rolf Guesten, Michael J. Kaufman

Abstract: We report the results of spectroscopic mapping observations carried out toward protostellar outflows in the BHR71, L1157, L1448, NGC 2071, and VLA 1623 molecular regions using the Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations, covering the 5.2 - 37 micron spectral region, provide detailed maps of the 8 lowest pure rotational lines of molecular hydrogen and of the… ▽ More We report the results of spectroscopic mapping observations carried out toward protostellar outflows in the BHR71, L1157, L1448, NGC 2071, and VLA 1623 molecular regions using the Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations, covering the 5.2 - 37 micron spectral region, provide detailed maps of the 8 lowest pure rotational lines of molecular hydrogen and of the [SI] 25.25 micron and [FeII] 26.0 micron fine structure lines. The molecular hydrogen lines, believed to account for a large fraction of the radiative cooling from warm molecular gas that has been heated by a non-dissociative shock, allow the energetics of the outflows to be elucidated. Within the regions mapped towards these 5 outflow sources, total H2 luminosities ranging from 0.02 to 0.75 L(solar) were inferred for the sum of the 8 lowest pure rotational transitions. By contrast, the much weaker [FeII] 26.0 micron fine structure transition traces faster, dissociative shocks; here, only a small fraction of the fast shock luminosity emerges as line radiation that can be detected with Spitzer/IRS. △ Less

Submitted 6 October, 2009; originally announced October 2009.

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

Journal ref: Astrophys.J.706:170-183,2009

Water, O2 and Ice in Molecular Clouds

Authors: David Hollenbach, Michael J. Kaufman, Edwin A. Bergin, Gary J. Melnick

Abstract: We model the temperature and chemical structure of molecular clouds as a function of depth into the cloud, assuming a cloud of constant density n illuminated by an external FUV (6 eV < E < 13.6 eV) flux G_0 (scaling factor in multiples of the local interstellar field). Extending previous photodissociation region models, we include the freezing of species, simple grain surface chemistry, and deso… ▽ More We model the temperature and chemical structure of molecular clouds as a function of depth into the cloud, assuming a cloud of constant density n illuminated by an external FUV (6 eV < E < 13.6 eV) flux G_0 (scaling factor in multiples of the local interstellar field). Extending previous photodissociation region models, we include the freezing of species, simple grain surface chemistry, and desorption (including FUV photodesorption) of ices. We also treat the opaque cloud interior with time-dependent chemistry. Here, under certain conditions, gas phase elemental oxygen freezes out as water ice and the elemental C/O abundance ratio can exceed unity, leading to complex carbon chemistry. Gas phase H2O and O2 peak in abundance at intermediate depth into the cloud, roughly A_V~3-8 from the surface, the depth proportional to ln(G_0/n). Closer to the surface, molecules are photodissociated. Deeper into the cloud, molecules freeze to grain surfaces. At intermediate depths photodissociation rates are attenuated by dust extinction, but photodesorption prevents total freezeout. For G_0 < 500, abundances of H2O and O2 peak at values ~10^(-7), producing columns ~10^(15) per cm^2, independent of G_0 and n. The peak abundances depend primarily on the product of the photodesorption yield of water ice and the grain surface area per H nucleus. At higher values of G_0, thermal desorption of O atoms from grains enhances the gas phase H2O peak abundance and column slightly, whereas the gas phase O2 peak abundance rises to ~10^(-5) and the column to ~2x10^(16) per cm^2. We present simple analytic equations for the abundances as a function of depth which clarify the dependence on parameters. The models are applied to observations of H2O, O2, and water ice in a number of sources, including B68, NGC 2024, and Rho Oph. △ Less

Submitted 2 March, 2009; v1 submitted 9 September, 2008; originally announced September 2008.

Comments: 70 pages including 17 figures

Journal ref: Astrophys.J.690:1497-1521,2009

Detection of Extended Hot Water in the Outflow from NGC 2071

Authors: Gary J. Melnick, Volker Tolls, David A. Neufeld, Yuan Yuan, Paule Sonnentrucker, Dan M. Watson, Edwin A. Bergin, Michael J. Kaufman

Abstract: We report the results of spectroscopic mapping observations carried out toward a ~1 min x 1 min region within the northern lobe of the outflow from NGC 2071 using the Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations covered the 5.2-37 um spectral region and have led to the detection of a number of ionic, atomic, and molecular lines, including fine-structure emission… ▽ More We report the results of spectroscopic mapping observations carried out toward a ~1 min x 1 min region within the northern lobe of the outflow from NGC 2071 using the Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations covered the 5.2-37 um spectral region and have led to the detection of a number of ionic, atomic, and molecular lines, including fine-structure emission of Si+, Fe+, S++, S, the S(0)-S(7) pure rotational lines of H2, the R(3) and R(4) transitions of HD, and at least 11 transitions of H2O. In addition, the 6.2, 7.4, 7.6, 7.9, 8.6 and 11.3 um PAH emission bands were also observed and several transitions of OH were tentatively detected. Most of the detected line transitions were strong enough to map including, for the first time, three transitions of hot H2O. We find that: (1) the water emission is extended; (2) the extended emission is aligned with the outflow; and, (3) the spatial distribution of the water emission generally follows that observed for H2. Based on the measured line intensities, we derive an HD abundance relative to H2 of 1.1-1.8 10^-5 and an H2O number density of 12-2 cm^3. The H2 density in the water-emitting region is not well constrained by our observations, but is likely between 3 10^4 and 10^6 cm^3, yielding an H2O abundance relative to H2 of between 2 10^-5 and 6 10^-4. Future observations planned for the Herschel Space Observatory should greatly improve the density estimate, and thus our knowledge of the H2O abundance, for the water-emitting regions reported here. Finally, we note a possible departure from the H2O ortho-to-para ratio of 3:1 expected for water formed in hot post-shocked gas, suggesting that a significant fraction of the water vapor we detect may arise from H2O sputtered from cold dust grains. △ Less

Submitted 5 May, 2008; originally announced May 2008.

Comments: 35 pages, 15 figures, 4 tables, accepted for publication in ApJ

Chemical Rates on Small Grains and PAHs: C^+ Recombination and H_2 Formation

Authors: Mark G. Wolfire, A. G. G. M. Tielens, David Hollenbach, Michael J. Kaufman

Abstract: We use observations of the CI, CII, HI, and H_2 column densities along lines of sight in the Galactic plane to determine the formation rate of H_2 on grains and to determine chemical reaction rates with Polycyclic Aromatic Hydrocarbons. Photodissociation region models are used to find the best fit parameters to the observed columns. We find the H_2 formation rate on grains has a low rate (R ~ 1… ▽ More We use observations of the CI, CII, HI, and H_2 column densities along lines of sight in the Galactic plane to determine the formation rate of H_2 on grains and to determine chemical reaction rates with Polycyclic Aromatic Hydrocarbons. Photodissociation region models are used to find the best fit parameters to the observed columns. We find the H_2 formation rate on grains has a low rate (R ~ 1 x 10^(-17) cm^(3) s^(-1)) along lines of sight with low column density (A_V < 0.25) and low molecular fraction (f_(H_2) < 10^(-4)). At higher column densities (0.25 < A_V <2.13), we find a rate of R ~ 3.5x10^(-17) cm^(3) s^(-1). The lower rate at low column densities could be the result of grain processing by interstellar shocks which may deplete the grain surface area or process the sites of H +H formation, thereby inhibiting H_2 production. Alternatively, the formation rate may be normal, and the low molecular fraction may be the result of lines of sight which graze larger clouds. Such lines of sight would have a reduced H_2 self-shielding compared to the line-of-sight column. We find the reaction C^+ +PAH^- --> C + PAH^0 is best fit with a rate 2.4 x 10^(-7) Φ_PAH T_2^(-0.5) cm^(3) s^(-1) with T_2= T/100 K and the reaction C^+ + PAH^0 --> C + PAH^+ is best fit with a rate 8.8x 10^(-9)Φ_PAH cm^(3) s^(-1). In high column density gas we find Φ_PAH ~ 0.4. In low column density gas, Φ_PAH is less well constrained with Φ_PAH ~ 0.2 - 0.4. △ Less

Submitted 2 March, 2008; originally announced March 2008.

Comments: 40 pages, 12 figure, accepted for publication in the Astrophysical Journal, June 10, 2008, vol 680 issue

SWAS Observations of Water in Molecular Outflows

Authors: J. Franklin, R. L. Snell, M. J. Kaufman, G. J. Melnick, D. A. Neufeld, D. J. Hollenbach, E. A. Bergin

Abstract: We present SWAS detections of the ground-state 1(10)-1(01) transition of o-H2O at 557 GHz in 18 molecular outflows. These results are combined with ground-based observations of the J=1-0 transitions of 12CO and 13CO obtained at the FCRAO and, for a subset of the outflows, data from ISO. Assuming the SWAS water line emission originates from the same gas traced by CO emission, we find that the out… ▽ More We present SWAS detections of the ground-state 1(10)-1(01) transition of o-H2O at 557 GHz in 18 molecular outflows. These results are combined with ground-based observations of the J=1-0 transitions of 12CO and 13CO obtained at the FCRAO and, for a subset of the outflows, data from ISO. Assuming the SWAS water line emission originates from the same gas traced by CO emission, we find that the outflowing gas in most outflows has an o-H2O abundance relative to H2 of between 10(-7) and 10(-6). Analysis of the water abundance as a function of outflow velocity reveals a strong dependence. The water abundance increases with velocity, and at the highest outflow velocities some outflows have relative o-H2O abundances of order 10(-4). However the mass of gas with such elevated water abundances represents less that 1% of the total outflow gas mass. The ISO LWS observations of high-J rotational lines of CO and the 179.5 micron transition of o-H2O provide evidence for a warmer outflow component than required to produce either the SWAS or FCRAO lines. The mass associated with the ISO emission is similar to that responsible for the highest velocity water emission detected by SWAS. The bulk of the outflowing gas has an abundance of o-H2O well below what would be expected if the gas has passed through a C-shock with shock velocities greater than 10 km/s. Gas-phase water can be depleted in the post-shock gas due to freeze-out onto grain mantles, however the rate of freeze-out is too slow to explain our results. Therefore we believe that only a small fraction of the outflowing molecular gas has passed through shocks strong enough to fully convert the gas-phase oxygen to water. This result has implications for the acceleration mechanism of the molecular gas in these outflows. △ Less

Submitted 13 November, 2007; originally announced November 2007.

Comments: Accepted for publication in Ap.J., 51 pages including 4 pages of figures

NICMOS Observations of Shocked H_2 in Orion

Authors: Sean W. J. Colgan, A. S. B. Schultz, M. J. Kaufman, E. F. Erickson, D. J. Hollenbach

Abstract: HST NICMOS narrowband images of the shocked molecular hydrogen emission in OMC-1 are analyzed to reveal new information on the BN/KL outflow. The outstanding morphological feature of this region is the array of molecular hydrogen ``fingers'' emanating from the general vicinity of IRc2 and the presence of several Herbig-Haro objects. The NICMOS images appear to resolve individual shock fronts. Th… ▽ More HST NICMOS narrowband images of the shocked molecular hydrogen emission in OMC-1 are analyzed to reveal new information on the BN/KL outflow. The outstanding morphological feature of this region is the array of molecular hydrogen ``fingers'' emanating from the general vicinity of IRc2 and the presence of several Herbig-Haro objects. The NICMOS images appear to resolve individual shock fronts. This work is a more quantitative and detailed analysis of our data from a previous paper (Schultz etal. 1999). Line strengths for the H_2 1--0 S(4) plus 2--1 S(6) lines at 1.89 micron are estimated from measurements with the Paschen_alpha continuum filter F190N at 1.90 micron, and continuum measurements at 1.66 and 2.15 micron. We compare the observed H_2 line strengths and ratios of the 1.89 micron and 2.12 micron 1--0 S(1) lines with models for molecular cloud shock waves. Most of the data cannot be fit by J-shocks, but are well matched by C-shocks with shock velocities in the range of 20--45 km/s and preshock densities of 10^{4} - 10^{6} cm^{-3}, similiar to values obtained in larger beam studies which averaged over many shocks. There is also some evidence that shocks with higher densities have lower velocities. △ Less

Submitted 18 September, 2007; originally announced September 2007.

Comments: 33 pages, 11 figures, 2 of those color

Spitzer IRS Observations of the Galactic Center: Shocked Gas in the Radio Arc Bubble

Authors: Janet P. Simpson, Sean W. J. Colgan, Angela S. Cotera, Edwin F. Erickson, David J. Hollenbach, Michael J. Kaufman, Robert H. Rubin

Abstract: We present Spitzer IRS spectra (R ~600, 10 - 38 micron) of 38 positions in the Galactic Center (GC), all at the same Galactic longitude and spanning plus/minus 0.3 degrees in latitude. Our positions include the Arches Cluster, the Arched Filaments, regions near the Quintuplet Cluster, the ``Bubble'' lying along the same line-of-sight as the molecular cloud G0.11-0.11, and the diffuse interstella… ▽ More We present Spitzer IRS spectra (R ~600, 10 - 38 micron) of 38 positions in the Galactic Center (GC), all at the same Galactic longitude and spanning plus/minus 0.3 degrees in latitude. Our positions include the Arches Cluster, the Arched Filaments, regions near the Quintuplet Cluster, the ``Bubble'' lying along the same line-of-sight as the molecular cloud G0.11-0.11, and the diffuse interstellar gas along the line-of-sight at higher Galactic latitudes. From measurements of the [O IV], [Ne II], [Ne III], [Si II], [S III], [S IV], [Fe II], [Fe III], and H_2 S(0), S(1), and S(2) lines we determine the gas excitation and ionic abundance ratios. The Ne/H and S/H abundance ratios are ~ 1.6 times that of the Orion Nebula. The main source of excitation is photoionization, with the Arches Cluster ionizing the Arched Filaments and the Quintuplet Cluster ionizing the gas nearby and at lower Galactic latitudes including the far side of the Bubble. In addition, strong shocks ionize gas to O^{+3} and destroy dust grains, releasing iron into the gas phase (Fe/H ~ 1.3 times 10^{-6} in the Arched Filaments and Fe/H ~ 8.8 times 10^{-6} in the Bubble). The shock effects are particularly noticeable in the center of the Bubble, but O$^{+3}$ is present in all positions. We suggest that the shocks are due to the winds from the Quintuplet Cluster Wolf-Rayet stars. On the other hand, the H_2 line ratios can be explained with multi-component models of warm molecular gas in photodissociation regions without the need for H_2 production in shocks. △ Less

Submitted 15 August, 2007; originally announced August 2007.

Comments: 51 pages, 17 figures To be published in the Astrophysical Journal

Spitzer spectral line mapping of supernova remnants: I. Basic data and principal component analysis

Authors: David A. Neufeld, David J. Hollenbach, Michael J. Kaufman, Ronald L. Snell, Gary J. Melnick, Edwin A. Bergin, Paule Sonnentrucker

Abstract: We report the results of spectroscopic mapping observations carried out toward small (1 x 1 arcmin) regions within the supernova remnants W44, W28, IC443, and 3C391 using the Infrared Spectrograph of the Spitzer Space Telescope. These observations, covering the 5.2 - 37 micron spectral region, have led to the detection of a total of 15 fine structure transitions of Ne+, Ne++, Si+, P+, S, S++, Cl… ▽ More We report the results of spectroscopic mapping observations carried out toward small (1 x 1 arcmin) regions within the supernova remnants W44, W28, IC443, and 3C391 using the Infrared Spectrograph of the Spitzer Space Telescope. These observations, covering the 5.2 - 37 micron spectral region, have led to the detection of a total of 15 fine structure transitions of Ne+, Ne++, Si+, P+, S, S++, Cl+, Fe+, and Fe++; the S(0) - S(7) pure rotational lines of molecular hydrogen; and the R(3) and R(4) transitions of hydrogen deuteride. In addition to these 25 spectral lines, the 6.2, 7.7, 8.6, 11.3 and 12.6 micron PAH emission bands were also observed. Most of the detected line transitions have proven strong enough to map in several sources, providing a comprehensive picture of the relative distribution of the various line emissions observable in the Spitzer/IRS bandpass. A principal component analysis of the spectral line maps reveals that the observed emission lines fall into five distinct groups, each of which may exhibit a distinct spatial distribution: (1) lines of S and H2 (J > 2); (2) the H2 S(0) line; (3) lines of ions with appearance potentials less than 13.6 eV; (4) lines of ions with appearance potentials greater than 13.6 eV, not including S++; (5) lines of S++. Lines of group (1) likely originate in molecular material subject to a slow, nondissociative shock that is driven by the overpressure within the supernova remnant, and lines in groups (3) - (5) are associated primarily with dissociative shock fronts with a range of (larger) shock velocities. The H2 S(0) line shows a low-density diffuse emission component, and - in some sources - a shock-excited component. △ Less

Submitted 17 April, 2007; originally announced April 2007.

Comments: 43 pages, including 21 figures. Accepted for publication in ApJ

Journal ref: Astrophys.J.664:890-908,2007

A PDR-Code Comparison Study

Authors: M. Roellig, N. P. Abel, T. Bell, F. Bensch, J. Black, G. J. Ferland, B. Jonkheid, I. Kamp, M. J. Kaufman, J. Le Bourlot, F. Le Petit, R. Meijerink, O. Morata, V. Ossenkopf, E. Roueff, G. Shaw, M. Spaans, A. Sternberg, J. Stutzki, W. -F. Thi, E. F. van Dishoeck, P. A. M. van Hoof, S. Viti, M. G. Wolfire

Abstract: We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been cr… ▽ More We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly. △ Less

Submitted 8 February, 2007; originally announced February 2007.

Comments: 28 pages, 16 figures, see all benchmark data under: http://www.ph1.uni-koeln.de/pdr-comparison

Journal ref: Astron.Astrophys.467:187,2007

Spitzer observations of hydrogen deuteride

Authors: David A. Neufeld, Joel D. Green, David J. Hollenbach, Paule Sonnentrucker, Gary J. Melnick, Edwin A. Bergin, Ronald L. Snell, William J. Forrest, Dan M. Watson, Michael J. Kaufman

Abstract: We report the detection of interstellar hydrogen deuteride (HD) toward the supernova remnant IC443, and the tentative detection of HD toward the Herbig Haro objects HH54 and HH7 and the star forming region GGD37 (Cepheus A West). Our detections are based upon spectral line mapping observations of the R(3) and R(4) rotational lines of HD, at rest wavelengths of 28.502 and 23.034 micron respective… ▽ More We report the detection of interstellar hydrogen deuteride (HD) toward the supernova remnant IC443, and the tentative detection of HD toward the Herbig Haro objects HH54 and HH7 and the star forming region GGD37 (Cepheus A West). Our detections are based upon spectral line mapping observations of the R(3) and R(4) rotational lines of HD, at rest wavelengths of 28.502 and 23.034 micron respectively, obtained using the Infrared Spectrograph onboard the Spitzer Space Telescope. The HD R(4)/R(3) line intensity ratio promises to be a valuable probe of the gas pressure in regions where it can be observed. The derived HD/H2 abundance ratios are 1.19(+0.35/-0.24)E-5, 1.80(+0.54/-0.32)E-5, and 1.41(+0.46/-0.33)E-5 respectively (68.3% confidence limits, based upon statistical errors alone) for IC443 (clump C), HH54, and HH7. If HD is the only significant reservoir of gas-phase deuterium in these sources, the inferred HD/H2 ratios are all consistent with a gas-phase elemental abundance [n(D)/n(H)](gas) ~ 7.5E-6, a factor 2 - 3 below the values obtained previously from observations of atomic deuterium in the local bubble and the Galactic halo. However, similarly low gas-phase deuterium abundances have been inferred previously for molecular gas clouds in the Orion region, and in atomic clouds along sight-lines within the Galactic disk to stars more distant than 500 pc from the Sun. △ Less

Submitted 26 June, 2006; originally announced June 2006.

Comments: 12 pages, including 5 figures. Accepted for publication in ApJ Letters

Journal ref: Astrophys.J.647:L33-L36,2006

First detection of [CII]158um at high redshift: vigorous star formation in the early universe

Authors: R. Maiolino, P. Cox, P. Caselli, A. Beelen, F. Bertoldi, C. L. Carilli, M. J. Kaufman, K. M. Menten, T. Nagao, A. Omont, A. Weiss, C. M. Walmsley, F. Walter

Abstract: We report the detection of the 2P_3/2 -> 2P_1/2 fine-structure line of C+ at 157.74 micron in SDSSJ114816.64+525150.3 (hereafter J1148+5251), the most distant known quasar, at z=6.42, using the IRAM 30-meter telescope. This is the first detection of the [CII] line at high redshift, and also the first detection in a Hyperluminous Infrared Galaxy (L_FIR > 10^13 Lsun). The [CII] line is detected at… ▽ More We report the detection of the 2P_3/2 -> 2P_1/2 fine-structure line of C+ at 157.74 micron in SDSSJ114816.64+525150.3 (hereafter J1148+5251), the most distant known quasar, at z=6.42, using the IRAM 30-meter telescope. This is the first detection of the [CII] line at high redshift, and also the first detection in a Hyperluminous Infrared Galaxy (L_FIR > 10^13 Lsun). The [CII] line is detected at a significance level of 8 sigma and has a luminosity of 4.4 x 10^9 Lsun. The L_[CII]/L_FIR ratio is 2 x 10^-4, about an order of magnitude smaller than observed in local normal galaxies and similar to the ratio observed in local Ultraluminous Infrared Galaxies. The [CII] line luminosity indicates that the host galaxy of this quasar is undergoing an intense burst of star formation with an estimated rate of ~3000 Msun/yr. The detection of C+ in SDSS J1148+5251 suggests a significant enrichment of metals at z ~ 6 (age of the universe ~870 Myr), although the data are consistent with a reduced carbon to oxygen ratio as expected from chemical evolutionary models of the early phases of galaxy formation. △ Less

Submitted 2 August, 2005; originally announced August 2005.

Comments: 5 pages, 2 figures, accepted by A&A Letters

High Resolution Observations of the Near-Infrared Emission from NGC 6822 Hubble V

Authors: Sungho Lee, Soojong Pak, Sang-Gak Lee, Christopher J. Davis, Michael J. Kaufman, Kenji Mochizuki, Daniel T. Jaffe

Abstract: We observed Hubble V, the brightest HII region complex in the dwarf irregular galaxy NGC 6822, at near-infrared (1.8--2.4 um) wavelengths using the Cooled Grating Spectrometer 4 (CGS4) at the United Kingdom Infra-Red Telescope (UKIRT). The line emission maps of Hubble V show the typical structure of a photo-dissociation region (PDR) where an ionized core, traced by compact He I emission (2.0587… ▽ More We observed Hubble V, the brightest HII region complex in the dwarf irregular galaxy NGC 6822, at near-infrared (1.8--2.4 um) wavelengths using the Cooled Grating Spectrometer 4 (CGS4) at the United Kingdom Infra-Red Telescope (UKIRT). The line emission maps of Hubble V show the typical structure of a photo-dissociation region (PDR) where an ionized core, traced by compact He I emission (2.0587 um) and Br_gamma emission (2.1661 um), is surrounded by an outer layer traced by molecular hydrogen (H_2) emission. The measured line ratios of H_2 2-1 S(1) (2.2477 um) / 1-0 S(1) (2.1218 um) from 0.2 to 0.6 and the un-shifted and un-resolved line profiles suggest that the H_2 emission originates purely from a PDR. By comparing the H_2 results with a PDR model, we conclude that Hubble V includes dense (10^{4.5} cm^{-3}) and warm PDRs. In this environment, most of the H_2 molecules are excited by far-UV photons (with a field strength of 10^{2-4} times that of the average interstellar field), although collisional processes de-excite H_2 and contribute significantly to the excitation of the first vibrational level. We expect that Hubble V is in the early stage of molecular cloud dissolution. △ Less

Submitted 7 June, 2005; v1 submitted 24 October, 2004; originally announced October 2004.

Comments: 9 pages, 3 figures, MNRAS in press, Replacement: much condensation in the length but small change in the conclusions

Journal ref: Mon.Not.Roy.Astron.Soc. 361 (2005) 1273-1280

The Production of HI in Photodissociation Regions and A Comparison with CO(1-0) Emission

Authors: Ronald J. Allen, Harold I. Heaton, Michael J. Kaufman

Abstract: The gas at the surfaces of molecular clouds in galaxies is heated and dissociated by photons from young stars both near and far. HI resulting from the dissociation of molecular hydrogen H2 emits hyperfine line emission at 21 cm, and warmed CO emits dipole rotational lines such as the 2.6 mm line of CO(1-0). We use previously developed models for photodissociation regions (PDRs) to compute the in… ▽ More The gas at the surfaces of molecular clouds in galaxies is heated and dissociated by photons from young stars both near and far. HI resulting from the dissociation of molecular hydrogen H2 emits hyperfine line emission at 21 cm, and warmed CO emits dipole rotational lines such as the 2.6 mm line of CO(1-0). We use previously developed models for photodissociation regions (PDRs) to compute the intensities of these HI and CO(1-0) lines as a function of the total volume density n in the cloud and the far ultraviolet flux G0 incident upon it and present the results in units familiar to observers. The intensities of these two lines behave differently with changing physical conditions in the PDR, and, taken together, the two lines can provide a ground-based radio astronomy diagnostic for determining n and G0 separately in distant molecular clouds. This diagnostic is particularly useful in the range Gzero <~ 100, 10 cm^{-3} <~ n <~ 10^5 cm^{-3}, which applies to a large fraction of the volume of the interstellar medium in galaxies. If the molecular cloud is located near discrete sources of far-UV (FUV) emission, the PDR-generated HI and CO(1-0) emission on the cloud surface can be more easily identified, appearing as layered ``blankets'' or ``blisters'' on the side of the cloud nearest to the FUV source. As an illustration, we consider the Galactic object G216 -2.5, i.e. ``Maddalena's Cloud'', which has been previously identified as a large PDR in the Galaxy. We determine that this cloud has n ~ 200 cm^{-3}, G0 ~ 0.8, consistent with other data. △ Less

Submitted 16 February, 2004; originally announced February 2004.

Comments: 13 Pages, 3 Figures. Accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J. 608 (2004) 306-313

Water Absorption From Line-of-Sight Clouds Toward W49A

Authors: Rene Plume, Michael J. Kaufman, David A. Neufeld, Ronald L. Snell, David J. Hollenbach, Paul F. Goldsmith, John Howe, Edwin A. Bergin, Gary J. Melnick, Frank Bensch

Abstract: We have observed 6 clouds along the line-of-sight toward W49A using the Submillimeter Wave Astronomy Satellite (SWAS) and several ground-based observatories. The ortho-H2O 1-0 and OH (1665 and 1667 MHz) transitions are observed in absorption, whereas the low-J CO, 13CO, and C18O lines, as well as the [CI] 1-0 transition, are seen in emission. By using both the o-H218O and o-H2O absorption lines,… ▽ More We have observed 6 clouds along the line-of-sight toward W49A using the Submillimeter Wave Astronomy Satellite (SWAS) and several ground-based observatories. The ortho-H2O 1-0 and OH (1665 and 1667 MHz) transitions are observed in absorption, whereas the low-J CO, 13CO, and C18O lines, as well as the [CI] 1-0 transition, are seen in emission. By using both the o-H218O and o-H2O absorption lines, we are able to constrain the column-averaged o-H_2O abundances in each line-of-sight cloud to within about an order of magnitude. Assuming the standard N(H2)/N(CO) ratio of 10^4, we find N(o-H2O)/N(H2) = 8.1 x 10^-8 - 4 x 10^-7 for three clouds with optically thin water lines. In three additional clouds, the H$_2$O lines are saturated so we have used observations of the H218O ground-state transition to find upper limits to the water abundance of 8.2x 10^-8 - 1.5x10^-6. We measure the OH abundance from the average of the 1665 and 1667 MHz observations and find N(OH)/N(H2) = 2.3x10^-7 - 1.1x10^-6. The o-H2O and OH abundances are similar to those determined for line-of-sight water absorption features towards W51 and Sgr B2 but are higher than those seen from water emission lines in molecular clouds. However, the clouds towards W49 have lower ratios of OH relative to H2O column densities than are predicted by simple models which assume that dissociative recombination is the primary formation pathway for OH and H2O. Building on the work of Neufeld et al. (2002), we present photo-chemistry models including additional chemical effects, which can also explain the observed OH and H2O column densities as well as the observed H2O/CO abundance ratios. △ Less

Submitted 22 January, 2004; originally announced January 2004.

Comments: 32 pages, 7 figures, To appear in ApJ April 10 issue

Journal ref: Astrophys.J. 605 (2004) 247-258

A Survey of 557 GHz Water Vapor Emission in the NGC 1333 Molecular Cloud

Authors: Edwin A. Bergin, Michael J. Kaufman, Gary J. Melnick, Ronald L. Snell, John E. Howe

Abstract: Using NASA's Submillimeter Wave Astronomy Satellite (SWAS) we have examined the production of water in quiescent and shocked molecular gas through a survey of the 556.936 GHz transition of ortho-H2O in the NGC 1333 molecular core. These observations reveal broad emission lines associated with the IRAS~2, IRAS~4, IRAS~7, and HH7-11 outflows. Towards 3 positions we detect narrow (~2-3 km/s) emissi… ▽ More Using NASA's Submillimeter Wave Astronomy Satellite (SWAS) we have examined the production of water in quiescent and shocked molecular gas through a survey of the 556.936 GHz transition of ortho-H2O in the NGC 1333 molecular core. These observations reveal broad emission lines associated with the IRAS~2, IRAS~4, IRAS~7, and HH7-11 outflows. Towards 3 positions we detect narrow (~2-3 km/s) emission lines clearly associated with the ambient gas. The SWAS observations, with a resolution of 4', are supplemented with observations from the Infrared Space Observatory (ISO) and by an unbiased survey of a 17' x 15' area, with 50" resolution, in the low-J transitions of CO, 13CO, C18O, N2H+, CH3OH, and SiO. Using these combined data sets, with consistent assumptions, we find beam-averaged o-H2O abundances of > 10^{-6} relative to H2O for all four outflows. A comparison of SWAS and ISO water data is consistent with non-dissociative shock models, provided the majority of the 557 GHz emission arises from cool post-shock material with enhanced abundances. In the ambient gas the o-H2O abundance is found to lie between 0.1-1 x 10^{-7} relative to H2 and is enhanced when compared to cold pre-stellar molecular cores. A comparison of the water emission with tracers of dense condensations and shock chemistry finds no clear correlation. However, the water emission appears to be associated with the presence of luminous external heating sources which power the reflection nebula and the photodissociation (PDR) region. Simple PDR models are capable of reproducing the water and high-J 13CO emission, suggesting that a PDR may account for the excitation of water in low density undepleted gas as suggested by Spaans & van Dishoeck (2001). △ Less

Submitted 9 September, 2002; originally announced September 2002.

Comments: 36 pages, 8 figures (2 color), accepted by ApJ

Journal ref: Astrophys.J. 582 (2003) 830-845

SWAS and Arecibo observations of H2O and OH in a diffuse cloud along the line-of-sight to W51

Authors: David A. Neufeld, Michael J. Kaufman, Paul F. Goldsmith, David J. Hollenbach, Rene Plume

Abstract: Observations of W51 with the Submillimeter Wave Astronomy Satellite (SWAS) have yielded the first detection of water vapor in a diffuse molecular cloud. The water vapor lies in a foreground cloud that gives rise to an absorption feature at an LSR velocity of 6 km/s. The inferred H2O column density is 2.5E+13 cm-2. Observations with the Arecibo radio telescope of hydroxyl molecules at ten positio… ▽ More Observations of W51 with the Submillimeter Wave Astronomy Satellite (SWAS) have yielded the first detection of water vapor in a diffuse molecular cloud. The water vapor lies in a foreground cloud that gives rise to an absorption feature at an LSR velocity of 6 km/s. The inferred H2O column density is 2.5E+13 cm-2. Observations with the Arecibo radio telescope of hydroxyl molecules at ten positions in W51 imply an OH column density of 8E+13 cm-2 in the same diffuse cloud. The observed H2O/OH ratio of ~ 0.3 is significantly larger than an upper limit derived previously from ultraviolet observations of the similar diffuse molecular cloud lying in front of HD 154368. The observed variation in H2O/OH likely points to the presence in one or both of these clouds of a warm (T > 400) gas component in which neutral-neutral reactions are important sources of OH and/or H2O. △ Less

Submitted 16 July, 2002; v1 submitted 16 July, 2002; originally announced July 2002.

Comments: 15 pages (AASTeX) including 4 (eps) figures. To appear in the Astrophysical Journal

Journal ref: Astrophys.J. 580 (2002) 278-284

ISO--LWS observations of the two nearby spiral galaxies: NGC6946 and NGC1313

Authors: A. Contursi, M. J. Kaufman, G. Helou, D. J. Hollenbach, J. Brauher, G. J. Stacey, D. A. Dale, S. Malhotra, M. Rubio, R. H. Rubin, S. D. Lord

Abstract: (Abridged) We present the analysis of the main FIR fine structure lines emission in NGC1313 and NGC6946. We calculate that a component probably associated with the diffuse disks contributes <~40% in N6946 and ~30 % in N1313 to the total [CII] emission. The main PDR physical parameters responsible for the neutral atomic gas emission in N1313 and N6946 do not significantly differ from what Malhotr… ▽ More (Abridged) We present the analysis of the main FIR fine structure lines emission in NGC1313 and NGC6946. We calculate that a component probably associated with the diffuse disks contributes <~40% in N6946 and ~30 % in N1313 to the total [CII] emission. The main PDR physical parameters responsible for the neutral atomic gas emission in N1313 and N6946 do not significantly differ from what Malhotra etal (2001) found by modelling the integrated emission of a sample of 60 normal galaxies,although there are evidences for a beam averaged contribution of a less active component inside NGC6946 higher than its contribution in the integrated emission of normal galaxies. CO and [CII] in N6946 are well correlated with a mean [CII]/CO ratio similar to that of the normal galaxies sample. In N1313 the [CII]/CO seems to systematically increase from the North to the south, along the S-shaped spiral arm, indicating much more inhomogeneous conditions than in N6946. HI and [CII] in N6946 are completely de-correlated, probably because they arise from different gas components. In N1313 we successfully detect two distinct gas components: a cirrus-like component where HI and [CII] are weakly correlated as observed in our Galaxy, and a component associated with dense PDRs completely de-correlated from HI as observed in N6946.Finally, we find that the HI residing in dense PDRs and presumably recently photo-dissociated, constitutes a few % of the total HI. In turn, this dense gas component produces most of the [CII] emission emitted by the atomic neutral medium. △ Less

Submitted 15 April, 2002; originally announced April 2002.

Comments: Latex, 100 pages, 11 Figures, 11 Tables. Accepted for publication in AJ

Far Infrared Spectroscopy of Normal Galaxies: Physical Conditions in the Interstellar Medium

Authors: S. Malhotra, M. J. Kaufman, D. Hollenbach, G. Helou, R. H. Rubin, J. Brauher, D. Dale, N. Y. Lu, S. Lord, G. Stacey, A. Contursi

Abstract: The most important cooling lines of the neutral interstellar medium (ISM) lie in the far-infrared (FIR). We present measurements by the Infrared Space Observatory Long Wavelength Spectrometer of seven lines from neutral and ionized ISM of 60 normal, star-forming galaxies. The galaxy sample spans a range in properties such as morphology, FIR colors (indicating dust temperature), and FIR/Blue rati… ▽ More The most important cooling lines of the neutral interstellar medium (ISM) lie in the far-infrared (FIR). We present measurements by the Infrared Space Observatory Long Wavelength Spectrometer of seven lines from neutral and ionized ISM of 60 normal, star-forming galaxies. The galaxy sample spans a range in properties such as morphology, FIR colors (indicating dust temperature), and FIR/Blue ratios (indicating star-formation activity and optical depth). In two-thirds of the galaxies in this sample, the [CII] line is proportional to FIR dust continuum. The other one-third show a smooth decline in [CII]/FIR with increasing F60/F100 and FIR/B, spanning a range of a factor of more than 50. Two galaxies, at the warm and active extreme of the range have [CII]/FIR < 2 \times 10^{-4} (3-sigma upper limit). This is due to increased positive grain charge in the warmer and more active galaxies, which leads to less efficient heating by photoelectrons from dust grains. The ratio of the two principal photodissociation region (PDR) cooling lines [CII]/[OI] shows a tight correlation with F60/F100, indicating that both gas and dust temperatures increase together. We derive a theoretical scaling between [NII] and [CII] from ionized gas and use it to separate [CII] emission from neutral PDRs and ionized gas. Comparison of PDR models of Kaufman et al. (1999) with observed ratios of (a) [OI]/[CII] and ([CII]+[OI])/FIR and (b) [OI]/FIR and F60/F100 yields far-UV flux G0 and gas density n. The derived G0 scales as n to the power 1.4. We interpret this correlation as arising from Stromgren sphere scalings if much of the line and continuum luminosity arises near star-forming regions. The differences in G0 and n may be due to differences in the physical properties of the star-forming clouds.(Short abstract) △ Less

Submitted 27 June, 2001; originally announced June 2001.

Comments: Accepted by the Astrophysical Journal

Journal ref: Astrophys.J.561:766-786,2001

Extended [C I] and ^{13}CO(5-4) Emission in M17SW

Authors: J. E. Howe, M. L. N. Ashby, E. A. Bergin, G. Chin, N. R. Erickson, P. F. Goldsmith, M. Harwit, D. J. Hollenbach, M. J. Kaufman, S. C. Kleiner, D. G. Koch, D. A. Neufeld, B. M. Patten, R. Plume, R. Schieder, R. L. Snell, J. R. Stauffer, V. Tolls, Z. Wang, G. Winnewisser, Y. F. Zhang, G. J. Melnick

Abstract: We mapped a 13 by 22 pc region in emission from 492 GHz [C I] and, for the first time, 551 GHz ^{13}CO(5-4) in the giant molecular cloud M17SW, using the Submillimeter Wave Astronomy Satellite. The morphologies of the [C I] and ^{13}CO emission are strikingly similar. The extent and intensity of the [C I] and ^{13}CO(5-4) emission is explained as arising from photodissociation regions on the sur… ▽ More We mapped a 13 by 22 pc region in emission from 492 GHz [C I] and, for the first time, 551 GHz ^{13}CO(5-4) in the giant molecular cloud M17SW, using the Submillimeter Wave Astronomy Satellite. The morphologies of the [C I] and ^{13}CO emission are strikingly similar. The extent and intensity of the [C I] and ^{13}CO(5-4) emission is explained as arising from photodissociation regions on the surfaces of embedded molecular clumps. Modeling of the ^{13}CO(5-4) emission in comparison to ^{13}CO(1-0) indicates a temperature gradient across the cloud, peaking to at least 63 K near the M17 ionization front and decreasing to at least 20 K at the western edge of the cloud. We see no correlation between gas density and column density. The beam-averaged column density of C I in the core is 1x10^{18} cm^-2, and the mean column density ratio N(C I)/N(CO) is about 0.4. The variations of N(C I)/N(CO) with position in M17SW indicate a similar clump size distribution throughout the cloud. △ Less

Submitted 19 October, 2000; originally announced October 2000.

Comments: 5 pages, 3 Postscript figures, uses aastex.cls, emulateapj5.sty (included), and apjfonts.sty (included)

Journal ref: ApJ (Letters), 539, L137 (2000)

Far Infrared and Submillimeter Emission from Galactic and Extragalactic Photo-Dissociation Regions

Authors: Michael J. Kaufman, Mark G. Wolfire, David J. Hollenbach, Michael L. Luhman

Abstract: Photodissociation Region (PDR) models are computed over a wide range of physical conditions, from those appropriate to giant molecular clouds illuminated by the interstellar radiation field to the conditions experienced by circumstellar disks very close to hot massive stars. These models use the most up-to-date values of atomic and molecular data, the most current chemical rate coefficients, and… ▽ More Photodissociation Region (PDR) models are computed over a wide range of physical conditions, from those appropriate to giant molecular clouds illuminated by the interstellar radiation field to the conditions experienced by circumstellar disks very close to hot massive stars. These models use the most up-to-date values of atomic and molecular data, the most current chemical rate coefficients, and the newest grain photoelectric heating rates which include treatments of small grains and large molecules. In addition, we examine the effects of metallicity and cloud extinction on the predicted line intensities. Results are presented for PDR models with densities over the range n=10^1-10^7 cm^-3 and for incident far-ultraviolet radiation fields over the range G_0=10^-0.5-10^6.5, for metallicities Z=1 and 0.1 times the local Galactic value, and for a range of PDR cloud sizes. We present line strength and/or line ratio plots for a variety of useful PDR diagnostics: [C II] 158 micron, [O I] 63 and 145 micron, [C I] 370 and 609 micron, CO J=1-0, J=2-1, J=3-2, J=6-5 and J=15-14, as well as the strength of the far-infrared continuum. These plots will be useful for the interpretation of Galactic and extragalactic far infrared and submillimeter spectra observable with ISO, SOFIA, SWAS, FIRST and other orbital and suborbital platforms. As examples, we apply our results to ISO and ground based observations of M82, NGC 278, and the Large Magellenic Cloud. △ Less

Submitted 19 July, 1999; originally announced July 1999.

Comments: 54 pages, 20 figures, accepted for publication in The Astrophysical Journal

Journal ref: Astrophys.J.527:795-813,1999

Thermal Water Vapor Emission from Shocked Regions in Orion

Authors: Martin Harwit, David A. Neufeld, Gary J. Melnick, Michael J. Kaufman

Abstract: Using the Long Wavelength Spectrometer (LWS) onboard the Infrared Space Observatory (ISO), we have observed thermal water vapor emission from a roughly circular field of view approximately 75 arc seconds in diameter centered on the Orion BN-KL region. The Fabry-Perot line strengths, line widths, and spectral line shifts observed in eight transitions between 71 and 125 microns show good agreement… ▽ More Using the Long Wavelength Spectrometer (LWS) onboard the Infrared Space Observatory (ISO), we have observed thermal water vapor emission from a roughly circular field of view approximately 75 arc seconds in diameter centered on the Orion BN-KL region. The Fabry-Perot line strengths, line widths, and spectral line shifts observed in eight transitions between 71 and 125 microns show good agreement with models of thermal emission arising from a molecular cloud subjected to a magnetohydrodynamic C-type shock. Both the breadth and the relative strengths of the observed lines argue for emission from a shock rather than from warm quiescent gas in the Orion core. Though one of the eight transitions appears anomalously strong, and may be subject to the effects of radiative pumping, the other seven indicate an H2O/H2 abundance ratio of order 5E-4, and a corresponding gas-phase oxygen-to-hydrogen abundance ratio of order 4E-4. Given current estimates of the interstellar, gas-phase, oxygen and carbon abundances in the solar vicinity, this value is consistent with theoretical shock models that predict the conversion into water of all the gas-phase oxygen that is not bound as CO. The overall cooling provided by rotational transitions of H2O in this region appears to be comparable to the cooling through rotational lines of CO, but is an order of magnitude lower than cooling through H2 emission. However, the model that best fits our observations shows cooling by H2O and CO dominant in that portion of the post-shock region where temperatures are below ~ 800 K and neither vibrational nor rotational radiative cooling by H2 is appreciable. △ Less

Submitted 26 February, 1998; originally announced February 1998.

Comments: 13 pages, 1 postscript figure, LaTeX, uses aaspp4.sty, accepted for publication in ApJ Letters