GalaxyFlow: Upsampling Hydrodynamical Simulations for Realistic Mock Stellar Catalogs

Cosmological N-body simulations of galaxies operate at the level of "star particles" with a mass resolution on the scale of thousands of solar masses. Turning these simulations into stellar mock catalogs requires "upsampling" the star particles into individual stars following the same phase-space density. In this paper, we introduce two new upsampling methods. First, we describe GalaxyFlow, a sophisticated upsampling method that utilizes normalizing flows to both estimate the stellar phase space density and sample from it. Second, we improve on existing upsamplers based on adaptive kernel density estimation, using maximum likelihood estimation to fine-tune the bandwidth for such algorithms in a way that improves both the density estimation accuracy and upsampling results. We demonstrate our upsampling techniques on a neighborhood of the Solar location in two simulated galaxies: Auriga 6 and h277. Both yield smooth stellar distributions that closely resemble the stellar densities seen in the Gaia DR3 catalog. Furthermore, we introduce a novel multi-model classifier test to compare the accuracy of different upsampling methods quantitatively. This test confirms that GalaxyFlow estimates the density of the underlying star particles more accurately than methods based on kernel density estimation, at the cost of being more computationally intensive.