From starlight to dark matter: a stochastic interpolation approach to map dark matter from stellar density

Abstract

The dark matter halo profile in galaxies holds key information about the nature of dark matter and galaxy formation. Constraining the dark matter profile of galaxies beyond the Milky Way traditionally requires expensive spectroscopic observations for kinematic information. In this paper, we explore a conditional generative model framework to map the dark matter profile of Milky Way-mass galaxies from stellar density profiles, more easily obtainable through large photometric imaging surveys. As a proof of concept, we train the model to learn a stochastic bridge between instrument systematics-free baryonic stellar distributions and underlying dark matter density maps from the DREAMS hydrodynamics simulation suite. We recover 2D dark matter density profiles with a typical accuracy of $\sim0.1$ dex ($\sim1.5$% of the truth in log scale). The stochastic sampling procedure provides uncertainty estimates of the predicted dark matter map, with typical values $\sim0.1$ dex. Out-of-domain tests with Milky Way-mass galaxies from IllustrisTNG and FIRE simulations show that, while the model can qualitatively be generalized to TNG50 galaxies from IllustrisTNG, the model is sensitive to the galaxy formation model, with $\sim0.2$-$0.4$ dex over-prediction for the inner profiles ($r\lesssim5$ kpc) of the FIRE test galaxies. Future work will explore training with additional suites of simulations and/or conditioning on additional information, such as multi-band images. Our results are a first step towards using generative models as a flexible, uncertainty-aware framework for turning forthcoming data from large imaging surveys into spatially resolved dark matter maps.