Source author record

Jacques Le Bourlot

Jacques Le Bourlot appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2025arXiv

PDRs4All XVII: Formation and excitation of HD in photodissociation regions. Application to the Orion Bar

The James Webb Space Telescope enabled the first detection of several rovibrational emission lines of HD in the Orion Bar, a prototypical photodissociation region. This provides an incentive to examine the physics of HD in dense and strong PDRs. Using the latest data available on HD excitation by collisional, radiative and chemical processes, our goal is to unveil HD formation and excitation processes in PDRs by comparing our state-of-the-art PDR model with observations made in the Orion Bar and discuss if and how HD can be used as a complementary tracer of physical parameters in the emitting region. We compute detailed PDR models, using an upgraded version of the Meudon PDR code, which are compared to NIRSpec data using excitation diagrams and synthetic emission spectra. The models predict that HD is mainly produced in the gas phase via the reaction D + H2 = H + HD at the front edge of the PDR and that the D/HD transition is located slightly closer to the edge than the H/H2 transition. Rovibrational levels are excited by UV pumping. In the observations, HD rovibrational emission is detected close to the H/H2 dissociation fronts of the Orion Bar and peaks where vibrationally excited H2 peaks, rather than at the maximum emission of pure rotational H2 levels. We derive an excitation temperature around Tex ~ 480 - 710 K. Due to high continuum in the Orion Bar, fringes lead to high noise levels beyond 15 $μ$m, no pure rotational lines of HD are detected. The comparison to PDR models shows that a range of thermal pressure P = (3-9)x10$^7$ K cm$^{-3}$ with no strong constraints on the intensity of the UV field are compatible with HD observations. This range of pressure is compatible with previous estimates from H2 observations with JWST. This is the first time that observations of HD emission lines in the near-infrared are used to put constraints on the thermal pressure in the PDR.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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]

preprint2020arXiv

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

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)

preprint2016arXiv

Efficient ortho-para conversion of H2 on interstellar grain surfaces

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.

preprint2016arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

Physical conditions in the central molecular zone inferred by H3+

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.

preprint2014arXiv

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

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.

preprint2012arXiv

Herschel observations of interstellar chloronium

Using the Herschel Space Observatory's Heterodyne Instrument for the Far-Infrared (HIFI), we have observed para-chloronium (H2Cl+) toward six sources in the Galaxy. We detected interstellar chloronium absorption in foreground molecular clouds along the sight-lines to the bright submillimeter continuum sources Sgr A (+50 km/s cloud) and W31C. Both the para-H2-35Cl+ and para-H2-37Cl+ isotopologues were detected, through observations of their 1(11)-0(00) transitions at rest frequencies of 485.42 and 484.23 GHz, respectively. For an assumed ortho-to-para ratio of 3, the observed optical depths imply that chloronium accounts for ~ 4 - 12% of chlorine nuclei in the gas phase. We detected interstellar chloronium emission from two sources in the Orion Molecular Cloud 1: the Orion Bar photodissociation region and the Orion South condensation. For an assumed ortho-to-para ratio of 3 for chloronium, the observed emission line fluxes imply total beam-averaged column densities of ~ 2.0E+13 cm-2 and ~ 1.2E+13 cm-2, respectively, for chloronium in these two sources. We obtained upper limits on the para-H2-35Cl+ line strengths toward H2 Peak 1 in the Orion Molecular cloud and toward the massive young star AFGL 2591. The chloronium abundances inferred in this study are typically at least a factor ~10 larger than the predictions of steady-state theoretical models for the chemistry of interstellar molecules containing chlorine. Several explanations for this discrepancy were investigated, but none has proven satisfactory, and thus the large observed abundances of chloronium remain puzzling.

preprint2012arXiv

Surface chemistry in the interstellar medium - I - H2 formation by Langmuir-Hinshelwood and Eley-Rideal mechanisms

H2 formation remains a major issue for the understanding of interstellar physics. We investigate H2 formation in the interstellar medium at the light of the most recent experimental and theoretical data. We implemented detailed H2 formation mechanisms on grains surface in the Meudon PDR code : i) Langmuir-Hinshelwood mechanism taking into account the contribution of the different sizes of dust grains in the diffusion processes and ii) the Eley-Rideal mechanism. We show that, thanks to these processes, H2 can be formed even in regions where dust temperature is larger than 25 K. We also show that formation by Eley-Rideal mechanism can be a significant source of heating of the gas. We derive line intensities for various astrophysical conditions. Such an approach results in an enhanced H2 formation rate compared to the standard formation determined by observations in absorption in the UV. We derive some H2 line intensities for isobaric and isochoric models.

preprint2011arXiv

Herschel Measurements of Molecular Oxygen in Orion

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

preprint2010arXiv

Excitation of H$_2$ in photodissociation regions as seen by Spitzer

We present spectroscopic observations obtained with the infrared Spitzer Space Telescope, which provide insight into the H$_2$ physics and gas energetics in photodissociation Regions (PDRs) of low to moderate far-ultraviolet (FUV) fields and densities. We analyze data on six well known Galactic PDRs (L1721, California, N7023E, Horsehead, rho Oph, N2023N), sampling a poorly explored range of excitation conditions ($χ\sim 5-10^3$), relevant to the bulk of molecular clouds in galaxies. Spitzer observations of H$_2$ rotational lines are complemented with H$_2$ data, including ro-vibrational line measurements, obtained with ground-based telescopes and ISO, to constrain the relative contributions of ultraviolet pumping and collisions to the H$_2$ excitation. The data analysis is supported by model calculations with the Meudon PDR code. The observed column densities of rotationally excited H$_2$ are observed to be much higher than PDR model predictions. In the lowest excitation PDRs, the discrepancy between the model and the data is about one order of magnitude for rotational levels $J \ge $3. We discuss whether an enhancement in the H$_2$ formation rate or a local increase in photoelectric heating, as proposed for brighter PDRs in former ISO studies, may improve the data-model comparison. We find that an enhancement in the H$_2$ formation rates reduces the discrepancy, but the models still fall short of the data. This large disagreement suggests that our understanding of the formation and excitation of H$_2$ and/or of PDRs energetics is still incomplete. We discuss several explanations, which could be further tested using the Herschel Space Telescope

preprint2006arXiv

A model for atomic and molecular interstellar gas: The Meudon PDR code

We present the revised ``Meudon'' model of Photon Dominated Region (PDR code), presently available on the web under the Gnu Public Licence at: http://aristote.obspm.fr/MIS. General organisation of the code is described down to a level that should allow most observers to use it as an interpretation tool with minimal help from our part. Two grids of models, one for low excitation diffuse clouds and one for dense highly illuminated clouds, are discussed, and some new results on PDR modelisation highlighted.