Showing 1–50 of 109 results for author: Wolfire, M

PDRs4All IX. Sulfur elemental abundance in the Orion Bar

Authors: Asunción Fuente, Evelyne Roueff, Franck Le Petit, Jacques Le Bourlot, Emeric Bron, Mark G. Wolfire, James F. Babb, Pei-Gen Yan, Takashi Onaka, John H. Black, Ilane Schroetter, Dries Van De Putte, Ameek Sidhu, Amélie Canin, Boris Trahin, Felipe Alarcón, Ryan Chown, Olga Kannavou, Olivier Berné, Emilie Habart, Els Peeters, Javier R. Goicoechea, Marion Zannese, Raphael Meshaka, Yoko Okada , et al. (9 additional authors not shown)

Abstract: One of the main problems in astrochemistry is determining the amount of sulfur in volatiles and refractories in the interstellar medium. The detection of the main sulfur reservoirs (icy H$_2$S and atomic gas) has been challenging, and estimates are based on the reliability of models to account for the abundances of species containing less than 1% of the total sulfur. The high sensitivity of the Ja… ▽ More One of the main problems in astrochemistry is determining the amount of sulfur in volatiles and refractories in the interstellar medium. The detection of the main sulfur reservoirs (icy H$_2$S and atomic gas) has been challenging, and estimates are based on the reliability of models to account for the abundances of species containing less than 1% of the total sulfur. The high sensitivity of the James Webb Space Telescope provides an unprecedented opportunity to estimate the sulfur abundance through the observation of the [S I] 25.249 $μ$m line. We used the [S III] 18.7 $μ$m, [S IV] 10.5 $μ$m, and [S l] 25.249 $μ$m lines to estimate the amount of sulfur in the ionized and molecular gas along the Orion Bar. For the theoretical part, we used an upgraded version of the Meudon photodissociation region (PDR) code to model the observations. New inelastic collision rates of neutral atomic sulfur with ortho- and para- molecular hydrogen were calculated to predict the line intensities. The [S III] 18.7 $μ$m and [S IV] 10.5 $μ$m lines are detected over the imaged region with a shallow increase (by a factor of 4) toward the HII region. We estimate a moderate sulfur depletion, by a factor of $\sim$2, in the ionized gas. The corrugated interface between the molecular and atomic phases gives rise to several edge-on dissociation fronts we refer to as DF1, DF2, and DF3. The [S l] 25.249 $μ$m line is only detected toward DF2 and DF3, the dissociation fronts located farthest from the HII region. The detailed modeling of DF3 using the Meudon PDR code shows that the emission of the [S l] 25.249 $μ$m line is coming from warm ($>$ 40 K) molecular gas located at A$_{\rm V}$ $\sim$ 1$-$5 mag from the ionization front. Moreover, the intensity of the [S l] 25.249 $μ$m line is only accounted for if we assume the presence of undepleted sulfur. △ Less

Submitted 4 June, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: 16 pages, 6 figures. Accepted for publication in Astronomy and Astrophysics

Modeling Ionized Gas in the Small Magellanic Cloud: The Wolf-Rayet Nebula N76

Authors: Elizabeth Tarantino, Alberto D. Bolatto, Rémy Indebetouw, Mónica Rubio, Karin M. Sandstrom, J. -D T. Smith, Daniel Stapleton, Mark Wolfire

Abstract: We present Cloudy modeling of infrared emission lines in the Wolf-Rayet (WR) nebula N76 caused by one of the most luminous and hottest WR stars in the low metallicity Small Magellanic Cloud. We use spatially resolved mid-infrared Spitzer/IRS and far-infrared Herschel/PACS spectroscopy to establish the physical conditions of the ionized gas. The spatially resolved distribution of the emission allow… ▽ More We present Cloudy modeling of infrared emission lines in the Wolf-Rayet (WR) nebula N76 caused by one of the most luminous and hottest WR stars in the low metallicity Small Magellanic Cloud. We use spatially resolved mid-infrared Spitzer/IRS and far-infrared Herschel/PACS spectroscopy to establish the physical conditions of the ionized gas. The spatially resolved distribution of the emission allows us to constrain properties much more accurately than using spatially integrated quantities. We construct models with a range of constant hydrogen densities between n$_H$ = 4 - 10 cm$^{-3}$ and a stellar wind-blown cavity of 10 pc which reproduces the intensity and shape of most ionized gas emission lines, including the high ionization lines [OIV] and [NeV], as well as [SIII], [SIV], [OIII], and [NeIII]. Our models suggest that the majority of [SiII] emission (91%) is produced at the edge of the HII region around the transition between ionized and atomic gas while very little of the [CII] emission (<5%) is associated with the ionized gas. The physical conditions of N76 are characterized by a hot HII region with a maximum electron temperature of T$_e$ ~ 24,000 K, electron densities that range from n$_e$ ~ 4 to 12 cm$^{-3}$, and high ionization parameters of log(U) ~ -1.15 to -1.77. By analyzing a low metallicty WR nebula with a single ionization source, this work gives valuable insights on the impact WR stars have to the galaxy-integrated ionized gas properties in nearby dwarf galaxies. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 27 pages, 10 figures, 2 tables. Accepted in the Astrophysical Journal

PDRs4All VIII: Mid-IR emission line inventory of the Orion Bar

Authors: Dries Van De Putte, Raphael Meshaka, Boris Trahin, Emilie Habart, Els Peeters, Olivier Berné, Felipe Alarcón, Amélie Canin, Ryan Chown, Ilane Schroetter, Ameek Sidhu, Christiaan Boersma, Emeric Bron, Emmanuel Dartois, Javier R. Goicoechea, Karl D. Gordon, Takashi Onaka, Alexander G. G. M. Tielens, Laurent Verstraete, Mark G. Wolfire, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Jan Cami, Sara Cuadrado , et al. (113 additional authors not shown)

Abstract: Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observat… ▽ More Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observational inventory of mid-IR emission lines, and spatially resolve the substructure of the PDR, with a mosaic cutting perpendicularly across the ionization front and three dissociation fronts. We extracted five spectra that represent the ionized, atomic, and molecular gas layers, and measured the most prominent gas emission lines. An initial analysis summarizes the physical conditions of the gas and the potential of these data. We identified around 100 lines, report an additional 18 lines that remain unidentified, and measured the line intensities and central wavelengths. The H I recombination lines originating from the ionized gas layer bordering the PDR, have intensity ratios that are well matched by emissivity coefficients from H recombination theory, but deviate up to 10% due contamination by He I lines. We report the observed emission lines of various ionization stages of Ne, P, S, Cl, Ar, Fe, and Ni, and show how certain line ratios vary between the five regions. We observe the pure-rotational H$_2$ lines in the vibrational ground state from 0-0 S(1) to 0-0 S(8), and in the first vibrationally excited state from 1-1 S(5) to 1-1 S(9). We derive H$_2$ excitation diagrams, and approximate the excitation with one thermal (~700 K) component representative of an average gas temperature, and one non-thermal component (~2700 K) probing the effect of UV pumping. We compare these results to an existing model for the Orion Bar PDR and highlight the differences with the observations. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 26 pages, 12 figures, 3 tables. Submitted to A&A, under review (1st revision)

Journal ref: A&A 687, A86 (2024)

A far-ultraviolet-driven photoevaporation flow observed in a protoplanetary disk

Authors: Olivier Berné, Emilie Habart, Els Peeters, Ilane Schroetter, Amélie Canin, Ameek Sidhu, Ryan Chown, Emeric Bron, Thomas J. Haworth, Pamela Klaassen, Boris Trahin, Dries Van De Putte, Felipe Alarcón, Marion Zannese, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Jan Cami, Sara Cuadrado, Emmanuel Dartois, Daniel Dicken, Meriem Elyajouri, Asunción Fuente, Javier R. Goicoechea , et al. (121 additional authors not shown)

Abstract: Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of… ▽ More Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modelling their kinematics and excitation allows us to constrain the physical conditions within the gas. We quantify the mass-loss rate induced by the FUV irradiation, finding it is sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant planet formation in the disk. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Journal ref: Science, 383, 6686, 2024

Formation of the Methyl Cation by Photochemistry in a Protoplanetary Disk

Authors: Olivier Berné, Marie-Aline Martin-Drumel, Ilane Schroetter, Javier R. Goicoechea, Ugo Jacovella, Bérenger Gans, Emmanuel Dartois, Laurent Coudert, Edwin Bergin, Felipe Alarcon, Jan Cami, Evelyne Roueff, John H. Black, Oskar Asvany, Emilie Habart, Els Peeters, Amelie Canin, Boris Trahin, Christine Joblin, Stephan Schlemmer, Sven Thorwirth, Jose Cernicharo, Maryvonne Gerin, Alexander Tielens, Marion Zannese , et al. (31 additional authors not shown)

Abstract: Forty years ago it was proposed that gas phase organic chemistry in the interstellar medium was initiated by the methyl cation CH3+, but hitherto it has not been observed outside the Solar System. Alternative routes involving processes on grain surfaces have been invoked. Here we report JWST observations of CH3+ in a protoplanetary disk in the Orion star forming region. We find that gas-phase orga… ▽ More Forty years ago it was proposed that gas phase organic chemistry in the interstellar medium was initiated by the methyl cation CH3+, but hitherto it has not been observed outside the Solar System. Alternative routes involving processes on grain surfaces have been invoked. Here we report JWST observations of CH3+ in a protoplanetary disk in the Orion star forming region. We find that gas-phase organic chemistry is activated by UV irradiation. △ Less

Submitted 6 January, 2024; originally announced January 2024.

Comments: Published in Nature

Journal ref: Nature 621, 56-59 (2023)

OH as a probe of the warm water cycle in planet-forming disks

Authors: Marion Zannese, Benoît Tabone, Emilie Habart, Javier R. Goicoechea, Alexandre Zanchet, Ewine F. van Dishoeck, Marc C. van Hemert, John H. Black, Alexander G. G. M. Tielens, A. Veselinova, P. G. Jambrina, M. Menendez, E. Verdasco, F. J. Aoiz, L. Gonzalez-Sanchez, Boris Trahin, Emmanuel Dartois, Olivier Berné, Els Peeters, Jinhua He, Ameek Sidhu, Ryan Chown, Ilane Schroetter, Dries Van De Putte, Amélie Canin , et al. (30 additional authors not shown)

Abstract: Water is a key ingredient for the emergence of life as we know it. Yet, its destruction and reformation in space remains unprobed in warm gas. Here, we detect the hydroxyl radical (OH) emission from a planet-forming disk exposed to external far-ultraviolet (FUV) radiation with the James Webb Space Telescope. The observations are confronted with the results of quantum dynamical calculations. The hi… ▽ More Water is a key ingredient for the emergence of life as we know it. Yet, its destruction and reformation in space remains unprobed in warm gas. Here, we detect the hydroxyl radical (OH) emission from a planet-forming disk exposed to external far-ultraviolet (FUV) radiation with the James Webb Space Telescope. The observations are confronted with the results of quantum dynamical calculations. The highly excited OH infrared rotational lines are the tell-tale signs of H2O destruction by FUV. The OH infrared ro-vibrational lines are attributed to chemical excitation via the key reaction O+H=OH+H which seeds the formation of water in the gas-phase. We infer that the equivalent of the Earth ocean's worth of water is destroyed per month and replenished. These results show that under warm and irradiated conditions water is destroyed and efficiently reformed via gas-phase reactions. This process, assisted by diffusive transport, could reduce the HDO/H2O ratio in the warm regions of planet-forming disks. △ Less

Submitted 22 December, 2023; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: Version submitted to Nature Astronomy

PDRs4All III: JWST's NIR spectroscopic view of the Orion Bar

Authors: Els Peeters, Emilie Habart, Olivier Berne, Ameek Sidhu, Ryan Chown, Dries Van De Putte, Boris Trahin, Ilane Schroetter, Amelie Canin, Felipe Alarcon, Bethany Schefter, Baria Khan, Sofia Pasquini, Alexander G. G. M. Tielens, Mark G. Wolfire, Emmanuel Dartois, Javier R. Goicoechea, Alexandros Maragkoudakis, Takashi Onaka, Marc W. Pound, Silvia Vicente, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma , et al. (113 additional authors not shown)

Abstract: (Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion… ▽ More (Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion Bar PDR as part of the PDRs4All JWST Early Release Science Program. The NIRSpec data reveal a forest of lines including, but not limited to, HeI, HI, and CI recombination lines, ionic lines, OI and NI fluorescence lines, Aromatic Infrared Bands (AIBs including aromatic CH, aliphatic CH, and their CD counterparts), CO2 ice, pure rotational and ro-vibrational lines from H2, and ro-vibrational lines HD, CO, and CH+, most of them detected for the first time towards a PDR. Their spatial distribution resolves the H and He ionisation structure in the Huygens region, gives insight into the geometry of the Bar, and confirms the large-scale stratification of PDRs. We observe numerous smaller scale structures whose typical size decreases with distance from Ori C and IR lines from CI, if solely arising from radiative recombination and cascade, reveal very high gas temperatures consistent with the hot irradiated surface of small-scale dense clumps deep inside the PDR. The H2 lines reveal multiple, prominent filaments which exhibit different characteristics. This leaves the impression of a "terraced" transition from the predominantly atomic surface region to the CO-rich molecular zone deeper in. This study showcases the discovery space created by JWST to further our understanding of the impact radiation from young stars has on their natal molecular cloud and proto-planetary disk, which touches on star- and planet formation as well as galaxy evolution. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 52 pages, 30 figures, submitted to A&A

Journal ref: A&A 685, A74 (2024)

Identifying physical structures in our Galaxy with Gaussian Mixture Models: An unsupervised machine learning technique

Authors: M. Tiwari, R. Kievit, S. Kabanovic, L. Bonne, F. Falasca, C. Guevara, R. Higgins, M. Justen, R. Karim, Ü. Kavak, C. Pabst, M. W. Pound, N. Schneider, R. Simon, J. Stutzki, M. Wolfire, A. G. G. M. Tielens

Abstract: We explore the potential of the Gaussian Mixture Model (GMM), an unsupervised machine learning method, to identify coherent physical structures in the ISM. The implementation we present can be used on any kind of spatially and spectrally resolved data set. We provide a step-by-step guide to use these models on different sources and data sets. Following the guide, we run the models on NGC 1977, RCW… ▽ More We explore the potential of the Gaussian Mixture Model (GMM), an unsupervised machine learning method, to identify coherent physical structures in the ISM. The implementation we present can be used on any kind of spatially and spectrally resolved data set. We provide a step-by-step guide to use these models on different sources and data sets. Following the guide, we run the models on NGC 1977, RCW 120 and RCW 49 using the [CII] 158 $μ$m mapping observations from the SOFIA telescope. We find that the models identified 6, 4 and 5 velocity coherent physical structures in NGC 1977, RCW 120 and RCW 49, respectively, which are validated by analysing the observed spectra towards these structures and by comparison to earlier findings. In this work we demonstrate that GMM is a powerful tool that can better automate the process of spatial and spectral analysis to interpret mapping observations. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 19 pages, 14 figures

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

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

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

Submitted 29 September, 2023; originally announced September 2023.

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

SOFIA FEEDBACK Survey: The Pillars of Creation in [C II] and Molecular Lines

Authors: Ramsey L. Karim, Marc W. Pound, Alexander G. G. M. Tielens, Maitraiyee Tiwari, Lars Bonne, Mark G. Wolfire, Nicola Schneider, Ümit Kavak, Lee G. Mundy, Robert Simon, Rolf Güsten, Jürgen Stutzki, Friedrich Wyrowski, Netty Honingh

Abstract: We investigate the physical structure and conditions of photodissociation regions (PDRs) and molecular gas within the Pillars of Creation in the Eagle Nebula using SOFIA FEEDBACK observations of the [C II] 158 micron line. These observations are velocity resolved to 0.5 km s$^{-1}$ and are analyzed alongside a collection of complimentary data with similar spatial and spectral resolution: the [O I]… ▽ More We investigate the physical structure and conditions of photodissociation regions (PDRs) and molecular gas within the Pillars of Creation in the Eagle Nebula using SOFIA FEEDBACK observations of the [C II] 158 micron line. These observations are velocity resolved to 0.5 km s$^{-1}$ and are analyzed alongside a collection of complimentary data with similar spatial and spectral resolution: the [O I] 63 micron line, also observed with SOFIA, and rotational lines of CO, HCN, HCO$^{+}$, CS, and N$_2$H$^{+}$. Using the superb spectral resolution of SOFIA, APEX, CARMA, and BIMA, we reveal the relationships between the warm PDR and cool molecular gas layers in context of the Pillars' kinematic structure. We assemble a geometric picture of the Pillars and their surroundings informed by illumination patterns and kinematic relationships and derive physical conditions in the PDRs associated with the Pillars. We estimate an average molecular gas density $n_{{\rm H}_2} \sim 1.3 \times 10^5$ cm$^{-3}$ and an average atomic gas density $n_{\rm H} \sim 1.8 \times 10^4$ cm$^{-3}$ and infer that the ionized, atomic, and molecular phases are in pressure equilibrium if the atomic gas is magnetically supported. We find pillar masses of 103, 78, 103, and 18 solar masses for P1a, P1b, P2, and P3 respectively, and evaporation times of $\sim$1-2 Myr. The dense clumps at the tops of the pillars are currently supported by the magnetic field. Our analysis suggests that ambipolar diffusion is rapid and these clumps are likely to collapse within their photoevaporation timescales. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: 42 pages, 16 figures. Accepted for publication in The Astronomical Journal

PDRs4All IV. An embarrassment of riches: Aromatic infrared bands in the Orion Bar

Authors: Ryan Chown, Ameek Sidhu, Els Peeters, Alexander G. G. M. Tielens, Jan Cami, Olivier Berné, Emilie Habart, Felipe Alarcón, Amélie Canin, Ilane Schroetter, Boris Trahin, Dries Van De Putte, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Emeric Bron, Sara Cuadrado, Emmanuel Dartois, Daniel Dicken, Meriem El-Yajouri, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Lina Issa , et al. (114 additional authors not shown)

Abstract: (Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $μ$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory o… ▽ More (Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $μ$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory of the AIBs found in the Orion Bar, along with mid-IR template spectra from five distinct regions in the Bar: the molecular PDR, the atomic PDR, and the HII region. We use JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program, PDRs4All (ID: 1288). We extract five template spectra to represent the morphology and environment of the Orion Bar PDR. The superb sensitivity and the spectral and spatial resolution of these JWST observations reveal many details of the AIB emission and enable an improved characterization of their detailed profile shapes and sub-components. While the spectra are dominated by the well-known AIBs at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 $μ$m, a wealth of weaker features and sub-components are present. We report trends in the widths and relative strengths of AIBs across the five template spectra. These trends yield valuable insight into the photochemical evolution of PAHs, such as the evolution responsible for the shift of 11.2 $μ$m AIB emission from class B$_{11.2}$ in the molecular PDR to class A$_{11.2}$ in the PDR surface layers. This photochemical evolution is driven by the increased importance of FUV processing in the PDR surface layers, resulting in a "weeding out" of the weakest links of the PAH family in these layers. For now, these JWST observations are consistent with a model in which the underlying PAH family is composed of a few species: the so-called 'grandPAHs'. △ Less

Submitted 5 September, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

Comments: 25 pages, 10 figures, to appear in A&A

Journal ref: A&A 685, A75 (2024)

PDRs4All II: JWST's NIR and MIR imaging view of the Orion Nebula

Authors: Emilie Habart, Els Peeters, Olivier Berné, Boris Trahin, Amélie Canin, Ryan Chown, Ameek Sidhu, Dries Van De Putte, Felipe Alarcón, Ilane Schroetter, Emmanuel Dartois, Sílvia Vicente, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Emeric Bron, Jan Cami, Sara Cuadrado, Daniel Dicken, Meriem Elyajouri, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Lina Issa , et al. (117 additional authors not shown)

Abstract: The JWST has captured the most detailed and sharpest infrared images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). We investigate the fundamental interaction of far-ultraviolet photons with molecular clouds. The transitions across the ionization front (IF), dissociation fron… ▽ More The JWST has captured the most detailed and sharpest infrared images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). We investigate the fundamental interaction of far-ultraviolet photons with molecular clouds. The transitions across the ionization front (IF), dissociation front (DF), and the molecular cloud are studied at high-angular resolution. These transitions are relevant to understanding the effects of radiative feedback from massive stars and the dominant physical and chemical processes that lead to the IR emission that JWST will detect in many Galactic and extragalactic environments. Due to the proximity of the Orion Nebula and the unprecedented angular resolution of JWST, these data reveal that the molecular cloud borders are hyper structured at small angular scales of 0.1-1" (0.0002-0.002 pc or 40-400 au at 414 pc). A diverse set of features are observed such as ridges, waves, globules and photoevaporated protoplanetary disks. At the PDR atomic to molecular transition, several bright features are detected that are associated with the highly irradiated surroundings of the dense molecular condensations and embedded young star. Toward the Orion Bar PDR, a highly sculpted interface is detected with sharp edges and density increases near the IF and DF. This was predicted by previous modeling studies, but the fronts were unresolved in most tracers. A complex, structured, and folded DF surface was traced by the H2 lines. This dataset was used to revisit the commonly adopted 2D PDR structure of the Orion Bar. JWST provides us with a complete view of the PDR, all the way from the PDR edge to the substructured dense region, and this allowed us to determine, in detail, where the emission of the atomic and molecular lines, aromatic bands, and dust originate. △ Less

Submitted 2 September, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

Journal ref: A&A 685, A73 (2024)

The PhotoDissociation Region Toolbox: Software and Models for Astrophysical Analysis

Authors: Marc W. Pound, Mark G. Wolfire

Abstract: The PhotoDissociation Region Toolbox provides comprehensive, easy-to-use, public software tools and models that enable an understanding of the interaction of the light of young, luminous, massive stars with the gas and dust in the Milky Way and in other galaxies. It consists of an open-source Python toolkit and photodissociation region models for analysis of infrared and millimeter/submillimeter l… ▽ More The PhotoDissociation Region Toolbox provides comprehensive, easy-to-use, public software tools and models that enable an understanding of the interaction of the light of young, luminous, massive stars with the gas and dust in the Milky Way and in other galaxies. It consists of an open-source Python toolkit and photodissociation region models for analysis of infrared and millimeter/submillimeter line and continuum observations obtained by ground-based and sub-orbital telescopes, and astrophysics space missions. Photodissociation regions (PDRs) include all of the neutral gas in the ISM where far-ultraviolet photons dominate the chemistry and/or heating. In regions of massive star formation, PDRs are created at the boundaries between the H II regions and neutral molecular cloud, as photons with energies 6 eV $ < h ν< $ 13.6 eV photodissociate molecules and photoionize metals. The gas is heated by photo-electrons from small grains and large molecules and cools mostly through far-infrared fine-structure lines like [O I] and [C II]. The models are created from state-of-the art PDR codes that includes molecular freeze-out; recent collision, chemical, and photo rates; new chemical pathways, such as for oxygen chemistry; and allow for both clumpy and uniform media. The models predict the emergent intensities of many spectral lines and FIR continuum. The tools find the best-fit models to the observations and provide insights into the physical conditions and chemical makeup of the gas and dust. The PDR Toolbox enables novel analysis of data from telescopes such as ISO, Spitzer, Herschel, STO, SOFIA, SWAS, APEX, ALMA, and JWST. △ Less

Submitted 14 October, 2022; originally announced October 2022.

Comments: 22 pages, 10 figures, includes code listings

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

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

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

Submitted 31 August, 2022; originally announced August 2022.

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

SOFIA FEEDBACK survey: PDR diagnostics of stellar feedback in different regions of RCW 49

Authors: M. Tiwari, M. Wolfire, M. W. Pound, E. Tarantino, R. Karim, L. Bonne, C. Buchbender, R. Güsten, C. Guevara, S. Kabanovic, Ü. Kavak, M. Mertens, N. Schneider, R. Simon, J. Stutzki, A. G. G. M. Tielens

Abstract: We quantified the effects of stellar feedback in RCW 49 by determining the physical conditions in different regions using the [CII] 158 $μ$m and [OI] 63 $μ$m observations from SOFIA, the $^{12}$CO (3-2) observations from APEX and the H$_2$ line observations from Spitzer telescopes. Large maps of RCW 49 were observed with the SOFIA and APEX telescopes, while the Spitzer observations were only avail… ▽ More We quantified the effects of stellar feedback in RCW 49 by determining the physical conditions in different regions using the [CII] 158 $μ$m and [OI] 63 $μ$m observations from SOFIA, the $^{12}$CO (3-2) observations from APEX and the H$_2$ line observations from Spitzer telescopes. Large maps of RCW 49 were observed with the SOFIA and APEX telescopes, while the Spitzer observations were only available towards three small areas. From our qualitative analysis, we found that the H$_2$ 0-0 S(2) emission line probes denser gas compared to the H$_2$ 0-0 S(1) line. In four regions ("northern cloud", "pillar", "ridge", and "shell"), we compared our observations with the updated PDR Toolbox models and derived the integrated far-ultraviolet flux between 6-13.6 eV ($G_{\rm 0}$), H nucleus density ($n$), temperatures and pressures. We found the ridge to have the highest $G_{\rm 0}$ (2.4 $\times$ 10$^3$ Habing units), while the northern cloud has the lowest $G_{\rm 0}$ (5 $\times$ 10$^2$ Habing units). This is a direct consequence of the location of these regions with respect to the Wd2 cluster. The ridge also has a high density (6.4 $\times$ 10$^3$ cm$^{-3}$), which is consistent with its ongoing star formation. Among the Spitzer positions, we found the one closest to the Wd2 cluster to be the densest, suggesting an early phase of star formation. Furthermore, the Spitzer position that overlaps with the shell was found to have the highest $G_{\rm 0}$ and we expect this to be a result of its proximity to an O9V star. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Comments: 18 pages, 12 figures

High angular resolution near-IR view of the Orion Bar revealed by Keck/NIRC2

Authors: Emilie Habart, Romane Le Gal, Carlos Alvarez, Els Peeters, Olivier Berné, Mark G. Wolfire, Javier R. Goicoechea, Thiébaut Schirmer, Emeric Bron, Markus Röllig

Abstract: Nearby Photo-Dissociation Regions (PDRs), where the gas and dust are heated by the far UV-irradiation emitted from stars, are ideal templates to study the main stellar feedback processes. With this study we aim to probe the detailed structures at the interfaces between ionized, atomic, and molecular gas in the Orion Bar. This nearby prototypical strongly irradiated PDR will be among the first targ… ▽ More Nearby Photo-Dissociation Regions (PDRs), where the gas and dust are heated by the far UV-irradiation emitted from stars, are ideal templates to study the main stellar feedback processes. With this study we aim to probe the detailed structures at the interfaces between ionized, atomic, and molecular gas in the Orion Bar. This nearby prototypical strongly irradiated PDR will be among the first targets of the James Webb Space Telescope (JWST) within the framework of the PDRs4All Early Release Science program. We employed the sub-arcsec resolution accessible with Keck-II NIRC2 and its adaptive optics system to obtain the most detailed and complete images, ever performed, of the vibrationally excited line H$_2$ 1-0 S(1) at 2.12~$μ$m, tracing the dissociation front, and the [FeII] and Br$γ$ lines, at 1.64 and 2.16~$μ$m respectively, tracing the ionization front. We obtained narrow-band filter images in these key gas line diagnostic over $\sim 40''$ at spatial scales of $\sim$0.1$''$ ($\sim$0.0002~pc or $\sim$40~AU at 414~pc). The Keck/NIRC2 observations spatially resolve a plethora of irradiated sub-structures such as ridges, filaments, globules and proplyds. A remarkable spatial coincidence between the H$_2$ 1-0 S(1) vibrational and HCO$^+$ J=4-3 rotational emission previously obtained with ALMA is observed. This likely indicates the intimate link between these two molecular species and highlights that in high pressure PDR the H/H$_2$ and C$^+$/C/CO transitions zones come closer as compared to a typical layered structure of a constant density PDR. This is in agreement with several previous studies that claimed that the Orion Bar edge is composed of very small, dense, highly irradiated PDRs at high thermal pressure immersed in a more diffuse environment. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Journal ref: A&A 673, A149 (2023)

Photodissociation and X-Ray Dominated Regions

Authors: Mark G. Wolfire, Livia Vallini, Mélanie Chevance

Abstract: The radiation from stars and active galactic nuclei (AGN) creates photodissociation regions (PDRs) and X-ray dominated regions (XDRs), where the chemistry or heating are dominated by far-ultraviolet (FUV) radiation or X-ray radiation, respectively. PDRs include a wide range of environments from the diffuse interstellar medium to dense star-forming regions. XDRs are found in the center of galaxies… ▽ More The radiation from stars and active galactic nuclei (AGN) creates photodissociation regions (PDRs) and X-ray dominated regions (XDRs), where the chemistry or heating are dominated by far-ultraviolet (FUV) radiation or X-ray radiation, respectively. PDRs include a wide range of environments from the diffuse interstellar medium to dense star-forming regions. XDRs are found in the center of galaxies hosting AGN, in protostellar disks, and in the vicinity of X-ray binaries. In this review, we describe the dominant thermal, chemical, and radiation transfer processes in PDRs and XDRs, as well as a brief description of models and their use to analyze observations. We then present recent results from Milky Way, nearby extragalactic, and high-redshift observations. Several important results are: $\bullet$ Velocity resolved PDR lines reveal the kinematics of the neutral atomic gas and provide constraints on the stellar feedback process. Their interpretation is, however, in dispute as observations suggest a prominent role for stellar winds while they are much less important in theoretical models. $\bullet$ A significant fraction of molecular mass resides in CO-dark gas especially in low-metallicity/highly irradiated environments. $\bullet$ The CO ladder and CI/CII ratios can determine if FUV or X-rays dominate the ISM heating of extragalactic sources. $\bullet$ With ALMA, PDR and XDR tracers are now routinely detected on galactic scales over cosmic time. This makes it possible to link the star formation history of the Universe to the evolution of the physical and chemical properties of the gas. △ Less

Submitted 15 February, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Comments: Submitted to Annual Reviews of Astronomy and Astrophysics, Volume 60, 2022

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

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

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

Submitted 10 February, 2022; originally announced February 2022.

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

PDRs4All: A JWST Early Release Science Program on radiative feedback from massive stars

Authors: Olivier Berné, Émilie Habart, Els Peeters, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Emeric Bron, Jan Cami, Stéphanie Cazaux, Emmanuel Dartois, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Yoko Okada, Takashi Onaka, Massimo Robberto, Markus Röllig, Alexander G. G. M. Tielens, Silvia Vicente, Mark G. Wolfire, Felipe Alarcon, C. Boersma, Ameélie Canin, Ryan Chown, Daniel Dicken , et al. (112 additional authors not shown)

Abstract: Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation… ▽ More Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: Submitted to PASP

[CII] 158$μ$m emission from Orion A. II. Photodissociation region physics

Authors: C. H. M. Pabst, J. R. Goicoechea, A. Hacar, D. Teyssier, O. Berné, M. G. Wolfire, R. D. Higgins, E. T. Chambers, S. Kabanovic, R. Güsten, J. Stutzki, C. Kramer, A. G. G. M. Tielens

Abstract: The [CII] 158$μ$m fine-structure line is the dominant cooling line of moderate-density photodissociation regions (PDRs) illuminated by moderately bright far-ultraviolet (FUV) radiation fields. We aim to understand the origin of [CII] emission and its relation to other tracers of gas and dust in PDRs. One focus is a study of the heating efficiency of interstellar gas as traced by the [CII] line to… ▽ More The [CII] 158$μ$m fine-structure line is the dominant cooling line of moderate-density photodissociation regions (PDRs) illuminated by moderately bright far-ultraviolet (FUV) radiation fields. We aim to understand the origin of [CII] emission and its relation to other tracers of gas and dust in PDRs. One focus is a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of the photoelectric heating of neutral gas by polycyclic aromatic hydrocarbon (PAH) molecules and very small grains. We make use of a one-square-degree map of velocity-resolved [CII] line emission toward the Orion Nebula complex, and split this out into the individual spatial components, the expanding Veil Shell, the surface of OMC4, and the PDRs associated with the compact HII region of M43 and the reflection nebula NGC 1977. We employed Herschel far-infrared photometric images to determine dust properties. Moreover, we compared with Spitzer mid-infrared photometry to trace hot dust and large molecules, and velocity-resolved IRAM 30m CO(2-1) observations of the molecular gas. The [CII] intensity is tightly correlated with PAH emission in the IRAC 8$μ$m band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is very different in the four subregions and is very sensitive to the detailed geometry. Constant-density PDR models are able to reproduce the observed [CII], CO(2-1), and integrated far-infrared (FIR) intensities. We observe strong variations in the photoelectric heating efficiency in the Veil Shell behind the Orion Bar and these variations are seemingly not related to the spectral properties of the PAHs. The [CII] emission from the Orion Nebula complex stems mainly from moderately illuminated PDR surfaces. Future observations with the James Webb Space Telescope can shine light on the PAH properties that may be linked to these variations. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Journal ref: A&A 658, A98 (2022)

Characterizing Spatial Variations of PAH Emission in the Reflection Nebula NGC 1333

Authors: C. Knight, E. Peeters, M. Wolfire, D. J. Stock

Abstract: Infrared emission features at 3.3, 6.2, 7.7, 8.6, and 11.2 $μ$m, attributed to polycyclic aromatic hydrocarbons, show variations in relative intensity, shape, and peak position. These variations depend on the physical conditions of the photodissociation region (PDR) in which strong PAH emission arises, but their relationship has yet to be fully quantified. We aim to better calibrate the response o… ▽ More Infrared emission features at 3.3, 6.2, 7.7, 8.6, and 11.2 $μ$m, attributed to polycyclic aromatic hydrocarbons, show variations in relative intensity, shape, and peak position. These variations depend on the physical conditions of the photodissociation region (PDR) in which strong PAH emission arises, but their relationship has yet to be fully quantified. We aim to better calibrate the response of PAH species to their environment using observations with matching apertures and spatial resolution. We present observations from the Field-Imaging Far-Infrared Line Spectrometer (FIFI-LS) on board the Stratospheric Observatory for Infrared Astronomy (SOFIA) of the gas cooling lines [OI] 63, 146 $μ$m and [CII] 158 $μ$m in the reflection nebula NGC 1333 and use archival dust continuum observations from the Photodetector Array Camera and Spectrometer (PACS) on board Herschel. We employ PDR modelling to derive the physical conditions and compare these with the characteristics of the PAH emission as observed with the Infrared Spectrometer (IRS) on board Spitzer. We find distinct spatial characteristics for the various PAH spectral components. We conclude that the ionic bands (6.2, 7.7, 8.6, and 11.0) and the 7-9 $μ$m emission are due to multiple PAH sub-populations and that the plateaus are distinct from the features perched on top. The 6-9 $μ$m PAH emission exhibit a significant change in behaviour between the irradiated PDR and diffuse outskirts, confirming these bands arise from multiple PAH sub-populations with different underlying molecular properties. We find multiple promising relationships between PAH ratios and the FUV radiation field strength but no clear correlations with the PAH ionization parameter. △ Less

Submitted 15 November, 2021; originally announced November 2021.

Comments: 22 pages, 21 figures, to be published in MNRAS

Observation and calibration strategies for large-scale multi-beam velocity-resolved mapping of the [CII] emission in the Orion molecular cloud

Authors: R. Higgins, S. Kabanovic, C. Pabst, D. Teyssier, J. R. Goicoechea, O. Berne, E. Chambers, M. Wolfire, S. Suri, C. Buchbender, Y. Okada, M. Mertens, A. Parikka, R. Aladro, H. Richter, R. Güsten, J. Stutzki, A. G. G. M. Tielens

Abstract: Context. The [CII] 158micron far-infrared fine-structure line is one of the dominant cooling lines of the star-forming interstellar medium (ISM). Hence [CII] emission originates in and thus can be used to trace a range of ISM processes. Velocity-resolved large-scale mapping of [CII] in star-forming regions provides a unique perspective of the kinematics of these regions and their interactions with… ▽ More Context. The [CII] 158micron far-infrared fine-structure line is one of the dominant cooling lines of the star-forming interstellar medium (ISM). Hence [CII] emission originates in and thus can be used to trace a range of ISM processes. Velocity-resolved large-scale mapping of [CII] in star-forming regions provides a unique perspective of the kinematics of these regions and their interactions with the exciting source of radiation. Aims. We explore the scientific applications of large-scale mapping of velocity-resolved [CII] observations. With the [CII] observations, we investigate the effect of stellar feedback on the ISM. We present the details of observation, calibration, and data reduction using a heterodyne array receiver mounted on an airborne observatory. Results. A square-degree [CII] map with a spectral resolution of 0.3 km/s is presented. The scientific potential of this data is summarized with discussion of mechanical and radiative stellar feedback, filament tracing using [CII], [CII] opacity effects, [CII] and carbon recombination lines, and [CII] interaction with the large molecular cloud. The data quality and calibration is discussed in detail, and new techniques are presented to mitigate the effects of unavoidable instrument deficiencies (e.g. baseline stability) and thus to improve the data quality. A comparison with a smaller [CII] map taken with the Herschel/Heterodyne Instrument for the Far-Infrared (HIFI) spectrometer is presented. △ Less

Submitted 1 July, 2021; v1 submitted 29 June, 2021; originally announced June 2021.

Journal ref: A&A 652, A77 (2021)

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

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

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

Submitted 1 June, 2021; originally announced June 2021.

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

[CII] $158\,μ\mathrm{m}$ line emission from Orion A. I. A template for extragalactic studies?

Authors: C. H. M. Pabst, A. Hacar, J. R. Goicoechea, D. Teyssier, O. Berné, M. G. Wolfire, R. D. Higgins, E. T. Chambers, S. Kabanovic, R. Güsten, J. Stutzki, C. Kramer, A. G. G. M. Tielens

Abstract: The [CII] $158\,μ\mathrm{m}$ fine-structure line is one of the dominant coolants of the neutral interstellar medium. It is hence one of the brightest far-infrared emission lines and can be observed not only in star-forming regions throughout the Galaxy, but also in the diffuse interstellar medium and in distant galaxies. [CII] line emission has been suggested to be a powerful tracer of star-format… ▽ More The [CII] $158\,μ\mathrm{m}$ fine-structure line is one of the dominant coolants of the neutral interstellar medium. It is hence one of the brightest far-infrared emission lines and can be observed not only in star-forming regions throughout the Galaxy, but also in the diffuse interstellar medium and in distant galaxies. [CII] line emission has been suggested to be a powerful tracer of star-formation. We aim to understand the origin of [CII] emission and its relation to other tracers of interstellar gas and dust. This includes a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of gas heating. We make use of a one-square-degree map of velocity-resolved [CII] line emission towards the Orion Nebula complex, including M43 and NGC 1977. The [CII] intensity is tightly correlated with PAH emission in the IRAC $8\,μ\mathrm{m}$ band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is affected by the detailed geometry of the region. We find particularly low [CII]-over-FIR intensity ratios towards large columns of (warm and cold) dust, which suggest the interpretation of the "[CII] deficit" in terms of a "FIR excess". A slight decrease in the FIR line-over-continuum intensity ratio can be attributed to a decreased heating efficiency of the gas. We find that, at the mapped spatial scales, predictions of the star-formation rate from [CII] emission underestimate the star-formation rate calculated from YSO counts in the Orion Nebula complex by an order of magnitude. [CII] emission from the Orion Nebula complex arises dominantly in the cloud surfaces, many viewed in edge-on geometry. [CII] emission from extended faint cloud surfaces may contribute significantly to the total [CII] emission on galactic scales. △ Less

Submitted 8 May, 2021; originally announced May 2021.

Journal ref: A&A 651, A111 (2021)

Characterizing the Multi-Phase Origin of [CII] Emission in M101 and NGC 6946 with Velocity Resolved Spectroscopy

Authors: Elizabeth Tarantino, Alberto D. Bolatto, Rodrigo Herrera-Camus, Andrew I. Harris, Mark Wolfire, Christof Buchbender, Kevin V. Croxall, Daniel A. Dale, Brent Groves, Rebecca C. Levy, Denise Riquelme, J. D. T. Smith, Jürgen Stutzki

Abstract: The [CII] fine-structure transition at 158 micron is frequently the brightest far-infrared line in galaxies. Due to its low ionization potential, C+ can trace the ionized, atomic, and molecular phases of the ISM. We present velocity resolved [CII] and [NII] pointed observations from SOFIA/GREAT on ~500 pc scales in the nearby galaxies M101 and NGC 6946 and investigate the multi-phase origin of [CI… ▽ More The [CII] fine-structure transition at 158 micron is frequently the brightest far-infrared line in galaxies. Due to its low ionization potential, C+ can trace the ionized, atomic, and molecular phases of the ISM. We present velocity resolved [CII] and [NII] pointed observations from SOFIA/GREAT on ~500 pc scales in the nearby galaxies M101 and NGC 6946 and investigate the multi-phase origin of [CII] emission over a range of environments. We show that ionized gas makes a negligible contribution to the [CII] emission in these positions using [NII] observations. We spectrally decompose the [CII] emission into components associated with the molecular and atomic phases using existing CO(2-1) and HI data and show that a peak signal-to-noise ratio of 10-15 is necessary for a reliable decomposition. In general, we find that in our pointings greater than or equal to 50% of the [CII] emission arises from the atomic phase, with no strong dependence on star formation rate, metallicity, or galactocentric radius. We do find a difference between pointings in these two galaxies, where locations in NGC 6946 tend to have larger fractions of [CII] emission associated with the molecular phase than in M101. We also find a weak but consistent trend for fainter [CII] emission to exhibit a larger contribution from the atomic medium. We compute the thermal pressure of the cold neutral medium through the [CII] cooling function and find log(P_th/k)=3.8-4.6 [K cm^-3], a value slightly higher than similar determinations, likely because our observations are biased towards star-forming regions. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 23 pages, 11 figures, 4 tables, accepted to The Astrophysical Journal

SOFIA FEEDBACK survey: exploring the dynamics of the stellar wind driven shell of RCW 49

Authors: M. Tiwari, R. Karim, M. W. Pound, M. Wolfire, A. Jacob, C. Buchbender, R. Güsten, C. Guevara, R. D. Higgins, S. Kabanovic, C. Pabst, O. Ricken, N. Schneider, R. Simon, J. Stutzki, A. G. G. M. Tielens

Abstract: We unveil the stellar wind driven shell of the luminous massive star-forming region of RCW 49 using SOFIA FEEDBACK observations of the [CII] 158 $μ$m line. The complementary dataset of the $^{12}$CO and $^{13}$CO J = 3 - 2 transitions is observed by the APEX telescope and probes the dense gas toward RCW 49. Using the spatial and spectral resolution provided by the SOFIA and APEX telescopes, we dis… ▽ More We unveil the stellar wind driven shell of the luminous massive star-forming region of RCW 49 using SOFIA FEEDBACK observations of the [CII] 158 $μ$m line. The complementary dataset of the $^{12}$CO and $^{13}$CO J = 3 - 2 transitions is observed by the APEX telescope and probes the dense gas toward RCW 49. Using the spatial and spectral resolution provided by the SOFIA and APEX telescopes, we disentangle the shell from a complex set of individual components of gas centered around RCW 49. We find that the shell of radius ~ 6 pc is expanding at a velocity of 13 km s$^{-1}$ toward the observer. Comparing our observed data with the ancillary data at X-Ray, infrared, sub-millimeter and radio wavelengths, we investigate the morphology of the region. The shell has a well defined eastern arc, while the western side is blown open and is venting plasma further into the west. Though the stellar cluster, which is ~ 2 Myr old gave rise to the shell, it only gained momentum relatively recently as we calculate the shell's expansion lifetime ~ 0.27 Myr, making the Wolf-Rayet star WR20a a likely candidate responsible for the shell's re-acceleration. △ Less

Submitted 9 April, 2021; originally announced April 2021.

Comments: 31 pages, 17 figures

The Case for Thermalization as a Contributor to the [CII] Deficit

Authors: Jessica Sutter, Daniel A. Dale, Karin Sandstrom, J. D. T. Smith, Alberto Bolatto, Mederic Boquien, Daniela Calzetti, Kevin V. Croxall, Ilse De Looze, Maude Galametz, Brent A. Groves, George Helou, Rodrigo Herrera-Camus, Leslie K. Hunt, Robert C. Kennicutt, Eric W. Pelligrini, Christine Wilson, Mark G. Wolfire

Abstract: The [CII] deficit, which describes the observed decrease in the ratio of [CII] 158 micron emission to continuum infrared emission in galaxies with high star formation surface densities, places a significant challenge to the interpretation of [CII] detections from across the observable universe. In an attempt to further decode the cause of the [CII] deficit, the [CII] and dust continuum emission fr… ▽ More The [CII] deficit, which describes the observed decrease in the ratio of [CII] 158 micron emission to continuum infrared emission in galaxies with high star formation surface densities, places a significant challenge to the interpretation of [CII] detections from across the observable universe. In an attempt to further decode the cause of the [CII] deficit, the [CII] and dust continuum emission from 18 Local Volume galaxies has been split based on conditions within the interstellar medium where it originated. This is completed using the Key Insights in Nearby Galaxies: a Far-Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) surveys and the wide-range of wavelength information, from UV to far-infrared emission lines, available for a selection of star-forming regions within these samples. By comparing these subdivided [CII] emissions to isolated infrared emission and other properties, we find that the thermalization (collisional de-excitation) of the [CII] line in HII regions plays a significant role in the deficit observed in our sample. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: To be published in MNRAS, 8 pages, 6 figures

FEEDBACK: a SOFIA Legacy Program to Study Stellar Feedback in Regions of Massive Star Formation

Authors: N. Schneider, R. Simon, C. Guevara, C. Buchbender, R. D. Higgins, Y. Okada, J. Stutzki, R. Guesten, L. D. Anderson, J. Bally, H. Beuther, L. Bonne, S. Bontemps, E. Chambers, T. Csengeri, U. U. Graf, A. Gusdorf, K. Jacobs, S. Kabanovic, R. Karim, M. Luisi, K. Menten, M. Mertens, B. Mookerjea, V. Ossenkopf-Okada , et al. (15 additional authors not shown)

Abstract: FEEDBACK is a SOFIA legacy program dedicated to study the interaction of massive stars with their environment. It performs a survey of 11 galactic high mass star forming regions in the 158 $μ$m (1.9 THz) line of CII and the 63 $μ$m (4.7 THz) line of OI. We employ the 14 pixel LFA and 7 pixel HFA upGREAT instrument to spectrally resolve (0.24 MHz) these FIR structure lines. With an observing time o… ▽ More FEEDBACK is a SOFIA legacy program dedicated to study the interaction of massive stars with their environment. It performs a survey of 11 galactic high mass star forming regions in the 158 $μ$m (1.9 THz) line of CII and the 63 $μ$m (4.7 THz) line of OI. We employ the 14 pixel LFA and 7 pixel HFA upGREAT instrument to spectrally resolve (0.24 MHz) these FIR structure lines. With an observing time of 96h, we will cover $\sim$6700 arcmin$^2$ at 14.1$''$ angular resolution for the CII line and 6.3$''$ for the OI line. The observations started in spring 2019 (Cycle 7). Our aim is to understand the dynamics in regions dominated by different feedback processes from massive stars such as stellar winds, thermal expansion, and radiation pressure, and to quantify the mechanical energy injection and radiative heating efficiency. The CII line provides the kinematics of the gas and is one of the dominant cooling lines of gas for low to moderate densities and UV fields. The OI line traces warm and high-density gas, excited in photodissociations regions with a strong UV field or by shocks. The source sample spans a broad range in stellar characteristics from single OB stars, to small groups of O stars, to rich young stellar clusters, to ministarburst complexes. It contains well-known targets such as Aquila, the Cygnus X region, M16, M17, NGC7538, NGC6334, Vela, and W43 as well as a selection of HII region bubbles, namely RCW49, RCW79, and RCW120. These CII maps, together with the less explored OI 63 $μ$m line, provide an outstanding database for the community. They will be made publically available and will trigger further studies and follow-up observations. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Journal ref: PASP 2020, Volume 132, Number 1016; https://iopscience.iop.org/article/10.1088/1538-3873/aba840

Herschel 158$ μ$m [CII] Observations of "CO-dark" Gas in the Perseus Giant Molecular Cloud

Authors: K. P. Hall, Snežana Stanimirović, Min-Young Lee, Mark Wolfire, Paul Goldsmith

Abstract: We present observations of velocity-resolved [CII] 158 $μ$m emission from both a dense and a more diffuse photodissociation region (PDR) in the Perseus giant molecular cloud (GMC) using the Heterodyne Instrument for the Far-Infrared onboard the Herschel Space Telescope. We detect [CII] emission from 80\% of the total positions, with a 95\% detection rate from the dense boundary region. The integra… ▽ More We present observations of velocity-resolved [CII] 158 $μ$m emission from both a dense and a more diffuse photodissociation region (PDR) in the Perseus giant molecular cloud (GMC) using the Heterodyne Instrument for the Far-Infrared onboard the Herschel Space Telescope. We detect [CII] emission from 80\% of the total positions, with a 95\% detection rate from the dense boundary region. The integrated intensity of the [CII] emission remains relatively constant across each boundary, despite the observed range in optical extinction between 1 mag and 10 mag. This flat profile indicates a constant heating and cooling rate within both regions observed. The integrated intensity of [CII] emission is reasonably well correlated with the neutral hydrogen (HI) column density, as well as total gas column density. This, in addition to the 80$'$ (7 pc) extent of the [CII] emission from cloud center, suggests that the HI envelope plays a dominant role in explaining the [CII] emission emanating from Perseus. We compare the [CII] and $^{12}$CO integrated intensities with predictions from a 1-D, two-sided slab PDR model and show that a simple core $+$ envelope, equilibrium model without an additional "CO-dark" H$_2$ component can reproduce observations well. Additional observations are needed to disentangle how much of the [CII] emission is associated with the "CO-dark" H$_2$ gas, as well as constrain spatial variations of the dust-to-gas ratio across Perseus. △ Less

Submitted 17 June, 2020; originally announced June 2020.

Molecular globules in the Veil bubble of Orion. IRAM 30m 12CO, 13CO, and C18O 2-1 expanded maps of Orion A

Authors: J. R. Goicoechea, C. H. M. Pabst, S. Kabanovic, M. G. Santa-Maria, N. Marcelino, A. G. G. M. Tielens, A. Hacar, O. Berne, C. Buchbender, S. Cuadrado, R. Higgins, C. Kramer, J. Stutzki, S. Suri, D. Teyssier, M. Wolfire

Abstract: Strong winds and ultraviolet (UV) radiation from O-type stars disrupt and ionize their molecular core birthplaces, sweeping up material into parsec-size shells. Owing to dissociation by starlight, the thinnest shells are expected to host low molecular abundances and therefore little star formation. Here, we expand previous maps taken with the IRAM 30m telescope and present square-degree 12CO and 1… ▽ More Strong winds and ultraviolet (UV) radiation from O-type stars disrupt and ionize their molecular core birthplaces, sweeping up material into parsec-size shells. Owing to dissociation by starlight, the thinnest shells are expected to host low molecular abundances and therefore little star formation. Here, we expand previous maps taken with the IRAM 30m telescope and present square-degree 12CO and 13CO (J=2-1) maps of the wind-driven "Veil bubble'' that surrounds the Trapezium cluster and its natal Orion molecular core (OMC). Although widespread and extended CO emission is largely absent from the Veil, we show that several CO "globules'' exist and are embedded in the [CII]158um-bright shell that confines the bubble. This includes the first detection of quiescent CO at negative LSR velocities in Orion. Given the harsh UV irradiation conditions in this translucent material, the detection of CO globules is surprising. These globules are small (R=7,100 AU), not massive (M=0.3M_Sun), and are moderately dense: n_ H=4x10^4 cm^-3 (median values). They are confined by the external pressure of the shell, P_ext/k~10^7 cm^-3 K, and are likely magnetically supported. They are either transient objects formed by instabilities or have detached from pre-existing molecular structures, sculpted by the passing shock associated with the expanding shell and by UV radiation from the Trapezium. Some represent the first stages in the formation of small pillars, others of isolated small globules. Although their masses do not suggest they will form stars, one globule matches the position of a known YSO. The lack of extended CO in the "Veil shell'' demonstrates that feedback from massive stars expels, agitates, and reprocesses most of the disrupted molecular cloud gas, thereby limiting the star-formation rate in the region. The presence of globules is a result of this feedback. △ Less

Submitted 5 May, 2020; v1 submitted 27 April, 2020; originally announced April 2020.

Comments: Accepted for publication in A&A. (Abridged abstract. English edited in this v2. Some figures have been bitmapped to lower resolution)

Journal ref: A&A 639, A1 (2020)

[CI](1-0) and [CI](2-1) in resolved local galaxies

Authors: Alison F. Crocker, Eric Pellegrini, J. -D. T. Smith, Bruce T. Draine, Christine D. Wilson, Mark Wolfire, Lee Armus, Elias Brinks, Daniel A. Dale, Brent Groves, Rodrigo Herrera-Camus, Leslie K. Hunt, Robert C. Kennicutt, Eric J. Murphy, Karin Sandstrom, Eva Schinnerer, Dimitra Rigopoulou, Erik Rosolowsky, Paul van der Werf

Abstract: We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4… ▽ More We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4-3) and [CI](2-1) lines; we hypothesize this is due to the similar upper level temperature of these two lines. We modeled the [CI] and CO line emission using large velocity gradient models combined with an empirical template. According to this modeling, the [CI](1-0) line is clearly dominated by the low-excitation component. We determine [CI] to molecular mass conversion factors for both the [CI](1-0) and [CI](2-1) lines, with mean values of $α_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $α_{\mathrm{[CI](2-1)}} = 34 $ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ with logarithmic root-mean-square spreads of 0.20 and 0.32 dex, respectively. The similar spread of $α_{\mathrm{[CI](1-0)}}$ to $α_{\mathrm{CO}}$ (derived using the CO(2-1) line) suggests that [CI](1-0) may be just as good a tracer of cold molecular gas as CO(2-1) in galaxies of this type. On the other hand, the wider spread of $α_{\mathrm{[CI](2-1)}}$ and the tight relation found between [CI](2-1) and CO(4-3) suggest that much of the [CI](2-1) emission may originate in warmer molecular gas. △ Less

Submitted 14 January, 2020; originally announced January 2020.

Comments: 15 pages, 10 figures, published in ApJ

Journal ref: The Astrophysical Journal 887 (2019) 105

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: The KINGFISH Sample

Authors: G. Aniano, B. T. Draine, L. K. Hunt, K. Sandstrom, D. Calzetti, R. C. Kennicutt, D. A. Dale, M. Galametz, K. D. Gordon, A. K. Leroy, J. -D. T. Smith, H. Roussel, M. Sauvage, F. Walter, L. Armus, A. D. Bolatto, M. Boquien, A. Crocker, I. De Looze, J. Donovan Meyer, G. Helou, J. Hinz, B. D. Johnson, J. Koda, A. Miller , et al. (8 additional authors not shown)

Abstract: Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starli… ▽ More Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: ApJ, accepted

Using [CII] 158 micron Emission From Isolated ISM Phases as a Star-Formation Rate Indicator

Authors: Jessica Sutter, Daniel A. Dale, Kevin V. Croxall, Eric W. Pelligrini, J. D. T. Smith, Phillip N. Appleton, Pedro Beirao, Alberto D. Bolatto, Daniela Calzetti, Alison Crocker, Ilse De Looze, Bruce Draine, Maud Galametz, Brent A. Groves, George Helou, Rodrigo Herrera-Camus, Leslie K. Hunt, Robert C. Kennicutt, Helene Roussel, Mark G. Wolfire

Abstract: The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that origi… ▽ More The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that originates from the ionized and neutral phases of the ISM. Sub-kiloparsec measurements of the [CII] 158 micron and [NII] 205 micron line in nearby galaxies have recently become available as part of the Key Insights in Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) programs. With the information from these two far-infrared lines along with the multi-wavelength suite of KINGFISH data, a calibration of the [CII] emission line as a star formation rate indicator and a better understanding of the [CII] deficit are pursued. [CII] emission is also compared to PAH emission in these regions to compare photoelectric heating from PAH molecules to cooling by [CII] in the neutral and ionized phases of the ISM. We find that the [CII] emission originating in the neutral phase of the ISM does not exhibit a deficit with respect to the infrared luminosity and is therefore preferred over the [CII] emission originating in the ionized phase of the ISM as a star formation rate indicator for the normal star-forming galaxies included in this sample. △ Less

Submitted 11 October, 2019; originally announced October 2019.

Comments: 19 pages, 8 figures, 8 tables

Chemical Abundances in a Turbulent Medium -- H$_2$, OH$^+$, H$_2$O$^+$, ArH$^+$

Authors: Shmuel Bialy, David Neufeld, Mark Wolfire, Amiel Sternberg, Blakesley Burkhart

Abstract: Supersonic turbulence results in strong density fluctuations in the interstellar medium (ISM), which have a profound effect on the chemical structure. Particularly useful probes of the diffuse ISM are the ArH$^+$, OH$^+$, H$_2$O$^+$ molecular ions, which are highly sensitive to fluctuations in the density and the H$_2$ abundance. We use isothermal magnetohydrodynamic (MHD) simulations of various s… ▽ More Supersonic turbulence results in strong density fluctuations in the interstellar medium (ISM), which have a profound effect on the chemical structure. Particularly useful probes of the diffuse ISM are the ArH$^+$, OH$^+$, H$_2$O$^+$ molecular ions, which are highly sensitive to fluctuations in the density and the H$_2$ abundance. We use isothermal magnetohydrodynamic (MHD) simulations of various sonic Mach numbers, $\mathcal{M}_s$, and density decorrelation scales, $y_{\rm dec}$, to model the turbulent density field. We post-process the simulations with chemical models and obtain the probability density functions (PDFs) for the H$_2$, ArH$^+$, OH$^+$ and H$_2$O$^+$ abundances. We find that the PDF dispersions increases with increasing $\mathcal{M}_s$ and $y_{\rm dec}$, as the magnitude of the density fluctuations increases, and as they become more coherent. Turbulence also affects the median abundances: when $\mathcal{M}_s$ and $y_{\rm dec}$ are high, low-density regions with low H$_2$ abundance become prevalent, resulting in an enhancement of ArH$^+$ compared to OH$^+$ and H$_2$O$^+$. We compare our models with Herschel observations. The large scatter in the observed abundances, as well as the high observed ArH$^+$/OH$^+$ and ArH$^+$/H$_2$O$^+$ ratios are naturally reproduced by our supersonic $(\mathcal{M}_s=4.5)$, large decorrelation scale $(y_{\rm dec}=0.8)$ model, supporting a scenario of a large-scale turbulence driving. The abundances also depend on the UV intensity, CR ionization rate, and the cloud column density, and the observed scatter may be influenced by fluctuations in these parameters. △ Less

Submitted 9 October, 2019; v1 submitted 26 September, 2019; originally announced September 2019.

Comments: Accepted for publication in the Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 885, Number 2, Pages 109, 2019

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

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

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

Submitted 13 August, 2019; originally announced August 2019.

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

Probing ISM Structure in Trumpler 14 & Carina I Using The Stratospheric Terahertz Observatory 2

Authors: Young Min Seo, Paul F. Goldsmith, Chris Walker, David J. Hollenbach, Mark G. Wolfire, Craig Kulesa, Volker Tolls, Pietro N. Bernasconi, Umit Kavak, Floris F. S. van der Tak, Russ Shipman, Jian Rong Gao, Alexander Tielens, Michael G. Burton, Harold Yorke, Erick Young, William L. Peters, Abram Young, Christopher Groppi, Kristina Davis, Jorge L. Pineda, William D. Langer, Jonathan H. Kawamura, Antony Stark, Gary Melnick , et al. (4 additional authors not shown)

Abstract: We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $μ$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a… ▽ More We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $μ$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a velocity resolution of 0.17 km s$^{-1}$. The observations cover a 0.25$^\circ$ by 0.28$^\circ$ area with central position {\it l} = 297.34$^\circ$, {\it b} = -0.60$^\circ$. The kinematics show that bright [C\,{\sc ii}] structures are spatially and spectrally correlated with the surfaces of CO clouds, tracing the photodissociation region and ionization front of each molecular cloud. Along 7 lines of sight that traverse Tr 14 into the dark ridge to the southwest, we find that the [C\,{\sc ii}] luminosity from the HII region is 3.7 times that from the PDR. In same los we find in the PDRs an average ratio of 1:4.1:5.6 for the mass in atomic gas:dark-CO gas: molecular gas traced by CO. Comparing multiple gas tracers including HI 21cm, [C\,{\sc ii}], CO, and radio recombination lines, we find that the HII regions of the Carina Nebula Complex are well-described as HII regions with one-side freely expanding towards us, consistent with the champagne model of ionized gas evolution. The dispersal of the GMC in this region is dominated by EUV photoevaporation; the dispersal timescale is 20-30 Myr. △ Less

Submitted 22 May, 2019; v1 submitted 22 March, 2019; originally announced March 2019.

Comments: ApJ accepted

Disruption of the Orion Molecular Core 1 by the stellar wind of the massive star $θ^1$ Ori C

Authors: C. Pabst, R. Higgins, J. R. Goicoechea, D. Teyssier, O. Berne, E. Chambers, M. Wolfire, S. T. Suri, R. Guesten, J. Stutzki, U. U. Graf, C. Risacher, A. G. G. M. Tielens

Abstract: Massive stars inject mechanical and radiative energy into the surrounding environment, which stirs it up, heats the gas, produces cloud and intercloud phases in the interstellar medium, and disrupts molecular clouds (the birth sites of new stars). Stellar winds, supernova explosions and ionization by ultraviolet photons control the lifetimes of molecular clouds. Theoretical studies predict that mo… ▽ More Massive stars inject mechanical and radiative energy into the surrounding environment, which stirs it up, heats the gas, produces cloud and intercloud phases in the interstellar medium, and disrupts molecular clouds (the birth sites of new stars). Stellar winds, supernova explosions and ionization by ultraviolet photons control the lifetimes of molecular clouds. Theoretical studies predict that momentum injection by radiation should dominate that by stellar winds, but this has been difficult to assess observationally. Velocity-resolved large-scale images in the fine-structure line of ionized carbon ([C II]) provide an observational diagnostic for the radiative energy input and the dynamics of the interstellar medium around massive stars. Here we report observations of a one-square-degree region (about 7 parsecs in diameter) of Orion molecular core -- the region nearest to Earth that exhibits massive-star formation -- at a resolution of 16 arcseconds (0.03 parsecs) in the [C II] line at 1.9 terahertz (158 micrometres). The results reveal that the stellar wind originating from the massive star $θ^{1}$ Orionis C has swept up the surrounding material to create a bubble roughly four parsecs in diameter with a 2,600-solar-mass shell, which is expanding at 13 kilometres per second. This finding demonstrates that the mechanical energy from the stellar wind is converted very efficiently into kinetic energy of the shell and causes more disruption of the Orion molecular core 1 than do photo-ionization and evaporation or future supernova explosions. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: Published in Nature

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

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

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

Submitted 14 August, 2018; originally announced August 2018.

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

Journal ref: ApJ 865 (2018) 127

Mapping low frequency carbon radio recombination lines towards Cassiopeia A at 340, 148, 54 and 43 MHz

Authors: Pedro Salas, J. B. Raymond Oonk, Reinout J. van Weeren, Mark G. Wolfire, Kimberly L. Emig, M. Carmen Toribio, Huub J. A. Röttgering, Alexander G. G. M. Tielens

Abstract: Quantitative understanding of the interstellar medium requires knowledge of its physical conditions. Low frequency carbon radio recombination lines (CRRLs) trace cold interstellar gas, and can be used to determine its physical conditions (e.g., electron temperature and density). In this work we present spatially resolved observations of the low frequency ($\leq390$ MHz) CRRLs centered around C… ▽ More Quantitative understanding of the interstellar medium requires knowledge of its physical conditions. Low frequency carbon radio recombination lines (CRRLs) trace cold interstellar gas, and can be used to determine its physical conditions (e.g., electron temperature and density). In this work we present spatially resolved observations of the low frequency ($\leq390$ MHz) CRRLs centered around C$268α$, C$357α$, C$494α$ and C$539α$ towards Cassiopeia A on scales of $\leq1.2$ pc. We compare the spatial distribution of CRRLs with other ISM tracers. This comparison reveals a spatial offset between the peak of the CRRLs and other tracers, which is very characteristic for photodissociation regions and that we take as evidence for CRRLs being preferentially detected from the surfaces of molecular clouds. Using the CRRLs we constrain the gas electron temperature and density. These constraints on the gas conditions suggest variations of less than a factor of two in pressure over $\sim1$ pc scales, and an average hydrogen density of $200$-$470$ cm$^{-3}$. From the electron temperature and density maps we also constrain the ionized carbon emission measure, column density and path length. Based on these, the hydrogen column density is larger than $10^{22}$ cm$^{-2}$, with a peak of $\sim4\times10^{22}$ cm$^{-2}$ towards the South of Cassiopeia A. Towards the southern peak the line of sight length is $\sim40$ pc over a $\sim2$ pc wide structure, which implies that the gas is a thin surface layer on a large (molecular) cloud that is only partially intersected by Cassiopeia A. These observations highlight the utility of CRRLs as tracers of low density extended HI and CO-dark gas halo's around molecular clouds. △ Less

Submitted 16 January, 2018; originally announced January 2018.

Comments: 18 pages, 15 figures, 4 tables, accepted by MNRAS

First Results from the $Herschel$ and ALMA Spectroscopic Surveys of the SMC: The Relationship Between [CII]-bright Gas and CO-bright Gas at Low Metallicity

Authors: Katherine E. Jameson, Alberto D. Bolatto, Mark Wolfire, Steven R. Warren, Rodrigo Herrera-Camus, Kevin Croxall, Eric Pellegrini, John-David Smith, Monica Rubio, Remy Indebetouw, Frank P. Israel, Margaret Meixner, Julia Roman-Duval, Jacco Th. van Loon, Erik Muller, Celia Verdugo, Hans Zinnecker, Yoko Okada

Abstract: The Small Magellanic Cloud (SMC) provides the only laboratory to study the structure of molecular gas at high resolution and low metallicity. We present results from the Herschel Spectroscopic Survey of the SMC (HS$^{3}$), which mapped the key far-IR cooling lines [CII], [OI], [NII], and [OIII] in five star-forming regions, and new ALMA 7m-array maps of $^{12}$CO and $^{13}$CO $(2-1)$ with coverag… ▽ More The Small Magellanic Cloud (SMC) provides the only laboratory to study the structure of molecular gas at high resolution and low metallicity. We present results from the Herschel Spectroscopic Survey of the SMC (HS$^{3}$), which mapped the key far-IR cooling lines [CII], [OI], [NII], and [OIII] in five star-forming regions, and new ALMA 7m-array maps of $^{12}$CO and $^{13}$CO $(2-1)$ with coverage overlapping four of the five HS$^{3}$ regions. We detect [CII] and [OI] throughout all of the regions mapped. The data allow us to compare the structure of the molecular clouds and surrounding photodissociation regions using $^{13}$CO, CO, [CII], and [OI] emission at $<10$" ($<3$ pc) scales. We estimate Av using far-IR thermal continuum emission from dust and find the CO/[CII] ratios reach the Milky Way value at high A$_{V}$ in the centers of the clouds and fall to $\sim{1/5-1/10}\times$ the Milky Way value in the outskirts, indicating the presence of translucent molecular gas not traced by bright CO emission. We estimate the amount of molecular gas traced by bright [CII] emission at low A$_{V}$ and bright CO emission at high A$_{V}$. We find that most of the molecular gas is at low A$_{V}$ and traced by bright [CII] emission, but that faint CO emission appears to extend to where we estimate the H$_{2}$-to-HI transition occurs. By converting our H$_{2}$ gas estimates to a CO-to-H$_{2}$ conversion factor ($X_{CO}$), we show that $X_{CO}$ is primarily a function of A$_{V}$, consistent with simulations and models of low metallicity molecular clouds. △ Less

Submitted 10 January, 2018; originally announced January 2018.

Comments: Accepted for publication in ApJ

[CII] emission from L1630 in the Orion B molecular cloud

Authors: C. H. M. Pabst, J. R. Goicoechea, D. Teyssier, O. Berné, B. B. Ochsendorf, M. G. Wolfire, R. D. Higgins, D. Riquelme, C. Risacher, J. Pety, F. Le Petit, E. Roueff, E. Bron, A. G. G. M. Tielens

Abstract: Observations towards L1630 in the Orion B molecular cloud, comprising the iconic Horsehead Nebula, allow us to study the interplay between stellar radiation and a molecular cloud under relatively benign conditions, that is, intermediate densities and an intermediate UV radiation field. Contrary to the well-studied Orion Molecular Cloud 1 (OMC1), which hosts much harsher conditions, L1630 has littl… ▽ More Observations towards L1630 in the Orion B molecular cloud, comprising the iconic Horsehead Nebula, allow us to study the interplay between stellar radiation and a molecular cloud under relatively benign conditions, that is, intermediate densities and an intermediate UV radiation field. Contrary to the well-studied Orion Molecular Cloud 1 (OMC1), which hosts much harsher conditions, L1630 has little star formation. We aim to relate the [CII] fine-structure line emission to the physical conditions predominant in L1630 and compare it to studies of OMC1. The [CII] $158\,μ\mathrm{m}$ emission from an area of $12' \times 17'$ in L1630 was observed using the upGREAT instrument onboard SOFIA. Of the [CII] emission from the mapped area 95%, $13\,L_{\odot}$, originates from the molecular cloud; the adjacent HII region contributes only 5%, that is, $1\,L_{\odot}$. From comparison with other data (CO (1-0)-line emission, far-infrared (FIR) continuum studies, emission from polycyclic aromatic hydrocarbons (PAHs)), we infer a gas density of the molecular cloud of $n_{\mathrm{H}}\sim 3\cdot 10^3\,\mathrm{cm^{-3}}$, with surface layers, including the Horsehead Nebula, having a density of up to $n_{\mathrm{H}}\sim 4\cdot 10^4\,\mathrm{cm^{-3}}$. The temperature of the surface gas is $T\sim 100\,\mathrm{K}$. The average [CII] cooling efficiency within the molecular cloud is $1.3\cdot 10^{-2}$. The fraction of the mass of the molecular cloud within the studied area that is traced by [CII] is only $8\%$. Our PDR models are able to reproduce the FIR-[CII] correlations and also the CO (1-0)-[CII] correlations. Finally, we compare our results on the heating efficiency of the gas with theoretical studies of photoelectric heating by PAHs, clusters of PAHs, and very small grains, and find the heating efficiency to be lower than theoretically predicted, a continuation of the trend set by other observations. △ Less

Submitted 19 July, 2017; originally announced July 2017.

Comments: 23 pages, 30 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 606, A29 (2017)

The Origins of [CII] Emission in Local Star-forming Galaxies

Authors: Kevin Croxall, J. D. T. Smith, Eric Pellegrini, Brent Groves, Alberto Bolatto, Rodrigo Herrera-Camus, Karin Sandstrom, Bruce T. Draine, Mark Wolfire, Lee Armus, Mederic Boquien, Bernhard Brandl, Daniel A. Dale, Maud Galametz, Leslie K. Hunt, Robert C. Kennicutt, Kathryn Kreckel, Dimitra Rigopoulou, Paul P. van der Werf, Christine D. Wilson

Abstract: The [CII] 158um fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, [CII] balances the heating, including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of [CII] emission remains unclear, because C+ can be found in multiple phases of the inters… ▽ More The [CII] 158um fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, [CII] balances the heating, including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of [CII] emission remains unclear, because C+ can be found in multiple phases of the interstellar medium. Here we measure the fractions of [CII] emission originating in the ionized and neutral gas phases of a sample of nearby galaxies. We use the [NII] 205um fine-structure line to trace the ionized medium, thereby eliminating the strong density dependence that exists in the ratio of [CII]/[NII] 122um. Using the FIR [CII] and [NII] emission detected by the KINGFISH and Beyond the Peak Herschel programs, we show that 60-80% of [CII] emission originates from neutral gas. We find that the fraction of [CII] originating in the neutral medium has a weak dependence on dust temperature and the surface density of star formation, and a stronger dependence on the gas-phase metallicity. In metal-rich environments, the relatively cooler ionized gas makes substantially larger contributions to total [CII] emission than at low abundance, contrary to prior expectations. Approximate calibrations of this metallicity trend are provided. △ Less

Submitted 14 July, 2017; originally announced July 2017.

Comments: 8 pages, accepted for publication in ApJ

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

Authors: David A. Neufeld, Mark G. Wolfire

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

Submitted 12 April, 2017; originally announced April 2017.

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

The PAH emission characteristics of the reflection nebula NGC2023

Authors: Els Peeters, Charles W. Bauschlicher, Jr., Louis J. Allamandola, Alexander G. G. M. Tielens, Alessandra Ricca, Mark G. Wolfire

Abstract: We present 5-20 micron spectral maps of the reflection nebula NGC2023 obtained with the Infrared Spectrograph SL and SH modes on board the Spitzer Space Telescope which reveal emission from polycyclic aromatic hydrocarbons (PAHs), C60, and H2 superposed on a dust continuum. We show that several PAH emission bands correlate with each other and exhibit distinct spatial distributions revealing a spat… ▽ More We present 5-20 micron spectral maps of the reflection nebula NGC2023 obtained with the Infrared Spectrograph SL and SH modes on board the Spitzer Space Telescope which reveal emission from polycyclic aromatic hydrocarbons (PAHs), C60, and H2 superposed on a dust continuum. We show that several PAH emission bands correlate with each other and exhibit distinct spatial distributions revealing a spatial sequence with distance from the illuminating star. We explore the distinct morphology of the 6.2, 7.7 and 8.6 micron PAH bands and find that at least two spatially distinct components contribute to the 7--9 micron PAH emission in NGC2023. We report that the PAH features behave independently of the underlying plateaus. We present spectra of compact oval PAHs ranging in size from C_66 to C_210, determined computationally using density functional theory, and investigate trends in the band positions and relative intensities as a function of PAH size, charge and geometry. Based on the NASA Ames PAH database, we discuss the 7--9 micron components in terms of band assignments and relative intensities. We assign the plateau emission to very small grains with possible contributions from PAH clusters and identify components in the 7--9 micron emission that likely originates in these structures. Based on the assignments and the observed spatial sequence, we discuss the photochemical evolution of the interstellar PAH family as they are more and more exposed to the radiation field of the central star in the evaporative flows associated with the PDRs in NGC2023. △ Less

Submitted 23 January, 2017; originally announced January 2017.

Comments: 42 pages, 25 figures, 2 animations, accepted for ApJ

Thermal Pressure in the Cold Neutral Medium of Nearby Galaxies

Authors: R. Herrera-Camus, A. Bolatto, M. Wolfire, E. Ostriker, B. Draine, A. Leroy, K. Sandstrom, L. Hunt, R. Kennicutt, D. Calzetti, J. D. Smith, K. Croxall, M. Galametz, I. de Looze, D. Dale, A. Crocker, B. Groves

Abstract: Dynamic and thermal processes regulate the structure of the multi-phase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $μ$m, HI, and CO data… ▽ More Dynamic and thermal processes regulate the structure of the multi-phase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $μ$m, HI, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of $\sim$1 kiloparsec. We find a distribution of thermal pressures in the $P_{\rm th}/k\sim10^3-10^5$ K cm$^{-3}$ range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of $P_{\rm th}/k\approx3600$ K cm$^{-3}$. These results are consistent with thermal pressure measurements using other observational methods. We find that $P_{\rm th}$ increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral medium could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about $\sim30$% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks. △ Less

Submitted 5 January, 2017; originally announced January 2017.

Comments: Accepted for publication in ApJ on Dec. 20th, 2016; 17 pages; 7 figures and 1 table

The Spatially Resolved [CII] Cooling Line Deficit in Galaxies

Authors: J. D. T. Smith, Kevin Croxall, Bruce Draine, Ilse De Looze, Karin Sandstrom, Lee Armus, Pedro Beirao, Alberto Bolatto, Mederic Boquien, Bernhard Brandl, Alison Crocker, Daniel A. Dale, Maud Galametz, Brent Groves, George Helou, Rodrigo Herrera-Camus, Leslie Hunt, Robert Kennicutt, Fabian Walter, Mark Wolfire

Abstract: We present [CII] 158um measurements from over 15,000 resolved regions within 54 nearby galaxies of the KINGFISH program to investigate the so-called [CII] "line cooling deficit" long known to occur in galaxies with different luminosities. The [CII]/TIR ratio ranges from above 1% to below 0.1% in the sample, with a mean value of 0.48+-0.21%. We find that the surface density of 24um emission dominat… ▽ More We present [CII] 158um measurements from over 15,000 resolved regions within 54 nearby galaxies of the KINGFISH program to investigate the so-called [CII] "line cooling deficit" long known to occur in galaxies with different luminosities. The [CII]/TIR ratio ranges from above 1% to below 0.1% in the sample, with a mean value of 0.48+-0.21%. We find that the surface density of 24um emission dominates this trend, with [CII]/TIR dropping as nuInu{24um} increases. Deviations from this overall decline are correlated with changes in the gas phase metal abundance, with higher metallicity associated with deeper deficits at a fixed surface brightness. We supplement the local sample with resolved [CII] measurements from nearby luminous infrared galaxies and high redshift sources from z=1.8-6.4, and find that star formation rate density drives a continuous trend of deepening [CII] deficit across six orders of magnitude in SFRD. The tightness of this correlation suggests that an approximate star formation rate density can be estimated directly from global measurements of [CII]/TIR, and a relation is provided to do so. Several low-luminosity AGN hosts in the sample show additional and significant central suppression of [CII]/TIR, but these deficit enhancements occur not in those AGN with the highest X-ray luminosities, but instead those with the highest central starlight intensities. Taken together, these results demonstrate that the [CII] cooling line deficit in galaxies likely arises from local physical phenomena in interstellar gas. △ Less

Submitted 4 November, 2016; originally announced November 2016.

Comments: 17 pages, 5 figures, accepted for publication in ApJ

A Simple and Accurate Network for Hydrogen and Carbon Chemistry in the ISM

Authors: Munan Gong, Eve C. Ostriker, Mark G. Wolfire

Abstract: Chemistry plays an important role in the interstellar medium (ISM), regulating heating and cooling of the gas, and determining abundances of molecular species that trace gas properties in observations. Although solving the time-dependent equations is necessary for accurate abundances and temperature in the dynamic ISM, a full chemical network is too computationally expensive to incorporate in nume… ▽ More Chemistry plays an important role in the interstellar medium (ISM), regulating heating and cooling of the gas, and determining abundances of molecular species that trace gas properties in observations. Although solving the time-dependent equations is necessary for accurate abundances and temperature in the dynamic ISM, a full chemical network is too computationally expensive to incorporate in numerical simulations. In this paper, we propose a new simplified chemical network for hydrogen and carbon chemistry in the atomic and molecular ISM. We compare results from our chemical network in detail with results from a full photo-dissociation region (PDR) code, and also with the Nelson & Langer (1999) (NL99) network previously adopted in the simulation literature. We show that our chemical network gives similar results to the PDR code in the equilibrium abundances of all species over a wide range of densities, temperature, and metallicities, whereas the NL99 network shows significant disagreement. Applying our network in 1D models, we find that the $\mathrm{CO}$-dominated regime delimits the coldest gas and that the corresponding temperature tracks the cosmic ray ionization rate in molecular clouds. We provide a simple fit for the locus of $\mathrm{CO}$ dominated regions as a function of gas density and column. We also compare with observations of diffuse and translucent clouds. We find that the $\mathrm{CO}$, $\mathrm{CHx}$ and $\mathrm{OHx}$ abundances are consistent with equilibrium predictions for densities $n=100-1000~\mathrm{cm^{-3}}$, but the predicted equilibrium $\mathrm{C}$ abundance is higher than observations, signaling the potential importance of non-equilibrium/dynamical effects. △ Less

Submitted 4 April, 2019; v1 submitted 27 October, 2016; originally announced October 2016.

Comments: A PDR code written in C++ and python based on this chemical network can be found at https://github.com/munan/pdr

Journal ref: ApJ, 843:38 (35pp), 2017 July 1

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

Authors: David A. Neufeld, Mark G. Wolfire

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

Submitted 1 July, 2016; originally announced July 2016.

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

The Ionized Gas in Nearby Galaxies as Traced by the [NII] 122 and 205 μm Transitions

Authors: R. Herrera-Camus, A. Bolatto, J. D. Smith, B. Draine, E. Pellegrini, M. Wolfire, K. Croxall, I. de Looze, D. Calzetti, R. Kennicutt, A. Crocker, L. Armus, P. van der Werf, K. Sandstrom, M. Galametz, B. Brandl, B. Groves, D. Rigopoulou, F. Walter, A. Leroy, M. Boquien, F. S. Tabatabaei, P. Beirao

Abstract: The [NII] 122 and 205 μm transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range,… ▽ More The [NII] 122 and 205 μm transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in the electron density within individual galaxies can be as a high as a factor of ~50, frequently with strong radial gradients. We find that $n_e$ increases as a function of infrared color, dust-weighted mean starlight intensity, and star formation rate surface density ($Σ_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between [NII] surface brightness and $Σ_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code. △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: 18 pages, 9 figures, accepted for publication in The Astrophysical Journal

Extended Carbon Line Emission in the Galaxy: Searching for Dark Molecular Gas along the G328 Sightline

Authors: Michael G. Burton, Michael C. B. Ashley, Catherine Braiding, Matthew Freeman, Craig Kulesa, Mark G. Wolfire, David J. Hollenbach, Gavin Rowell, James Lau

Abstract: We present spectral data cubes of the [CI] 809GHz, 12CO 115GHz, 13CO 110GHz and HI 1.4GHz line emission from an 1 square degree region along the l = 328° (G328) sightline in the Galactic Plane. Emission arises principally from gas in three spiral arm crossings along the sight line. The distribution of the emission in the CO and [CI] lines is found to be similar, with the [CI] slightly more extende… ▽ More We present spectral data cubes of the [CI] 809GHz, 12CO 115GHz, 13CO 110GHz and HI 1.4GHz line emission from an 1 square degree region along the l = 328° (G328) sightline in the Galactic Plane. Emission arises principally from gas in three spiral arm crossings along the sight line. The distribution of the emission in the CO and [CI] lines is found to be similar, with the [CI] slightly more extended, and both are enveloped in extensive HI. Spectral line ratios per voxel in the data cubes are found to be similar across the entire extent of the Galaxy. However, towards the edges of the molecular clouds the [CI]/13CO and 12CO/13CO line ratios rise by ~50%, and the [CI]/HI ratio falls by ~10$%. We attribute this to these sightlines passing predominantly through the surfaces of photodissociation regions (PDRs), where the carbon is found mainly as C or C+, while the H2 is mostly molecular, and the proportion of atomic gas also increases. We undertake modelling of the PDR emission from low density molecular clouds excited by average interstellar radiation fields and cosmic-ray ionization to quantify this comparison, finding that depletion of sulfur and reduced PAH abundance is needed to match line fluxes and ratios. Roughly one-third of the molecular gas along the sightline is found to be associated with this surface region, where the carbon is largely not to be found in CO. ~10% of the atomic hydrogen along the sightline is cold gas within PDRs. △ Less

Submitted 19 August, 2015; originally announced August 2015.

Comments: 66 Pages, 21 Figures (2 in Appendix), 4 Tables, Accepted for Publication in The Astrophysical Journal, 13 August 2015, Ref ApJ/519068/ART

Journal ref: 2015 Astrophysical Journal, 811, #13