skip to main content

Title: Robust clustering of the local Milky Way stellar kinematic substructures with Gaia eDR3

Understanding local stellar kinematic substructures in the solar neighbourhood helps build a complete picture of the formation of the Milky Way, as well as an empirical phase space distribution of dark matter that would inform detection experiments. We apply the clustering algorithm hdbscan on the Gaia early third data release to identify a list of stable clusters in velocity space and action-angle space by taking into account the measurement uncertainties and studying the stability of the clustering results. We find 1405 (497) stars in 23 (6) robust clusters in velocity space (action-angle space) that are consistently not associated with noise. We discuss the kinematic properties of these structures and study whether many of the small clusters belong to a similar larger cluster based on their chemical abundances. They are attributed to the known structures: the Gaia Sausage-Enceladus, the Helmi Stream, and globular cluster NGC 3201 are found in both spaces, while NGC 104 and the thick disc (Sequoia) are identified in velocity space (action-angle space). Although we do not identify any new structures, we find that the hdbscan member selection of already known structures is unstable to input kinematics of the stars when resampled within their uncertainties. We therefore present the stable subset of local kinematic structures, which are consistently identified by the clustering algorithm, and emphasize the need to take into account error propagation during both the manual and automated identification of stellar structures, both for existing ones as well as future discoveries.

more » « less
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 2623-2648
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The Milky Way has accreted many ultra-faint dwarf galaxies (UFDs), and stars from these galaxies can be found throughout our Galaxy today. Studying these stars provides insight into galaxy formation and early chemical enrichment, but identifying them is difficult. Clustering stellar dynamics in 4D phase space (E,Lz,Jr,Jz) is one method of identifying accreted structure that is currently being utilized in the search for accreted UFDs. We produce 32 simulated stellar halos using particle tagging with the Caterpillar simulation suite and thoroughly test the abilities of different clustering algorithms to recover tidally disrupted UFD remnants. We perform over 10,000 clustering runs, testing seven clustering algorithms, roughly twenty hyperparameter choices per algorithm, and six different types of data sets each with up to 32 simulated samples. Of the seven algorithms, HDBSCAN most consistently balances UFD recovery rates and cluster realness rates. We find that, even in highly idealized cases, the vast majority of clusters found by clustering algorithms do not correspond to real accreted UFD remnants and we can generally only recover 6% of UFDs remnants at best. These results focus exclusively on groups of stars from UFDs, which have weak dynamic signatures compared to the background of other stars. The recoverable UFD remnants are those that accreted recently,zaccretion≲ 0.5. Based on these results, we make recommendations to help guide the search for dynamically linked clusters of UFD stars in observational data. We find that real clusters generally have higher median energy andJr, providing a way to help identify real versus fake clusters. We also recommend incorporating chemical tagging as a way to improve clustering results.

    more » « less
  2. Abstract A number of studies based on the data collected by the Hubble Space Telescope (HST) GO-13297 program “HST Legacy Survey of Galactic Globular Clusters: Shedding UV Light on Their Populations and Formation” have investigated the photometric properties of a large sample of Galactic globular clusters and revolutionized our understanding of their stellar populations. In this paper, we expand upon previous studies by focusing our attention on the stellar clusters’ internal kinematics. We computed proper motions for stars in 56 globular clusters and one open cluster by combining the GO-13297 images with archival HST data. The astrophotometric catalogs released with this paper represent the most complete and homogeneous collection of proper motions of stars in the cores of stellar clusters to date, and expand the information provided by the current (and future) Gaia data releases to much fainter stars and into the crowded central regions. We also census the general kinematic properties of stellar clusters by computing the velocity dispersion and anisotropy radial profiles of their bright members. We study the dependence on concentration and relaxation time, and derive dynamical distances. Finally, we present an in-depth kinematic analysis of the globular cluster NGC 5904. 
    more » « less

    As an introduction of a kinematic survey of Magellanic Cloud (MC) star clusters, we report on the dynamical masses and mass-to-light ratios (M/L) of NGC 419 (Small Magellanic Cloud) and NGC 1846 (Large Magellanic Cloud). We have obtained more than one hundred high-resolution stellar spectra in and around each cluster using the multi-object spectrograph M2FS on the Magellan/Clay Telescope. Line-of-sight velocities and positions of the stars observed in each cluster were used as input to an expectation-maximization algorithm used to estimate cluster membership probabilities, resulting in samples of 46 and 52 likely members (PM ≥ 50 per cent) in NGC 419 and NGC 1846, respectively. This process employed single-mass King models constrained by the structural parameters of the clusters and provided self-consistent dynamical mass estimates for both clusters. Our best-fitting results show that NGC 419 has a projected central velocity dispersion of $2.44^{+0.37}_{-0.21}$ km s−1, corresponding to a total mass of $7.6^{+2.5}_{-1.3}\times 10^4\ {\rm M}_{\odot }$ and V-band M/L ratio of $0.22^{+0.08}_{-0.05}$ in solar units. For NGC 1846, the corresponding results are $2.04^{+0.28}_{-0.24}$ km s−1, $5.4^{+1.5}_{-1.4}\times 10^4\ {\rm M}_{\odot }$, and $0.32^{+0.11}_{-0.11}$. The mean metallicities of NGC 419 and NGC 1846 are found to be $\rm [Fe/H]=-0.84\pm 0.19$ and −0.70 ± 0.08, respectively, based on the spectra of likely cluster members. We find marginal statistical evidence of rotation in both clusters, though in neither cluster does rotation alter our mass estimates significantly. We critically compare our findings with those of previous kinematic studies of these two clusters in order to evaluate the consistency of our observational results and analytic tools.

    more » « less
  4. Abstract

    Orbital characteristics based on Gaia Early Data Release 3 astrometric parameters are analyzed for ∼4000 metal-poor stars ([Fe/H] ≤ −0.8) compiled from the Best and Brightest survey. Selected as metal-poor candidates based on broadband near- and far-IR photometry, 43% of these stars had medium-resolution (1200 ≲R≲ 2000) validation spectra obtained over a 7 yr campaign from 2014 to 2020 with a variety of telescopes. The remaining stars were chosen based on photometric metallicity determinations from the Huang et al. recalibration of the Sky Mapper Southern Survey. Dynamical clusters of these stars are obtained from the orbital energy and cylindrical actions using theHDBSCANunsupervised learning algorithm. We identify 52 dynamically tagged groups (DTGs) with between five and 21 members; 18 DTGs have at least 10 member stars. Milky Way (MW) substructures such as Gaia-Sausage-Enceladus, the Metal-Weak Thick-Disk, Thamnos, the Splashed Disk, and the Helmi Stream are identified. Associations with MW globular clusters are determined for eight DTGs; no recognized MW dwarf galaxies were associated with any of our DTGs. Previously identified dynamical groups are also associated with our DTGs, with emphasis placed on their structural determination and possible new identifications. Chemically peculiar stars are identified as members of several DTGs, with six DTGs that are associated withr-process-enhanced stars. We demonstrate that the mean carbon andα-element abundances of our DTGs are correlated with their mean metallicity in an understandable manner. Similarly, we find that the mean metallicity, carbon, andα-element abundances are separable into different regions of the mean rotational-velocity space.

    more » « less
  5. Abstract Gaia DR2 has provided an unprecedented wealth of information about the positions and motions of stars in our Galaxy, and has highlighted the degree of disequilibria in the disc. As we collect data over a wider area of the disc it becomes increasingly appealing to start analysing stellar actions and angles, which specifically label orbit space, instead of their current phase space location. Conceptually, while $\bar{x}$ and $\bar{v}$ tell us about the potential and local interactions, grouping in action puts together stars that have similar frequencies and hence similar responses to dynamical effects occurring over several orbits. Grouping in actions and angles refines this further to isolate stars which are travelling together through space and hence have shared histories. Mixing these coordinate systems can confuse the interpretation. For example, it has been suggested that by moving stars to their guiding radius, the Milky Way spiral structure is visible as ridge-like overdensities in the Gaia data (Khoperskov et al. 2020). However, in this work, we show that these features are in fact the known kinematic moving groups, both in the Lz − φ and the vR − vφ planes. Using simulations we show how this distinction will become even more important as we move to a global view of the Milky Way. As an example, we show that the radial velocity wave seen in the Galactic disc in Gaia and APOGEE should become stronger in the action-angle frame, and that it can be reproduced by transient spiral structure. 
    more » « less