Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighboring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational potential to multiple stellar streams simultaneously by maximizing the clustering of the stream stars in action space. We apply this technique to members of the GD-1, Pal 5, Orphan and Helmi streams, exploiting both the individual and combined data sets. We describe the Galactic potential with a St\"ackel model, and vary up to five parameters simultaneously. We find that we can only constrain the enclosed mass, and that the strongest constraints come from the GD-1 and Pal 5 streams whose combined data set yields $M(< 20\ \mathrm{kpc}) = 3.47^{+0.95}_{-1.44} \times 10^{11} \ M_{\odot}$. When including the Orphan and Helmi stream in the data set, the mass uncertainty increases to $M(< 20\ \mathrm{kpc}) = 3.12^{+5.69}_{-1.07} \times 10^{11} \ M_{\odot}$, indicative of the hidden systematic errors in fits to individual streams. These systematics are likely due to insufficiencies in our St\"ackel model of the potential and to the limited phase space explored by the data set of 4 streams, so the larger uncertainty of the combined fit is the more robust measure of the actual uncertainty in the Milky Way potential.

The Gaia mission has provided the largest catalogue ever of sources with tangential velocity information. Using this catalogue for dynamical studies is, however, difficult because most of the stars lack line-of-sight velocity measurements. Recently, we presented a selection of $\sim 10^7$ halo stars with accurate distances that have been selected based on their photometry and proper motions. Using this sample, we model the tail of the velocity distribution in the stellar halo, locally and as a function of distance. Our goal is to measure the escape velocity, and herewith to constrain the mass of our Galaxy. We fit the tail of the velocity distribution with a power-law distribution, a commonly used approach first established by Leonard & Tremaine (1990). For the first time ever we use tangential velocities measured accurately for an unprecedented number of halo stars to estimate the escape velocity. For the solar neighbourhood we determine a very precise lower limit to the escape velocity: $497^{+40}_{-24}~{\rm km/s}$. Our best guess is that this value is underestimated by $\sim 10\%$. For a spherical NFW halo in a Milky Way potential, this translates into a lower limit on the mass of $M_{200} = 0.67^{+0.30}_{-0.15}\cdot10^{12}~{\rm M}_\odot$ and concentration parameter $c = 15.0^{2.6}_{-2.3}$. The behaviour of the escape velocity towards the inner Galaxy follows expectations from Milky Way dynamical models, but curiously we find that its value appears to increase with distance beyond the solar radius. This behaviour is due to a change in the shape of the velocity distribution, and could be related to the presence of velocity clumps. A tentative analysis of the escape velocity as a function of $(R,z)$ shows that it decreases too slowly compared to what is expected for a spherical halo using standard values for the characteristic parameters describing the galactic disc.

Circa 10 billion years ago the Milky Way merged with a massive satellite, Gaia-Enceladus. To gain insight into the properties of its debris we analyse in detail the suite of simulations from Villalobos & Helmi (2008), which includes an experiment that produces a good match to the kinematics of nearby halo stars inferred from Gaia data. We compare the kinematic distributions of stellar particles in the simulations and study the distribution of debris in orbital angular momentum, eccentricity and energy, and its relation to the mass-loss history of the simulated satellite. We confirm that Gaia-Enceladus probably fell in on a retrograde, 30$^\circ$ inclination orbit. We find that while 75% of the debris in our preferred simulation has large eccentricity ($> 0.8$), roughly 9% has eccentricity smaller than 0.6. Star particles lost early have large retrograde motions, and a subset of these have low eccentricity. Such stars would be expected to have lower metallicities as they stem from the outskirts of the satellite, and hence naively they could be confused with debris associated with a separate system. These considerations seem to apply to some of the stars from the postulated Sequoia galaxy. When a massive discy galaxy merges, it leaves behind debris with a complex phase-space structure, a large range of orbital properties, and a range of chemical abundances. Observationally, this results in substructures with very different properties, which can be misinterpreted as implying independent progeny. Detailed chemical abundances of large samples of stars and tailored hydrodynamical simulations are critical to resolving such conundrums.

The Gaia mission has provided the largest ever astrometric chart of the Milky Way. Using it to map the Galactic halo is helpful for disentangling its merger history. The identification of halo stars in Gaia DR2 with reliable distance estimates requires special methods because such stars are typically farther away and scarce. We apply the reduced proper motion (RPM) method to identify halo main sequence stars on the basis of Gaia photometry and proper motions. Using the colour-absolute-magnitude relation for this type of stars, we calculate photometric distances. Our selection results in a set of $\sim10^7$ tentative main sequence halo stars with typical distance uncertainties of $7\%$ and with median velocity errors of 20 km/s. The median distance of our sample is $\sim 4.4$ kpc, with the faintest stars located at $\sim 16$ kpc. The spatial distribution of the stars in our sample is centrally concentrated. Visual inspection of the mean velocities of stars on the sky reveals large-scale patterns as well as clear imprints of the GD-1 stream and tentative hints of the Jhelum and Leiptr streams. Incompleteness and selection effects limit our ability to interpret the patterns reliably as well as to identify new substructures. We define a pseudo-velocity space by setting to zero the line-of-sight velocities of our sample stars. In this space, we recover several known structures such as the footprint of Gaia-Enceladus (i.e. the Gaia-Sausage) as well as the Helmi streams and some other retrograde substructures (Sequoia, Thamnos). We show that the two-point velocity correlation function reveals significant clustering on scales smaller than 100 km/s, of similar amplitude as found for the 6D Gaia halo sample. This clustering indicates the presence of nearby streams that are predominantly phase-mixed.

Aims. Several kinematic and chemical substructures have been recently found amongst Milky Way halo stars with retrograde motions. It is currently unclear how these various structures are related to each other. This Letter aims to shed light on this issue. Methods. We explore the retrograde halo with an augmented version of the Gaia DR2 RVS sample, extended with data from three large spectroscopic surveys, namely RAVE, APOGEE and LAMOST. In this dataset, we identify several structures using the HDBSCAN clustering algorithm. We discuss their properties and possible links using all the available chemical and dynamical information. Results. In concordance with previous work, we find that stars with [Fe/H] $<-1$ have more retrograde motions than those with [Fe/H] $>-1$. The retrograde halo contains a mixture of debris from objects like Gaia-Enceladus, Sequoia, and even the chemically defined thick-disc. We find that the Sequoia has a smaller range in orbital energies than previously suggested and is confined to high-energy. Sequoia could be a small galaxy in itself, but since it overlaps both in integrals-of-motion space and chemical abundance space with the less bound debris of Gaia-Enceladus, its nature cannot be fully settled yet. In the low-energy part of the halo we find evidence for at least one more distinct structure: Thamnos. Stars in Thamnos are on low inclination, mildly eccentric retrograde orbits, moving at $v_{\phi}\approx-150$ km/s, and are chemically distinct from the other structures. Conclusions. Even with the excellent Gaia DR2 data it remains challenging to piece together all the fragments found in the retrograde halo. At this point, we are very much in need of large datasets with high-quality high-resolution spectra and tailored high-resolution hydrodynamical simulations of galaxy mergers.

$Context$. The assembly history experienced by the Milky Way is currently being unveiled thanks to the data provided by the $Gaia$ mission. It is likely that the globular cluster system of our Galaxy has followed a similarly intricate formation path. $Aims$. To constrain this formation path, we explore the link between the globular clusters and the known merging events that the Milky Way has experienced. $Methods$. To this end, we combined the kinematic information provided by $Gaia$ for almost all Galactic clusters, with the largest sample of cluster ages available after carefully correcting for systematic errors. To identify clusters with a common origin we analysed their dynamical properties, particularly in the space of integrals of motion. $Results$. We find that about 40% of the clusters likely formed in situ. A similarly large fraction, 35%, appear to be possibly associated to known merger events, in particular to $Gaia$-Enceladus (19%), the Sagittarius dwarf galaxy (5%), the progenitor of the Helmi streams (6%), and to the Sequoia galaxy (5%), although some uncertainty remains due to the degree of overlap in their dynamical characteristics. Of the remaining clusters, 16% are tentatively associated to a group with high binding energy, while the rest are all on loosely bound orbits and likely have a more heterogeneous origin. The resulting age-metallicity relations are remarkably tight and differ in their detailed properties depending on the progenitor, providing further confidence on the associations made. $Conclusions$. We provide a table listing the likely associations. Improved kinematic data by future Gaia data releases and especially a larger, systematic error-free sample of cluster ages would help to further solidify our conclusions.

The halo of the Milky Way has long been hypothesized to harbour significant amounts of merger debris. This view has been supported over more than a decade by wide-field photometric surveys which have revealed the outer halo to be lumpy. The recent release of Gaia DR2 is allowing us to establish that mergers also have been important and possibly built up the majority of the inner halo. In this work we focus on the Helmi streams, a group of streams crossing the Solar vicinity and known for almost two decades. We characterize their properties and relevance for the build-up of the Milky Way's halo. We identify new members of the Helmi streams in an unprecedented dataset with full phase-space information combining Gaia DR2, and the APOGEE DR2, RAVE DR5 and LAMOST DR4 spectroscopic surveys. Based on the orbital properties of the stars, we find new stream members up to a distance of 5 kpc from the Sun, which we characterize using photometry and metallicity information. We also perform N-body experiments to constrain the time of accretion and properties of the progenitor of the streams. We find nearly 600 new members of the Helmi streams. Their HR diagram reveals a broad age range, from approximately 11 to 13 Gyr, while their metallicity distribution goes from $\sim$ 2.3 to $\sim$1.0, and peaks at [Fe/H] $\sim$1.5. These findings confirm that the streams originate in a dwarf galaxy. Furthermore, we find 7 globular clusters to be likely associated, and which follow a well-defined age-metallicity sequence whose properties suggest a relatively massive progenitor object. Our N-body simulations favour a system with a stellar mass of $\sim 10^8\,\mathrm{M}_\odot$ accreted $5 - 8$ Gyr ago. The debris from the Helmi streams is an important donor to the MilkyWay halo, contributing approximately 15\% of its mass in field stars and 10\% of its globular clusters.

The assembly process of our Galaxy can be retrieved using the motions and chemistry of individual stars. Chemo-dynamical studies of the nearby halo have long hinted at the presence of multiple components such as streams, clumps, duality and correlations between the stars' chemical abundances and orbital parameters. More recently, the analysis of two large stellar surveys have revealed the presence of a well-populated chemical elemental abundance sequence, of two distinct sequences in the colour-magnitude diagram, and of a prominent slightly retrograde kinematic structure all in the nearby halo, which may trace an important accretion event experienced by the Galaxy. Here report an analysis of the kinematics, chemistry, age and spatial distribution of stars in a relatively large volume around the Sun that are mainly linked to two major Galactic components, the thick disk and the stellar halo. We demonstrate that the inner halo is dominated by debris from an object which at infall was slightly more massive than the Small Magellanic Cloud, and which we refer to as Gaia-Enceladus. The stars originating in Gaia-Enceladus cover nearly the full sky, their motions reveal the presence of streams and slightly retrograde and elongated trajectories. Hundreds of RR Lyrae stars and thirteen globular clusters following a consistent age-metallicity relation can be associated to Gaia-Enceladus on the basis of their orbits. With an estimated 4:1 mass-ratio, the merger with Gaia-Enceladus must have led to the dynamical heating of the precursor of the Galactic thick disk and therefore contributed to the formation of this component approximately 10 Gyr ago. These findings are in line with simulations of galaxy formation, which predict that the inner stellar halo should be dominated by debris from just a few massive progenitors.

The early Universe presented a star formation environment that was almost devoid of heavy elements. The lowest metallicity stars thus provide a unique window into the earliest Galactic stages, but are exceedingly rare and difficult to find. Here we present the discovery of an ultra-metal-poor star, Pristine_221.8781+9.7844, using narrow-band Ca H&K photometry from the Pristine survey. Follow-up medium and high-resolution spectroscopy confirms the ultra-metal-poor nature of Pristine_221.8781+9.7844 ([Fe/H] = -4.66 +/- 0.13 in 1D LTE) with an enhancement of 0.3-0.4 dex in alpha-elements relative to Fe, and an unusually low carbon abundance. We derive an upper limit of A(C) = 5.6, well below typical A(C) values for such ultra metal-poor stars. This makes Pristine_221.8781+9.7844 one of the most metal-poor stars; in fact, it is very similar to the most metal-poor star known (SDSS J102915+172927). The existence of a class of ultra metal-poor stars with low(er) carbon abundances suggest that there must have been several formation channels in the early Universe through which long-lived, low-mass stars were formed.

We explore the phase-space structure of nearby halo stars identified kinematically from Gaia DR2 data. We focus on their distribution in velocity and in "integrals of motion" space as well as on their photometric properties. Our sample of stars selected to be moving at a relative velocity of at least 210 km/s with respect to the Local Standard of Rest, contains an important contribution from the low rotational velocity tail of the disk(s). The $V_R$-distribution of these stars depicts a small asymmetry similar to that seen for the faster rotating thin disk stars near the Sun. We also identify a prominent, slightly retrograde "blob", which traces the metal-poor halo main sequence reported by Gaia Collaboration et al. (2018d). We also find many small clumps especially noticeable in the tails of the velocity distribution of the stars in our sample. Their HR diagrams disclose narrow sequences characteristic of simple stellar populations. This stream-frosting confirms predictions from cosmological simulations, namely that substructure is most apparent amongst the fastest moving stars, typically reflecting more recent accretion events.

There is long tradition extending more than a century on the identification of moving groups in the Solar neighborhood. However, with the advent of large kinematic surveys, and especially of the upcoming Gaia data releases there is a need for more sophisticated and automated substructure finders. We analyze the TGASxRAVE dataset to identify moving groups in the Galactic disk. These groups of stars may then be used to map dynamical and star formation processes in the vicinity of the Sun. We use the ROCKSTAR algorithm, a "friends-of-friends"-like substructure finder in 6D phase-space, and analyze the Hertzsprung-Russell diagrams of the groups identified. We find 125 moving groups within 300 pc of the Sun, containing on average 50 stars, and with 3D velocity dispersions smaller than 10 km/s. Most of these groups were previously unknown. Our photometric analysis allows us to isolate a subsample of 30 statistically significant groups likely composed of stars that were born together.

We model the time evolution of gaps in tidal streams caused by the impact of a dark matter subhalo, while both orbit a spherical gravitational potential. To this end, we make use of the simple behaviour of orbits in action-angle space. A gap effectively results from the divergence of two nearby orbits whose initial phase-space separation is, for very cold thin streams, largely given by the impulse induced by the subhalo. We find that in a spherical potential the size of a gap increases linearly with time for sufficiently long timescales. We have derived an analytic expression that shows how the growth rate depends on the mass of the perturbing subhalo, its scale and its relative velocity with respect to the stream. We have verified these scalings using N-body simulations and find excellent agreement. For example, a subhalo of mass 10^8 Msun directly impacting a very cold thin stream on an inclined orbit can induce a gap that may reach a size of several tens of kpc after a few Gyr. The gap size fluctuates importantly with phase on the orbit, and it is largest close to pericentre. This indicates that it may not be fully straightforward to invert the spectrum of gaps present in a stream to recover the mass spectrum of the subhalos.