Source author record

Pierre Gratier

Pierre Gratier 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

8works
6topics
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

8 published item(s)

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

Revised mass-radius relationships for water-rich rocky planets more irradiated than the runaway greenhouse limit

Mass-radius relationships for water-rich rocky planets are usually calculated assuming most water is present in condensed (either liquid or solid) form. Planet density estimates are then compared to these mass-radius relationships, even when these planets are more irradiated than the runaway greenhouse irradiation limit (around 1.1~times the insolation at Earth for planets orbiting a Sun-like star), for which water has been shown to be unstable in condensed form and would instead form a thick H2O-dominated atmosphere. Here we use the LMD Generic numerical climate model to derive new mass-radius relationships appropriate for water-rich rocky planets that are more irradiated than the runaway greenhouse irradiation limit, meaning planets endowed with a steam, water-dominated atmosphere. For a given water-to-rock mass ratio, these new mass-radius relationships lead to planet bulk densities much lower than calculated when water is assumed to be in condensed form. In other words, using traditional mass-radius relationships for planets that are more irradiated than the runaway greenhouse irradiation limit tends to dramatically overestimate -- possibly by several orders of magnitude -- their bulk water content. In particular, this result applies to TRAPPIST-1 b, c, and d, which can accommodate a water mass fraction of at most 2, 0.3 and 0.08 %, respectively, assuming planetary core with a terrestrial composition. In addition, we show that significant changes of mass-radius relationships (between planets less and more irradiated than the runaway greenhouse limit) can be used to remove bulk composition degeneracies in multiplanetary systems such as TRAPPIST-1. Finally, we provide an empirical formula for the H2O steam atmosphere thickness which can be used to construct mass-radius relationships for any water-rich, rocky planet more irradiated than the runaway greenhouse irradiation threshold.

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

Detection of CH3SH in protostar IRAS 16293-2422

The nature of the main sulphur reservoir in star forming regions is a long standing mystery. The observed abundance of sulphur-bearing species in dense clouds is only about 0.1 per cent of the same quantity in diffuse clouds. Therefore, the main sulphur species in star forming regions of the interstellar medium are still unknown. IRAS 16293-2422 is one of the regions where production of S-bearing species is favourable due to its conditions which allows the evaporation of ice mantles. We carried out observations in the 3 mm band towards the solar type protostar IRAS 16293-2422 with the IRAM 30m telescope. We observed a single frequency setup with the EMIR heterodyne 3 mm receiver with an Lower Inner (LI) tuning frequency of 89.98 GHz. Several lines of the complex sulphur species CH3SH were detected. Observed abundances are compared with simulations using the NAUTILUS gas-grain chemical model. Modelling results suggest that CH3SH has the constant abundance of 4e-9 (compared to H2) for radii lower than 200 AU and is mostly formed on the surfaces. Detection of CH3SH indicates that there may be several new families of S-bearing molecules (which could form starting from CH3SH) which have not been detected or looked for yet.

preprint2013arXiv

An observational investigation of the identity of B11244 ($l$-C$_3$H$^+$/C$_3$H$^-$)

Pety et al. (2012) have reported the detection of eight transitions of a closed-shell, linear molecule (B11244) in observations toward the Horsehead PDR, which they attribute to the $l$-C$_3$H$^+$ cation. Recent high-level \textit{ab initio} calculations have called this assignment into question; the anionic C$_3$H$^-$ molecule has been suggested as a more likely candidate. Here, we examine observations of the Horsehead PDR, Sgr B2(N), TMC-1, and IRC+10216 in the context of both $l$-C$_3$H$^+$ and C$_3$H$^-$. We find no observational evidence of $K_a = 1$ lines which should be present were the carrier indeed C$_3$H$^-$. Additionally, we find a strong anti-correlation between the presence of known molecular anions and B11244 in these regions. Finally, we discuss the formation and destruction chemistry of C$_3$H$^-$ in the context of the physical conditions in the regions. Based on these results, we conclude there is little evidence to support the carrier is C$_3$H$^-$.

preprint2011arXiv

The CO-to-H2 Conversion Factor From Infrared Dust Emission Across the Local Group

We estimate the conversion factor relating CO emission to H2 mass, alpha_CO, in five Local Group galaxies that span approximately an order of magnitude in metallicity - M31, M 33, the Large Magellanic Cloud (LMC), NGC 6822, and the Small Magellanic Cloud (SMC). We model the dust mass along the line of sight from infrared (IR) emission and then solve for the alpha_CO that best allows a single gas-to-dust ratio (delta_GDR) to describe each system. This approach remains sensitive to CO-dark envelopes of H2 surrounding molecular clouds. In M 31, M 33, and the LMC we find alpha_CO \approx 3-9 M_sun pc^-2 (K km s^-1)^-1, consistent with the Milky Way value within the uncertainties. The two lowest metallicity galaxies in our sample, NGC 6822 and the SMC (12 + log(O/H) \approx 8.2 and 8.0), exhibit a much higher alpha_CO. Our best estimates are α_NGC6822 \approx 30 M_sun/pc^-2 (K km s^-1)^-1 and α_SMC \approx 70 M_sun/pc^-2 (K km s-1)-1. These results are consistent with the conversion factor becoming CO a strong function of metallicity around 12 + log(O/H) \sim 8.4 - 8.2. We favor an interpretation where decreased dust-shielding leads to the dominance of CO-free envelopes around molecular clouds below this metallicity.

preprint2010arXiv

Molecular Cloud Formation and the Star Formation Efficiency in M~33

Does star formation proceed in the same way in large spirals such as the Milky Way and in smaller chemically younger galaxies? Earlier work suggests a more rapid transformation of H$_2$ into stars in these objects but (1) a doubt remains about the validity of the H$_2$ mass estimates and (2) there is currently no explanation for why star formation should be more efficient. M~33, a local group spiral with a mass $\sim 10$\% and a metallicity half that of the Galaxy, represents a first step towards the metal poor Dwarf Galaxies. We have searched for molecular clouds in the outer disk of M~33 and present here a set of detections of both $^{12}$CO and $^{13}$CO, including the only detections (for both lines) beyond the R$_{25}$ radius in a subsolar metallicity galaxy. The spatial resolution enables mass estimates for the clouds and thus a measure of the $N({\rm H}_2) / I_{\rm CO}$ ratio, which in turn enables a more reliable calculation of the H$_2$ mass. Our estimate for the outer disk of M~33 is $N({\rm H}_2) / I_{\rm CO(1-0)} \sim 5 \times 10^{20} \,{\rm cm^{-2}/(K{\rm \ km \ s^{-1}})}$ with an estimated uncertainty of a factor $\le 2$. While the $^{12/13}$CO line ratios do not provide a reliable measure of $N({\rm H}_2) / I_{\rm CO}$, the values we find are slightly greater than Galactic and corroborate a somewhat higher $N({\rm H}_2) / I_{\rm CO}$ value. Comparing the CO observations with other tracers of the interstellar medium, no reliable means of predicting where CO would be detected was identified. In particular, CO detections were often not directly on local HI or FIR or H$α$ peaks, although generally in regions with FIR emission and high HI column density. The results presented here provide support for the quicker transformation of H$_2$ into stars in M~33 than in large local universe spirals.