Spatially Coherent 3D Distributions of HI and CO in the Milky Way

TL;DR

This paper presents a novel Bayesian approach using 3D Gaussian processes to reconstruct the 3D spatial distributions of HI and CO gases in the Milky Way, addressing previous limitations in resolution, optical depth assumptions, and velocity field modeling.

Contribution

It introduces a new method that models correlations in the interstellar medium and accounts for absorption effects, resulting in high-fidelity, spatially coherent 3D maps of Galactic gas components.

Findings

Maps confirm Galactic disc warping and flaring.

Good agreement with 3D dust maps on scales >400 pc.

Enhanced resolution and accuracy over previous models.

Abstract

The spatial distribution of the gaseous components of the Milky Way is of great importance for a number of different fields, for example, Galactic structure, star formation, and cosmic rays. However, obtaining distance information to gaseous clouds in the interstellar medium from Doppler-shifted line emission is notoriously difficult given our vantage point in the Galaxy. It requires spatial knowledge of gas velocities and generally suffers from distance ambiguities. Previous works often assumed the optically thin limit, had a fixed velocity field, and lacked resolution overall. We aim to overcome these issues and improve previous reconstructions of the gaseous constituents of the interstellar medium of the Galaxy. We used 3D Gaussian processes to model correlations in the interstellar medium, including correlations between different lines of sight, and enforce a spatially coherent structure in the prior. For modelling the transport of radiation from the emitting gas to us as observers, we took absorption effects into account. A special numerical grid ensures that there is high resolution nearby. We inferred the spatial distributions of atomic hydrogen (HI), carbon monoxide (CO), their emission line widths, and the Galactic velocity field in a joint Bayesian inference. We further constrained these fields with complementary data from Galactic masers and young stellar object clusters. Our main result consists of a set of samples that implicitly contain statistical uncertainties. The resulting maps are spatially coherent and reproduce the data with high fidelity. We confirm previous findings regarding the warping and flaring of the Galactic disc. A comparison with 3D dust maps reveals a good agreement on scales larger than approximately 400 pc. While our results are not free of artefacts, they present a big step forward in obtaining high-quality 3D maps of the interstellar medium.