Data from the Gaia satellite show that the solar neighbourhood of the Milky Way's stellar halo is imprinted with substructure from several accretion events. Evidence of these events is found in "the Shards", stars clustering with high significance in both action space and metallicity. Stars in the Shards share a common origin, likely as ancient satellite galaxies of the Milky Way, so will be embedded in dark matter (DM) counterparts. These "Dark Shards" contain two substantial streams (S1 and S2), as well as several retrograde, prograde and lower energy substructures. The retrograde stream S1 has a very high Earth-frame speed of $\sim 550$ km s$^{-1}$ while S2 moves on a prograde, but highly polar orbit and enhances the peak of the speed distribution at around $300$ km s$^{-1}$. The presence of the Dark Shards locally leads to modifications of many to the fundamental properties of experimental DM signals. The~S2 stream in particular gives rise to an array of effects in searches for axions and in the time dependence of nuclear recoils: shifting the peak day, inducing nonsinusoidal distortions, and increasing the importance of the gravitational focusing of DM by the Sun. Dark Shards additionally bring new features for directional signals, while also enhancing the DM flux towards Cygnus.

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 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 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 have distinct dynamical properties. The very retrograde substructure Rg5 found previously by Myeong et al. is also retrieved. We report 6 new DTGs with highly retrograde orbits, 2 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 $\alpha$-poor) and the Splashed Disk, sent up into eccentric halo orbits from the encounter (and is $\alpha$-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] $\sim -$3. This finding indicates that its progenitor might be an ultra-faint dwarf galaxy that has experienced $r$-process enrichment from neutron star mergers.

The Gaia Sausage is the major accretion event that built the stellar halo of the Milky Way galaxy. Here, we provide dynamical and chemical evidence for a second substantial accretion episode, distinct from the Gaia Sausage. The Sequoia Event provided the bulk of the high energy retrograde stars in the stellar halo, as well as the recently discovered globular cluster FSR 1758. There are up to 6 further globular clusters, including $\omega$~Centauri, as well as many of the retrograde substructures in Myeong et al. (2018), associated with the progenitor dwarf galaxy, named the Sequoia. The stellar mass in the Sequoia galaxy is $\sim 5 \times 10^{7} M_\odot$, whilst the total mass is $\sim 10^{10} M_\odot$, as judged from abundance matching or from the total sum of the globular cluster mass. Although clearly less massive than the Sausage, the Sequoia has a distinct chemo-dynamical signature. The strongly retrograde Sequoia stars have a typical eccentricity of $\sim0.6$, whereas the Sausage stars have no clear net rotation and move on predominantly radial orbits. On average, the Sequoia stars have lower metallicity by $\sim 0.3$ dex and higher abundance ratios as compared to the Sausage. We conjecture that the Sausage and the Sequoia galaxies may have been associated and accreted at a comparable epoch.

The recently discovered S1 stream passes through the Solar neighbourhood on a low inclination, counter-rotating orbit. The progenitor of S1 is a dwarf galaxy with a total mass comparable to the present-day Fornax dwarf spheroidal, so the stream is expected to have a significant DM component. We compute the effects of the S1 stream on WIMP and axion detectors as a function of the density of its unmeasured dark component. In WIMP detectors the S1 stream supplies more high energy nuclear recoils so will marginally improve DM detection prospects. We find that even if S1 comprises less than 10% of the local density, multi-ton xenon WIMP detectors can distinguish the S1 stream from the bulk halo in the relatively narrow mass range between 5 and 25 GeV. In directional WIMP detectors such as CYGNUS, S1 increases DM detection prospects more substantially since it enhances the anisotropy of the WIMP signal. Finally, we show that axion haloscopes possess by far the greatest potential sensitivity to the S1 stream. Once the axion mass has been discovered, the distinctive velocity distribution of S1 can easily be extracted from the axion power spectrum.

We use the SDSS-Gaia catalogue to search for substructure in the stellar halo. The sample comprises 62\,133 halo stars with full phase space coordinates and extends out to heliocentric distances of $\sim 10$ kpc. As actions are conserved under slow changes of the potential, they permit identification of groups of stars with a common accretion history. We devise a method to identify halo substructures based on their clustering in action space, using metallicity as a secondary check. This is validated against smooth models and numerical constructed stellar halos from the Aquarius simulations. We identify 21 substructures in the SDSS-Gaia catalogue, including 7 high significance, high energy and retrograde ones. We investigate whether the retrograde substructures may be material stripped off the atypical globular cluster $\omega$~Centauri. Using a simple model of the accretion of the progenitor of the $\omega$~Centauri, we tentatively argue for the possible association of up to 5 of our new substructures (labelled Rg1, Rg3, Rg4, Rg6 and Rg7) with this event. This sets a minimum mass of $5 \times 10^8 M_\odot$ for the progenitor, so as to bring $\omega$~Centauri to its current location in action -- energy space. Our proposal can be tested by high resolution spectroscopy of the candidates to look for the unusual abundance patterns possessed by $\omega$~Centauri stars.

The Gaia Sausage is an elongated structure in velocity space discovered by Belokurov et al. (2018) using the kinematics of metal-rich halo stars. It was created by a massive dwarf galaxy ($\sim 5 \times 10^{10} M_\odot$) on a strongly radial orbit that merged with the Milky Way at a redshift $z\lesssim 3$. We search forthe associated Sausage Globular Clusters by analysing the structure of 91 Milky Way globular clusters (GCs) in action space using the Gaia Data Release 2 catalogue, complemented with Hubble Space Telescope proper motions. There is a characteristic energy $E_{\rm crit}$ which separates the in situ objects, such as the bulge/disc clusters, from the accreted objects, such as the young halo clusters. There are 15 old halo GCs that have $E > E_{\rm crit}$. Eight of the high energy, old halo GCs are strongly clumped in azimuthal and vertical action, yet strung out like beads on a chain at extreme radial action. They are very radially anisotropic ($\beta \sim 0.95$) and move on orbits that are all highly eccentric ($e \gtrsim 0.80$). They also form a track in the age-metallicity plane distinct from the bulk of the Milky Way GCs and compatible with a dwarf spheroidal origin. These properties are consistent with GCs associated with the merger event that gave rise to the Gaia Sausage.

We analyse the structure of the local stellar halo of the Milky Way using $\sim$ 60000 stars with full phase space coordinates extracted from the SDSS--{\it Gaia} catalogue. We display stars in action space as a function of metallicity in a realistic axisymmetric potential for the Milky Way Galaxy. The metal-rich population is more distended towards high radial action $J_R$ as compared to azimuthal or vertical action, $J_\phi$ or $J_z$. It has a mild prograde rotation $(\langle v_\phi \rangle \approx 25$ km s$^{-1}$), is radially anisotropic and highly flattened with axis ratio $q \approx 0.6 - 0.7$. The metal-poor population is more evenly distributed in all three actions. It has larger prograde rotation $(\langle v_\phi \rangle \approx 50$ km s$^{-1}$), a mild radial anisotropy and a roundish morphology ($q\approx 0.9$). We identify two further components of the halo in action space. There is a high energy, retrograde component that is only present in the metal-rich stars. This is suggestive of an origin in a retrograde encounter, possibly the one that created the stripped dwarf galaxy nucleus, $\omega$Centauri. Also visible as a distinct entity in action space is a resonant component, which is flattened and prograde. It extends over a range of metallicities down to [Fe/H] $\approx -3$. It has a net outward radial velocity $\langle v_R \rangle \approx 12$ km s$^{-1}$ within the Solar circle at $|z| <3.5$ kpc. The existence of resonant stars at such extremely low metallicities has not been seen before.

We use the SDSS-Gaia Catalogue to identify six new pieces of halo substructure. SDSS-Gaia is an astrometric catalogue that exploits SDSS data release 9 to provide first epoch photometry for objects in the Gaia source catalogue. We use a version of the catalogue containing $245\,316$ stars with all phase space coordinates within a heliocentric distance of $\sim 10$ kpc. We devise a method to assess the significance of halo substructures based on their clustering in velocity space. The two most substantial structures are multiple wraps of a stream which has undergone considerable phase mixing (S1, with 94 members) and a kinematically cold stream (S2, with 61 members). The member stars of S1 have a median position of ($X,Y,Z$) = ($8.12, -0.22, 2.75$) kpc and a median metallicity of [Fe/H] $= -1.78$. The stars of S2 have median coordinates ($X,Y,Z$) = ($8.66, 0.30, 0.77$) kpc and a median metallicity of [Fe/H] $= -1.91$. They lie in velocity space close to some of the stars in the stream reported by Helmi et al. (1999). By modelling, we estimate that both structures had progenitors with virial masses $\approx 10^{10} M_\odot$ and infall times $\gtrsim 9$ Gyr ago. Using abundance matching, these correspond to stellar masses between $10^6$ and $10^7 M_\odot$. These are somewhat larger than the masses inferred through the mass-metallicity relation by factors of 5 to 15. Additionally, we identify two further substructures (S3 and S4 with 55 and 40 members) and two clusters or moving groups (C1 and C2 with 24 and 12) members. In all 6 cases, clustering in kinematics is found to correspond to clustering in both configuration space and metallicity, adding credence to the reliability of our detections.

We identify a halo substructure in the Tycho Gaia Astrometric Solution (TGAS) dataset, cross-matched with the RAVE-on data release. After quality cuts, the stars with large radial action ($J_R > 800$ kms$^{-1}$ kpc) are extracted. A subset of these stars is clustered in longitude and velocity and can be selected with further cuts. The 14 stars are centered on $(X,Y,Z) \approx (9.0,-1.0,-0.6)$ kpc and form a coherently moving structure in the halo with median $(v_R,v_\phi,v_z) = (167.33,0.86,-94.85)$ kms$^{-1}$. They are all metal-poor giants with median [Fe/H] $=-0.83$. To guard against the effects of distance errors, we compute spectrophotometric distances for the 8 out of 14 stars where this is possible. We find that 6 of the stars are still comoving. These 6 stars also have a much tighter [Fe/H] distribution $\sim -0.7$ with one exception ([Fe/H] = -2.12). We conclude that the existence of the comoving cluster is stable against changes in distance estimation and conjecture that this is the dissolving remnant of a long-ago accreted globular cluster.

We report the discovery of tidal tails around the two outer halo globular clusters, Eridanus and Palomar 15, based on $gi$-band images obtained with DECam at the CTIO 4-m Blanco Telescope. The tidal tails are among the most remote stellar streams presently known in the Milky Way halo. Cluster members have been determined from the color-magnitude diagrams and used to establish the radial density profiles, which show, in both cases, a strong departure in the outer regions from the best-fit King profile. Spatial density maps reveal tidal tails stretching out on opposite sides of both clusters, extending over a length of $\sim$760 pc for Eridanus and $\sim$1160 pc for Palomar 15. The great circle projected from the Palomar 15 tidal tails encompasses the Galactic Center, while that for Eridanus passes close to four dwarf satellite galaxies, one of which (Sculptor) is at a comparable distance to that of Eridanus.