Showing 1–50 of 77 results for author: Petit, F L

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

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)

Bias versus variance when fitting multi-species molecular lines with a non-LTE radiative transfer model

Authors: Antoine Roueff, Jérôme Pety, Maryvonne Gerin, Léontine Ségal, Javier Goicoechea, Harvey Liszt, Pierre Gratier, Ivana Bešlić, Lucas Einig, M. Gaudel, Jan Orkisz, Pierre Palud, Miriam Santa-Maria, Victor de Souza Magalhaes, Antoine Zakardjian, Sebastien Bardeau, Emeric E. Bron, Pierre Chainais, Simon Coudé, Karine Demyk, Viviana Guzman Veloso, Annie Hughes, David Languignon, François Levrier, Dariusz C Lis , et al. (6 additional authors not shown)

Abstract: Robust radiative transfer techniques are requisite for efficiently extracting the physical and chemical information from molecular rotational lines.We study several hypotheses that enable robust estimations of the column densities and physical conditions when fitting one or two transitions per molecular species. We study the extent to which simplifying assumptions aimed at reducing the complexity… ▽ More Robust radiative transfer techniques are requisite for efficiently extracting the physical and chemical information from molecular rotational lines.We study several hypotheses that enable robust estimations of the column densities and physical conditions when fitting one or two transitions per molecular species. We study the extent to which simplifying assumptions aimed at reducing the complexity of the problem introduce estimation biases and how to detect them.We focus on the CO and HCO+ isotopologues and analyze maps of a 50 square arcminutes field. We used the RADEX escape probability model to solve the statistical equilibrium equations and compute the emerging line profiles, assuming that all species coexist. Depending on the considered set of species, we also fixed the abundance ratio between some species and explored different values. We proposed a maximum likelihood estimator to infer the physical conditions and considered the effect of both the thermal noise and calibration uncertainty. We analyzed any potential biases induced by model misspecifications by comparing the results on the actual data for several sets of species and confirmed with Monte Carlo simulations. The variance of the estimations and the efficiency of the estimator were studied based on the Cram{é}r-Rao lower bound.Column densities can be estimated with 30% accuracy, while the best estimations of the volume density are found to be within a factor of two. Under the chosen model framework, the peak 12CO(1--0) is useful for constraining the kinetic temperature. The thermal pressure is better and more robustly estimated than the volume density and kinetic temperature separately. Analyzing CO and HCO+ isotopologues and fitting the full line profile are recommended practices with respect to detecting possible biases.Combining a non-local thermodynamic equilibrium model with a rigorous analysis of the accuracy allows us to obtain an efficient estimator and identify where the model is misspecified. We note that other combinations of molecular lines could be studied in the future. △ Less

Submitted 29 March, 2024; originally announced March 2024.

Comments: Astronomy and Astrophysics - A\&A, In press

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

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)

HCN emission from translucent gas and UV-illuminated cloud edges revealed by wide-field IRAM 30m maps of Orion B GMC: Revisiting its role as tracer of the dense gas reservoir for star formation

Authors: M. G. Santa-Maria, J. R. Goicoechea, J. Pety, M. Gerin, J. H. Orkisz, F. Le Petit, L. Einig, P. Palud, V. de Souza Magalhaes, I. Bešlić, L. Segal, S. Bardeau, E. Bron, P. Chainais, J. Chanussot, P. Gratier, V. V. Guzmán, A. Hughes, D. Languignon, F. Levrier, D. C. Lis, H. S. Liszt, J. Le Bourlot, Y. Oya, K. Öberg , et al. (6 additional authors not shown)

Abstract: We present 5 deg^2 (~250 pc^2) HCN, HNC, HCO+, and CO J=1-0 maps of the Orion B GMC, complemented with existing wide-field [CI] 492 GHz maps, as well as new pointed observations of rotationally excited HCN, HNC, H13CN, and HN13C lines. We detect anomalous HCN J=1-0 hyperfine structure line emission almost everywhere in the cloud. About 70% of the total HCN J=1-0 luminosity arises from gas at A_V <… ▽ More We present 5 deg^2 (~250 pc^2) HCN, HNC, HCO+, and CO J=1-0 maps of the Orion B GMC, complemented with existing wide-field [CI] 492 GHz maps, as well as new pointed observations of rotationally excited HCN, HNC, H13CN, and HN13C lines. We detect anomalous HCN J=1-0 hyperfine structure line emission almost everywhere in the cloud. About 70% of the total HCN J=1-0 luminosity arises from gas at A_V < 8 mag. The HCN/CO J=1-0 line intensity ratio shows a bimodal behavior with an inflection point at A_V < 3 mag typical of translucent gas and UV-illuminated cloud edges. We find that most of the HCN J=1-0 emission arises from extended gas with n(H2) ~< 10^4 cm^-3, even lower density gas if the ionization fraction is > 10^-5 and electron excitation dominates. This result explains the low-A_V branch of the HCN/CO J=1-0 intensity ratio distribution. Indeed, the highest HCN/CO ratios (~0.1) at A_V < 3 mag correspond to regions of high [CI] 492 GHz/CO J=1-0 intensity ratios (>1) characteristic of low-density PDRs. Enhanced FUV radiation favors the formation and excitation of HCN on large scales, not only in dense star-forming clumps. The low surface brightness HCN and HCO+ J=1-0 emission scale with I_FIR (a proxy of the stellar FUV radiation field) in a similar way. Together with CO J=1-0, these lines respond to increasing I_FIR up to G0~20. On the other hand, the bright HCN J=1-0 emission from dense gas in star-forming clumps weakly responds to I_FIR once the FUV radiation field becomes too intense (G0>1500). The different power law scalings (produced by different chemistries, densities, and line excitation regimes) in a single but spatially resolved GMC resemble the variety of Kennicutt-Schmidt law indexes found in galaxy averages. As a corollary for extragalactic studies, we conclude that high HCN/CO J=1-0 line intensity ratios do not always imply the presence of dense gas. △ Less

Submitted 18 September, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: accepted for publication in A&A. 24 pages, 18 figures, plus Appendix. Abridged Abstract

Journal ref: A&A 679, A4 (2023)

Neural network-based emulation of interstellar medium models

Authors: Pierre Palud, Lucas Einig, Franck Le Petit, Emeric Bron, Pierre Chainais, Jocelyn Chanussot, Jérôme Pety, Pierre-Antoine Thouvenin, David Languignon, Ivana Bešlić, Miriam G. Santa-Maria, Jan H. Orkisz, Léontine E. Ségal, Antoine Zakardjian, Sébastien Bardeau, Maryvonne Gerin, Javier R. Goicoechea, Pierre Gratier, Viviana V. Guzman, Annie Hughes, François Levrier, Harvey S. Liszt, Jacques Le Bourlot, Antoine Roueff, Albrecht Sievers

Abstract: The interpretation of observations of atomic and molecular tracers in the galactic and extragalactic interstellar medium (ISM) requires comparisons with state-of-the-art astrophysical models to infer some physical conditions. Usually, ISM models are too time-consuming for such inference procedures, as they call for numerous model evaluations. As a result, they are often replaced by an interpolatio… ▽ More The interpretation of observations of atomic and molecular tracers in the galactic and extragalactic interstellar medium (ISM) requires comparisons with state-of-the-art astrophysical models to infer some physical conditions. Usually, ISM models are too time-consuming for such inference procedures, as they call for numerous model evaluations. As a result, they are often replaced by an interpolation of a grid of precomputed models. We propose a new general method to derive faster, lighter, and more accurate approximations of the model from a grid of precomputed models. These emulators are defined with artificial neural networks (ANNs) designed and trained to address the specificities inherent in ISM models. Indeed, such models often predict many observables (e.g., line intensities) from just a few input physical parameters and can yield outliers due to numerical instabilities or physical bistabilities. We propose applying five strategies to address these characteristics: 1) an outlier removal procedure; 2) a clustering method that yields homogeneous subsets of lines that are simpler to predict with different ANNs; 3) a dimension reduction technique that enables to adequately size the network architecture; 4) the physical inputs are augmented with a polynomial transform to ease the learning of nonlinearities; and 5) a dense architecture to ease the learning of simple relations. We compare the proposed ANNs with standard classes of interpolation methods to emulate the Meudon PDR code, a representative ISM numerical model. Combinations of the proposed strategies outperform all interpolation methods by a factor of 2 on the average error, reaching 4.5% on the Meudon PDR code. These networks are also 1000 times faster than accurate interpolation methods and require ten to forty times less memory. This work will enable efficient inferences on wide-field multiline observations of the ISM. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Journal ref: A&A 678, A198 (2023)

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)

Deep learning denoising by dimension reduction: Application to the ORION-B line cubes

Authors: Lucas Einig, Jérôme Pety, Antoine Roueff, Paul Vandame, Jocelyn Chanussot, Maryvonne Gerin, Jan H. Orkisz, Pierre Palud, Miriam Garcia Santa-Maria, Victor de Souza Magalhaes, Ivana Bešlić, Sébastien Bardeau, Emeric E. Bron, Pierre Chainais, Javier R Goicoechea, Pierre Gratier, Viviana Guzman Veloso, Annie Hughes, Jouni Kainulainen, David Languignon, Rosine Lallement, François Levrier, Dariuscz C. Lis, Harvey Liszt, Jacques Le Bourlot , et al. (7 additional authors not shown)

Abstract: Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR… ▽ More Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR) are major challenges for consistent analysis and interpretation.Aims. We search for a denoising method of the low SNR regions of the studied data cubes that would allow to recover the low SNR emission without distorting the signals with high SNR.Methods. We perform an in-depth data analysis of the 13 CO and C 17 O (1 -- 0) data cubes obtained as part of the ORION-B large program performed at the IRAM 30m telescope. We analyse the statistical properties of the noise and the evolution of the correlation of the signal in a given frequency channel with that of the adjacent channels. This allows us to propose significant improvements of typical autoassociative neural networks, often used to denoise hyperspectral Earth remote sensing data. Applying this method to the 13 CO (1 -- 0) cube, we compare the denoised data with those derived with the multiple Gaussian fitting algorithm ROHSA, considered as the state of the art procedure for data line cubes.Results. The nature of astronomical spectral data cubes is distinct from that of the hyperspectral data usually studied in the Earth remote sensing literature because the observed intensities become statistically independent beyond a short channel separation. This lack of redundancy in data has led us to adapt the method, notably by taking into account the sparsity of the signal along the spectral axis. The application of the proposed algorithm leads to an increase of the SNR in voxels with weak signal, while preserving the spectral shape of the data in high SNR voxels.Conclusions. The proposed algorithm that combines a detailed analysis of the noise statistics with an innovative autoencoder architecture is a promising path to denoise radio-astronomy line data cubes. In the future, exploring whether a better use of the spatial correlations of the noise may further improve the denoising performances seems a promising avenue. In addition, △ Less

Submitted 24 July, 2023; originally announced July 2023.

Journal ref: A&A 677, A158 (2023)

The extremely sharp transition between molecular and ionized gas in the Horsehead nebula

Authors: C. Hernández-Vera, V. V. Guzmán, J. R. Goicoechea, V. Maillard, J. Pety, F. Le Petit, M. Gerin, E. Bron, E. Roueff, A. Abergel, T. Schirmer, J. Carpenter, P. Gratier, K. Gordon, K. Misselt

Abstract: (Abridged) Massive stars can determine the evolution of molecular clouds with their strong ultraviolet (UV) radiation fields. Moreover, UV radiation is relevant in setting the thermal gas pressure in star-forming clouds, whose influence can extend from the rims of molecular clouds to entire star-forming galaxies. Probing the fundamental structure of nearby molecular clouds is therefore crucial to… ▽ More (Abridged) Massive stars can determine the evolution of molecular clouds with their strong ultraviolet (UV) radiation fields. Moreover, UV radiation is relevant in setting the thermal gas pressure in star-forming clouds, whose influence can extend from the rims of molecular clouds to entire star-forming galaxies. Probing the fundamental structure of nearby molecular clouds is therefore crucial to understand how massive stars shape their surrounding medium and how fast molecular clouds are destroyed, specifically at their UV-illuminated edges, where models predict an intermediate zone of neutral atomic gas between the molecular cloud and the surrounding ionized gas whose size is directly related to the exposed physical conditions. We present the highest angular resolution (~$0.5$", corresponding to $207$ au) and velocity-resolved images of the molecular gas emission in the Horsehead nebula, using CO J=3-2 and HCO$^+$ J=4-3 observations with ALMA. We find that CO and HCO$^+$ are present at the edge of the cloud, very close to the ionization (H$^+$/H) and dissociation fronts (H/H$_2$), suggesting a very thin layer of neutral atomic gas (<$650$ au) and a small amount of CO-dark gas ($A_V=0.006-0.26$ mag) for stellar UV illumination conditions typical of molecular clouds in the Milky Way. The new ALMA observations reveal a web of molecular gas filaments with an estimated thermal gas pressure of $P_{\mathrm{th}} = (2.3 - 4.0) \times 10^6$ K cm$^{-3}$, and the presence of a steep density gradient at the cloud edge that can be well explained by stationary isobaric PDR models with pressures consistent with our estimations. However, in the HII region and PDR interface, we find $P_{\mathrm{th,PDR}} > P_{\mathrm{th,HII}}$, suggesting the gas is slightly compressed. Therefore, dynamical effects cannot be completely ruled out and even higher angular observations will be needed to unveil their role. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 15 pages, 1 table, 9 figures; Accepted for publication in A&A

Journal ref: A&A 677, A152 (2023)

Understanding the temperatures of H3+ and H2 in diffuse interstellar sightlines

Authors: Jacques Le Bourlot, Evelyne Roueff, Franck Le Petit, Florian Kehrein, Annika Oetjens, Holger Kreckel

Abstract: The triatomic hydrogen ion H3+ is one of the most important species for the gas phase chemistry of the interstellar medium. Observations of H3+ are used to constrain important physical and chemical parameters of interstellar environments. However, the temperatures inferred from the two lowest rotational states of H3+ in diffuse lines of sight - typically the only ones observable - appear consisten… ▽ More The triatomic hydrogen ion H3+ is one of the most important species for the gas phase chemistry of the interstellar medium. Observations of H3+ are used to constrain important physical and chemical parameters of interstellar environments. However, the temperatures inferred from the two lowest rotational states of H3+ in diffuse lines of sight - typically the only ones observable - appear consistently lower than the temperatures derived from H2 observations in the same sightlines. All previous attempts at modelling the temperatures of H3+ in the diffuse interstellar medium failed to reproduce the observational results. Here we present new studies, comparing an independent master equation for H3+ level populations to results from the Meudon PDR code for photon dominated regions. We show that the populations of the lowest rotational states of H3+ are strongly affected by the formation reaction and that H3+ ions experience incomplete thermalisation before their destruction by free electrons. Furthermore, we find that for quantitative analysis more than two levels of H3+ have to be considered and that it is crucial to include radiative transitions as well as collisions with H2. Our models of typical diffuse interstellar sightlines show very good agreement with observational data, and thus they may finally resolve the perceived temperature difference attributed to these two fundamental species. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Gas kinematics around filamentary structures in the Orion B cloud

Authors: Mathilde Gaudel, Jan H. Orkisz, Maryvonne Gerin, Jérôme Pety, Antoine Roueff, Antoine Marchal, François Levrier, Marc-Antoine Miville-Deschênes, Javier R. Goicoechea, Evelyne Roueff, Franck Le Petit, Victor de Souza Magalhaes, Pierre Palud, Miriam G. Santa-Maria, Maxime Vono, Sébastien Bardeau, Emeric Bron, Pierre Chainais, Jocelyn Chanussot, Pierre Gratier, Viviana Guzman, Annie Hughes, Jouni Kainulainen, David Languignon, Jacques Le Bourlot , et al. (5 additional authors not shown)

Abstract: Understanding the initial properties of star-forming material and how they affect the star formation process is key. From an observational point of view, the feedback from young high-mass stars on future star formation properties is still poorly constrained. In the framework of the IRAM 30m ORION-B large program, we obtained observations of the translucent and moderately dense gas, which we used t… ▽ More Understanding the initial properties of star-forming material and how they affect the star formation process is key. From an observational point of view, the feedback from young high-mass stars on future star formation properties is still poorly constrained. In the framework of the IRAM 30m ORION-B large program, we obtained observations of the translucent and moderately dense gas, which we used to analyze the kinematics over a field of 5 deg^2 around the filamentary structures. We used the ROHSA algorithm to decompose and de-noise the C18O(1-0) and 13CO(1-0) signals by taking the spatial coherence of the emission into account. We produced gas column density and mean velocity maps to estimate the relative orientation of their spatial gradients. We identified three cloud velocity layers at different systemic velocities and extracted the filaments in each velocity layer. The filaments are preferentially located in regions of low centroid velocity gradients. By comparing the relative orientation between the column density and velocity gradients of each layer from the ORION-B observations and synthetic observations from 3D kinematic toy models, we distinguish two types of behavior in the dynamics around filaments: (i) radial flows perpendicular to the filament axis that can be either inflows (increasing the filament mass) or outflows and (ii) longitudinal flows along the filament axis. The former case is seen in the Orion B data, while the latter is not identified. We have also identified asymmetrical flow patterns, usually associated with filaments located at the edge of an HII region. This is the first observational study to highlight feedback from HII regions on filament formation and, thus, on star formation in the Orion B cloud. This simple statistical method can be used for any molecular cloud to obtain coherent information on the kinematics. △ Less

Submitted 25 November, 2022; originally announced November 2022.

Comments: 45 pages. Abridged abstract. Accepted by Astronomy and Astrophysics

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

OH mid-infrared emission as a diagnostic of H$_2$O UV photodissociation. II. Application to interstellar PDRs

Authors: Marion Zannese, Benoît Tabone, Emilie Habart, Franck Le Petit, Ewine F. van Dishoeck, Emeric Bron

Abstract: Water photodissociation in the 114 - 144 nm UV range forms excited OH which emits at mid-infrared wavelengths via highly excited rotational lines. These lines have only been detected with Spitzer in several proto-planetary disks and shocks. Previous studies have shown they are a unique diagnostic for water photodissociation. Thanks to its high sensitivity and angular resolution, the James Webb Spa… ▽ More Water photodissociation in the 114 - 144 nm UV range forms excited OH which emits at mid-infrared wavelengths via highly excited rotational lines. These lines have only been detected with Spitzer in several proto-planetary disks and shocks. Previous studies have shown they are a unique diagnostic for water photodissociation. Thanks to its high sensitivity and angular resolution, the James Webb Space Telescope (JWST) could be able to detect them in other environments such as interstellar Photo-Dissociation Regions (PDRs). In order to predict the emerging spectrum of OH, we use the Meudon PDR Code to compute the thermal and chemical structure of PDRs. The influence of thermal pressure ($P_{\rm th}/k$ = $n_{\rm H} T_{\rm K}$) and UV field strength on the integrated intensities, as well as their detectability with the JWST are studied in details. OH mid-IR emission is predicted to originate very close to the H$^0$/H$_2$ transition and is directly proportional to the column density of water photodissociated in that layer. Because neutral gas-phase formation of water requires relatively high temperatures ($T_{\rm K} \gtrsim 300~$K), the resulting OH mid-IR lines are primarily correlated with the temperature at this position, and are therefore brighter in regions with high pressure. This implies that these lines are predicted to be only detectable in strongly irradiated PDRs ($G_0^{\rm incident}$ $>$ 10$^3$) with high thermal pressure ($P_{\rm th}/k$ $\gtrsim$ 5$\times$10$^7$ K cm$^{-3}$). In the latter case, OH mid-IR lines are less dependent on the strength of the incident UV field. The detection in PDRs like the Orion bar, which should be possible, is also investigated. To conclude, OH mid-IR lines observable by JWST are a promising diagnostics for dense and strongly irradiated PDRs. △ Less

Submitted 26 October, 2022; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: Accepted on the 17/10/2022 by Astronomy & Astrophysics (A&A)

Journal ref: A&A 671, A41 (2023)

Multi-line Observations, Models, and Data Needed to Understand the Nature of UV-irradiated Interstellar Matter

Authors: Javier R. Goicoechea, Sara Cuadrado, Franck Le Petit

Abstract: Far-ultraviolet photons from OB-type massive stars regulate the heating, ionization, and chemistry of much of the neutral interstellar gas in star-forming galaxies. The interaction of FUV radiation and interstellar matter takes place in environments broadly known as photodissociation regions (PDRs). PDR line diagnostics are the smoking gun of the radiative feedback from massive stars. Improving ou… ▽ More Far-ultraviolet photons from OB-type massive stars regulate the heating, ionization, and chemistry of much of the neutral interstellar gas in star-forming galaxies. The interaction of FUV radiation and interstellar matter takes place in environments broadly known as photodissociation regions (PDRs). PDR line diagnostics are the smoking gun of the radiative feedback from massive stars. Improving our understanding of stellar feedback in the ISM requires quantifying the energy budget, gas dynamics, and chemical composition of PDR environments. This goal demands astronomical instrumentation able to deliver multi-line spectroscopic images of the ISM (of the Milky Way and nearby galaxies). It also requires interdisciplinary collaborations to obtain the rate coefficients and cross sections of the many microphysical processes that occur in the ISM and that are included in models such as the Meudon PDR code. △ Less

Submitted 21 July, 2022; originally announced July 2022.

Comments: To appear in "The European Physical Journal Web of Conferences". Invited paper at "Multi-line Diagnostics of the Interstellar Medium", Nice 2022

From the Circumnuclear Disk in the Galactic Center to thick, obscuring tori of AGNs -- Modelling the molecular emission of a parsec-scale torus as found in NGC1068

Authors: B. Vollmer, R. I. Davies, P. Gratier, Th. Lizee, M. Imanishi, J. F. Gallimore, C. M. V. Impellizzeri, S. Garcia-Burillo, F. Le Petit

Abstract: The accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimetre interferometry within the inner ~20pc of the active galaxy NGC1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions a… ▽ More The accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimetre interferometry within the inner ~20pc of the active galaxy NGC1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions and determine the physical properties of their gas distribution, we present a modelling effort that includes (i) a simple dynamical simulations involving partially inelastic collisions between disk gas clouds, (ii) an analytical model of a turbulent clumpy gas disk calibrated by the dynamical model and observations, (iii) local turbulent and cosmic ray gas heating and cooling via H2O, H2, and CO emission, and (iv) determination of the molecular abundances. We also consider photodissociation regions (PDR) where gas is directly illuminated by the central engine. We compare the resulting model datacubes of the CO, HCN, HCO+, and CS brightness temperatures to available observations. In both cases the kinematics can be explained by one or two clouds colliding with a pre-existing ring, in a prograde sense for the CND and retrograde for NGC1068. And, with only dense disk clouds, the line fluxes can be reproduced to within a factor of about two. To avoid self-absorption of the intercloud medium, turbulent heating at the largest scales, comparable to the disk height, has to be decreased by a factor of 50-200. Our models indicate that turbulent mechanical energy input is the dominant gas heating mechanism within the thick gas disks. In N1068, while the bulk of the AGN X-ray radiation is absorbed in a layer of Compton-thick gas inside the dust sublimation radius, the optical/UV radiation may enhance the molecular line emission from photodissociation regions by ~50% at the inner edge of the gas ring. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: 36 pages, 32 figures, accepted for publication in A&A

Journal ref: A&A 665, A102 (2022)

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

Dynamical effects of the radiative stellar feedback on the H I-to-H2 transition

Authors: Vincent Maillard, Emeric Bron, Franck Le Petit

Abstract: The atomic-to-molecular hydrogen (H/H2) transition has been extensively studied as it controls the fraction of gas in a molecular state in an interstellar cloud. This fraction is linked to star-formation by the Schmidt-Kennicutt law. While theoretical estimates of the column density of the H I layer have been proposed for static photodissociation regions (PDRs), Herschel and well-resolved ALMA (At… ▽ More The atomic-to-molecular hydrogen (H/H2) transition has been extensively studied as it controls the fraction of gas in a molecular state in an interstellar cloud. This fraction is linked to star-formation by the Schmidt-Kennicutt law. While theoretical estimates of the column density of the H I layer have been proposed for static photodissociation regions (PDRs), Herschel and well-resolved ALMA (Atacama Large Millimeter Array) observations have revealed dynamical effects in star forming regions, caused by the process of photoevaporation. We extend the analytic study of the H/H2 transition to include the effects of the propagation of the ionization front, in particular in the presence of photoevaporation at the walls of blister H II regions, and we find its consequences on the total atomic hydrogen column density at the surface of clouds in the presence of an ultraviolet field, and on the properties of the H/H2 transition. We solved semi-analytically the differential equation giving the H2 column density profile by taking into account H2 formation on grains, H2 photodissociation, and the ionization front propagation dynamics modeled as advection of the gas through the ionization front. Taking this advection into account reduces the width of the atomic region compared to static models. The atomic region may disappear if the ionization front velocity exceeds a certain value, leading the H/H2 transition and the ionization front to merge. For both dissociated and merged configurations, we provide analytical expressions to determine the total H I column density. Our results take the metallicity into account. Finally, we compared our results to observations of PDRs illuminated by O-stars, for which we conclude that the dynamical effects are strong, especially for low-excitation PDRs. △ Less

Submitted 15 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: Accepted by Astronomy & Astrophysics on September, 11th 2021. 18 pages, 16 figures

Journal ref: A&A 656, A65 (2021)

Influence of the nano-grain depletion in photon-dominated regions: Application to the gas physics and chemistry in the Horsehead

Authors: T. Schirmer, E. Habart, N. Ysard, E. Bron, J. Le Bourlot, L. Verstraete, A. Abergel, A. P. Jones, E. Roueff, F. Le Petit

Abstract: The large disparity in physical conditions from the diffuse interstellar medium (ISM) to denser clouds such as photon-dominated regions (PDRs) triggers an evolution of the dust properties (i.e. composition, size, and shape). The gas physics and chemistry are tightly connected to these dust properties and are therefore affected by dust evolution and especially the nano-grain depletion in the outer… ▽ More The large disparity in physical conditions from the diffuse interstellar medium (ISM) to denser clouds such as photon-dominated regions (PDRs) triggers an evolution of the dust properties (i.e. composition, size, and shape). The gas physics and chemistry are tightly connected to these dust properties and are therefore affected by dust evolution and especially the nano-grain depletion in the outer irradiated part of PDRs. We highlight the influence of nano-grain depletion on the gas physics and chemistry in the Horsehead nebula, a prototypical PDR. We used a model for atomic and molecular gas in PDRs, the Meudon PDR code, using diffuse ISM-like dust and Horsehead-like dust to study the influence of nano-grain depletion on the gas physics and chemistry, focusing on the impact on photoelectric heating and H2 formation and, therefore, on the H2 gas lines. We find that nano-grain depletion in the Horsehead strongly affects gas heating through the photoelectric effect and thus the gas temperature and the H2 formation, hence the H -> H2 position. Consequently, the first four pure rotational lines of H2 (e.g. 0-0 S(0), S(1), S(2), and S(3)) vary by a factor of 2 to 14. The 0-0 S(3) line that is often underestimated in models is underestimated even more when taking nano-grain depletion into account due to the decrease in gas heating through the photoelectric effect. This strongly suggests that our understanding of the excitation of H2 and/or of heating processes in the Horsehead, and more generally in PDRs, is still incomplete. Nano-grain depletion in the outer part of the Horsehead has a strong influence on several gas tracers that will be prominent in JWST observations of irradiated clouds. We therefore need to take this depletion into account in order to improve our understanding of the Horsehead, and more generally PDRs, and to contribute to the optimal scientific return of the mission. △ Less

Submitted 16 April, 2021; originally announced April 2021.

Comments: Accepted in A&A, comments welcome

Journal ref: A&A 649, A148 (2021)

A molecular wind blows out of the Kuiper belt

Authors: Quentin Kral, J. E. Pringle, Aurélie Guilbert-Lepoutre, Luca Matrà, Julianne I. Moses, Emmanuel Lellouch, Mark C. Wyatt, Nicolas Biver, Dominique Bockelée-Morvan, Amy Bonsor, Franck Le Petit, G. Randall Gladstone

Abstract: Gas has been detected in many exoplanetary systems ($>$10 Myr), thought to be released in the destruction of volatile-rich planetesimals orbiting in exo-Kuiper belts. In this letter, we aim to explore whether gas is also expected in the Kuiper belt (KB) in our Solar System. To quantify the gas release in our Solar System, we use models for gas release that have been applied to extrasolar planetary… ▽ More Gas has been detected in many exoplanetary systems ($>$10 Myr), thought to be released in the destruction of volatile-rich planetesimals orbiting in exo-Kuiper belts. In this letter, we aim to explore whether gas is also expected in the Kuiper belt (KB) in our Solar System. To quantify the gas release in our Solar System, we use models for gas release that have been applied to extrasolar planetary systems, as well as a physical model that accounts for gas released due to the progressive internal warming of large planetesimals. We find that only bodies larger than about 4 km can still contain CO ice after 4.6 Gyr of evolution. This finding may provide a clue as to why Jupiter-family comets, thought to originate in the Kuiper belt, are deficient in CO compared to Oort-clouds comets. We predict that gas is still produced in the KB right now at a rate of $2 \times 10^{-8}$ M$_\oplus$/Myr for CO and orders of magnitude more when the Sun was younger. Once released, the gas is quickly pushed out by the Solar wind. Therefore, we predict a gas wind in our Solar System starting at the KB location and extending far beyond with regards to the heliosphere with a current total CO mass of $\sim 2 \times 10^{-12}$ M$_\oplus$. We also predict the existence of a slightly more massive atomic gas wind made of carbon and oxygen (neutral and ionized) with a mass of $\sim 10^{-11}$ M$_\oplus$. We predict that gas is currently present in our Solar System beyond the Kuiper belt and that although it cannot be detected with current instrumentation, it could be observed in the future with an in situ mission using an instrument similar to Alice on New Horizons with larger detectors. Our model of gas release due to slow heating may also work for exoplanetary systems and provide the first real physical mechanism for the gas observations. △ Less

Submitted 9 September, 2021; v1 submitted 2 April, 2021; originally announced April 2021.

Comments: accepted for publication as a Letter to the editor in A&A; abstract shortened; 15 pages, 5 figures, 4 tables

Journal ref: A&A 653, L11 (2021)

Quantitative inference of the $H_2$ column densities from 3 mm molecular emission: A case study towards Orion B

Authors: Pierre Gratier, Jérôme Pety, Emeric Bron, Antoine Roueff, Jan H. Orkisz, Maryvonne Gerin, Victor de Souza Magalhaes, Mathilde Gaudel, Maxime Vono, Sébastien Bardeau, Jocelyn Chanussot, Pierre Chainais, Javier R. Goicoechea, Viviana V. Guzmán, Annie Hughes, Jouni Kainulainen, David Languignon, Jacques Le Bourlot, Franck Le Petit, François Levrier, Harvey Liszt, Nicolas Peretto, Evelyne Roueff, Albrecht Sievers

Abstract: Molecular hydrogen being unobservable in cold molecular clouds, the column density measurements of molecular gas currently rely either on dust emission observation in the far-IR or on star counting. (Sub-)millimeter observations of numerous trace molecules are effective from ground based telescopes, but the relationships between the emission of one molecular line and the H2 column density (NH2) is… ▽ More Molecular hydrogen being unobservable in cold molecular clouds, the column density measurements of molecular gas currently rely either on dust emission observation in the far-IR or on star counting. (Sub-)millimeter observations of numerous trace molecules are effective from ground based telescopes, but the relationships between the emission of one molecular line and the H2 column density (NH2) is non-linear and sensitive to excitation conditions, optical depths, abundance variations due to the underlying physico-chemistry. We aim to use multi-molecule line emission to infer NH2 from radio observations. We propose a data-driven approach to determine NH2 from radio molecular line observations. We use supervised machine learning methods (Random Forests) on wide-field hyperspectral IRAM-30m observations of the Orion B molecular cloud to train a predictor of NH2, using a limited set of molecular lines as input, and the Herschel-based dust-derived NH2 as ground truth output. For conditions similar to the Orion B molecular cloud, we obtain predictions of NH2 within a typical factor of 1.2 from the Herschel-based estimates. An analysis of the contributions of the different lines to the predictions show that the most important lines are $^{13}$CO(1-0), $^{12}$CO(1-0), C$^{18}$O(1-0), and HCO$^+$(1-0). A detailed analysis distinguishing between diffuse, translucent, filamentary, and dense core conditions show that the importance of these four lines depends on the regime, and that it is recommended to add the N$_2$H$^+$(1-0) and CH$_3$OH(20-10) lines for the prediction of NH2 in dense core conditions. This article opens a promising avenue to directly infer important physical parameters from the molecular line emission in the millimeter domain. The next step will be to try to infer several parameters simultaneously (e.g., NH2 and far-UV illumination field) to further test the method. [Abridged] △ Less

Submitted 31 August, 2020; originally announced August 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 645, A27 (2021)

Tracers of the ionization fraction in dense and translucent gas: I. Automated exploitation of massive astrochemical model grids

Authors: Emeric Bron, Evelyne Roueff, Maryvonne Gerin, Jérôme Pety, Pierre Gratier, Franck Le Petit, Viviana Guzman, Jan H. Orkisz, Victor de Souza Magalhaes, Mathilde Gaudel, Maxime Vono, Sébastien Bardeau, Pierre Chainais, Javier R. Goicoechea, Annie Hughes, Jouni Kainulainen, David Languignon, Jacques Le Bourlot, François Levrier, Harvey Liszt, Karin Öberg, Nicolas Peretto, Antoine Roueff, Albrecht Sievers

Abstract: The ionization fraction plays a key role in the physics and chemistry of the neutral interstellar medium, from controlling the coupling of the gas to the magnetic field to allowing fast ion-neutral reactions that drive interstellar chemistry. Most estimations of the ionization fraction have relied on deuterated species such as DCO+, whose detection is limited to dense cores representing an extreme… ▽ More The ionization fraction plays a key role in the physics and chemistry of the neutral interstellar medium, from controlling the coupling of the gas to the magnetic field to allowing fast ion-neutral reactions that drive interstellar chemistry. Most estimations of the ionization fraction have relied on deuterated species such as DCO+, whose detection is limited to dense cores representing an extremely small fraction of the volume of the giant molecular clouds they are part of. As large field-of-view hyperspectral maps become available, new tracers may be found. We search for the best observable tracers of the ionization fraction based on a grid of astrochemical models. We build grids of models that sample randomly a large space of physical conditions (unobservable quantities such as gas density, temperature, etc.) and compute the corresponding observables (line intensities, column densities) and the ionization fraction. We estimate the predictive power of each potential tracer by training a Random Forest model to predict the ionization fraction from that tracer, based on these model grids. In both translucent medium and cold dense medium conditions, several observable tracers with very good predictive power for the ionization fraction are found. Several tracers in cold dense medium conditions are found to be better and more widely applicable than the traditional DCO+/HCO+ ratio. We also provide simpler analytical fits for estimating the ionization fraction from the best tracers, and for estimating the associated uncertainties. We discuss the limitations of the present study and select a few recommended tracers in both types of conditions. The method presented here is very general and can be applied to the measurement of any other quantity of interest (cosmic ray flux, elemental abundances, etc.) from any type of model (PDR models, time-dependent chemical models, etc.). (abridged) △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 29 pages, 15 figures. Accepted for publication in A&A

Journal ref: A&A 645, A28 (2021)

C18O, 13CO, and 12CO abundances and excitation temperatures in the Orion B molecular cloud: An analysis of the precision achievable when modeling spectral line within the Local Thermodynamic Equilibrium approximation

Authors: Antoine Roueff, Maryvonne Gerin, Pierre Gratier, Francois Levrier, Jerome Pety, Mathilde Gaudel, Javier R. Goicoechea, Jan H. Orkisz, Victor de Souza Magalhaes, Maxime Vono, Sebastien Bardeau, Emeric Bron, Jocelyn Chanussot, Pierre Chainais, Viviana V. Guzman, Annie Hughes, Jouni Kainulainen, David Languignon, Jacques Le Bourlot, Franck Le Petit, Harvey S. Liszt, Antoine Marchal, Marc-Antoine Miville-Deschenes, Nicolas Peretto, Evelyne Roueff , et al. (1 additional authors not shown)

Abstract: CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the C… ▽ More CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the Cramer Rao Bound (CRB) technique to analyze and estimate the precision of the physical parameters in the framework of local-thermodynamic-equilibrium excitation and radiative transfer with an additive white Gaussian noise. We propose a maximum likelihood estimator to infer the physical conditions from the 1-0 and 2-1 transitions of CO isotopologues. Simulations show that this estimator is unbiased and efficient for a common range of excitation temperatures and column densities (Tex > 6 K, N > 1e14 - 1e15 cm-2). Contrary to the general assumptions, the different CO isotopologues have distinct excitation temperatures, and the line intensity ratios between different isotopologues do not accurately reflect the column density ratios. We find mean fractional abundances that are consistent with previous determinations towards other molecular clouds. However, significant local deviations are inferred, not only in regions exposed to UV radiation field but also in shielded regions. These deviations result from the competition between selective photodissociation, chemical fractionation, and depletion on grain surfaces. We observe that the velocity dispersion of the C18O emission is 10% smaller than that of 13CO. The substantial gain resulting from the simultaneous analysis of two different rotational transitions of the same species is rigorously quantified. The CRB technique is a promising avenue for analyzing the estimation of physical parameters from the fit of spectral lines. △ Less

Submitted 17 May, 2020; originally announced May 2020.

Comments: 27 pages, 23 PDF figures. Accepted for publication in A&A. Uses aa latex macro

Journal ref: A&A 645, A26 (2021)

The role of C/O in nitrile astrochemistry in PDRs and planet-forming disks

Authors: Romane Le Gal, Madison T. Brady, Karin I. Öberg, Evelyne Roueff, Franck Le Petit

Abstract: Complex nitriles, such as HC3N, and CH3CN, are observed in a wide variety of astrophysical environments, including at relatively high abundances in photon-dominated regions (PDR) and the UV exposed atmospheres of planet-forming disks. The latter have been inferred to be oxygen-poor, suggesting that these observations may be explained by organic chemistry in C-rich environments. In this study we fi… ▽ More Complex nitriles, such as HC3N, and CH3CN, are observed in a wide variety of astrophysical environments, including at relatively high abundances in photon-dominated regions (PDR) and the UV exposed atmospheres of planet-forming disks. The latter have been inferred to be oxygen-poor, suggesting that these observations may be explained by organic chemistry in C-rich environments. In this study we first explore if the PDR complex nitrile observations can be explained by gas-phase PDR chemistry alone if the elemental C/O ratio is elevated. In the case of the Horsehead PDR, we find that gas-phase chemistry with C/O $\gtrsim$ 0.9 can indeed explain the observed nitrile abundances, increasing predicted abundances by several orders of magnitude compared to standard C/O assumptions. We also find that the nitrile abundances are sensitive to the cosmic ray ionization treatment, and provide constraints on the branching ratios between CH3CN and CH3NC productions. In a fiducial disk model, an elevated C/O ratio increases the CH3CN and HC3N productions by more than an order of magnitude, bringing abundance predictions within an order of magnitude to what has been inferred from observations. The C/O ratio appears to be a key variable in predicting and interpreting complex organic molecule abundances in photon-dominated regions across a range of scales. △ Less

Submitted 3 October, 2019; originally announced October 2019.

Comments: 15 pages, 7 figures, 2 tables, accepted for publication in ApJ

ESA Voyage 2050 white paper: A complete census of the gas phases in and around galaxies, far-UV spectropolarimetry as a prime tool for understanding galaxy evolution and star formation

Authors: V. Lebouteiller, C. Gry. H. Yan, P. Richter, B. Godard, E. B. Jenkins, D. Welty, N. Lehner, P. Guillard, J. Roman-Duval, E Roueff, F. Leone, D. Kunth, J. C. Howk, P. Boissé, F. Boulanger, E. Bron, B. James, J. Le Bourlot, F. Le Petit, M. Pieri, V. Valdivia

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

Submitted 6 September, 2019; originally announced September 2019.

Comments: This ESA white paper draws from and expands on the following white papers previously submitted for the Astro 2020 Decadal Survey, Gry et al. (2019), Lebouteiller et al. (2019), and Yan et al. (2019). The original ESA white paper was reformatted as a 2-column paper for the present document. (arXiv admin note: substantial text overlap with arXiv:1903.07295, arXiv:1903.06250, arXiv:1903.11065)

The water line emission and ortho-to-para ratio in the Orion Bar photon-dominated region

Authors: T. Putaud, X. Michaut, F. Le Petit, E. Roueff, D. C. Lis

Abstract: A very low ortho-to-para ratio (OPR) of 0.1-0.5 was previously reported in the Orion Bar photon-dominated region (PDR), based on observations of two optically thin $\mathrm{H_2^{18}O}$ lines which were analyzed by using a single-slab large velocity gradient model. The corresponding spin temperature does not coincide with the kinetic temperature of the molecular gas in this UV-illuminated region. T… ▽ More A very low ortho-to-para ratio (OPR) of 0.1-0.5 was previously reported in the Orion Bar photon-dominated region (PDR), based on observations of two optically thin $\mathrm{H_2^{18}O}$ lines which were analyzed by using a single-slab large velocity gradient model. The corresponding spin temperature does not coincide with the kinetic temperature of the molecular gas in this UV-illuminated region. This was interpreted as an indication of water molecules being formed on cold icy grains which were subsequently released by UV photodesorption. A more complete set of water observations in the Orion Bar, including seven $\mathrm{H_2^{16}O}$ lines and one $\mathrm{H_2^{18}O}$ line, carried out using Herschel/HIFI instrument, was reanalyzed using the Meudon PDR code to derive gas-phase water abundance and the OPR, taking into account the steep density and temperature gradients present in the region. The model line intensities are in good agreement with the observations assuming that water molecules formed with an OPR corresponding to thermal equilibrium conditions at the local kinetic temperature of the gas and when solely considering gas-phase chemistry and water gas-grain exchanges through adsorption and desorption. Gas-phase water is predicted to arise from a region deep into the cloud, corresponding to a visual extinction of $A_{\mathrm{V}} \sim 9$, with a $\mathrm{H_2^{16}O}$ fractional abundance of $\sim 2\times 10^{-7}$ and column density of $(1.4 \pm 0.8) \times 10^{15}$ cm$^{-2}$ for a total cloud depth of $A_{\mathrm{V}}=15$. A line-of-sight average ortho-to-para ratio of $2.8 \pm 0.2$ is derived. The observational data are consistent with a nuclear spin isomer repartition corresponding to the thermal equilibrium at a temperature of $36 \pm 2$ K, much higher than the spin temperature previously reported for this region and close to the gas kinetic temperature in the water-emitting gas. △ Less

Submitted 7 October, 2019; v1 submitted 1 August, 2019; originally announced August 2019.

Journal ref: A&A 632, A8 (2019)

Radiative and mechanical feedback into the molecular gas in the Large Magellanic Cloud. II. 30 Doradus

Authors: Min-Young Lee, Suzanne Madden, Franck Le Petit, Antoine Gusdorf, Pierre Lesaffre, Ronin Wu, Vianney Lebouteiller, Frederic Galliano, Melanie Chevance

Abstract: With an aim of probing the physical conditions and excitation mechanisms of warm molecular gas in individual star-forming regions, we performed Herschel SPIRE FTS observations of 30 Doradus in the LMC. In our FTS observations, important FIR cooling lines in the ISM, including CO J=4-3 to 13-12, [CI] 370 micron, and [NII] 205 micron, were clearly detected. In combination with ground-based CO data,… ▽ More With an aim of probing the physical conditions and excitation mechanisms of warm molecular gas in individual star-forming regions, we performed Herschel SPIRE FTS observations of 30 Doradus in the LMC. In our FTS observations, important FIR cooling lines in the ISM, including CO J=4-3 to 13-12, [CI] 370 micron, and [NII] 205 micron, were clearly detected. In combination with ground-based CO data, we then constructed CO spectral line energy distributions (SLEDs) on 10 pc scales over a 60 pc x 60 pc area and found that the shape of the observed CO SLEDs considerably changes across 30 Doradus, e.g., the peak transition varies from J=6-5 to 10-9, while the slope characterized by the high-to-intermediate J ratio ranges from 0.4 to 1.8. To examine the source(s) of these variations in CO transitions, we analyzed the CO observations, along with [CII] 158 micron, [CI] 370 micron, [OI] 145 micron, H2 0-0 S(3), and FIR luminosity data, using state-of-the-art models of PDRs and shocks. Our detailed modeling showed that the observed CO emission likely originates from highly-compressed (thermal pressure ~ 1e7-1e9 K cm-3) clumps on 0.7-2 pc scales, which could be produced by either UV photons (UV radiation field ~ 1e3-1e5 Mathis fields) or low-velocity C-type shocks (pre-shock medium density ~ 1e4-1e6 cm-3 and shock velocity ~ 5-10 km s-1). Considering the stellar content in 30 Doradus, however, we tentatively excluded the stellar origin of CO excitation and concluded that low-velocity shocks driven by kpc scale processes (e.g., interaction between the Milky Way and the Magellanic Clouds) are likely the dominant source of heating for CO. The shocked CO-bright medium was then found to be warm (temperature ~ 100-500 K) and surrounded by a UV-regulated low pressure component (a few (1e4-1e5) K cm-3) that is bright in [CII] 158 micron, [CI] 370 micron, [OI] 145 micron, and FIR dust continuum emission. △ Less

Submitted 24 May, 2019; originally announced May 2019.

Comments: Accepted for publication in A&A

Journal ref: A&A 628, A113 (2019)

Vacuum ultraviolet photodesorption and photofragmentation of formaldehyde-containing ices

Authors: Géraldine Féraud, Mathieu Bertin, Claire Romanzin, Rémi Dupuy, Franck Le Petit, Evelyne Roueff, Laurent Philippe, Xavier Michaut, Pascal Jeseck, Jean-Hugues Fillion

Abstract: Non-thermal desorption from icy grains containing H$_2$CO has been invoked to explain the observed H$_2$CO gas phase abundances in ProtoPlanetary Disks (PPDs) and Photon Dominated Regions (PDRs). Photodesorption is thought to play a key role, however no absolute measurement of the photodesorption from H$_2$CO ices were performed up to now, so that a default value is used in the current astrophysic… ▽ More Non-thermal desorption from icy grains containing H$_2$CO has been invoked to explain the observed H$_2$CO gas phase abundances in ProtoPlanetary Disks (PPDs) and Photon Dominated Regions (PDRs). Photodesorption is thought to play a key role, however no absolute measurement of the photodesorption from H$_2$CO ices were performed up to now, so that a default value is used in the current astrophysical models. As photodesorption yields differ from one molecule to the other, it is crucial to experimentally investigate photodesorption from H$_2$CO ices. We measured absolute wavelength-resolved photodesorption yields from pure H$_2$CO ices, H$_2$CO on top of a CO ice (H$_2$CO/CO), and H$_2$CO mixed with CO ice (H$_2$CO:CO) irradiated in the Vacuum UltraViolet (VUV) range (7-13.6~eV). Photodesorption from a pure H$_2$CO ice releases H$_2$CO in the gas phase, but also fragments, such as CO and H$_2$. Energy-resolved photodesorption spectra, coupled with InfraRed (IR) and Temperature Programmed Desorption (TPD) diagnostics, showed the important role played by photodissociation and allowed to discuss photodesorption mechanisms. For the release of H$_2$CO in the gas phase, they include Desorption Induced by Electronic Transitions (DIET), indirect DIET through CO-induced desorption of H$_2$CO and photochemical desorption. We found that H$_2$CO photodesorbs with an average efficiency of $\sim 4-10 \times 10^{-4}$ molecule/photon, in various astrophysical environments. H$_2$CO and CO photodesorption yields and photodesorption mechanisms, involving photofragmentation of H$_2$CO, can be implemented in astrochemical codes. The effects of photodesorption on gas/solid abundances of H$_2$CO and all linked species from CO to Complex Organic Molecules (COMs), and on the H$_2$CO snowline location, are now on the verge of being unravelled. △ Less

Submitted 20 May, 2019; originally announced May 2019.

Comments: ACS Earth and Space Chemistry; Complex Organic Molecules (COMs) in Star-Forming Regions special issue

Oxygen fractionation in dense molecular clouds

Authors: Jean-Christophe Loison, Valentine Wakelam, Pierre Gratier, Kevin M. Hickson, Aurore Bacmann, Marcelino Agùndez, Nuria Marcelino, José Cernicharo, Viviana Guzman, Maryvonne Gerin, Javier R. Goicoechea, Evelyne Roueff, Franck Le Petit, Jérome Pety, Asunción Fuente, Pablo Riviere-Marichalar

Abstract: We have developed the first gas-grain chemical model for oxygen fractionation (also including sulphur fractionation) in dense molecular clouds, demonstrating that gas-phase chemistry generates variable oxygen fractionation levels, with a particularly strong effect for NO, SO, O2, and SO2. This large effect is due to the efficiency of the neutral 18O + NO, 18O + SO, and 18O + O2 exchange reactions.… ▽ More We have developed the first gas-grain chemical model for oxygen fractionation (also including sulphur fractionation) in dense molecular clouds, demonstrating that gas-phase chemistry generates variable oxygen fractionation levels, with a particularly strong effect for NO, SO, O2, and SO2. This large effect is due to the efficiency of the neutral 18O + NO, 18O + SO, and 18O + O2 exchange reactions. The modeling results were compared to new and existing observed isotopic ratios in a selection of cold cores. The good agreement between model and observations requires that the gas-phase abundance of neutral oxygen atoms is large in the observed regions. The S16O/S18O ratio is predicted to vary substantially over time showing that it can be used as a sensitive chemical proxy for matter evolution in dense molecular clouds. △ Less

Submitted 23 February, 2019; originally announced February 2019.

Comments: 22 pages, 5 figures

A dynamically young, gravitationally stable network of filaments in Orion B

Authors: Jan H. Orkisz, Nicolas Peretto, Jérôme Pety, Maryvonne Gerin, François Levrier, Emeric Bron, Sébastien Bardeau, Javier R. Goicoechea, Pierre Gratier, Viviana V. Guzmán, Annie Hughes, David Languignon, Franck Le Petit, Harvey S. Liszt, Karin Öberg, Evelyne Roueff, Albrecht Sievers, Pascal Tremblin

Abstract: Filaments are a key step on the path that leads from molecular clouds to star formation. However, their characteristics are heavily debated, and the exact processes that lead to their formation and fragmentation into dense cores still remain to be fully understood. We aim at characterising the mass, kinematics, and stability against gravitational collapse of a statistically significant sample of f… ▽ More Filaments are a key step on the path that leads from molecular clouds to star formation. However, their characteristics are heavily debated, and the exact processes that lead to their formation and fragmentation into dense cores still remain to be fully understood. We aim at characterising the mass, kinematics, and stability against gravitational collapse of a statistically significant sample of filaments in the Orion B molecular cloud, which is renown for its very low star formation efficiency. We characterise the gas column densities and kinematics over a field of 1.9 deg$^2$, using C$^{18}$O(J=1-0) data from the IRAM-30m large programme ORION-B. Using two different Hessian-based filters, we extract and compare two filamentary networks, each containing over 100 filaments. Independent of the extraction method, the filaments have widths of 0.12$\pm$0.04 pc, and show a wide range of linear (1 - 100 $M_{\odot}$pc$^{-1}$) and volume densities (2.10$^3$ - 2.10$^5$ cm$^{-3}$). Compared to previous studies, the filament population is dominated by low-density, thermally sub-critical structures, suggesting that most of the identified filaments are not collapsing to form stars. In fact, only ~1% of the Orion B cloud mass covered by our observations can be found in super-critical, star-forming filaments, explaining the low star formation efficiency of the region. The velocity profiles observed across the filaments show quiescence in the centre, and coherency in the plane of the sky, despite being mostly supersonic. The filaments in Orion B apparently belong to a continuum which contains a few elements comparable to already studied star-forming filaments as well as many lower-density, gravitationally unbound structures. This comprehensive study of the Orion B filaments shows that the mass fraction in super-critical filaments is a key factor in determining star formation efficiency. △ Less

Submitted 11 March, 2019; v1 submitted 6 February, 2019; originally announced February 2019.

Comments: 22 pages, 15 figures. Accepted by A&A

Journal ref: A&A 624, A113 (2019)

The chemistry of disks around T Tauri and Herbig Ae/Be stars

Authors: Marcelino Agundez, Evelyne Roueff, Franck Le Petit, Jacques Le Bourlot

Abstract: Infrared and (sub-)mm observations of disks around T Tauri and Herbig Ae/Be stars point to a chemical differentiation between both types of disks, with a lower detection rate of molecules in disks around hotter stars. To investigate the potential underlying causes we perform a comparative study of the chemistry of T Tauri and Herbig Ae/Be disks, using a model that pays special attention to photoch… ▽ More Infrared and (sub-)mm observations of disks around T Tauri and Herbig Ae/Be stars point to a chemical differentiation between both types of disks, with a lower detection rate of molecules in disks around hotter stars. To investigate the potential underlying causes we perform a comparative study of the chemistry of T Tauri and Herbig Ae/Be disks, using a model that pays special attention to photochemistry. The warmer disk temperatures and higher ultraviolet flux of Herbig stars compared to T Tauri stars induce some differences in the disk chemistry. In the hot inner regions, H2O, and simple organic molecules like C2H2, HCN, and CH4 are predicted to be very abundant in T Tauri disks and even more in Herbig Ae/Be disks, in contrast with infrared observations that find a much lower detection rate of water and simple organics toward disks around hotter stars. In the outer regions, the model indicates that the molecules typically observed in disks, like HCN, CN, C2H, H2CO, CS, SO, and HCO+, do not have drastic abundance differences between T Tauri and Herbig Ae disks. Some species produced under the action of photochemistry, like C2H and CN, are predicted to have slightly lower abundances around Herbig Ae stars due to a narrowing of the photochemically active layer. Observations indeed suggest that these radicals are somewhat less abundant in Herbig Ae disks, although in any case the inferred abundance differences are small, of a factor of a few at most. A clear chemical differentiation between both types of disks concerns ices, which are expected to be more abundant in Herbig Ae disks. The global chemical behavior of T Tauri and Herbig Ae/Be disks is quite similar. The main differences are driven by the warmer temperatures of the latter, which result in a larger reservoir or water and simple organics in the inner regions and a lower mass of ices in the outer disk. △ Less

Submitted 26 March, 2018; originally announced March 2018.

Comments: Accepted for publication in A&A. Abstract abridged

Journal ref: A&A 616, A19 (2018)

Structure of photodissociation fronts in star-forming regions revealed by observations of high-J CO emission lines with Herschel

Authors: C. Joblin, E. Bron, C. Pinto, P. Pilleri, F. Le Petit, M. Gerin, J. Le Bourlot, A. Fuente, O. Berne, J. R. Goicoechea, E. Habart, M. Koehler, D. Teyssier, Z. Nagy, J. Montillaud, C. Vastel, J. Cernicharo, M. Roellig, V. Ossenkopf-Okada, E. A. Bergin

Abstract: In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J… ▽ More In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J lines that trace warm molecular gas at PDR interfaces. We observed the 12CO and 13CO ladders in two prototypical PDRs, the Orion Bar and NGC 7023 NW using the instruments onboard Herschel. We also considered line emission from key species in the gas cooling of PDRs (C+, O, H2) and other tracers of PDR edges such as OH and CH+. All the intensities are collected from Herschel observations, the literature and the Spitzer archive and are analyzed using the Meudon PDR code. A grid of models was run to explore the parameter space of only two parameters: thermal gas pressure and a global scaling factor that corrects for approximations in the assumed geometry. We conclude that the emission in the high-J CO lines, which were observed up to Jup=23 in the Orion Bar (Jup=19 in NGC7023), can only originate from small structures of typical thickness of a few 1e-3 pc and at high thermal pressures (Pth~1e8 K cm-3). Compiling data from the literature, we found that the gas thermal pressure increases with the intensity of the UV radiation field given by G0, following a trend in line with recent simulations of the photoevaporation of illuminated edges of molecular clouds. This relation can help rationalising the analysis of high-J CO emission in massive star formation and provides an observational constraint for models that study stellar feedback on molecular clouds. △ Less

Submitted 15 April, 2018; v1 submitted 11 January, 2018; originally announced January 2018.

Comments: Accepted for publication in A&A; 20 pages, 11 figures

Journal ref: A&A 615, A129 (2018)

Constraining physical conditions for the PDR of Trumpler 14 in the Carina Nebula

Authors: Ronin Wu, Emeric Bron, Takashi Onaka, Franck Le Petit, Frederic Galliano, David Languignon, Tomohiko Nakamura, Yoko Okada

Abstract: We investigate the physical conditions of the CO gas near the young star cluster, Trumpler 14 of the Carina Nebula. The observations presented in this work are taken with the Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) onboard the Herschel Space Observatory. Our field of view covers the edge of a cavity carved by Trumpler 14 about… ▽ More We investigate the physical conditions of the CO gas near the young star cluster, Trumpler 14 of the Carina Nebula. The observations presented in this work are taken with the Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) onboard the Herschel Space Observatory. Our field of view covers the edge of a cavity carved by Trumpler 14 about $1\,\mathrm{Myr}$ ago and marks the transition from HII regions to photo-dissociation regions. With the state-of-the-art Meudon PDR code, we successfully derive the physical conditions, which include the thermal pressure ($P$) and the scaling factor of radiation fields ($G_{\mathrm{UV}}$), from the observed CO spectral line energy distributions~(SLEDs) in the observed region. The derived $G_{\mathrm{UV}}$ values generally show an excellent agreement with the UV radiation fields created by nearby OB-stars and thus confirm that the main excitation source of the observed CO emission are the UV-photons provided by the massive stars. The derived thermal pressure is between $0.5-3\,\times\,10^{8}\,\mathrm{K\,cm^{-3}}$ with the highest values found along the ionization front in Car I-E region facing Trumpler 14, hinting that the cloud structure is similar to the recent observations of the Orion Bar. Comparing the derived thermal pressure with the radiation fields, we report the first observationally-derived and spatially-resolved $P \sim 2\times10^4\,G_{\mathrm{UV}}$ relationship. As direct comparisons of the modeling results to the observed $^{13}\mathrm{CO}$, [OI] $63\,\mathrm{μm}$, and [CII] $158\,\mathrm{μm}$ intensities are not straightforward, we urge the readers to be cautious when constraining the physical conditions of PDRs with combinations of $^{12}\mathrm{CO}$, $^{13}\mathrm{CO}$, [CI], [OI] $63\,\mathrm{μm}$, and [CII] $158\,\mathrm{μm}$ observations. △ Less

Submitted 27 February, 2018; v1 submitted 5 January, 2018; originally announced January 2018.

Comments: Accepted for publication by A&A, comments are welcome

Journal ref: A&A 618, A53 (2018)

The ESO Diffuse Interstellar Bands Large Exploration Survey: EDIBLES I. Project description, survey sample and quality assessment

Authors: Nick Cox, Jan Cami, Amin Farhang, Jonathan Smoker, Ana Monreal-Ibero, Rosine Lallement, Peter Sarre, Charlotte Marshall, Keith Smith, Christopher Evans, Pierre Royer, Harold Linnartz, Martin Cordiner, Christine Joblin, Jacco van Loon, Bernard Foing, Neil Bhatt, Emeric Bron, Meriem Elyajouri, Alex de Koter, Pascale Ehrenfreund, Atefeh Javadi, Lex Kaper, Habib Khosroshadi, Mike Laverick , et al. (5 additional authors not shown)

Abstract: The carriers of the diffuse interstellar bands (DIBs) are largely unidentified molecules ubiquitously present in the interstellar medium (ISM). After decades of study, two strong and possibly three weak near-infrared DIBs have recently been attributed to the C60+ fullerene based on observational and laboratory measurements. There is great promise for the identification of the over 400 other known… ▽ More The carriers of the diffuse interstellar bands (DIBs) are largely unidentified molecules ubiquitously present in the interstellar medium (ISM). After decades of study, two strong and possibly three weak near-infrared DIBs have recently been attributed to the C60+ fullerene based on observational and laboratory measurements. There is great promise for the identification of the over 400 other known DIBs, as this result could provide chemical hints towards other possible carriers. In an effort to systematically study the properties of the DIB carriers, we have initiated a new large-scale observational survey: the ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES). The main objective is to build on and extend existing DIB surveys to make a major step forward in characterising the physical and chemical conditions for a statistically significant sample of interstellar lines-of-sight, with the goal to reverse-engineer key molecular properties of the DIB carriers. EDIBLES is a filler Large Programme using the Ultraviolet and Visual Echelle Spectrograph at the Very Large Telescope at Paranal, Chile. It is designed to provide an observationally unbiased view of the presence and behaviour of the DIBs towards early-spectral-type stars whose lines-of-sight probe the diffuse-to-translucent ISM. Such a complete dataset will provide a deep census of the atomic and molecular content, physical conditions, chemical abundances and elemental depletion levels for each sightline. Achieving these goals requires a homogeneous set of high-quality data in terms of resolution (R ~ 70000 -- 100000), sensitivity (S/N up to 1000 per resolution element), and spectral coverage (305--1042 nm), as well as a large sample size (100+ sightlines). In this first paper the goals, objectives and methodology of the EDIBLES programme are described and an initial assessment of the data is provided. △ Less

Submitted 4 August, 2017; originally announced August 2017.

Comments: 29 pages (including 14 pages appendix), 15 figures, accepted by A&A

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

[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)

Herschel survey and modelling of externally-illuminated photoevaporating protoplanetary disks

Authors: Jason Champion, Olivier Berné, Silvia Vicente, Inga Kamp, Franck Le Petit, Antoine Gusdorf, Christine Joblin, Javier R. Goicoechea

Abstract: Protoplanetary disks undergo substantial mass-loss by photoevaporation, a mechanism which is crucial to their dynamical evolution. However, the processes regulating the gas energetics have not been well constrained by observations so far. We aim at studying the processes involved in disk photoevaporation when it is driven by far-UV photons. We present a unique Herschel survey and new ALMA observat… ▽ More Protoplanetary disks undergo substantial mass-loss by photoevaporation, a mechanism which is crucial to their dynamical evolution. However, the processes regulating the gas energetics have not been well constrained by observations so far. We aim at studying the processes involved in disk photoevaporation when it is driven by far-UV photons. We present a unique Herschel survey and new ALMA observations of four externally-illuminated photoevaporating disks (a.k.a. proplyds). For the analysis of these data, we developed a 1D model of the photodissociation region (PDR) of a proplyd, based on the Meudon PDR code and computed the far infrared line emission. We successfully reproduce most of the observations and derive key physical parameters, i.e. densities at the disk surface of about $10^{6}$ cm$^{-3}$ and local gas temperatures of about 1000 K. Our modelling suggests that all studied disks are found in a transitional regime resulting from the interplay between several heating and cooling processes that we identify. These differ from those dominating in classical PDRs, i.e. grain photo-electric effect and cooling by [OI] and [CII] FIR lines. This energetic regime is associated to an equilibrium dynamical point of the photoevaporation flow: the mass-loss rate is self-regulated to set the envelope column density at a value that maintains the temperature at the disk surface around 1000 K. From our best-fit models, we estimate mass-loss rates - of the order of $10^{-7}$ $\mathrm{M}_\odot$/yr - that are in agreement with earlier spectroscopic observation of ionised gas tracers. This holds only if we assume an evaporation flow launched from the disk surface at sound speed (supercritical regime). We have identified the energetic regime regulating FUV-photoevaporation in proplyds. This regime could be implemented into models of the dynamical evolution of protoplanetary disks. △ Less

Submitted 1 February, 2017; originally announced February 2017.

Comments: 17 pages (6 more with appendices), 11 figures, accepted for publication in A&A

Journal ref: A&A 604, A69 (2017)

Dissecting the molecular structure of the Orion B cloud: Insight from Principal Component Analysis

Authors: Pierre Gratier, Emeric Bron, Maryvonne Gerin, Jérôme Pety, Viviana V. Guzman, Jan Orkisz, Sébastien Bardeau, Javier R. Goicoechea, Franck Le Petit, Harvey Liszt, Karin Öberg, Nicolas Peretto, Evelyne Roueff, Albrech Sievers, Pascal Tremblin

Abstract: Context. The combination of wideband receivers and spectrometers currently available in (sub-)millimeter observatories deliver wide- field hyperspectral imaging of the interstellar medium. Tens of spectral lines can be observed over degree wide fields in about fifty hours. This wealth of data calls for restating the physical questions about the interstellar medium in statistical terms. Aims. We ai… ▽ More Context. The combination of wideband receivers and spectrometers currently available in (sub-)millimeter observatories deliver wide- field hyperspectral imaging of the interstellar medium. Tens of spectral lines can be observed over degree wide fields in about fifty hours. This wealth of data calls for restating the physical questions about the interstellar medium in statistical terms. Aims. We aim at gaining information on the physical structure of the interstellar medium from a statistical analysis of many lines from different species over a large field of view, without requiring detailed radiative transfer or astrochemical modeling. Methods. We coupled a nonlinear rescaling of the data with one of the simplest multivariate analysis methods, namely the Principal Component Analysis, to decompose the observed signal into components that we interpret first qualitatively and then quantitatively based on our deep knowledge of the observed region and of the astrochemistry at play. Results. We identify 3 principal components, linear compositions of line brightness temperatures, that are correlated at various levels with the column density, the volume density and the UV radiation field. Conclusions. When sampling a sufficiently diverse mixture of physical parameters, it is possible to decompose the molecular emission in order to gain physical insight on the observed interstellar medium. This opens a new avenue for future studies of the interstellar medium. △ Less

Submitted 16 January, 2017; originally announced January 2017.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 599, A100 (2017)

Turbulence and star formation efficiency in molecular clouds: solenoidal versus compressive motions in Orion B

Authors: Jan H. Orkisz, Jérôme Pety, Maryvonne Gerin, Emeric Bron, Viviana V. Guzmán, Sébastien Bardeau, Javier R. Goicoechea, Pierre Gratier, Franck Le Petit, François Levrier, Harvey Liszt, Karin Öberg, Nicolas Peretto, Evelyne Roueff, Albrecht Sievers, Pascal Tremblin

Abstract: The nature of turbulence in molecular clouds is one of the key parameters that control star formation efficiency: compressive motions, as opposed to solenoidal motions, can trigger the collapse of cores, or mark the expansion of Hii regions. We try to observationally derive the fractions of momentum density ($ρv$) contained in the solenoidal and compressive modes of turbulence in the Orion B molec… ▽ More The nature of turbulence in molecular clouds is one of the key parameters that control star formation efficiency: compressive motions, as opposed to solenoidal motions, can trigger the collapse of cores, or mark the expansion of Hii regions. We try to observationally derive the fractions of momentum density ($ρv$) contained in the solenoidal and compressive modes of turbulence in the Orion B molecular cloud and relate these fractions to the star formation efficiency in the cloud. The implementation of a statistical method developed by Brunt & Federrath (2014), applied to a $^{13}$CO(J=1-0) datacube obtained with the IRAM-30m telescope, allows us to retrieve 3-dimensional quantities from the projected quantities provided by the observations, yielding an estimate of the compressive versus solenoidal ratio in various regions of the cloud. Despite the Orion B molecular cloud being highly supersonic (mean Mach number $\sim$ 6), the fractions of motion in each mode diverge significantly from equipartition. The cloud's motions are on average mostly solenoidal (excess > 8 % with respect to equipartition), which is consistent with its low star formation rate. On the other hand, the motions around the main star-forming regions (NGC 2023 and NGC 2024) prove to be strongly compressive. We have successfully applied to observational data a method that was so far only tested on simulations, and have shown that there can be a strong intra-cloud variability of the compressive and solenoidal fractions, these fractions being in turn related to the star formation efficiency. This opens a new possibility for star-formation diagnostics in galactic molecular clouds. △ Less

Submitted 4 January, 2017; originally announced January 2017.

Comments: accepted by A&A

Journal ref: A&A 599, A99 (2017)

The anatomy of the Orion B Giant Molecular Cloud: A local template for studies of nearby galaxies

Authors: Jérôme Pety, Viviana V. Guzmán, Jan H. Orkisz, Harvey S. Liszt, Maryvonne Gerin, Emeric Bron, Sébastien Bardeau, Javier R. Goicoechea, Pierre Gratier, Franck Le Petit, François Levrier, Karin I. Oberg, Evelyne Roueff, Albrecht Sievers

Abstract: We aim to develop the Orion B Giant Molecular Cloud (GMC) as a local template for interpreting extra-galactic molecular line observations. We use the wide-band receiver at the IRAM-30m to spatially and spectrally resolve the Orion B GMC. The observations cover almost 1 square degree at 26" resolution with a bandwidth of 32 GHz from 84 to 116 GHz in only two tunings. Among the mapped spectral lines… ▽ More We aim to develop the Orion B Giant Molecular Cloud (GMC) as a local template for interpreting extra-galactic molecular line observations. We use the wide-band receiver at the IRAM-30m to spatially and spectrally resolve the Orion B GMC. The observations cover almost 1 square degree at 26" resolution with a bandwidth of 32 GHz from 84 to 116 GHz in only two tunings. Among the mapped spectral lines are the 12CO, 13CO, C18O, C17O, HCN, HNC, 12CN, CCH, HCO+, N2H+ (1-0), and 12CS, 32SO, SiO, c-C3H2, CH3OH (2-1) transitions. We introduce the molecular anatomy of the Orion B GMC, including relations between line intensities and gas column density or far-UV radiation fields, and correlations between selected line and line ratios. We also obtain a dust-traced gas mass that is less than about one third the CO-traced mass, using the standard Xco conversion factor. The presence of overluminous CO can be traced back to the dependence of the CO intensity on UV illumination. In fact, while most lines show some dependence on the UV radiation field, CN and CCH are the most sensitive. Moreover dense cloud cores are almost exclusively traced by N2H+. Other traditional high density tracers, such as HCN (1-0), are also easily detected in extended translucent regions at a typical density of about 500 H2 cm-3. In general, we find no straightforward relation between line critical density and the fraction of the line luminosity coming from dense gas regions. Our initial findings demonstrate that the relations between line (ratio) intensities and environment in GMCs are more complicated than often assumed. Sensitivity (i.e., the molecular column density), excitation, and above all chemistry contribute to the observed line intensity distributions. They must be considered together when developing the next generation of extra-galactic molecular line diagnostics of mass, density, temperature and radiation field. △ Less

Submitted 12 November, 2016; originally announced November 2016.

Comments: 32 pages, 29 PDF figures. Accepted for publication in A&A. Uses aa and natbib latex macros

Journal ref: A&A 599, A98 (2017)

Radiative and mechanical feedback into the molecular gas in the Large Magellanic Cloud. I. N159W

Authors: Min-Young Lee, Suzanne Madden, Vianney Lebouteiller, Antoine Gusdorf, Benjamin Godard, Ronin Wu, Maud Galametz, Diane Cormier, Franck Le Petit, Evelyne Roueff, Emeric Bron, Lynn Carlson, Melanie Chevance, Yasuo Fukui, Frederic Galliano, Sacha Hony, Annie Hughes, Remy Indebetouw, Franck Israel, Akiko Kawamura, Jacques Le Bourlot, Pierre Lesaffre, Margaret Meixner, Erik Muller, Omnarayani Nayak , et al. (3 additional authors not shown)

Abstract: We present Herschel SPIRE Fourier Transform Spectrometer (FTS) observations of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In our observations, a number of far-infrared cooling lines including CO(4-3) to CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected. With an aim of investigating the physical conditions and excitation processes of molecu… ▽ More We present Herschel SPIRE Fourier Transform Spectrometer (FTS) observations of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In our observations, a number of far-infrared cooling lines including CO(4-3) to CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected. With an aim of investigating the physical conditions and excitation processes of molecular gas, we first construct CO spectral line energy distributions (SLEDs) on 10 pc scales by combining the FTS CO transitions with ground-based low-J CO data and analyze the observed CO SLEDs using non-LTE radiative transfer models. We find that the CO-traced molecular gas in N159W is warm (kinetic temperature of 153-754 K) and moderately dense (H2 number density of (1.1-4.5)e3 cm-3). To assess the impact of the energetic processes in the interstellar medium on the physical conditions of the CO-emitting gas, we then compare the observed CO line intensities with the models of photodissociation regions (PDRs) and shocks. We first constrain the properties of PDRs by modelling Herschel observations of [OI] 145, [CII] 158, and [CI] 370 micron fine-structure lines and find that the constrained PDR components emit very weak CO emission. X-rays and cosmic-rays are also found to provide a negligible contribution to the CO emission, essentially ruling out ionizing sources (ultraviolet photons, X-rays, and cosmic-rays) as the dominant heating source for CO in N159W. On the other hand, mechanical heating by low-velocity C-type shocks with ~10 km/s appears sufficient enough to reproduce the observed warm CO. △ Less

Submitted 14 June, 2016; originally announced June 2016.

Comments: accepted for publication in A&A

Journal ref: A&A 596, A85 (2016)

A milestone toward understanding PDR properties in the extreme environment of LMC-30Dor

Authors: M. Chevance, S. C. Madden, V. Lebouteiller, B. Godard, D. Cormier, F. Galliano, S. Hony, R. Indebetouw, J. Le Bourlot, M. Y. Lee, F. Le Petit, E. Pellegrini, E. Roueff, R. Wu

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

Submitted 11 March, 2016; originally announced March 2016.

Comments: 20 pages, 23 figures, accepted in A&A

Efficient ortho-para conversion of H2 on interstellar grain surfaces

Authors: Emeric Bron, Franck Le Petit, Jacques Le Bourlot

Abstract: Context: Fast surface conversion between ortho- and para-H2 has been observed in laboratory studies, and this mechanism has been proposed to play a role in the control of the ortho-para ratio in the interstellar medium. Observations of rotational lines of H2 in Photo-Dissociation Regions (PDRs) have indeed found significantly lower ortho-para ratios than expected at equilibrium. The mechanisms con… ▽ More Context: Fast surface conversion between ortho- and para-H2 has been observed in laboratory studies, and this mechanism has been proposed to play a role in the control of the ortho-para ratio in the interstellar medium. Observations of rotational lines of H2 in Photo-Dissociation Regions (PDRs) have indeed found significantly lower ortho-para ratios than expected at equilibrium. The mechanisms controlling the balance of the ortho-para ratio in the interstellar medium thus remain incompletely understood, while this ratio can affect the thermodynamical properties of the gas (equation of state, cooling function). Aims: We aim to build an accurate model of ortho-para conversion on dust surfaces based on the most recent experimental and theoretical results, and to validate it by comparison to observations of H2 rotational lines in PDRs. Methods: We propose a statistical model of ortho-para conversion on dust grains with fluctuating dust temperatures, based on a master equation approach. This computation is then coupled to full PDR models and compared to PDR observations. Results: We show that the observations of rotational H2 lines indicate a high conversion efficiency on dust grains, and that this high efficiency can be accounted for if taking dust temperature fluctuations into account with our statistical model of surface conversion. Simpler models neglecting the dust temperature fluctuations do not reach the high efficiency deduced from the observations. Moreover, this high efficiency induced by dust temperature fluctuations is quite insensitive to the values of microphysical parameters of the model. Conclusions: Ortho-para conversion on grains is thus an efficient mechanism in most astrophysical conditions that can play a significant role in controlling the ortho-para ratio. △ Less

Submitted 17 January, 2016; originally announced January 2016.

Comments: Accepted in Astronomy & Astrophysics

Journal ref: A&A 588, A27 (2016)

HI-to-H2 transitions in the Perseus Molecular Cloud

Authors: Shmuel Bialy, Amiel Sternberg, Min-Young Lee, Franck Le Petit, Evelyne Roueff

Abstract: We apply the Sternberg et al. (2014) theoretical model to analyze HI and H2 observations in the Perseus molecular cloud. We constrain the physical properties of the HI shielding envelopes and the nature of the HI-to-H2 transitions. Our analysis (Bialy et al. 2015) implies that in addition to cold neutral gas (CNM), less dense thermally-unstable gas (UNM) significantly contributes to the shielding… ▽ More We apply the Sternberg et al. (2014) theoretical model to analyze HI and H2 observations in the Perseus molecular cloud. We constrain the physical properties of the HI shielding envelopes and the nature of the HI-to-H2 transitions. Our analysis (Bialy et al. 2015) implies that in addition to cold neutral gas (CNM), less dense thermally-unstable gas (UNM) significantly contributes to the shielding of the H2 cores in Perseus. △ Less

Submitted 9 November, 2015; originally announced November 2015.

Comments: 4 pages, 1 figure. To appear in proceedings of "The 6th Zermatt ISM Symposium: Conditions and Impact of Star Formation From Lab to Space", eds. R. Simon and M. Rollig

Physical conditions in the central molecular zone inferred by H3+

Authors: Franck Le Petit, Maxime Ruaud, Emeric Bron, Benjamin Godard, Evelyne Roueff, David Languignon, Jacques Le Bourlot

Abstract: The H3+ molecule has been detected in many lines of sight within the central molecular zone (CMZ) with exceptionally large column densities and unusual excitation properties compared to diffuse local clouds. The detection of the (3,3) metastable level has been suggested to be the signature of warm and diffuse gas in the CMZ. We use the Meudon PDR code to re-examine the relationship between the col… ▽ More The H3+ molecule has been detected in many lines of sight within the central molecular zone (CMZ) with exceptionally large column densities and unusual excitation properties compared to diffuse local clouds. The detection of the (3,3) metastable level has been suggested to be the signature of warm and diffuse gas in the CMZ. We use the Meudon PDR code to re-examine the relationship between the column density of H3+ and the cosmic-ray ionization rate, $ζ$, up to large values of $ζ$. We study the impact of the various mechanisms that can excite H3+ in its metastable state. We produce grids of PDR models exploring different parameters ($ζ$, size of clouds, metallicity) and infer the physical conditions that best match the observations toward ten lines of sight in the CMZ. For one of them, Herschel observations of HF, OH+, H2O+, and H3O+ can be used as additional constraints. We check that the results found for H3+ also account for the observations of these molecules. We find that the linear relationship between N(H3+) and $ζ$ only holds up to a certain value of the cosmic-ray ionization rate, which depends on the proton density. A value $ζ\sim 1 - 11 \times 10^{-14}$ s$^{-1}$ explains both the large observed H3+ column density and its excitation in the metastable level (3,3) in the CMZ. It also reproduces N(OH+), N(H2O+) and N(H3O+) detected toward Sgr B2(N). We confirm that the CMZ probed by H3+ is diffuse, nH $\lesssim$ 100 cm-3 and warm, T $\sim$ 212-505 K. This warm medium is due to cosmic-ray heating. We also find that the diffuse component probed by H3+ must fill a large fraction of the CMZ. Finally, we suggest the warm gas in the CMZ enables efficient H2 formation via chemisorption sites as in PDRs. This contributes to enhance the abundance of H3+ in this high cosmic-ray flux environment. △ Less

Submitted 8 October, 2015; originally announced October 2015.

Comments: accepted in A&A

Journal ref: A&A 585, A105 (2016)

IVOA recommendation: Parameter Description Language Version 1.0

Authors: Carlo Maria Zwolf, Paul Harrison, Julian Garrido, Jose Enrique Ruiz, Franck Le Petit

Abstract: This document discusses the definition of the Parameter Description Language (PDL). In this language parameters are described in a rigorous data model. With no loss of generality, we will represent this data model using XML. It intends to be a expressive language for self-descriptive web services exposing the semantic nature of input and output parameters, as well as all necessary complex constrai… ▽ More This document discusses the definition of the Parameter Description Language (PDL). In this language parameters are described in a rigorous data model. With no loss of generality, we will represent this data model using XML. It intends to be a expressive language for self-descriptive web services exposing the semantic nature of input and output parameters, as well as all necessary complex constraints. PDL is a step forward towards true web services interoperability. △ Less

Submitted 28 September, 2015; originally announced September 2015.

HI-to-H2 Transitions in the Perseus Molecular Cloud

Authors: Shmuel Bialy, Amiel Sternberg, Min-Young Lee, Franck Le Petit, Evelyne Roueff

Abstract: We use the Sternberg et al. (2014) theory for interstellar atomic to molecular (HI-to-H$_2$) conversion to analyze HI-to-H$_2$ transitions in five (low-mass) star-forming and dark regions in the Perseus molecular cloud, B1, B1E, B5, IC348, and NGC1333. The observed HI mass surface densities of 6.3 to 9.2 M$_{\odot}$ pc$^{-2}$ are consistent with HI-to-H$_2$ transitions dominated by HI-dust shieldi… ▽ More We use the Sternberg et al. (2014) theory for interstellar atomic to molecular (HI-to-H$_2$) conversion to analyze HI-to-H$_2$ transitions in five (low-mass) star-forming and dark regions in the Perseus molecular cloud, B1, B1E, B5, IC348, and NGC1333. The observed HI mass surface densities of 6.3 to 9.2 M$_{\odot}$ pc$^{-2}$ are consistent with HI-to-H$_2$ transitions dominated by HI-dust shielding in predominantly atomic envelopes. For each source, we constrain the dimensionless parameter $αG$, and the ratio $I_{\rm UV}/n$, of the FUV intensity to hydrogen gas density. We find $αG$ values from 5.0 to 26.1, implying characteristic atomic hydrogen densities 11.8 to 1.8 cm$^{-3}$, for $I_{\rm UV} \approx 1$ appropriate for Perseus. Our analysis implies that the dusty HI shielding layers are probably multiphased, with thermally unstable UNM gas in addition to cold CNM within the 21 cm kinematic radius. △ Less

Submitted 19 July, 2015; v1 submitted 22 May, 2015; originally announced May 2015.

Comments: 5 pages, 2 Figures. Minor improvements suggested by the referee. Accepted for publication in the Astrophysical Journal

Journal ref: ApJ, 809, 122 (2015)

$Herschel$ SPIRE-FTS observations of RCW 120

Authors: J. A. Rodón, A. Zavagno, J. -P. Baluteau, E. Habart, M. Köhler, J. Le Bourlot, F. Le Petit, A. Abergel

Abstract: The expansion of Galactic HII regions can trigger the formation of a new generation of stars. However, little is know about the physical conditions that prevail in these regions. We study the physical conditions that prevail in specific zones towards expanding HII regions that trace representative media such as the photodissociation region, the ionized region, and condensations with and without on… ▽ More The expansion of Galactic HII regions can trigger the formation of a new generation of stars. However, little is know about the physical conditions that prevail in these regions. We study the physical conditions that prevail in specific zones towards expanding HII regions that trace representative media such as the photodissociation region, the ionized region, and condensations with and without ongoing star formation. We use the SPIRE Fourier Transform Spectrometer (FTS) on board $Herschel$ to observe the HII region RCW 120. Continuum and lines are observed in the $190-670\,μ$m range. Line intensities and line ratios are obtained and used as physical diagnostics of the gas. We used the Meudon PDR code and the RADEX code to derive the gas density and the radiation field at nine distinct positions including the PDR surface and regions with and without star-formation activity. For the different regions we detect the atomic lines [NII] at $205\,μ$m and [CI] at $370$ and $609\,μ$m, the $^{12}{\rm CO}$ ladder between the $J=4$ and $J=13$ levels and the $^{13}{\rm CO}$ ladder between the $J=5$ and $J=14$ levels, as well as CH$ ^{+} $ in absorption. We find gas temperatures in the range $45-250\,$K for densities of $10^4-10^6\,{\rm cm}^{-3}$, and a high column density on the order of $N_{\rm H}\sim10^{22}\,{\rm cm}^{-2}$ that is in agreement with dust analysis. The ubiquitousness of the atomic and CH$ ^{+} $ emission suggests the presence of a low-density PDR throughout RCW 120. High-excitation lines of CO indicate the presence of irradiated dense structures or small dense clumps containing young stellar objects, while we also find a less dense medium ($N_{\rm H}\sim10^{20}\,{\rm cm}^{-2}$) with high temperatures ($80-200\,$K). △ Less

Submitted 24 April, 2015; originally announced April 2015.

Comments: 11 pages, 11 figures, accepted by A&A

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

Surface chemistry in the Interstellar Medium II. $\mathrm{H}_2$ formation on dust with random temperature fluctuations

Authors: Emeric Bron, Jacques Le Bourlot, Franck Le Petit

Abstract: The $\mathrm{H}_2$ formation on grains is known to be sensitive to dust temperature, which is also known to fluctuate for small grain sizes due to photon absorption. We aim at exploring the consequences of simultaneous fluctuations of the dust temperature and the adsorbed H-atom population on the $\mathrm{H}_2$ formation rate under the full range of astrophysically relevant UV intensities and gas… ▽ More The $\mathrm{H}_2$ formation on grains is known to be sensitive to dust temperature, which is also known to fluctuate for small grain sizes due to photon absorption. We aim at exploring the consequences of simultaneous fluctuations of the dust temperature and the adsorbed H-atom population on the $\mathrm{H}_2$ formation rate under the full range of astrophysically relevant UV intensities and gas conditions. The master equation approach is generalized to coupled fluctuations in both the grain's temperature and its surface population and solved numerically. The resolution can be simplified in the case of the Eley-Rideal mechanism, allowing a fast computation. For the Langmuir-Hinshelwood mechanism, it remains computationally expensive, and accurate approximations are constructed. We find the Langmuir-Hinshelwood mechanism to become an efficient formation mechanism in unshielded photon dominated region (PDR) edge conditions when taking those fluctuations into account, despite hot average dust temperatures. It reaches an importance comparable to the Eley-Rideal mechanism. However, we show that a simpler rate equation treatment gives qualitatively correct observable results in full cloud simulations under most astrophysically relevant conditions. Typical differences are a factor of 2-3 on the intensities of the $\mathrm{H}_2$ $v=0$ lines. We also find that rare fluctuations in cloud cores are sufficient to significantly reduce the formation efficiency. Our detailed analysis confirms that the usual approximations used in numerical models are adequate when interpreting observations, but a more sophisticated statistical analysis is required if one is interested in the details of surface processes. △ Less

Submitted 16 July, 2014; originally announced July 2014.

Comments: 21 pages, 28 figures, accepted in A&A

Journal ref: A&A 569, A100 (2014)

Chemical complexity in the Horsehead photodissociation region

Authors: Viviana V. Guzmán, Jérôme Pety, Pierre Gratier, Javier R. Goicoechea, Maryvonne Gerin, Evelyne Roueff, Franck Le Petit, Jacques Le Bourlot

Abstract: The interstellar medium is known to be chemically complex. Organic molecules with up to 11 atoms have been detected in the interstellar medium, and are believed to be formed on the ices around dust grains. The ices can be released into the gas-phase either through thermal desorption, when a newly formed star heats the medium around it and completely evaporates the ices; or through non-thermal deso… ▽ More The interstellar medium is known to be chemically complex. Organic molecules with up to 11 atoms have been detected in the interstellar medium, and are believed to be formed on the ices around dust grains. The ices can be released into the gas-phase either through thermal desorption, when a newly formed star heats the medium around it and completely evaporates the ices; or through non-thermal desorption mechanisms, such as photodesorption, when a single far-UV photon releases only a few molecules from the ices. The first one dominates in hot cores, hot corinos and strongly UV-illuminated PDRs, while the second one dominates in colder regions, such as low UV-field PDRs. This is the case of the Horsehead were dust temperatures are ~20-30K, and therefore offers a clean environment to investigate what is the role of photodesorption. We have carried-out an unbiased spectral line survey at 3, 2 and 1mm with the IRAM-30m telescope in the Horsehead nebula, with an unprecedented combination of bandwidth high spectral resolution and sensitivity. Two positions were observed: the warm PDR and a cold condensation shielded from the UV field (dense core), located just behind the PDR edge. We summarize our recently published results from this survey and present the first detection of the complex organic molecules HCOOH, CH2CO, CH3CHO and CH3CCH in a PDR. These species together with CH3CN present enhanced abundances in the PDR compared to the dense core. This suggests that photodesorption is an efficient mechanism to release complex molecules into the gas-phase in far-UV illuminated regions. △ Less

Submitted 30 April, 2014; originally announced April 2014.

Comments: 15 pages, 7 figures, 7 tables, Accepted in Faraday discussions 168

HI-to-H2 Transitions and H I Column Densities in Galaxy Star-Forming Regions

Authors: Amiel Sternberg, Franck Le Petit, Evelyne Roueff, Jacques Le Bourlot

Abstract: We present new analytic theory and radiative transfer computations for the atomic to molecular (HI-to-H2) transitions, and the build-up of atomic-hydrogen (HI) gas columns, in optically thick interstellar clouds, irradiated by far-ultraviolet photodissociating radiation fields. We derive analytic expressions for the total HI column densities for (1D) planar slabs, for beamed or isotropic radiation… ▽ More We present new analytic theory and radiative transfer computations for the atomic to molecular (HI-to-H2) transitions, and the build-up of atomic-hydrogen (HI) gas columns, in optically thick interstellar clouds, irradiated by far-ultraviolet photodissociating radiation fields. We derive analytic expressions for the total HI column densities for (1D) planar slabs, for beamed or isotropic radiation fields, from the weak- to strong-field limits, for gradual or sharp atomic to molecular transitions, and for arbitrary metallicity. Our expressions may be used to evaluate the HI column densities as functions of the radiation field intensity and the H2-dust-limited dissociation flux, the hydrogen gas density, and the metallicity-dependent H2 formation rate-coefficient and far-UV dust-grain absorption cross-section. We make the distinction between "HI-dust" and "H2-dust" opacity, and we present computations for the "universal H2-dust-limited effective dissociation bandwidth". We validate our analytic formulae with Meudon PDR code computations for the HI-to-H2 density profiles, and total HI column densities. We show that our general 1D formulae predict HI columns and H2 mass fractions that are essentially identical to those found in more complicated (and approximate) spherical (shell/core) models. We apply our theory to compute H2 mass fractions and star-formation thresholds for individual clouds in self-regulated galaxy disks, for a wide range of metallicities. Our formulae for the HI columns and H2 mass fractions may be incorporated into hydrodynamics simulations for galaxy evolution. △ Less

Submitted 27 February, 2015; v1 submitted 20 April, 2014; originally announced April 2014.

Comments: 33 pages, 15 figures, published in ApJS