-
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.
-
The magnetic field in the Flame nebula
Authors:
Ivana Bešlić,
Simon Coudé,
Dariusz C. Lis,
Maryvonne Gerin,
Paul F. Goldsmith,
Jerome Pety,
Antoine Roueff,
Karine Demyk,
Charles D. Dowell,
Lucas Einig,
Javier R. Goicoechea,
Francois Levrier,
Jan Orkisz,
Nicolas Peretto,
Miriam G. Santa-Maria,
Nathalie Ysard,
Antoine Zakardjian
Abstract:
Star formation is essential in galaxy evolution and the cycling of matter. The support of interstellar clouds against gravitational collapse by magnetic (B-) fields has been proposed to explain the low observed star formation efficiency in galaxies and the Milky Way. Despite the Planck satellite providing a 5-15' all-sky map of the B-field geometry in the diffuse interstellar medium, higher spatia…
▽ More
Star formation is essential in galaxy evolution and the cycling of matter. The support of interstellar clouds against gravitational collapse by magnetic (B-) fields has been proposed to explain the low observed star formation efficiency in galaxies and the Milky Way. Despite the Planck satellite providing a 5-15' all-sky map of the B-field geometry in the diffuse interstellar medium, higher spatial resolution observations are required to understand the transition from diffuse gas to gravitationally unstable filaments. NGC 2024, the Flame Nebula, in the nearby Orion B molecular cloud, contains a young, expanding HII region and a dense filament that harbors embedded protostellar objects. Therefore, NGC 2024 is an excellent opportunity to study the role of B-fields in the formation, evolution, and collapse of filaments, as well as the dynamics and effects of young HII regions on the surrounding molecular gas. We combine new 154 and 216 micron dust polarization measurements carried out using the HAWC+ instrument aboard SOFIA with molecular line observations of 12CN(1-0) and HCO+(1-0) from the IRAM 30-meter telescope to determine the B-field geometry and to estimate the plane of the sky magnetic field strength across the NGC 2024. The HAWC+ observations show an ordered B-field geometry in NGC 2024 that follows the morphology of the expanding HII region and the direction of the main filament. The derived plane of the sky B-field strength is moderate, ranging from 30 to 80 micro G. The strongest B-field is found at the northern-west edge of the HII region, characterized by the highest gas densities and molecular line widths. In contrast, the weakest field is found toward the filament in NGC 2024. The B-field has a non-negligible influence on the gas stability at the edges of the expanding HII shell (gas impacted by the stellar feedback) and the filament (site of the current star formation).
△ Less
Submitted 7 February, 2024; v1 submitted 30 January, 2024;
originally announced January 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.
-
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.
-
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.
-
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.
-
Bird's eye view of molecular clouds in the Milky Way: II. Cloud kinematics from sub-pc to kpc scales
Authors:
Andri Spilker,
Jouni Kainulainen,
Jan Orkisz
Abstract:
The kinematics of molecular gas are crucial for setting the stage for star formation. One key question related to the kinematic properties of gas is how they depend on the spatial scale. We aim to describe the CO spectra, velocity dispersions, and especially the linewidth-size relation, of molecular gas from cloud (parsec-) scales to kiloparsec scales in a complete region within the Milky Way disk…
▽ More
The kinematics of molecular gas are crucial for setting the stage for star formation. One key question related to the kinematic properties of gas is how they depend on the spatial scale. We aim to describe the CO spectra, velocity dispersions, and especially the linewidth-size relation, of molecular gas from cloud (parsec-) scales to kiloparsec scales in a complete region within the Milky Way disk. We utilise the census of molecular clouds within 2 kpc from our earlier work, together with CO emission data for them from the literature. We study the kinematics and the Larson's relations for the sample of individual clouds. We also mimic a face-on view of the Milky Way and analyse the kinematics of the clouds within apertures of 0.25-2 kpc in size. In this way, we describe the scale-dependency of the CO gas kinematics and Larson's relations. We describe the spectra of CO gas at cloud scales and in apertures between 0.25-2 kpc in our survey area. The spectra within the apertures are relatively symmetric but show non-Gaussian high-velocity wings. At cloud-scales, our sample shows a linewidth-size relation σ_v=1.5*R^{0.3\pm0.1} with a large scatter. The mass-size relation in the sample of clouds is M_{CO}= 794*R^{1.5\pm0.5}. The relations are also present for the apertures at kpc-scales. A suggestive dependency on galactic environment is seen, with apertures closer to the Galactic centre and the Sagittarius spiral arm having slightly higher velocity dispersions. We explore the possible effect of a diffuse component in the survey area, and find that such a component would widen the CO spectra and could flatten the linewidth-size relation. Understanding the nature of the possible diffuse CO component and its effects on observations is crucial for connecting Galactic and extragalactic data.
△ Less
Submitted 13 September, 2022;
originally announced September 2022.
-
The effect of viewing angle on the Kennicutt-Schmidt relation of the local molecular clouds
Authors:
Jouni Kainulainen,
Sara Rezaei Kh.,
Andri Spilker,
Jan Orkisz
Abstract:
The Gaia data give us an unprecedented view to the 3-dimensional (3D) structure of molecular clouds in the Solar neighbourhood. We study how the projected areas and masses of clouds, and consequently the Kennicutt-Schmidt relation (KS-relation), depend on the viewing angle. We derive the probability distributions of the projected areas and masses for nine clouds within 400 pc from the Sun using 3D…
▽ More
The Gaia data give us an unprecedented view to the 3-dimensional (3D) structure of molecular clouds in the Solar neighbourhood. We study how the projected areas and masses of clouds, and consequently the Kennicutt-Schmidt relation (KS-relation), depend on the viewing angle. We derive the probability distributions of the projected areas and masses for nine clouds within 400 pc from the Sun using 3D dust distribution data from the literature. We find that the viewing angle can have a dramatic effect on the observed areas and masses of individual clouds. The joint probability distributions of the areas and masses are strongly correlated, relatively flat, and can show multiple peaks. The typical ranges and 50% quartiles of the distributions are roughly 100-200% and 20-80% of the median value, respectively, making viewing angle effects important for all individual clouds. The threshold value used to define the cloud areas is also important; our analysis suggests that the clouds become more anisotropic for smaller thresholds (larger scales). On average, the areas and masses of the plane-of-the-sky and face-on projections agree, albeit with a large scatter. This suggests that sample averages of areas and masses are relatively free of viewing angle effects, which is important to facilitate comparisons of extragalactic and galactic data. Ultimately, our results demonstrate that a cloud's location in the KS-relation is affected by viewing angle in a non-trivial manner. However, the KS-relation of our sample as a whole is not strongly affected by these effects, because the co-variance of the areas and masses causes the observed mean column density to remain relatively constant.
△ Less
Submitted 15 November, 2021;
originally announced November 2021.
-
Bird's eye view of molecular clouds in the Milky Way: I. Column density and star formation from sub-pc to kpc scales
Authors:
Andri Spilker,
Jouni Kainulainen,
Jan Orkisz
Abstract:
Describing how the properties of the interstellar medium combine across size-scales is crucial for understanding star formation scaling laws and connecting Galactic and extragalactic data of molecular clouds. We describe how the statistical structure of clouds, and its connection to star formation, changes from sub-parsec to kiloparsec scales in a complete region within the Milky Way disk. We buil…
▽ More
Describing how the properties of the interstellar medium combine across size-scales is crucial for understanding star formation scaling laws and connecting Galactic and extragalactic data of molecular clouds. We describe how the statistical structure of clouds, and its connection to star formation, changes from sub-parsec to kiloparsec scales in a complete region within the Milky Way disk. We build a census of molecular clouds within 2 kpc from the Sun using literature. We examine the dust-based column density probability distributions (N-PDFs) of the clouds and their relation to star formation traced by young stellar objects (YSOs). We then examine our survey region from the outside, within apertures of varying sizes, and describe how the N-PDFs and their relation to star formation changes with the size-scale. The N-PDFs of the clouds are not well described by any single simple model; use of any single model may bias the interpretation of the N-PDFs. The top-heaviness of the N-PDFs correlates with star formation activity, and the correlation changes with Galactic environment (spiral-/inter-arm regions). We find that the density contrast of clouds may be more intimately linked to star formation than the dense gas mass fraction. The aperture-averaged N-PDFs vary with the size-scale and are more top-heavy for larger apertures. The top-heaviness of the aperture N-PDFs correlates with star formation activity up to roughly 0.5 kpc size-scale, depending on the environment. Our results suggest that the relations between cloud structure and star formation are environment specific and best captured by relative quantities (e.g., the density contrast). Finally, we show how the density structures of individual clouds give rise to a kpc-scale Kennicutt-Schmidt relationship as a combination of sampling effects and blending of different galactic environments.
△ Less
Submitted 2 September, 2021; v1 submitted 10 August, 2021;
originally announced August 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.
-
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.
-
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.
-
Autonomous Gaussian decomposition of the Galactic Ring Survey. I. Global statistics and properties of the 13CO emission data
Authors:
M. Riener,
J. Kainulainen,
H. Beuther,
J. D. Henshaw,
J. H. Orkisz,
Y. Wang
Abstract:
The analysis of large molecular line surveys of the Galactic plane is essential for our understanding of the gas kinematics on Galactic scales, in particular its link with the formation and evolution of dense structures in the interstellar medium. An approximation of the emission peaks with Gaussian functions allows for an efficient and straightforward extraction of useful physical information con…
▽ More
The analysis of large molecular line surveys of the Galactic plane is essential for our understanding of the gas kinematics on Galactic scales, in particular its link with the formation and evolution of dense structures in the interstellar medium. An approximation of the emission peaks with Gaussian functions allows for an efficient and straightforward extraction of useful physical information contained in the shape and Doppler-shifted frequency of the emission lines contained in these enormous data sets. In this work we present an overview and first results of a Gaussian decomposition of the entire Galactic Ring Survey (GRS) 13CO (1-0) data that consists of about 2.3 million spectra. We performed the decomposition with the fully automated GaussPy+ algorithm and fitted about 4.6 million Gaussian components to the GRS spectra. We discuss the statistics of the fit components and relations between the fitted intensities, velocity centroids, and velocity dispersions. We find that the magnitude of the velocity dispersion values increase toward the inner Galaxy and around the Galactic midplane, which we speculate is partly due to the influence of the Galactic bar and regions with higher non-thermal motions located in the midplane, respectively. We also use our decomposition results to infer global properties of the gas emission and find that the number of fit components used per spectrum is indicative for the amount of structure along the line of sight. We find that the emission lines from regions located on the far side of the Galaxy show increased velocity dispersion values, likely due to beam averaging effects. We demonstrate how this trend has the potential to aid in characterising Galactic structure by disentangling emission that is belonging to the nearby Aquila Rift molecular cloud from emission that is more likely associated with the Perseus and Outer spiral arms.
△ Less
Submitted 15 November, 2019;
originally announced November 2019.
-
GaussPy+: A fully automated Gaussian decomposition package for emission line spectra
Authors:
Manuel Riener,
Jouni Kainulainen,
Jonathan D. Henshaw,
Jan H. Orkisz,
Claire E. Murray,
Henrik Beuther
Abstract:
Our understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity structure of emission line observations. One approach to study the velocity structure is to decompose the spectra into individual velocity components; this leads to a description of the dataset that is significantly reduced in complexity. However, this decomposition requires full automati…
▽ More
Our understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity structure of emission line observations. One approach to study the velocity structure is to decompose the spectra into individual velocity components; this leads to a description of the dataset that is significantly reduced in complexity. However, this decomposition requires full automation lest it becomes prohibitive for large datasets, such as Galactic plane surveys. We developed GaussPy+, a fully automated Gaussian decomposition package that can be applied to emission line datasets, especially large surveys of HI and isotopologues of CO. We built our package upon the existing GaussPy algorithm and significantly improved its performance for noisy data. New functionalities of GaussPy+ include: i) automated preparatory steps, such as an accurate noise estimation, which can also be used as standalone applications; ii) an improved fitting routine; iii) an automated spatial refitting routine that can add spatial coherence to the decomposition results by refitting spectra based on neighbouring fit solutions. We thoroughly tested the performance of GaussPy+ on synthetic spectra and a test field from the Galactic Ring Survey. We found that GaussPy+ can deal with cases of complex emission and even low to moderate signal-to-noise values.
△ Less
Submitted 25 June, 2019;
originally announced June 2019.
-
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.
-
Exploring the inner parsecs of Active Galactic Nuclei using near-infrared high resolution polarimetric simulations with MontAGN
Authors:
Lucas Grosset,
Daniel Rouan,
Damien Gratadour,
Didier Pelat,
Jan Orkisz,
Frédéric Marin,
René Goosmann
Abstract:
Aims: In this paper, we aim to constrain the properties of dust structures in the central first parsecs of Active Galactic Nuclei (AGN). Our goal is to study the required optical depth and composition of different dusty and ionised structures.
Methods: We developed a radiative transfer code, MontAGN, optimised for polarimetric observations in the infrared. With both this code and STOKES, designe…
▽ More
Aims: In this paper, we aim to constrain the properties of dust structures in the central first parsecs of Active Galactic Nuclei (AGN). Our goal is to study the required optical depth and composition of different dusty and ionised structures.
Methods: We developed a radiative transfer code, MontAGN, optimised for polarimetric observations in the infrared. With both this code and STOKES, designed to be relevant from the hard X-ray band to near-infrared wavelengths, we investigate the polarisation emerging from a characteristic model of the AGN environment. For that purpose, we compare predictions of our models with previous infrared observations of NGC 1068, and try to reproduce several key polarisation patterns revealed by polarisation mapping.
Results: We constrain the required dust structures as well as their densities. More precisely, we find out that the electron density inside the ionisation cone is about $2.0 \times 10^9$ m$^{-3}$ . With structures constituted of spherical grains of constant density, we also highlight that the torus should be thicker than 20 in term of K band optical depth to block direct light from the centre. It should also have a stratification in density, with a lesser dense outer rim with an optical depth at 2.2 μm typically between 0.8 and 4 for observing the double scattering effect previously proposed.
Conclusions: We bring constraints on the dust structures in the inner parsecs of an AGN model supposed to describe NGC 1068. When compared to observations, this leads to optical depth of at least 20 in Ks band for the torus of NGC 1068, corresponding to $τ_V \approx 170$, which is within the range of current estimation based on observations. In the future, we will improve our study by including non uniform dust structures and aligned elongated grains to constrain other possible interpretations of the observations.
△ Less
Submitted 6 December, 2017;
originally announced December 2017.
-
Clustering the Orion B giant molecular cloud based on its molecular emission
Authors:
Emeric Bron,
Chloé Daudon,
Jérôme Pety,
François Levrier,
Maryvonne Gerin,
Pierre Gratier,
Jan H. Orkisz,
Viviana Guzman,
Sébastien Bardeau,
Javier R. Goicoechea,
Harvey Liszt,
Karin Öberg,
Nicolas Peretto,
Albrecht Sievers,
Pascal Tremblin
Abstract:
Previous attempts at segmenting molecular line maps of molecular clouds have focused on using position-position-velocity data cubes of a single line to separate the spatial components of the cloud. In contrast, wide field spectral imaging with large spectral bandwidth in the (sub)mm domain now allows to combine multiple molecular tracers to understand the different physical and chemical phases tha…
▽ More
Previous attempts at segmenting molecular line maps of molecular clouds have focused on using position-position-velocity data cubes of a single line to separate the spatial components of the cloud. In contrast, wide field spectral imaging with large spectral bandwidth in the (sub)mm domain now allows to combine multiple molecular tracers to understand the different physical and chemical phases that constitute giant molecular clouds. We aim at using multiple tracers (sensitive to different physical processes) to segment a molecular cloud into physically/chemically similar regions (rather than spatially connected components). We use a machine learning clustering method (the Meanshift algorithm) to cluster pixels with similar molecular emission, ignoring spatial information. Simple radiative transfer models are used to interpret the astrophysical information uncovered by the clustering. A clustering analysis based only on the J=1-0 lines of 12CO, 13CO and C18O reveals distinct density/column density regimes (nH~100, 500, and >1000 cm-3), closely related to the usual definitions of diffuse, translucent and high-column-density regions. Adding two UV-sensitive tracers, the (1-0) lines of HCO+ and CN, allows us to distinguish two clearly distinct chemical regimes, characteristic of UV-illuminated and UV-shielded gas. The UV-illuminated regime shows overbright HCO+ and CN emission, which we relate to photochemical enrichment. We also find a tail of high CN/HCO+ intensity ratio in UV-illuminated regions. Finer distinctions in density classes (nH~7E3, and 4E4 cm-3) for the densest regions are also identified, likely related to the higher critical density of the CN and HCO+ (1-0) lines. The association of simultaneous multi-line, wide-field mapping and powerful machine learning methods such as the Meanshift algorithm reveals how to decode the complex information available in molecular tracers.
△ Less
Submitted 19 October, 2017;
originally announced October 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.
-
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.
-
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.