-
Spatially Coherent 3D Distributions of HI and CO in the Milky Way
Authors:
Laurin Söding,
Gordian Edenhofer,
Torsten A. Enßlin,
Philipp Frank,
Ralf Kissmann,
Vo Hong Minh Phan,
Andrés Ramírez,
Hanieh Zhandinejad,
Philipp Mertsch
Abstract:
The spatial distribution of the gaseous components of the Milky Way is of great importance for a number of different fields, e.g. 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 unique vantage point in the Galaxy. It requires precise knowled…
▽ More
The spatial distribution of the gaseous components of the Milky Way is of great importance for a number of different fields, e.g. 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 unique vantage point in the Galaxy. It requires precise knowledge of gas velocities and generally suffers from distance ambiguities.
Previous works often assumed the optically thin limit (no absorption), a fixed velocity field, and lack resolution overall. We aim to overcome these issues and improve previous reconstructions of the gaseous constituents of the interstellar medium of the Galaxy.
We use 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 take absorption effects into account. A special numerical grid ensures high resolution nearby. We infer the spatial distributions of HI, CO, their emission line-widths, and the Galactic velocity field in a joint Bayesian inference. We further constrain 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.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
The influence of the 3D Galactic gas structure on cosmic-ray transport and gamma-ray emission
Authors:
Andrés Ramírez,
Gordian Edenhofer,
Torsten A. Enßlin,
Philipp Frank,
Philipp Mertsch,
Vo Hong Minh Phan,
Laurin Söding,
Hanieh Zhandinejad,
Ralf Kissmann
Abstract:
Cosmic rays (CRs) play a major role in the dynamics of the interstellar medium (ISM). Their interactions and transport ionize, heat, and push the ISM thereby coupling different regions of it. The spatial distribution of CRs depends on the distribution of their sources as well as the ISM constituents they interact with, such as gas, starlight, and magnetic fields. Particularly, gas interacts closel…
▽ More
Cosmic rays (CRs) play a major role in the dynamics of the interstellar medium (ISM). Their interactions and transport ionize, heat, and push the ISM thereby coupling different regions of it. The spatial distribution of CRs depends on the distribution of their sources as well as the ISM constituents they interact with, such as gas, starlight, and magnetic fields. Particularly, gas interacts closely with CRs, influencing CR fluxes and gamma -ray emission. We illustrate the influence of 3D gas structures on CR transport and gamma -ray emission. We use the PICARD code and multiple samples of recent 3D reconstructions of the HI and H$_2$ Galactic gas constituents to investigate the impact on the transport of CRs and emission of gamma -rays. We find the necessary transport parameters to reproduce local measurements of CR fluxes, and see that they depend on the local distribution of gas density and structure. The distribution of CR fluxes exhibits energy-dependent structures that vary for all CR species due to their corresponding loss processes. Regions of enhanced secondary (primary) species are spatially correlated (anti-correlated) with the gas density. We observe a high sensitivity of the gamma -ray emission on the contrast of gas structures, as those determine the 3D spatial distributions of hadronic interactions and bremsstrahlung. We find that corresponding gas-induced structures in the distribution of CR electrons are also visible in Inverse Compton (IC) emission. Due to the aforementioned sensitivity, the analysis of CR data for CR sources and transport parameters requires the usage of accurate 3D gas maps.
△ Less
Submitted 2 July, 2024;
originally announced July 2024.
-
The "C": The large Chameleon-Musca-Coalsack cloud
Authors:
Gordian Edenhofer,
João Alves,
Catherine Zucker,
Laura Posch,
Torsten A. Enßlin
Abstract:
Recent advancements in 3D dust mapping have transformed our understanding of the Milky Way's local interstellar medium, enabling us to explore its structure in three spatial dimensions for the first time. In this Letter, we use the most recent 3D dust map by Edenhofer et al. to study the well-known Chameleon, Musca, and Coalsack cloud complexes, located about 200 pc from the Sun. We find that thes…
▽ More
Recent advancements in 3D dust mapping have transformed our understanding of the Milky Way's local interstellar medium, enabling us to explore its structure in three spatial dimensions for the first time. In this Letter, we use the most recent 3D dust map by Edenhofer et al. to study the well-known Chameleon, Musca, and Coalsack cloud complexes, located about 200 pc from the Sun. We find that these three complexes are not isolated but rather connect to form a surprisingly well-defined half-ring, constituting a single C-shaped cloud with a radius of about 50 pc, a thickness of about 45 pc, and a total mass of about $5 \times 10^{4} \mathrm{M}_{\odot}$, or $9 \times 10^{4} \mathrm{M}_{\odot}$ if including everything in the vicinity of the C-shaped cloud. Despite the absence of an evident feedback source at its center, the dynamics of young stellar clusters associated with the C structure suggest that a single supernova explosion about 4 Myr to 10 Myr ago likely shaped this structure. Our findings support a single origin story for these cloud complexes, suggesting that they were formed by feedback-driven gas compression, and offer new insights into the processes that govern the birth of star-forming clouds in feedback-dominated regions, such as the Scorpius-Centaurus association.
△ Less
Submitted 29 June, 2024; v1 submitted 15 April, 2024;
originally announced April 2024.
-
JAXbind: Bind any function to JAX
Authors:
Jakob Roth,
Martin Reinecke,
Gordian Edenhofer
Abstract:
JAX is widely used in machine learning and scientific computing, the latter of which often relies on existing high-performance code that we would ideally like to incorporate into JAX. Reimplementing the existing code in JAX is often impractical and the existing interface in JAX for binding custom code either limits the user to a single Jacobian product or requires deep knowledge of JAX and its C++…
▽ More
JAX is widely used in machine learning and scientific computing, the latter of which often relies on existing high-performance code that we would ideally like to incorporate into JAX. Reimplementing the existing code in JAX is often impractical and the existing interface in JAX for binding custom code either limits the user to a single Jacobian product or requires deep knowledge of JAX and its C++ backend for general Jacobian products. With JAXbind we drastically reduce the effort required to bind custom functions implemented in other programming languages with full support for Jacobian-vector products and vector-Jacobian products to JAX. Specifically, JAXbind provides an easy-to-use Python interface for defining custom, so-called JAX primitives. Via JAXbind, any function callable from Python can be exposed as a JAX primitive. JAXbind allows a user to interface the JAX function transformation engine with custom derivatives and batching rules, enabling all JAX transformations for the custom primitive.
△ Less
Submitted 27 June, 2024; v1 submitted 13 March, 2024;
originally announced March 2024.
-
Non-parametric Bayesian reconstruction of Galactic magnetic fields using Information Field Theory: The inclusion of line-of-sight information in ultra-high energy cosmic ray backtracking
Authors:
Alexandros Tsouros,
Abhijit B. Bendre,
Gordian Edenhofer,
Torsten Enßlin,
Philipp Frank,
Michalis Mastorakis,
Vasiliki Pavlidou
Abstract:
(abridged) Ultra-high energy cosmic rays (UHECRs) are extremely energetic charged particles with energies surpassing $10^{18}$ eV. Their sources remain elusive, obscured by deflections caused by the Galactic magnetic field (GMF). This challenge is further complicated by our limited understanding of the three-dimensional structure of the GMF, as current GMF observations consist primarily of quantit…
▽ More
(abridged) Ultra-high energy cosmic rays (UHECRs) are extremely energetic charged particles with energies surpassing $10^{18}$ eV. Their sources remain elusive, obscured by deflections caused by the Galactic magnetic field (GMF). This challenge is further complicated by our limited understanding of the three-dimensional structure of the GMF, as current GMF observations consist primarily of quantities integrated along the line-of-sight (LOS). Nevertheless, data from upcoming stellar polarisation surveys along with Gaia's stellar parallax data are expected to yield local GMF measurements.. In this work, we employ methods of Bayesian statistical inference in order to sample the posterior distribution of the GMF within part of the Galaxy. By assuming a known rigidity and arrival direction of an UHECR, we backtrack its trajectory through various GMF configurations drawn from the posterior distribution. Our objective is to rigorously evaluate our algorithm's performance in scenarios that closely mirror the setting of expected future applications. In pursuit of this, we condition the posterior to synthetic integrated LOS measurements of the GMF, in addition to synthetic local POS-component measurements. In this proof of concept work, we assume the ground truth to be a magnetic field produced by a dynamo simulation of the Galactic ISM. Our results demonstrate that for all locations of the observed arrival direction on the POS, our algorithm is able to substantially update our knowledge on the original arrival direction of UHECRs with rigidity $E/Z = 5 \times 10^{19}$ eV, even in the case of complete absence of LOS information. If integrated data is included in the inference, then the regions of the celestial sphere where the maximum error occurs diminishes greatly. Even in those regions the maximum error is diminished by a factor of about $3$ in the specific setting studied.
△ Less
Submitted 8 March, 2024;
originally announced March 2024.
-
The Local Bubble is a Local Chimney: A New Model from 3D Dust Mapping
Authors:
Theo J. O'Neill,
Catherine Zucker,
Alyssa A. Goodman,
Gordian Edenhofer
Abstract:
Leveraging a high-resolution 3D dust map of the solar neighborhood from Edenhofer et al. (2024), we derive a new 3D model for the dust-traced surface of the Local Bubble, the supernova-driven cavity surrounding the Sun. We find that the surface of the Local Bubble is highly irregular in shape, with its peak extinction surface falling at an average distance of 170 pc from the Sun (spanning 70-600+…
▽ More
Leveraging a high-resolution 3D dust map of the solar neighborhood from Edenhofer et al. (2024), we derive a new 3D model for the dust-traced surface of the Local Bubble, the supernova-driven cavity surrounding the Sun. We find that the surface of the Local Bubble is highly irregular in shape, with its peak extinction surface falling at an average distance of 170 pc from the Sun (spanning 70-600+ pc) with a typical thickness of 35 pc and a total dust-traced mass of $(6.0 \pm 0.7) \times 10^5 \ \rm{M}_{\odot}$. The Local Bubble displays an extension in the Galactic Northern hemisphere that is morphologically consistent with representing a "Local Chimney." We argue this chimney was likely created by the "bursting" of this supernova-driven superbubble, leading to the funneling of interstellar medium ejecta into the lower Galactic halo. We find that many well-known dust features and molecular clouds fall on the surface of the Local Bubble and that several tunnels to other adjacent cavities in the interstellar medium may be present. Our new, parsec-resolution view of the Local Bubble may be used to inform future analysis of the evolution of nearby gas and young stars, the investigation of direct links between the solar neighborhood and the Milky Way's lower halo, and numerous other applications.
△ Less
Submitted 11 July, 2024; v1 submitted 7 March, 2024;
originally announced March 2024.
-
Re-Envisioning Numerical Information Field Theory (NIFTy.re): A Library for Gaussian Processes and Variational Inference
Authors:
Gordian Edenhofer,
Philipp Frank,
Jakob Roth,
Reimar H. Leike,
Massin Guerdi,
Lukas I. Scheel-Platz,
Matteo Guardiani,
Vincent Eberle,
Margret Westerkamp,
Torsten A. Enßlin
Abstract:
Imaging is the process of transforming noisy, incomplete data into a space that humans can interpret. NIFTy is a Bayesian framework for imaging and has already successfully been applied to many fields in astrophysics. Previous design decisions held the performance and the development of methods in NIFTy back. We present a rewrite of NIFTy, coined NIFTy.re, which reworks the modeling principle, ext…
▽ More
Imaging is the process of transforming noisy, incomplete data into a space that humans can interpret. NIFTy is a Bayesian framework for imaging and has already successfully been applied to many fields in astrophysics. Previous design decisions held the performance and the development of methods in NIFTy back. We present a rewrite of NIFTy, coined NIFTy.re, which reworks the modeling principle, extends the inference strategies, and outsources much of the heavy lifting to JAX. The rewrite dramatically accelerates models written in NIFTy, lays the foundation for new types of inference machineries, improves maintainability, and enables interoperability between NIFTy and the JAX machine learning ecosystem.
△ Less
Submitted 15 June, 2024; v1 submitted 26 February, 2024;
originally announced February 2024.
-
High-resolution Pan-STARRS and SMA observations of IRAS 23077+6707: A giant edge-on protoplanetary disk
Authors:
Kristina Monsch,
Joshua B. Lovell,
Ciprian T. Berghea,
Gordian Edenhofer,
Garrett K. Keating,
Sean M. Andrews,
Ammar Bayyari,
Jeremy J. Drake,
David J. Wilner
Abstract:
We present resolved images of IRAS 23077+6707 ("Dracula's Chivito") in 1.3 mm/225 GHz thermal dust and CO gas emission with the Submillimeter Array (SMA) and optical (0.5-$0.8\,μ\mathrm{m}$) scattered light with the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). The Pan-STARRS data show a bipolar distribution of optically scattering dust that is characteristic for disks observe…
▽ More
We present resolved images of IRAS 23077+6707 ("Dracula's Chivito") in 1.3 mm/225 GHz thermal dust and CO gas emission with the Submillimeter Array (SMA) and optical (0.5-$0.8\,μ\mathrm{m}$) scattered light with the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). The Pan-STARRS data show a bipolar distribution of optically scattering dust that is characteristic for disks observed at high inclinations. Its scattered light emission spans ${\sim}14''$, with two highly asymmetric filaments extending along the upper bounds of each nebula by ${\sim}9''$. The SMA data measure 1.3 mm continuum dust as well as $^{12}$CO, $^{13}$CO and C$^{18}$O $J$=2$-$1 line emission over $12''$-$14''$ extents, with the gas presenting the typical morphology of a disk in Keplerian rotation, in both position-velocity space and in each CO line spectrum. IRAS 23077+6707 has no reported distance estimate, but if it is located in the Cepheus star-forming region (180-800 pc), it would have a radius spanning thousands of au. Taken together, we infer IRAS 23077+6707 to be a giant and gas-rich edge-on protoplanetary disk, which to our knowledge is the largest in extent so far discovered.
△ Less
Submitted 14 May, 2024; v1 submitted 2 February, 2024;
originally announced February 2024.
-
A parsec-scale Galactic 3D dust map out to 1.25 kpc from the Sun
Authors:
Gordian Edenhofer,
Catherine Zucker,
Philipp Frank,
Andrew K. Saydjari,
Joshua S. Speagle,
Douglas Finkbeiner,
Torsten Enßlin
Abstract:
High-resolution 3D maps of interstellar dust are critical for probing the underlying physics shaping the structure of the interstellar medium, and for foreground correction of astrophysical observations affected by dust. We aim to construct a new 3D map of the spatial distribution of interstellar dust extinction out to a distance of 1.25 kpc from the Sun. We leveraged distance and extinction estim…
▽ More
High-resolution 3D maps of interstellar dust are critical for probing the underlying physics shaping the structure of the interstellar medium, and for foreground correction of astrophysical observations affected by dust. We aim to construct a new 3D map of the spatial distribution of interstellar dust extinction out to a distance of 1.25 kpc from the Sun. We leveraged distance and extinction estimates to 54 million nearby stars derived from the Gaia BP/RP spectra. Using the stellar distance and extinction information, we inferred the spatial distribution of dust extinction. We modeled the logarithmic dust extinction with a Gaussian process in a spherical coordinate system via iterative charted refinement and a correlation kernel inferred in previous work. In total, our posterior has over 661 million degrees of freedom. We probed the posterior distribution using the variational inference method MGVI. Our 3D dust map has an angular resolution of up to 14' (Nside = 256), and we achieve parsec-scale distance resolution, sampling the dust in 516 logarithmically spaced distance bins spanning 69 pc to 1250 pc. We generated 12 samples from the variational posterior of the 3D dust distribution and release the samples alongside the mean 3D dust map and its corresponding uncertainty. Our map resolves the internal structure of hundreds of molecular clouds in the solar neighborhood and will be broadly useful for studies of star formation, Galactic structure, and young stellar populations. It is available for download in a variety of coordinate systems online and can also be queried via the publicly available dustmaps Python package.
△ Less
Submitted 26 February, 2024; v1 submitted 2 August, 2023;
originally announced August 2023.
-
Reconstructing Galactic magnetic fields from local measurements for backtracking ultra-high-energy cosmic rays
Authors:
Alexandros Tsouros,
Gordian Edenhofer,
Torsten Enßlin,
Michalis Mastorakis,
Vasiliki Pavlidou
Abstract:
(abridged) Ultra-high energy cosmic rays (UHECRs) are highly energetic charged particles with energies exceeding $10^{18}$ eV. Identifying their sources and production mechanism can provide insight into many open questions in astrophysics and high energy physics. However, the Galactic magnetic field (GMF) deflects UHECRs, and the high uncertainties in our current understanding of the $3$-dimension…
▽ More
(abridged) Ultra-high energy cosmic rays (UHECRs) are highly energetic charged particles with energies exceeding $10^{18}$ eV. Identifying their sources and production mechanism can provide insight into many open questions in astrophysics and high energy physics. However, the Galactic magnetic field (GMF) deflects UHECRs, and the high uncertainties in our current understanding of the $3$-dimensional structure of the GMF does not permit us to accurately determine their true arrival direction on the plane of the sky (PoS). This difficulty arises from the fact that currently all GMF observations are integrated along the line-of-sight (LoS). Upcoming stellar optopolarimetric surveys as well as Gaia data on stellar parallaxes, are expected to provide local measurements of the GMF in the near future. In this paper, we evaluate the reconstruction of the GMF in a limited region of the Galaxy given sparse and local GMF measurements within that region, through Bayesian inference using principles of Information Field Theory. We backtrack UHECRs through GMF configurations drawn from the posterior to improve our knowledge of their true arrival directions. We show that, for a weakly turbulent GMF, it is possible to correct for its effect on the observed arrival direction of UHECRs to within $\sim 3^\circ$. For completely turbulent fields, we show that our procedure can still be used to significantly improve our knowledge on the true arrival direction of UHECRs.
△ Less
Submitted 2 October, 2023; v1 submitted 17 March, 2023;
originally announced March 2023.
-
The Musca molecular cloud: The perfect "filament" is still a sheet
Authors:
A. Tritsis,
F. Bouzelou,
R. Skalidis,
K. Tassis,
T. Enßlin,
G. Edenhofer
Abstract:
The true 3-dimensional (3D) morphology of the Musca molecular cloud is a topic that has received significant attention lately. Given that Musca does not exhibit intense star-formation activity, unveiling its shape has the potential of also revealing crucial information regarding the physics that dictates the formation of the first generation of stars within molecular clouds. Here, we revisit the s…
▽ More
The true 3-dimensional (3D) morphology of the Musca molecular cloud is a topic that has received significant attention lately. Given that Musca does not exhibit intense star-formation activity, unveiling its shape has the potential of also revealing crucial information regarding the physics that dictates the formation of the first generation of stars within molecular clouds. Here, we revisit the shape of Musca and we present a comprehensive array of evidence pointing towards a shape that is extended along the line-of-sight dimension: (a) 3D maps of differential extinction; (b) new non-local thermodynamic equilibrium radiative transfer simulations of CO rotational transitions from a sheet-like, magnetically-dominated simulated cloud; (c) an effective/critical density analysis of available CO observations; (d) indirect consequences that a filamentary structure would have had, from a theoretical star-formation perspective. We conclude that the full collection of observational evidence strongly suggests that Musca has a sheet-like geometry.
△ Less
Submitted 8 June, 2022;
originally announced June 2022.
-
The Galactic 3D large-scale dust distribution via Gaussian process regression on spherical coordinates
Authors:
R. H. Leike,
G. Edenhofer,
J. Knollmüller,
C. Alig,
P. Frank,
T. A. Enßlin
Abstract:
Knowing the Galactic 3D dust distribution is relevant for understanding many processes in the interstellar medium and for correcting many astronomical observations for dust absorption and emission. Here, we aim for a 3D reconstruction of the Galactic dust distribution with an increase in the number of meaningful resolution elements by orders of magnitude with respect to previous reconstructions, w…
▽ More
Knowing the Galactic 3D dust distribution is relevant for understanding many processes in the interstellar medium and for correcting many astronomical observations for dust absorption and emission. Here, we aim for a 3D reconstruction of the Galactic dust distribution with an increase in the number of meaningful resolution elements by orders of magnitude with respect to previous reconstructions, while taking advantage of the dust's spatial correlations to inform the dust map. We use iterative grid refinement to define a log-normal process in spherical coordinates. This log-normal process assumes a fixed correlation structure, which was inferred in an earlier reconstruction of Galactic dust. Our map is informed through 111 Million data points, combining data of PANSTARRS, 2MASS, Gaia DR2 and ALLWISE. The log-normal process is discretized to 122 Billion degrees of freedom, a factor of 400 more than our previous map. We derive the most probable posterior map and an uncertainty estimate using natural gradient descent and the Fisher-Laplace approximation. The dust reconstruction covers a quarter of the volume of our Galaxy, with a maximum coordinate distance of $16\,\text{kpc}$, and meaningful information can be found up to at distances of $4\,$kpc, still improving upon our earlier map by a factor of 5 in maximal distance, of $900$ in volume, and of about eighteen in angular grid resolution. Unfortunately, the maximum posterior approach chosen to make the reconstruction computational affordable introduces artifacts and reduces the accuracy of our uncertainty estimate. Despite of the apparent limitations of the presented 3D dust map, a good part of the reconstructed structures are confirmed by independent maser observations. Thus, the map is a step towards reliable 3D Galactic cartography and already can serve for a number of tasks, if used with care.
△ Less
Submitted 25 April, 2022;
originally announced April 2022.
-
On the Three-Dimensional Structure of Local Molecular Clouds
Authors:
Catherine Zucker,
Alyssa Goodman,
João Alves,
Shmuel Bialy,
Eric W. Koch,
Joshua S. Speagle,
Michael M. Foley,
Douglas Finkbeiner,
Reimar Leike,
Torsten Enßlin,
Joshua E. G. Peek,
Gordian Edenhofer
Abstract:
We leverage the 1 pc spatial resolution of the Leike et al. 2020 3D dust map to characterize the three-dimensional structure of nearby molecular clouds ($d \lesssim 400$ pc). We start by "skeletonizing" the clouds in 3D volume density space to determine their "spines," which we project on the sky to constrain cloud distances with $\approx 1\%$ uncertainty. For each cloud, we determine an average r…
▽ More
We leverage the 1 pc spatial resolution of the Leike et al. 2020 3D dust map to characterize the three-dimensional structure of nearby molecular clouds ($d \lesssim 400$ pc). We start by "skeletonizing" the clouds in 3D volume density space to determine their "spines," which we project on the sky to constrain cloud distances with $\approx 1\%$ uncertainty. For each cloud, we determine an average radial volume density profile around its 3D spine and fit the profiles using Gaussian and Plummer functions. The radial volume density profiles are well-described by a two-component Gaussian function, consistent with clouds having broad, lower-density outer envelopes and narrow, higher-density inner layers. The ratio of the outer to inner envelope widths is $\approx 3:1$. We hypothesize that these two components may be tracing a transition between atomic and diffuse molecular gas or between the unstable and cold neutral medium. Plummer-like models can also provide a good fit, with molecular clouds exhibiting shallow power-law wings with density, $n$, falling off like $n^{-2}$ at large radii. Using Bayesian model selection, we find that parameterizing the clouds' profiles using a single Gaussian is disfavored. We compare our results with 2D dust extinction maps, finding that the 3D dust recovers the total cloud mass from integrated approaches with fidelity, deviating only at higher levels of extinction ($A_V \gtrsim 2 - 3$ mag). The 3D cloud structure described here will enable comparisons with synthetic clouds generated in simulations, offering unprecedented insight into the origins and fates of molecular clouds in the interstellar medium.
△ Less
Submitted 27 September, 2021; v1 submitted 20 September, 2021;
originally announced September 2021.