| 1. |
- Bron, E., et al.
(author)
-
Tracers of the ionization fraction in dense and translucent gas: I. Automated exploitation of massive astrochemical model grids
- 2021
-
In: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 645
-
Journal article (peer-reviewed)abstract
- Context. The ionization fraction in the neutral interstellar medium (ISM) plays a key role in the physics and chemistry of the ISM, 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 that they are part of. As large field-of-view hyperspectral maps become available, new tracers may be found. The growth of observational datasets is paralleled by the growth of massive modeling datasets and new methods need to be devised to exploit the wealth of information they contain. Aims. We search for the best observable tracers of the ionization fraction based on a grid of astrochemical models, with the broader aim of finding a general automated method applicable to searching for tracers of any unobservable quantity based on grids of models. Methods. We built grids of models that randomly sample a large range of physical conditions (unobservable quantities such as gas density, temperature, elemental abundances, etc.) and computed the corresponding observables (line intensities, column densities) and the ionization fraction. We estimated 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. Results. In both translucent medium and cold dense medium conditions, we found several observable tracers with very good predictive power for the ionization fraction. Many 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. Conclusions. 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.).
|
|
| 2. |
- Einig, Lucas, et al.
(author)
-
Deep learning denoising by dimension reduction: Application to the ORION-B line cubes
- 2023
-
In: Astronomy and Astrophysics. - 0004-6361 .- 1432-0746. ; 677
-
Journal article (peer-reviewed)abstract
- Context. The availability of large bandwidth receivers for millimeter radio telescopes allows for the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain a lot of information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with an inhomogenous signal-to-noise ratio (S/N) are major challenges for consistent analysis and interpretation. Aims. We searched for a denoising method of the low S/N regions of the studied data cubes that would allow the low S/N emission to be recovered without distorting the signals with a high S/N. Methods. We performed an in-depth data analysis of the 13CO and C17O (1-0) data cubes obtained as part of the ORION-B large program performed at the IRAM 30 m telescope. We analyzed 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 has allowed us to propose significant improvements of typical autoassociative neural networks, often used to denoise hyperspectral Earth remote sensing data. Applying this method to the 13CO (1-0) cube, we were able to 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 in the S/N in voxels with a weak signal, while preserving the spectral shape of the data in high S/N 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 to be a promising avenue. In addition, dealing with the multiplicative noise associated with the calibration uncertainty at high S/N would also be beneficial for such large data cubes.
|
|
| 3. |
- Gaudel, Mathilde, et al.
(author)
-
Gas kinematics around filamentary structures in the Orion B cloud
- 2023
-
In: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 670
-
Journal article (peer-reviewed)abstract
- Context. 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. Aims. In the framework of the IRAM 30m ORION-B large program, we obtained observations of the translucent (2 ≤ AV < 6 mag) and moderately dense gas (6 ≤ AV < 15 mag), which we used to analyze the kinematics over a field of 5 deg2 around the filamentary structures. Methods. We used the Regularized Optimization for Hyper-Spectral Analysis (ROHSA) algorithm to decompose and de-noise the C 18 O(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. Results. 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 H II region. Conclusions. This is the first observational study to highlight feedback from H II 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.
|
|
| 4. |
- Gratier, Pierre, et al.
(author)
-
Quantitative inference of the H2 column densities from 3mm molecular emission: case study towards Orion B
- 2021
-
In: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 645
-
Journal article (peer-reviewed)abstract
- Context. Based on the finding that molecular hydrogen is unobservable in cold molecular clouds, the column density measurements of molecular gas currently rely either on dust emission observation in the far-infrared, which requires space telescopes, or on star counting, which is limited in angular resolution by the stellar density. The (sub)millimeter observations of numerous trace molecules can be effective using ground-based telescopes, but the relationship between the emission of one molecular line and the H-2 column density is non-linear and sensitive to excitation conditions, optical depths, and abundance variations due to the underlying physico- chemistry.Aims. We aim to use multi-molecule line emission to infer the H-2 molecular column density from radio observations.Methods. We propose a data-driven approach to determine the H-2 gas column densities from radio molecular line observations. We use supervised machine-learning methods (random forest) on wide-field hyperspectral IRAM-30m observations of the Orion B molecular cloud to train a predictor of the H-2 column density, using a limited set of molecular lines between 72 and 116 GHz as input, and the Herschel-based dust-derived column densities as "ground truth" output.Results. For conditions similar to those of the Orion B molecular cloud, we obtained predictions of the H-2 column density within a typical factor of 1.2 from the Herschel-based column density estimates. A global analysis of the contributions of the different lines to the predictions show that the most important lines are (CO)-C-13(1-0), (CO)-C-12(1-0), (CO)-O-18(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 that the N2H+(1-0) and CH3OH(2(0)-1(0)) lines be added for the prediction of the H-2 column density in dense core conditions.Conclusions. This article opens a promising avenue for advancing direct inferencing of important physical parameters from the molecular line emission in the millimeter domain. The next step will be to attempt to infer several parameters simultaneously (e.g., the column density and far-UV illumination field) to further test the method.
|
|
| 5. |
- Roueff, Antoine, et al.
(author)
-
Bias versus variance when fitting multi-species molecular lines with a non-LTE radiative transfer model: Application to the estimation of the gas temperature and volume density
- 2024
-
In: Astronomy and Astrophysics. - 0004-6361 .- 1432-0746. ; 686
-
Journal article (peer-reviewed)abstract
- Context. Robust radiative transfer techniques are requisite for efficiently extracting the physical and chemical information from molecular rotational lines. Aims. 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. Methods. 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. Results. 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. Conclusions. 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.
|
|
| 6. |
- Roueff, Antoine, et al.
(author)
-
C18O, 13CO, and 12CO abundances and excitation temperatures in the Orion B molecular cloud: Analysis of the achievable precision in modeling spectral lines within the approximation of the local thermodynamic equilibrium
- 2021
-
In: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 645
-
Journal article (peer-reviewed)abstract
- Context. CO isotopologue transitions are routinely observed in molecular clouds for the purpose of probing the column density of the gas and the elemental ratios of carbon and oxygen, in addition to tracing the kinematics of the environment. Aims. Our study is aimed 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, which includes IC 434, NGC 2023, and the Horsehead pillar. Methods. We used the Cramer Rao bound (CRB) technique to analyze and estimate the precision of the physical parameters in the framework of local-thermodynamic-equilibrium (LTE) excitation and radiative transfer with added 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 proves efficient for a common range of excitation temperatures and column densities (Tex > 6 K, N > 1014-1015 cm-2). Results. Contrary to general assumptions, the various 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 the 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. Conclusions. The CRB technique is a promising avenue for analyzing the estimation of physical parameters from the fit of spectral lines. Future works will generalize its application to non-LTE excitation and radiative transfer methods.
|
|
| 7. |
- Santa-Maria, Miriam G., et al.
(author)
-
HCN emission from translucent gas and UV-illuminated cloud edges revealed by wide-field IRAM 30 m maps of the Orion B GMC: Revisiting its role as a tracer of the dense gas reservoir for star formation
- 2023
-
In: Astronomy and Astrophysics. - 0004-6361 .- 1432-0746. ; 679
-
Journal article (peer-reviewed)abstract
- Context. Massive stars form within dense clumps inside giant molecular clouds (GMCs). Finding appropriate chemical tracers of the dense gas (n(H2) > several 104 cm-3 or AV > 8 mag) and linking their line luminosity with the star formation rate is of critical importance. Aims. Our aim is to determine the origin and physical conditions of the HCN-emitting gas and study their relation to those of other molecules. Methods. In the context of the IRAM 30m ORION-B large program, we present 5 deg2 (~250 pc2) 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 compare the observed HCN line intensities with radiative transfer models including line overlap effects and electron excitation. Furthermore, we study the HCN/HNC isomeric abundance ratio with updated photochemical models. Results. We spectroscopically resolve the HCN J = 1-0 hyperfine structure (HFS) components (and partially resolved J = 2-1 and 3-2 components). We detect anomalous HFS line intensity (and line width) ratios almost everywhere in the cloud. About 70% of the total HCN J = 1-0 luminosity, L′(HCN J = 1-0) = 110 K km s-1 pc-2, arises from AV < 8 mag. The HCN/CO J = 1-0 line intensity ratio, widely used as a tracer of the dense gas fraction, shows a bimodal behavior with an inflection point at AV < 3 mag typical of translucent gas and illuminated cloud edges. We find that most of the HCN J = 1-0 emission arises from extended gas with n(H2) < 104 cm-3, and even lower density gas if the ionization fraction is χe ≥10-5 and electron excitation dominates. This result contrasts with the prevailing view of HCN J = 1-0 emission as a tracer of dense gas and explains the low-AV branch of the HCN/CO J = 1-0 intensity ratio distribution. Indeed, the highest HCN/CO ratios (~ 0.1) at AV < 3 mag correspond to regions of high [CI] 492 GHz/CO J = 1-0 intensity ratios (>1) characteristic of low-density photodissociation regions. The low surface brightness (≲ 1 K km s-1) and extended HCN and HCO+ J = 1-0 emission scale with IFIR -a proxy of the stellar far-ultraviolet (FUV) radiation field -in a similar way. Together with CO J = 1-0, these lines respond to increasing IFIR up to G0 ≅ 20. On the other hand, the bright HCN J = 1-0 emission (> 6 K km s-1) from dense gas in star-forming clumps weakly responds to IFIR once the FUV field becomes too intense (G0 > 1500). In contrast, HNC J = 1-0 and [CI] 492 GHz lines weakly respond to IFIR for all G0. 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. Conclusions. Given the widespread and extended nature of the [CI] 492 GHz emission, as well as its spatial correlation with that of HCO+, HCN, and 13CO J = 1-0 lines (in this order), we argue that the edges of GMCs are porous to FUV radiation from nearby massive stars. Enhanced FUV radiation favors the formation and excitation of HCN on large scales, not only in dense star-forming clumps, and it leads to a relatively low value of the dense gas mass to total luminosity ratio, α (HCN) = 29 M⊙ /(K km s-1pc2) in Orion B. 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, which may be better traced by HNC than by HCN.
|
|
