We use the second Gaia data release in combination with the catalog of Sanders & Das (2018) to dissect the Milky Way disc in phase-space and relative ages. We confirm and report the existence of ...multiple velocity moving groups at low azimuthal velocities and angular momenta, below Arcturus, regularly separated by \(\sim~18-20\,\rm{km s^{-1}}\) in azimuthal velocity. Such features were predicted to exist more than ten years ago from the reaction of the Milky Way to a perturbation in the disc undergoing phase-mixing. These structures appear slightly younger than their phase-space surroundings, arguing against an extra-galactic origin. We also identify in relative age many of the classical ridges in the plane of azimuthal velocity vs. Galactocentric radius, traditionally associated with resonance features. These ridges are also younger than their phase-space surroundings in accordance with predictions from recent state-of-the-art cosmological hydrodynamical simulations of Milky Way-like galaxies. We study the response of dynamically young and old stellar disc populations to resonances from an analytic model of a large bar, which, remarkably, qualitatively reproduces the trends seen in the data. Our results re-inforce the idea that the Galactic disc is currently being shaped by both internal and external perturbations, and that, while absolute isochrone ages have to be taken with great care, exploring the dynamical structure of the disc in stellar ages, especially with future asteroseismic data, will provide much stronger constraints than metallicity/abundance trends alone.
The virial masses of ultra-diffuse galaxies (UDGs) have been estimated using the kinematics and abundance of their globular cluster populations, leading to disparate results. Some studies conclude ...that UDGs reside in massive dark matter halos while others, controversially, argue for the existence of UDGs with no dark matter at all. Here we show that these results arise because the uncertainties of these mass estimates have been substantially underestimated. Indeed, applying the same procedure to the well-studied Fornax dwarf spheroidal would conclude that it has an "overmassive" dark halo or, alternatively, that it lacks dark matter. We corroborate our argument with self-consistent mocks of tracers in cosmological halos, showing that masses from samples with \(5 < N < 10\) tracers (assuming no measurement errors) are uncertain by at least an order of magnitude. Finally, we estimate masses of UDGs with HST imaging in Coma and show that their recent mass measurements (with adequate uncertainties) are in agreement with that of other dwarfs, such as Fornax.. We also provide bias and scatter factors for a range of sample sizes and measurement errors, of wider applicability.
Mergers and tidal interactions between massive galaxies and their dwarf satellites are a fundamental prediction of the Lambda-Cold Dark Matter cosmology. These events are thought to provide important ...observational diagnostics of nonlinear structure formation. Stellar streams in the Milky Way and Andromeda are spectacular evidence for ongoing satellite disruption. However, constructing a statistically meaningful sample of tidal streams beyond the Local Group has proven a daunting observational challenge, and the full potential for deepening our understanding of galaxy assembly using stellar streams has yet to be realised. Here we introduce the Stellar Stream Legacy Survey, a systematic imaging survey of tidal features associated with dwarf galaxy accretion around a sample of ~3100 nearby galaxies within z~0.02, including about 940 Milky Way analogues. Our survey exploits public deep imaging data from the DESI Legacy Imaging Surveys, which reach surface brightness as faint as ~29 mag/arcsec^2 in the r band. As a proof of concept of our survey, we report the detection and broad-band photometry of 24 new stellar streams in the local Universe. We discuss how these observations can yield new constraints on galaxy formation theory through comparison to mock observations from cosmological galaxy simulations. These tests will probe the present-day mass assembly rate of galaxies, the stellar populations and orbits of satellites, the growth of stellar halos and the resilience of stellar disks to satellite bombardment.
Signatures of vertical disequilibrium have been observed across the Milky Way's disk. These signatures manifest locally as unmixed phase-spirals in \(z\)--\(v_z\) space ("snails-in-phase") and ...globally as nonzero mean \(z\) and \(v_z\) which wraps around as a physical spiral across the \(x\)--\(y\) plane ("snails-in-space"). We explore the connection between these local and global spirals through the example of a satellite perturbing a test-particle Milky Way (MW)-like disk. We anticipate our results to broadly apply to any vertical perturbation. Using a \(z\)--\(v_z\) asymmetry metric we demonstrate that in test-particle simulations: (a) multiple local phase-spiral morphologies appear when stars are binned by azimuthal action \(J_\phi\), excited by a single event (in our case, a satellite disk-crossing); (b) these distinct phase-spirals are traced back to distinct disk locations; and (c) they are excited at distinct times. Thus, local phase-spirals offer a global view of the MW's perturbation history from multiple perspectives. Using a toy model for a Sagittarius (Sgr)-like satellite crossing the disk, we show that the full interaction takes place on timescales comparable to orbital periods of disk stars within \(R \lesssim 10\) kpc. Hence such perturbations have widespread influence which peaks in distinct regions of the disk at different times. This leads us to examine the ongoing MW-Sgr interaction. While Sgr has not yet crossed the disk (currently, \(z_{Sgr} \approx -6\) kpc, \(v_{z,Sgr} \approx 210\) km/s), we demonstrate that the peak of the impact has already passed. Sgr's pull over the past 150 Myr creates a global \(v_z\) signature with amplitude \(\propto M_{Sgr}\), which might be detectable in future spectroscopic surveys.
Indications of disequilibrium throughout the Milky Way (MW) highlight the need for compact,flexible, non-parametric descriptions of phase--space distributions of galaxies. We present a new ...representation of the current Dark Matter (DM) distribution and potential derived from N-body simulations of the Milky Way and Large Magellanic Cloud (LMC) system using Basis Function Expansions (BFEs). We incorporate methods to maximize the physical signal in the representation. As a result, the simulations of \(10^8\) DM particles representing the MW--LMC system can be described by 354 coefficients. We find that the LMC induces asymmetric perturbations (odd l, m) to the MW's halo, which are not well-described by oblate, prolate, or triaxial halos. Furthermore, the energy in high-order even modes (l,m \(\geq\) 2) is similar to average triaxial halos found in cosmological simulations. As such, the response of the MW's halo to the LMC must be accounted for in order to recover the imprints of its assembly history. The LMC causes the outer halo (\(\geq\) 30 kpc) to shift from the disk center of mass (COM) by \(\sim\)15-25 kpc at present day, manifesting as a dipole in the BFE and in the radial velocities of halo stars. The shift depends on the LMC's infall mass, the distortion of the LMC's halo and the MW halo response. Within 30 kpc, halo tracers are expected to orbit the COM of the MW's disk, regardless of LMC infall mass. The LMC's halo is also distorted by MW tides, we discuss the implications for its mass loss and the subsequent effects on current Magellanic satellites.
Using Gaia DR2, we trace the Anticenter Stream (ACS) in various stellar populations across the sky and find that it is kinematically and spatially decoupled from the Monoceros Ring. Using stars from ...{\sc lamost} and {\sc segue}, we show that the ACS is systematically more metal-poor than Monoceros by \(0.1\) dex with indications of a narrower metallicity spread. Furthermore, the ACS is predominantly populated of old stars (\(\sim 10\,\rm{Gyr}\)), whereas Monoceros has a pronounced tail of younger stars (\(6-10\, \rm{Gyr}\)) as revealed by their cumulative age distributions. Put togehter, all of this evidence support predictions from simulations of the interaction of the Sagittarius dwarf with the Milky Way, which argue that the Anticenter Stream (ACS) is the remains of a tidal tail of the Galaxy excited during Sgr's first pericentric passage after it crossed the virial radius, whereas Monoceros consists of the composite stellar populations excited during the more extended phases of the interaction. We suggest that the ACS can be used to constrain the Galactic potential, particularly its flattening, setting strong limits on the existence of a dark disc. Importantly, the ACS can be viewed as a stand-alone fossil of the chemical enrichment history of the Galactic disc.
Several works have recently applied Jeans modelling to Gaia-based datasets to
infer the circular velocity curve for the Milky Way. Such works have
consistently found evidence for a continuous decline ...in the rotation curve
beyond $\sim$15kpc possibly indicative of a light dark matter halo. We used
Gaia DR3 RVS data, supplemented with Bayesian distances to determine the radial
variation of the second moments of the velocity distribution for stars close to
the Galactic plane. We have used these profiles to determine the rotation curve
using the Jeans equations under the assumption of axisymmetry and explored how
they vary with azimuth and above and below the Galactic disk plane. We have
applied the same methodology to an N-body simulation of a Milky Way-like galaxy
impacted by a satellite akin the Sagittarius dwarf and to the Auriga suite of
cosmological simulations. We reveal evidence of disequilibrium and deviations
from axisymmetry closer in. We find that the second moment of $V_R$ flattens
out at $R \gtrsim 12.5$kpc, and that the second moment of $V_{\phi}$ is
different above and below the plane for $R \gtrsim 11$kpc. The simulations
indicate that these features are typical of galaxies that have been perturbed
by external satellites. They also suggest that the difference between the true
circular velocity curve and that inferred from Jeans equations can be as high
as 15$\%$, but is likely of order 10$\%$ for the Milky Way. This is of larger
amplitude than the systematics inherent to Jeans equations. However, if the
density of the tracer population were truncated at large radii, the erroneous
conclusion of a steeply declining rotation curve can be reached. We find that
steady-state axisymmetric Jeans modelling becomes less robust at large radii,
indicating that particular caution is needed when interpreting the rotation
curve inferred in those regions.
Stars aligned in thin stream-like features (feathers), with widths of \(\delta\sim1-10^{\circ}\) and lengths as large as \(\Delta l\sim180^{\circ}\), have been observed towards the Anticenter of our ...Galaxy and their properties mapped in abundances and phase-space. We study their origin by analysing similar features arising in an N-body simulation of a Galactic disc interacting with a Sagittarius-like dwarf spheroidal galaxy (Sgr). By following the orbits of the particles identified as contributing to feathers backwards in time, we trace their excitation to one of Sgr's previous pericentric passages. These particles initially span a large range of phase-angles but a tight range of radii, suggesting they provide a probe of populations in distinct annuli in the outer Galactic disc. The structures are long lived and persist after multiple passages on timescales of \(\sim4 \,\rm{Gyrs}\). On the sky, they exhibit oscillatory motion that can be traced with a single orbit mapped over much of their full length and with amplitudes and gradients similar to those observed. We demonstrate how these properties of feathers may be exploited to measure the potential, its flattening, as well as infer the strength of recent potential perturbations.
Satellite galaxies are predicted to generate gravitational density wakes as they orbit within the dark matter (DM) halos of their hosts, causing their orbits to decay over time. The recent infall of ...the Milky Way's (MW) most massive satellite galaxy, the Large Magellanic Cloud (LMC), affords us the unique opportunity to study this process in action. In this work, we present high-resolution (\(m_{dm} = 4 \times 10^4 M_{\odot}\) ) N-body simulations of the MW-LMC interaction over the past 2 Gyr. We quantify the impact of the LMC's passage on the density and kinematics of the MW's DM halo and the observability of these structures in the MW's stellar halo. The LMC is found to generate pronounced Local and Global wakes in both the DM and stellar halos, leads to both local overdensities and distinct kinematic patterns that should be observable with ongoing and future surveys. Specifically, the Global Wake will result in redshifted radial velocities of stars in the North and blueshifts in the South, at distances larger than 45 kpc. The Local Wake traces the orbital path of the LMC through the halo (50-200 kpc), resulting in a stellar overdensity with a distinct, tangential kinematic pattern that persists to the present day. The detection of the MW's halo response will constrain: the infall mass of the LMC and its orbital trajectory, the mass of the MW, and it may inform us about the nature of the dark matter particle itself.
The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer ...properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter GC19) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself as well as the LMC-induced DM wake) and small scales (due to substructure). Using the GC19 simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components, and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs Bullock & Johnston (2005) suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream -- the most recent, massive accretion event in the MW apart from the LMC -- and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC's infall.