The dual origin of the Galactic thick disc and halo from the gas-rich Gaia-Enceladus Sausage merger

We analyse a set of cosmological magnetohydrodynamic simulations of the formation of Milky Way-mass galaxies identified to have a prominent radially anisotropic stellar halo component similar to the so-called 'Gaia Sausage' found in the Gaia data. We examine the effects of the progenitor of the Sausage (the Gaia-Enceladus Sausage, GES) on the formation of major galactic components analogous to the Galactic thick disc and inner stellar halo. We find that the GES merger is likely to have been gas-rich and contribute 10-50 per cent of gas to a merger-induced centrally concentrated starburst that results in the rapid formation of a compact, rotationally supported thick disc that occupies the typical chemical thick disc region of chemical abundance space. We find evidence that gas-rich mergers heated the proto-disc of the Galaxy, scattering stars on to less-circular orbits such that their rotation velocity and metallicity positively correlate, thus contributing an additional component that connects the Galactic thick disc to the inner stellar halo. We demonstrate that the level of kinematic heating of the proto-galaxy correlates with the kinematic state of the population before the merger, the progenitor mass, and orbital eccentricity of the merger. Furthermore, we show that the mass and time of the merger can be accurately inferred from local stars on counter-rotating orbits

The effects of dynamical substructure on Milky Way mass estimates from the high velocity tail of the local stellar halo

We investigate the impact of dynamical streams and substructure on estimates of the local escape speed and total mass of Milky Way-mass galaxies from modelling the high velocity tail of local halo stars. We use a suite of high-resolution, magneto-hydrodynamical cosmological zoom-in simulations, which resolve phase space substructure in local volumes around solar-like positions. We show that phase space structure varies significantly between positions in individual galaxies and across the suite. Substructure populates the high velocity tail unevenly and leads to discrepancies in the mass estimates. We show that a combination of streams, sample noise and truncation of the high velocity tail below the escape speed leads to a distribution of mass estimates with a median that falls below the true value by $\sim 20 \%$, and a spread of a factor of 2 across the suite. Correcting for these biases, we derive a revised value for the Milky Way mass presented in Deason et al. of $1.29 ^{+0.37}_{-0.47} \times 10^{12}$ $\rm M_{\odot}$.Comment: Re-submitted to MNRAS Letters after minor revisio

The effects of dynamical substructure on Milky Way mass estimates from the high-velocity tail of the local stellar halo

We investigate the impact of dynamical streams and substructure on estimates of the local escape speed and total mass of Milky-Way-mass galaxies from modelling the high-velocity tail of local halo stars. We use a suite of high-resolution magnetohydrodynamical cosmological zoom-in simulations that resolve phase space substructure in local volumes around solar-like positions. We show that phase space structure varies significantly between positions in individual galaxies and across the suite. Substructure populates the high-velocity tail unevenly and leads to discrepancies in the mass estimates. We show that a combination of streams, sample noise, and truncation of the high-velocity tail below the escape speed leads to a distribution of mass estimates with a median that falls below the true value by ∼20 per cent ∼20 per cent ⁠, and a spread of a factor of 2 across the suite. Correcting for these biases, we derive a revised value for the Milky Way mass presented in Deason et al. of 1.29 +0.37 −0.47 × 10 12 M ⊙ 1.29−0.47+0.37×1012M⊙ ⁠

HALO7D III: Chemical Abundances of Milky Way Halo Stars from Medium Resolution Spectra

The Halo Assembly in Lambda Cold Dark Matter: Observations in 7 Dimensions (HALO7D) survey measures the kinematics and chemical properties of stars in the Milky Way (MW) stellar halo to learn about the formation of our Galaxy. HALO7D consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of MW halo main sequence turn-off (MSTO) stars in the four CANDELS fields. HALO7D consists of deep pencil beams, making it complementary to other contemporary wide-field surveys. We present the [Fe/H] and [alpha/Fe] abundances for 113 HALO7D stars in the Galactocentric radial range of $\sim10-40$ kpc. Using the full 7D chemodynamical data (3D positions, 3D velocities, and abundances) of HALO7D, we measure the velocity anisotropy, $\beta$, of the halo velocity ellipsoid for each field and for different metallicity-binned subsamples. We find that two of the four fields have stars on very radial orbits while the remaining two have stars on more isotropic orbits. Separating the stars into high, mid, and low [Fe/H] bins at $-2.2$ dex and $-1.1$ dex for each field separately, we find differences in the anisotropies between the fields and between the bins; some fields appear dominated by radial orbits in all bins while other fields show variation between the [Fe/H] bins. These chemodynamical differences are evidence that the HALO7D fields have different fractional contributions from the progenitors that built up the MW stellar halo. Our results highlight the additional information that is available on smaller spatial scales when compared to results from a spherical average of the stellar halo.Comment: 32 pages, 15 figure

The Pisces Plume and the Magellanic wake

Using RR Lyrae stars in the Gaia Data Release 2 and Pan-STARRS1 we study the properties of the Pisces overdensity, a diffuse substructure in the outer halo of the Milky Way. We show that along the line of sight, Pisces appears as a broad and long plume of stars stretching from 40 to 110 kpc with a steep distance gradient. On the sky Pisces’s elongated shape is aligned with the Magellanic Stream. Using follow-up VLT FORS2 spectroscopy, we have measured the velocity distribution of the Pisces candidate member stars and have shown it to be as broad as that of the Galactic halo but offset to negative velocities. Using a suite of numerical simulations, we demonstrate that the structure has many properties in common with the predicted behaviour of the Magellanic wake, i.e. the Galactic halo overdensity induced by the infall of the Magellanic Clouds

Simulating cosmological substructure in the solar neighbourhood

We explore the predictive power of cosmological, hydrodynamical simulations for stellar phase-space substructure and velocity correlations with the AURIGA simulations and AURIGAIA mock Gaia catalogues. We show that at the solar circle the AURIGA simulations commonly host phase-space structures in the stellar component that have constant orbital energies and arise from accreted subhaloes. These structures can persist for a few Gyr, even after coherent streams in position space have been erased. We also explore velocity two-point correlation functions and find this diagnostic is not deterministic for particular clustering patterns in phase space. Finally, we explore these structure diagnostics with the AURIGAIA catalogues and show that current catalogues have the ability to recover some structures in phase space but careful consideration is required to separate physical structures from numerical structures arising from catalogue generation methods

Aurigaia: mock Gaia DR2 stellar catalogues from the Auriga cosmological simulations

We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic ΛCDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with V 20 deg. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. (2015) that preserves the phase-space distribution of the model stars. These publicly available catalogues contain five-parameter astrometry, radial velocities, multiband photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that (i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A - and B -dwarf stars up to radii of ∼15 kpc, and (ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to (i) develop and test analysis methods for the Gaia DR2 data, (ii) gauge the limitations and biases of the data, and (iii) interpret the data in the light of theoretical predictions from realistic ab initio simulations of galaxy formation in the ΛCDM cosmological model

Beta Dips in the Gaia Era: Simulation Predictions of the Galactic Velocity Anisotropy Parameter (β) for Stellar Halos

The velocity anisotropy parameter, β, is a measure of the kinematic state of orbits in the stellar halo, which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for β as a function of radius from three suites of simulations, including accretion-only and cosmological hydrodynamic simulations. We find that the two types of simulations are consistent and predict strong radial anisotropy ($\langle \beta \rangle \sim 0.7$) for Galactocentric radii greater than 10 kpc. Previous observations of β for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r ~ 20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in β. We find that dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z ~ 1 can result in a present-day low (isotropic to tangential) value of β over a broad range of radii and angles. While all of these mechanisms are plausible drivers for the β dip observed in the MW, each mechanism in the simulations has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of β(r) is required for measuring the mass of the MW halo, we note that significant transient dips in β could cause an overestimate of the halo's mass when using spherical Jeans equation modeling