Stellar ejecta gradually enrich the gas out of which subsequent stars form, making the least chemically enriched stellar systems direct fossils of structures formed in the early Universe1. Although a few hundred stars with metal content below 1,000th of the solar iron content are known in the Galaxy2–4, none of them inhabit globular clusters, some of the oldest known stellar structures. These show metal content of at least approximately 0.2% of the solar metallicity ([Fe / H] ≳ − 2.7). This metallicity floor appears universal5,6, and it has been proposed that protogalaxies that merged into the galaxies we observe today were simply not massive enough to form clusters that survived to the present day7. Here we report observations of a stellar stream, C-19, whose metallicity is less than 0.05% of the solar metallicity ([Fe/H]=−3.38±0.06(statistical)±0.20(systematic)). The low metallicity dispersion and the chemical abundances of the C-19 stars show that this stream is the tidal remnant of the most metal-poor globular cluster ever discovered, and is significantly below the purported metallicity floor: clusters with significantly lower metallicities than observed today existed in the past and contributed their stars to the Milky Way halo. 

The density of dark matter near the Sun, ρ_{DM, ⊙}, is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. We assess the impact of these assumptions on determinations of ρ_{DM, ⊙} by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia–Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρ_{DM, ⊙} from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ∼1.5 in the recovered surface mass density ∼2 kpc above the disc plane, implying a fractional error on ρ_{DM, ⊙} of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρ_{DM, ⊙}, consistent with the true value within our 95 per cent confidence intervals for realistic 20 per cent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 per cent to ρ_{DM, ⊙} (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.

We report the first evidence of an ongoing collision between two star clusters in our Galaxy, namely IC 4665 and Collinder 350. These are open clusters located at a distance of ∼330 pc from the Sun and ∼100 pc above the Galactic plane, and they both have prograde motions with only a small difference in their velocities (Collinder 350 moves ∼5kms−1∼5kms−1 faster than IC 4665); as inferred from ESA/Gaia based catalogue. Interestingly, the two clusters are physically separated by only ∼36 pc in space; a distance that is smaller than the sum of their respective radii. Furthermore, the clusters exhibit signatures of elongated stellar density distributions, and we also detect an onset of an inter-cluster stellar bridge. Moreover, the orbit analysis suggests that the younger cluster IC 4665 (age = 53 Myr) must have formed at a distance >500 pc away from Collinder 350 (age = 617 Myr). These findings together imply that the two clusters do not represent merging of two objects in a binary system; rather, what we are witnessing is an actual collision between two independently formed star clusters. This collision phenomenon provides a unique opportunity to explore new aspects of formation and evolution theory of star clusters.

The central density profiles in dwarf galaxy halos depend strongly on the nature of dark matter (DM). Recently, in Malhan et al. we employed N-body simulations to show that the cuspy cold DM subhalos predicted by cosmological simulations can be differentiated from cored subhalos using the properties of accreted globular cluster (GC) streams since these GCs experience tidal stripping within their parent halos prior to accretion onto the Milky Way. We previously found that clusters that are accreted within cuspy subhalos produce streams with larger physical widths and higher dispersions in line-of-sight velocity and angular momentum than streams that are accreted within cored subhalos. Here, we use the same suite of simulations to demonstrate that the dispersion in the tangential velocities of streams () is also sensitive to the central DM density profiles of their parent dwarfs and GCs that they were accreted from; cuspy subhalos produce streams with larger  than those accreted inside cored subhalos. Using Gaia EDR3 observations of multiple GC streams we compare their  values with simulations. The measured  values are consistent with both an "in situ" origin and with accretion inside cored subhalos of M ∼ 10^8–9 M_⊙ (or very low-mass cuspy subhalos of mass ∼10⁸ M_⊙). Despite the large current uncertainties in , we find a low probability that any of the progenitor GCs were accreted from cuspy subhalos of M ≳ 10 ⁹M_⊙. The uncertainties on Gaia tangential velocity measurements are expected to decrease in future and will allow for stronger constraints on subhalo DM density profiles.

We present a new spectroscopic study of the dwarf galaxy Boötes I (Boo I) with data from the Anglo-Australian Telescope and its AAOmega spectrograph together with the Two Degree Field multi-object system. We observed 36 high-probability Boo I stars selected using Gaia Early Data Release 3 proper motions and photometric metallicities from the Pristine survey. Out of those, 27 are found to be Boo I stars, resulting in an excellent success rate of 75 per cent at finding new members. Our analysis uses a new pipeline developed to estimate radial velocities and equivalent widths of the calcium triplet lines from Gaussian and Voigt line profile fits. The metallicities of 16 members are derived, including 3 extremely metal-poor stars ([Fe/H] < −3.0), which translates into a success rate of 25 per cent at finding them with the combination of Pristine and Gaia. Using the large spatial extent of our new members that spans up to 4.1 half-light radii and spectroscopy from the literature, we find a systemic velocity gradient of 0.40 ± 0.10 km s⁻¹ arcmin⁻¹ and a small but resolved metallicity gradient of −0.008 ± 0.003 dex arcmin⁻¹. Finally, we show that Boo I is more elongated than previously thought with an ellipticity of ϵ = 0.68 ± 0.15. Its velocity and metallicity gradients as well as its elongation suggest that Boo I may have been affected by tides, a result supported by direct dynamical modelling.

We use the photometric metallicities provided by the panoramic Pristine survey to study the veracity and derive the metallicities of the numerous stellar streams found by the application of the STREAMFINDER algorithm to the Gaia Early Data Release 3 data. All 26 streams present in Pristine show a clear metallicity distribution function, which provides an independent check of the reality of these structures, supporting the reliability of STREAMFINDER in finding streams and the power of Pristine to measure precise metallicities. We further present six candidate structures with coherent phase-space and metallicity signals that are very likely streams. The majority of studied streams are very metal-poor (14 structures with [Fe/H] < −2.0) and include three systems with [Fe/H] < −2.9 (C-11, C-19, and C-20). These streams could be the closest debris of low-luminosity dwarf galaxies or may have originated from globular clusters of significantly lower metallicity than any known current Milky Way globular cluster. Our study shows that the promise of the Gaia data for Galactic Archeology studies can be substantially strengthened by quality photometric metallicities, allowing us to peer back into the earliest epochs of the formation of our Galaxy and its stellar halo constituents.

We present maps of the stellar streams detected in the Gaia Data Release 2 (DR2) and Early Data Release 3 (EDR3) catalogs using the STREAMFINDER algorithm. We also report the spectroscopic follow-up of the brighter DR2 stream members obtained with the high-resolution CFHT/ESPaDOnS and VLT/UVES spectrographs as well as with the medium-resolution NTT/EFOSC2 spectrograph. Two new stellar streams that do not have a clear progenitor are detected in DR2 (named Hrid and Gunnthra), and seven are detected in EDR3 (named Gaia-6 to Gaia-12). Several candidate streams are also identified. The software also finds very long tidal tails associated with the 15 globular clusters: NGC 288, NGC 1261, NGC 1851, NGC 2298, NGC 2808, NGC 3201, M68, omega Cen, NGC 5466, Palomar 5, M5, NGC 6101, M92, NGC 6397, and NGC 7089. These stellar streams will be used in subsequent contributions in this series to chart the properties of the Galactic acceleration field on similar to 100 pc to similar to 100 kpc scales.

Stellar streams produced from dwarf galaxies provide direct evidence of the hierarchical formation of the Milky Way. Here, we present the first comprehensive study of the LMS-1 stellar stream, that we detect by searching for wide streams in the Gaia EDR3 data set using the STREAMFINDER algorithm. This stream was recently discovered by Yuan et al. We detect LMS-1 as a 60° long stream to the north of the galactic bulge, at a distance of ∼20 kpc from the Sun, together with additional components that suggest that the overall stream is completely wrapped around the inner Galaxy. Using spectroscopic measurements from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, the Sloan Digital Sky Survey, and the Apache Point Observatory Galactic Evolution Experiment, we infer that the stream is very metal-poor (〈[Fe/H]〉 = −2.1) with a significant metallicity dispersion (σ_[Fe/H] = 0.4), and it possesses a large radial velocity dispersion (σ_v = 20 ± 4 km s⁻¹). These estimates together imply that LMS-1 is a dwarf galaxy stream. The orbit of LMS-1 is close to polar, with an inclination of 75° to the galactic plane. Both the orbit and metallicity of LMS-1 are remarkably similar to the globular clusters NGC 5053, NGC 5024, and the stellar stream Indus. These findings make LMS-1 an important contributor to the stellar population of the inner Milky Way halo.

The steepness of the central density profiles of dark matter (DM) in low-mass galaxy haloes (e.g. dwarf galaxies) is a powerful probe of the nature of DM. We propose a novel scheme to probe the inner profiles of galaxy subhaloes using stellar streams. We show that the present-day morphological and dynamical properties of accreted globular cluster (GC) streams - those produced from tidal stripping of GCs that initially evolved within satellite galaxies and later merged with the Milky Way (MW) - are sensitive to the central DM density profile and mass of their parent satellites. GCs that accrete within cuspy cold dark matter (CDM) subhaloes produce streams that are physically wider and dynamically hotter than streams that accrete inside cored subhaloes. A first comparison of MW streams 'GD-1' and 'Jhelum' (likely of accreted GC origin) with our simulations indicates a preference for cored subhaloes. If these results hold up in future data, the implication is that either the DM cusps were erased by baryonic feedback, or their subhaloes naturally possessed cored density profiles implying particle physics models beyond CDM. Moreover, accreted GC streams are highly structured and exhibit complex morphological features (e.g. parallel structures and 'spurs'). This implies that the accretion scenario can naturally explain the recently observed peculiarities in some of the MW streams. We also propose a novel mechanism for forming 'gaps' in stellar streams when the remnant of the parent subhalo (which hosted the GC) later passes through the GC stream. This encounter can last a longer time (and have more of an impact) than the random encounters with DM subhaloes previously considered, because the GC stream and its parent subhalo are on similar orbits with small relative velocities. Current and future surveys of the MW halo will uncover numerous faint stellar streams and provide the data needed to substantiate our preliminary tests with this new probe of DM.

We present a new spectroscopic study of the faint Milky Way satellite Sagittarius II. Using multiobject spectroscopy from the Fibre Large Array Multi-Element Spectrograph, we supplement the data set of Longeard et al. with 47 newly observed stars, 19 of which are identified as members of the satellite. These additional member stars are used to put tighter constraints on the dynamics and the metallicity properties of the system. We find a low velocity dispersion of sigma(SgrII)(v) = 1.7 +/- 0.5 km s(-1), in agreement with the dispersion of Milky Way globular clusters of similar luminosity. We confirm the very metal-poor nature of the satellite ([Fe/H](spectro)(SgrII) = -2.23 +/- 0.07) and find that the metallicity dispersion of Sgr II is not resolved, reaching only 0.20 at the 95 per cent confidence limit. No star with a metallicity below -2.5 is confidently detected. Therefore, despite the unusually large size of the system (r(h) = 35.5(-1.2)(-1.4) pc), we conclude that Sgr II is an old and metal-poor globular cluster of the Milky Way.