We search for dynamical substructures in the LAMOST DR3 very metal-poor (VMP) star catalog. After cross-matching with Gaia DR2, there are ∼ 3300 VMP stars with available high-quality astrometric information that have halo-like kinematics. We apply a method based on self-organizing maps, StarGO, to find groups clustered in the 4D space of orbital energy and angular momentum. We identify 57 dynamically tagged groups, which we label DTG-1 to DTG-57. Most of them belong to existing massive substructures in the nearby halo, such as the Gaia Sausage or Sequoia. The stream identified by Helmi et al. is recovered, but the two disjoint portions of the substructure appear to have distinct dynamical properties. The very retrograde substructure Rg5 found previously by Myeong et al. is also retrieved. We report six new DTGs with highly retrograde orbits, two with very prograde orbits, and 12 with polar orbits. By mapping other datasets (APOGEE halo stars, and catalogs of r-process-enhanced and CEMP stars) onto the trained neuron map, we can associate stars with detailed chemical abundances to the DTGs, and look for associations with chemically peculiar stars. The highly eccentric Gaia Sausage groups contain representatives both of debris from the satellite itself (which is α-poor) and the Splashed Disk, sent up into eccentric halo orbits from the encounter (and is α-rich). The new prograde substructures also appear to be associated with the Splashed Disk. The DTGs belonging to the Gaia Sausage host two relatively metal-rich r-II stars and six CEMP stars in different sub-classes, consistent with the idea that the Gaia Sausage progenitor is a massive dwarf galaxy. Rg5 is dynamically associated with two highly r-process-enhanced stars with [Fe/H] ∼ −3. This finding indicates that its progenitor might be an ultra-faint dwarf galaxy that has experienced r-process enrichment from neutron star mergers.

We derive dynamical parameters for a large sample of 446 r-process-enhanced (RPE) metal-poor stars in the halo and disk systems of the Milky Way, based on data releases from the R-Process Alliance, supplemented by additional literature samples. This sample represents more than a 10-fold increase in size relative to that previously considered by Roederer et al. and, by design, covers a larger range of r-process-element enrichment levels. We test a number of clustering analysis methods on the derived orbital energies and other dynamical parameters for this sample, ultimately deciding on application of the HDBSCAN algorithm, which obtains 30 individual chemodynamically tagged groups (CDTGs); 21 contain between 3 and 5 stars, and 9 contain between 6 and 12 stars. Even though the clustering was performed solely on the basis of their dynamical properties, the stars in these CDTGs exhibit statistically significant similarities in their metallicity ([Fe/H]), carbonicity ([C/Fe]), and neutron-capture element ratios ([Sr/Fe], [Ba/Fe], and [Eu/Fe]). These results demonstrate that the RPE stars in these CDTGs have likely experienced common chemical-evolution histories, presumably in their parent satellite galaxies or globular clusters, prior to being disrupted into the Milky Way’s halo. We also confirm the previous claim that the orbits of the RPE stars preferentially exhibit pericentric distances that are substantially lower than the present distances of surviving ultrafaint dwarf and canonical dwarf spheroidal galaxies, consistent with the disruption hypothesis. The derived dynamical parameters for several of our CDTGs indicate their association with previously known substructures, dynamically tagged groups, and RPE groups.