Issue |
A&A
Volume 686, June 2024
|
|
---|---|---|
Article Number | A255 | |
Number of page(s) | 37 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/202449148 | |
Published online | 20 June 2024 |
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
1
Université de Toulon, Aix Marseille Univ, CNRS, IM2NP,
Toulon,
France
e-mail: antoine.roueff@univ-tln.fr
2
IRAM,
300 rue de la Piscine,
38406
Saint-Martin-d’Hères,
France
3
LERMA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités,
75014
Paris,
France
4
Instituto de Física Fundamental (CSIC).
Calle Serrano 121,
28006,
Madrid,
Spain
5
National Radio Astronomy Observatory,
520 Edgemont Road,
Charlottesville,
VA
22903,
USA
6
Laboratoire d’Astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N,
Allée Geoffroy Saint-Hilaire,
33615
Pessac,
France
7
Univ. Grenoble Alpes, Inria, CNRS, Grenoble INP, GIPSA-Lab,
Grenoble
38000,
France
8
Chalmers University of Technology, Department of Space, Earth and Environment,
412 93
Gothenburg,
Sweden
9
Univ. Lille, CNRS, Centrale Lille,
UMR 9189 - CRIStAL,
59651
Villeneuve d’Ascq,
France
10
LERMA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités,
92190
Meudon,
France
11
Institut de Recherche en Astrophysique et Planétologie (IRAP), Université Paul Sabatier,
Toulouse cedex 4,
France
12
Instituto de Astrofísica, Pontificia Universidad Católica de Chile,
Av. Vicuña Mackenna 4860,
7820436 Macul, Santiago,
Chile
13
Laboratoire de Physique de l’Ecole normale supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université de Paris, Sorbonne Paris Cité,
Paris,
France
14
Jet Propulsion Laboratory, California Institute of Technology,
4800 Oak Grove Drive,
Pasadena,
CA
91109,
USA
15
Department of Earth, Environment, and Physics, Worcester State University,
Worcester,
MA
01602,
USA
16
Harvard-Smithsonian Center for Astrophysics,
60 Garden Street,
Cambridge,
MA,
02138,
USA
17
School of Physics and Astronomy, Cardiff University,
Queen’s buildings,
Cardiff
CF24 3AA,
UK
Received:
2
January
2024
Accepted:
19
February
2024
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.
Key words: line: profiles / radiative transfer / methods: data analysis / methods: statistical / ISM: clouds / ISM: general
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.