Comparing simulated Milky Way satellite galaxies with observations using unsupervised clustering

We develop a new analysis method that allows us to compare multidimensional observables to a theoretical model. The method is based on unsupervised clustering algorithms which assign the observational and simulated data to clusters in high dimensionality. From the clustering result, a goodness of fi...

Full description

Saved in:
Bibliographic Details
Main Authors: Chen, Li-Hsin (Author) , Hartwig, Tilman (Author) , Klessen, Ralf S. (Author) , Glover, Simon (Author)
Format: Article (Journal)
Language:English
Published: 2022 November 10
In: Monthly notices of the Royal Astronomical Society
Year: 2022, Volume: 517, Issue: 4, Pages: 6140-6149
ISSN:1365-2966
DOI:10.1093/mnras/stac2897
Online Access:Verlag, Volltext: https://doi.org/10.1093/mnras/stac2897
Verlag, lizenzpflichtig, Volltext: https://academic.oup.com/mnras/article/517/4/6140/6754737?login=true
Get full text
Author Notes:Li-Hsin Chen, Tilman Hartwig, Ralf S. Klessen and Simon C.O. Glover
Description
Summary:We develop a new analysis method that allows us to compare multidimensional observables to a theoretical model. The method is based on unsupervised clustering algorithms which assign the observational and simulated data to clusters in high dimensionality. From the clustering result, a goodness of fit (the p-value) is determined with the Fisher-Freeman-Halton test. We first show that this approach is robust for 2D Gaussian distributions. We then apply the method to the observed MW satellites and simulated satellites from the fiducial model of our semi-analytic code a-sloth . We use the following five observables of the galaxies in the analysis: stellar mass, virial mass, heliocentric distance, mean stellar metallicity [Fe/H], and stellar metallicity dispersion σ[Fe/H]. A low p-value returned from the analysis tells us that our a-sloth fiducial model does not reproduce the mean stellar metallicity of the observed MW satellites well. We implement an ad hoc improvement to the physical model and show that the number of dark matter merger trees which have a p-value > 0.01 increases from 3 to 6. This method can be extended to data with higher dimensionality easily. We plan to further improve the physical model in a-sloth using this method to study elemental abundances of stars in the observed MW satellites.
Item Description:Gesehen am 19.01.2023
Physical Description:Online Resource
ISSN:1365-2966
DOI:10.1093/mnras/stac2897