We study the late-time evolution of the central regions of two Milky Way-like simulations of galaxies formed in a cosmological context, one hosting a fast bar and the other a slow one. We find that bar length, R_b, measurements fluctuate on a dynamical timescale by up to 100%, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15%, correlating with R_b. The Tremaine-Weinberg-method estimates of the bars' instantaneous pattern speeds show variations around the mean of up to ~20%, typically anti-correlating with the bar length and strength. Through power spectrum analyses, we establish that these bar pulsations, with a period in the range ~60-200 Myr, result from its interaction with multiple spiral modes, which are coupled with the bar. Because of the presence of odd spiral modes, the two bar halves typically do not connect at exactly the same time to a spiral arm, and their individual lengths can be significantly offset. We estimated that in about 50% of bar measurements in Milky Way-mass external galaxies, the bar lengths of SBab type galaxies are overestimated by ~15% and those of SBbc types by ~55%. Consequently, bars longer than their corotation radius reported in the literature, dubbed "ultra-fast bars", may simply correspond to the largest biases. Given that the Scutum-Centaurus arm is likely connected to the near half of the Milky Way bar, recent direct measurements may be overestimating its length by 1-1.5 kpc, while its present pattern speed may be 5-10 km/s/kpc smaller than its time-averaged value.

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 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.

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.

We use Gaia DR2 astrometric and line-of-sight velocity information combined with two sets of distances obtained with a Bayesian inference method to study the 3D velocity distribution in the Milky Way disc. We search for variations in all Galactocentric cylindrical velocity components ($V_{\phi}$, $V_R$ and $V_z$) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent work showing that bulk vertical motions in the $R\text{-}z$ plane are consistent with a combination of breathing and bending modes. In the $x\text{-}y$ plane, we show that, although the amplitudes change, the structure produced by these modes is mostly invariant as a function of distance from the plane. Comparing to two different Galactic disc models, we demonstrate that the observed patterns can drastically change in short time intervals, showing the complexity of understanding the origin of vertical perturbations. A strong radial $V_R$ gradient was identified in the inner disc, transitioning smoothly from $16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $30^\circ<\phi<45^\circ$ ahead of the Sun-Galactic centre line, to $-16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $-45^\circ<\phi<-30^\circ$ lagging the solar azimuth. We use a simulation with no significant recent mergers to show that exactly the opposite trend is expected from a barred potential, but overestimated distances can flip this trend to match the data. Alternatively, using an $N$-body simulation of the Sagittarius dwarf-Milky Way interaction, we demonstrate that a major recent perturbation is necessary to reproduce the observations. Such an impact may have strongly perturbed the existing bar or even triggered its formation in the last $1\text{-}2$ Gyr.

Orbital properties of stars, computed from their six-dimensional phase space measurements and an assumed Galactic potential, are used to understand the structure and evolution of the Galaxy. Stellar actions, computed from orbits, have the attractive quality of being invariant under certain assumptions and are therefore used as quantitative labels of a star's orbit. We report a subtle but important systematic error that is induced in the actions as a consequence of local midplane variations expected for the Milky Way. This error is difficult to model because it is non-Gaussian and bimodal, with neither mode peaking on the null value. An offset in the vertical position of the Galactic midplane of $\sim15\,\text{pc}$ for a thin disk-like orbit or $\sim 120\,\text{pc}$ for a thick disk-like orbit induces a $25\%$ systematic error in the vertical action $J_z$. In FIRE simulations of Milky Way-mass galaxies, these variations are on the order of $\sim100\,\text{pc}$ at the solar circle. From observations of the mean vertical velocity variation of $\sim5\text{--}10\,\text{km}\,\text{s}^{-1}$ with radius, we estimate that the Milky Way midplane variations are $\sim60\text{--}170\,\text{pc}$, consistent with three-dimensional dust maps. Action calculations and orbit integrations, which assume the global and local midplanes are identical, are likely to include this induced error, depending on the volume considered. Variation in the local standard of rest or distance to the Galactic center causes similar issues. The variation of the midplane must be taken into account when performing dynamical analysis across the large regions of the disk accessible to Gaia and future missions.

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.

Recent studies have shown that the passage of a massive satellite through the disk of a spiral galaxy can induce vertical wobbles in the disk and produce features such as in-plane rings and phase-space spirals. Here we analyze a high-resolution N-body simulation of a live stellar disk perturbed by the recent passage of a massive dwarf galaxy that induces such perturbations. We study the implications of the phase-space structures for the estimate of the matter density through traditional Jeans modelling. The dwarf satellite excites rapid time-variations in the potential, leading to a significant bias of the local matter surface density determined through such a method. In particular, while the Jeans modelling gives reasonable estimates in the most overdense regions of the disk, we show that it tends to overestimate the dynamical surface density in underdense regions. In these regions, the phase-space spiral is indeed more marked in the surface of section of height vs. vertical velocity. This prediction can be verified with future Gaia data releases. Our finding is highly relevant for future attempts at determining the dynamical surface density of the outer Milky Way disk as a function of radius. The outer disk of the Milky Way is indeed heavily perturbed, and \textit{Gaia} DR2 data have clearly shown that such phase-space perturbations are even present locally. While our results show that traditional Jeans modelling should give reliable results in overdense regions of the disk, the important biases in underdense regions call for the development of non-equilibrium methods to estimate the dynamical matter density locally and in the outer disk.

In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potential and then explore the perturbative effects of the Large Magellanic Cloud (LMC). We find that the LMC can produce precisely the observed motion-track mismatch and we therefore use the Orphan stream to measure the mass of the Cloud. We simultaneously fit the Milky Way and LMC potentials and infer that a total LMC mass of $1.38^{+0.27}_{-0.24} \times10^{11}\,\rm{M_\odot}$ is required to bend the Orphan Stream, showing for the first time that the LMC has a large and measurable effect on structures orbiting the Milky Way. This has far-reaching consequences for any technique which assumes that tracers are orbiting a static Milky Way. Furthermore, we measure the Milky Way mass within 50 kpc to be $3.80^{+0.14}_{-0.11}\times10^{11} M_\odot$. Finally, we use these results to predict that, due to the reflex motion of the Milky Way in response to the LMC, the outskirts of the Milky Way's stellar halo should exhibit a bulk, upwards motion.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics.

Studying our Galaxy, the Milky Way (MW), gives us a close-up view of the interplay between cosmology, dark matter, and galaxy formation. In the next decade our understanding of the MW's dynamics, stellar populations, and structure will undergo a revolution thanks to planned and proposed astrometric, spectroscopic and photometric surveys, building on recent advances by the Gaia astrometric survey. Together, these new efforts will measure three-dimensional positions and velocities and numerous chemical abundances for stars to the MW's edge and well into the Local Group, leading to a complete multidimensional view of our Galaxy. Studies of the multidimensional Milky Way beyond the Gaia frontier---from the edge of the Galactic disk to the edge of our Galaxy's dark matter halo---will unlock new scientific advances across astrophysics, from constraints on dark matter to insights into galaxy formation.

We analyse an N-body simulation of the interaction of the Milky Way (MW) with a Sagittarius-like dSph (Sgr), looking for signatures which may be attributed to its orbital history in the phase space volume around the Sun in light of $Gaia$ DR2 discoveries. The repeated impacts of Sgr excite coupled vertical and radial oscillations in the disc which qualitatively, and to a large degree quantitatively are able to reproduce many features in the 6D $Gaia$ DR2 samples, from the median $V_{R},V_{\phi}, V_{z} $ velocity maps to the local $\delta\rho(v_{z},z)$ phase-space spiral which is a manifestation of the global disc response to coupled oscillations within a given volume. The patterns in the large-scale velocity field are well described by tightly wound spirals and vertical corrugations excited from Sgr's impacts. We show that the last pericentric passage of Sgr resets the formation of the local present-day $\delta\rho(v_z,z)$ spiral and situate its formation around 500-800 Myr. As expected $\delta\rho(vz,z)$ grows in size and decreases in woundedness as a function of radius in both the $Gaia$ DR2 data and simulations. This is the first N-body model able to explain so many of the features in the data on different scales. We demonstrate how to use the full extent of the Galactic disc to date perturbations dating from Myr to Gyr, probe the underlying potential and constrain the mass-loss history of Sgr. $\delta\rho(vz,z)$ looks the same in all stellar populations age bins down to the youngest ages which rules out a bar buckling origin.

We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of the large number of RR Lyrae stars in the OS, we extract accurate distances and proper motions across the entire stretch of the tidal debris studied. As delineated by the GDR2 RR Lyrae, the stream exhibits two prominent twists in its shape on the sky which are accompanied by changes in the tangential motion. We complement the RR Lyrae maps with those created using GDR2 Red Giants and the DECam Legacy Survey Main Sequence Turn-Off stars. The behavior of the OS track on the sky is consistent across all three tracers employed. We detect a strong non-zero motion in the across-stream direction for a substantial portion of the stream. Such a misalignment between the debris track and the streaming velocity cannot be reproduced in a static gravitational potential and signals an interaction with a massive perturber.

Dark matter-baryon scaling relations in galaxies are important in order to constrain galaxy formation models. Here, we provide a modern quantitative assessment of those relations, by modelling the rotation curves of galaxies from the Spitzer Photometry and Accurate Rotation Curves (SPARC) database with the Einasto dark halo model. We focus in particular on the comparison between the original SPARC parameters, with constant mass-to-light ratios for bulges and disks, and the parameters for which galaxies follow the tightest radial acceleration relation. We show that fits are improved in the second case, and that the pure halo scaling relations also become tighter. We report that the density at the radius where the slope is -2 is strongly anticorrelated to this radius, and to the Einasto index. The latter is close to unity for a large number of galaxies, indicative of large cores. In terms of dark matter-baryon scalings, we focus on relations between the core properties and the extent of the baryonic component, which are relevant to the cusp-core transformation process. We report a positive correlation between the core size of halos with small Einasto index and the stellar disk scale-length, as well as between the averaged dark matter density within 2 kpc and the baryon-induced rotational velocity at that radius. This finding is related to the consequence of the radial acceleration relation on the diversity of rotation curve shapes, quantified by the rotational velocity at 2 kpc. While a tight radial acceleration relation slightly decreases the observed diversity compared to the original SPARC parameters, the diversity of baryon-induced accelerations at 2 kpc is sufficient to induce a large diversity, incompatible with current hydrodynamical simulations of galaxy formation, while maintaining a tight radial acceleration relation.

We present the kinematics of Blue Straggler (BS) stars identified in the Canada-France-Imaging-Survey (CFIS), covering 4000 deg$^2$ on the sky in the $u$-band. The BS sample, characterised through CFIS and Pan-STARRS photometry, has been kinematically decomposed into putative halo and disc populations after cross-matching with Gaia astrometry and SDSS/SEGUE/LAMOST spectroscopy. This decomposition clearly reveals the strong flaring of the outer Milky Way disc. In particular, we show that we can detect this flaring up to a vertical height of $|Z| \simeq 8$ kpc at a Galactocentric distance of $R\sim $27 kpc. While some small level of flaring is expected for extended discs built up by radial migration, we demonstrate that this observed very strong flaring of the Milky Way disc is more likely consistent with it being dynamically heated by the repeated passage of the Sagittarius dwarf spheroidal galaxy through the midplane.

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.

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.

Our Galaxy is thought to have undergone an active evolutionary history dominated by star formation, the accretion of cold gas, and, in particular, mergers up to 10 gigayear ago. The stellar halo reveals rich fossil evidence of these interactions in the form of stellar streams, substructures, and chemically distinct stellar components. The impact of dwarf galaxy mergers on the content and morphology of the Galactic disk is still being explored. Recent studies have identified kinematically distinct stellar substructures and moving groups, which may have extragalactic origin. However, there is mounting evidence that stellar overdensities at the outer disk/halo interface could have been caused by the interaction of a dwarf galaxy with the disk. Here we report detailed spectroscopic analysis of 14 stars drawn from two stellar overdensities, each lying about 5 kiloparsecs above and below the Galactic plane - locations suggestive of association with the stellar halo. However, we find that the chemical compositions of these stars are almost identical, both within and between these groups, and closely match the abundance patterns of the Milky Way disk stars. This study hence provides compelling evidence that these stars originate from the disk and the overdensities they are part of were created by tidal interactions of the disk with passing or merging dwarf galaxies.

The Monoceros Ring (also known as the Galactic Anticenter Stellar Structure) and A13 are stellar overdensities at estimated heliocentric distances of $d \sim 11$ kpc and 15 kpc observed at low Galactic latitudes towards the anticenter of our Galaxy. While these overdensities were initially thought to be remnants of a tidally-disrupted satellite galaxy, an alternate scenario is that they are composed of stars from the Milky Way (MW) disk kicked out to their current location due to interactions between a satellite galaxy and the disk. To test this scenario, we study the stellar populations of the Monoceros Ring and A13 by measuring the number of RR Lyrae and M giant stars associated with these overdensities. We obtain low-resolution spectroscopy for RR Lyrae stars in the two structures and measure radial velocities to compare with previously measured velocities for M giant stars in the regions of the Monoceros Ring and A13, to assess the fraction of RR Lyrae to M giant stars ($f_{RR:MG}$) in A13 and Mon/GASS. We perform velocity modeling on 153 RR Lyrae stars (116 in the Monoceros Ring and 37 in A13) and find that both structures have very low $f_{RR:MG}$. The results support a scenario in which stars in A13 and Mon/GASS formed in the MW disk. We discuss a possible association between Mon/GASS, A13, and the Triangulum-Andromeda overdensity based on their similar velocity distributions and $f_{RR:MG}$.

We present N-body simulations of a Sagittarius-like dwarf spheroidal galaxy (Sgr) that follow its orbit about the Milky Way (MW) since its first crossing of the Galaxy's virial radius to the present day. As Sgr orbits around the MW, it excites vertical oscillations, corrugating and flaring the Galactic stellar disc. These responses can be understood by a two-phase picture in which the interaction is first dominated by torques from the wake excited by Sgr in the MW dark halo before transitioning to tides from Sgr's direct impact on the disc at late times. We show for the first time that a massive Sgr model simultaneously reproduces the locations and motions of arc-like over densities, such as the Monoceros Ring and the Triangulum Andromeda stellar clouds, that have been observed at the extremities of the disc, while also satisfying the solar neighbourhood constraints on the vertical structure and streaming motions of the disc. In additional simulations, we include the Large Magellanic Cloud (LMC) self consistently with Sgr. The LMC introduces coupling through constructive and destructive interference, but no new corrugations. In our models, the excitation of the current structure of the outer disk can be traced to interactions as far back as 6-7 Gyr ago (corresponding to $z\leq1$). Given the apparently quiescent accretion history of the MW over this timescale, this places Sgr as the main culprit behind the vertical oscillations of the disc and the last major accretion event for the Galaxy with the capacity to modulate its chemodynamical structure.

Deep photometric surveys of the Milky Way have revealed diffuse structures encircling our Galaxy far beyond the "classical" limits of the stellar disk. This paper reviews results from our own and other observational programs, which together suggest that, despite their extreme positions, the stars in these structures were formed in our Galactic disk. Mounting evidence from recent observations and simulations implies kinematic connections between several of these distinct structures. This suggests the existence of collective disk oscillations that can plausibly be traced all the way to asymmetries seen in the stellar velocity distribution around the Sun. There are multiple interesting implications of these findings: they promise new perspectives on the process of disk heating, they provide direct evidence for a stellar halo formation mechanism in addition to the accretion and disruption of satellite galaxies, and, they motivate searches of current and near-future surveys to trace these oscillations across the Galaxy. Such maps could be used as dynamical diagnostics in the emerging field of "Galactoseismology", which promises to model the history of interactions between the Milky Way and its entourage of satellites, as well examine the density of our dark matter halo. As sensitivity to very low surface brightness features around external galaxies increases, many more examples of such disk oscillations will likely be identified. Statistical samples of such features not only encode detailed information about interaction rates and mergers, but also about long sought-after dark matter halo densities and shapes. Models for the Milky Way's own Galactoseismic history will therefore serve as a critical foundation for studying the weak dynamical interactions of galaxies across the universe.

We use high-resolution $N$-body simulations to study the effect of a galactic disc on the dynamical evolution of dark matter substructures with orbits and structural parameters extracted from the Aquarius A-2 merger tree (Springel et al. 2008). Satellites are modelled as equilibrium $N$-body realizations of generalized Hernquist profiles with $2\times10^6$ particles and injected in the analytical evolving host potential at $z_\mathrm{infall}$, defined by the peak of their mass evolution. We select all substructures with $M_{200}(z_\mathrm{infall})\geq 10^8\,\mathrm{M_\odot}$ and first pericentric distances $r_p<r_{200}$. Motivated by observations of Milky Way dwarf spheroidal galaxies, we also explore satellite models with cored dark matter profiles with a fixed core size $r_c=0.8\,a_s$ where $a_s$ is the Hernquist scale radius. We find that models with cuspy satellites have twice as many surviving substructures at $z=0$ than their cored counterparts, and four times as many if we only consider those on orbits with $r_p\lesssim0.1\,r_{200}$. For a given profile, adding an evolving disc potential reduces the number of surviving substructures further by a factor of $\lesssim2$ for satellites on orbits that penetrate the disc ($r_p\lesssim 20\,\mathrm{kpc}$). For large $r_p$, where tidal forces and the effect of the disc become negligible, the number of satellites per pericentre bin converges to similar values for all four models.

We present N-body and hydrodynamical simulations of the response of the Milky Way's baryonic disc to the presence of the Large Magellanic Cloud during a first infall scenario. For a fiducial galactic model reproducing the gross properties of the Galaxy, we explore a set of six initial conditions for the LMC of varying mass which all evolve to fit the measured constraints on its current position and velocity with respect to the Galactic Center. We find that the LMC can produce strong disturbances - warping of the stellar and gaseous discs - in the Galaxy, without violating constraints from the phase-space distribution of stars in the Solar Neighbourhood. All models correctly reproduce the phases of the warp and its anti-symmetrical shape about the disc's mid-plane. If the warp is due to the LMC alone, then the largest mass model is favoured ($2.5\times10^{11}\,\rm{M_{\odot}}$). Still, some quantitative discrepancies remain, including deficits in height of $\Delta Z=0.7 \, \rm{kpc}$ at $R=22 \,\rm{kpc}$ and $\Delta Z=0.7\, \rm{kpc}$ at $R=16 \,\rm{kpc}$. This suggests that even higher infall masses for the LMC's halo are allowed by the data. A comparison with the vertical perturbations induced by a heavy Sagittarius dSph model ($10^{11}\,\rm{M_{\odot}}$) suggest that positive interference with the LMC warp is expected at $R=16 \, \rm{kpc}$. We conclude that the vertical structure of the Galactic disc beyond the Solar Neighbourhood may jointly be shaped by its most massive satellites. As such, the current structure of the MW suggests we are seeing the process of disc heating by satellite interactions in action.

The Triangulum-Andromeda stellar clouds (TriAnd1 and TriAnd2) are a pair of concentric ring- or shell-like over-densities at large $R$ ($\approx$ 30 kpc) and $Z$ ($\approx$ -10 kpc) in the Galactic halo that are thought to have been formed from the accretion and disruption of a satellite galaxy. This paper critically re-examines this formation scenario by comparing the number ratio of RR Lyrae to M giant stars associated with the TriAnd clouds with other structures in the Galaxy. The current data suggest a stellar population for these over-densities ($f_{\rm RR:MG} < 0.38$ at 95% confidence) quite unlike any of the known satellites of the Milky Way ($f_{\rm RR:MG} \approx 0.5$ for the very largest and $f_{\rm RR:MG} >>1$ for the smaller satellites) and more like the population of stars born in the much deeper potential well inhabited by the Galactic disk ($f_{\rm RR:MG} < 0.01$). N-body simulations of a Milky-Way-like galaxy perturbed by the impact of a dwarf galaxy demonstrate that, in the right circumstances, concentric rings propagating outwards from that Galactic disk can plausibly produce similar over-densities. These results provide dramatic support for the recent proposal by Xu et al. (2015) that, rather than stars accreted from other galaxies, the TriAnd clouds could represent stars kicked-out from our own disk. If so, these would be the first populations of disk stars to be found in the Galactic halo and a clear signature of the importance of this second formation mechanism for stellar halos more generally. Moreover, their existence at the very extremities of the disk places strong constraints on the nature of the interaction that formed them.

We run controlled N-body experiments to study the evolution of the dark matter (DM) halo profiles of dwarf galaxies driven by the accretion of DM substructures. Our initial conditions assume that supernova feedback erases the primordial DM cusps of haloes with $10^{9}-10^{10} \rm{M_{\odot}}$ at $z=0$. The orbits and masses of the infalling substructures are borrowed from the {\it Aquarius} simulations. Our experiments show that a fraction of haloes that undergo 1:3 down to 1:30 mergers are susceptible to reform a DM cusp by $z\approx 0$. Cusp regrowth is driven the accretion of DM substructures that are dense enough reach the central regions of the main halo before being tidally disrupted. The infall of substructures with a mass ratio above 1:6 on the mean of the reported mass-concentration relation systematically lead to cusp regrowth. Between 1:6 to 1:8, and 1:8 to 1:30 substructures need to be located one and two-sigma above the mean, respectively. The merging timescales of these dense, low-mass substructures is relatively long $(5-11 \rm{Gyrs})$, which may pose a timescale problem for the longevity of DM cores in dwarf galaxies. These results suggest that a certain level of scatter in the central density slopes of galactic haloes acted-on by feedback is to be expected given the stochastic mass accretion histories of low-mass haloes and the diverse star formation histories observed in the Local Group dwarves.

We present cosmological N-body resimulations of the assembly of the Brightest Cluster Galaxies (BCGs) in rich clusters. At $z=2$ we populate dark matter subhalos with self-gravitating stellar systems whose abundance and structure match observed high-redshift galaxies. By $z=0$, mergers have built much larger galaxies at cluster centre. Their dark matter density profiles are shallower than in corresponding dark-matter-only simulations, but their total mass density profiles (stars + dark matter) are quite similar. Differences are found only at radii where the effects of central black holes may be significant. Dark matter density slopes shallower than $\gamma=1.0$ occur for $r/r_{200} < 0.015$, close to the half-light radii of the BCGs. Our experiments support earlier suggestions that NFW-like profiles are an attractor for the hierarchical growth of structure in collisionless systems -- total mass density profiles asymptote to the solution found in dark-matter-only simulations over the radial range where mergers produce significant mixing between stars and dark matter. Simulated dark matter fractions are substantially higher in BCGs than in field ellipticals, reaching 80\% within the half-light radius. We also estimate that supermassive black hole mergers should create BCG cores as large as $r_{c}\sim 3 \,\mathrm{kpc}$. The good agreement of all these properties with recent observational studies of BCG structure suggests that dissipational processes have not played a dominant role in the assembly of the observed systems.

We generate stellar distribution functions (DFs) in triaxial haloes in order to examine the reliability of slopes $\Gamma\equiv \Delta {\rm log} M / \Delta {\rm log} r$ inferred by applying mass estimators of the form $M\propto R_e\sigma^2$ (i.e. assuming spherical symmetry, where $R_e$ and $\sigma$ are luminous effective radius and global velocity dispersion, respectively) to two stellar sub-populations independently tracing the same gravitational potential. The DFs take the form $f(E)$, are dynamically stable, and are generated within triaxial potentials corresponding directly to subhaloes formed in cosmological dark-matter-only simulations of Milky Way and galaxy cluster haloes. Additionally, we consider the effect of different tracer number density profiles (cuspy and cored) on the inferred slopes of mass profiles. For the isotropic DFs considered here, we find that halo triaxiality tends to introduce an anti-correlation between $R_e$ and $\sigma$ when estimated for a variety of viewing angles. The net effect is a negligible contribution to the systematic error associated with the slope of the mass profile, which continues to be dominated by a bias toward greater overestimation of masses for more-concentrated tracer populations. We demonstrate that simple mass estimates for two distinct tracer populations can give reliable (and cosmologically meaningful) lower limits for $\Gamma$, irrespective of the degree of triaxiality or shape of the tracer number density profile.

We study the formation and evolution of Brightest Cluster Galaxies starting from a $z=2$ population of quiescent ellipticals and following them to $z=0$. To this end, we use a suite of nine high-resolution dark matter-only simulations of galaxy clusters in a $\Lambda$CDM universe. We develop a scheme in which simulation particles are weighted to generate realistic and dynamically stable stellar density profiles at $z=2$. Our initial conditions assign a stellar mass to every identified dark halo as expected from abundance matching; assuming there exists a one-to-one relation between the visible properties of galaxies and their host haloes. We set the sizes of the luminous components according to the observed relations for $z\sim2$ massive quiescent galaxies. We study the evolution of the mass-size relation, the fate of satellite galaxies and the mass aggregation of the cluster central. From $z=2$, these galaxies grow on average in size by a factor 5 to 10 of and in mass by 2 to 3. The stellar mass growth rate of the simulated BCGs in our sample is of 1.9 in the range $0.3<z<1.0$ consistent with observations, and of 1.5 in the range $0.0<z<0.3$. Furthermore the satellite galaxies evolve to the present day mass-size relation by $z=0$. Assuming passively evolving stellar populations, we present surface brightness profiles for our cluster centrals which resemble those observed for the cDs in similar mass clusters both at $z=0$ and at $z=1$. This demonstrates that the $\Lambda$CDM cosmology does indeed predict minor and major mergers to occur in galaxy clusters with the frequency and mass ratio distribution required to explain the observed growth in size of passive galaxies since $z=2$. Our experiment shows that Brightest Cluster Galaxies can form through dissipationless mergers of quiescent massive $z=2$ galaxies, without substantial additional star formation.

We build a sample of distant (D > 80 kpc) stellar halo stars with measured radial velocities. Faint ($20 < g <22$) candidate blue horizontal branch (BHB) stars were selected using the deep, but wide, multi-epoch Sloan Digital Sky Survey photometry. Follow-up spectroscopy for these A-type stars was performed using the VLT-FORS2 instrument. We classify stars according to their Balmer line profiles, and find 7 are bona fide BHB stars and 31 are blue stragglers (BS). Owing to the magnitude range of our sample, even the intrinsically fainter BS stars can reach out to D ~ 90 kpc. We complement this sample of A-type stars with intrinsically brighter, intermediate-age, asymptotic giant branch stars. A set of 4 distant cool carbon stars is compiled from the literature and we perform spectroscopic follow-up on a further 4 N-type carbon stars using the WHT-ISIS instrument. Altogether, this provides us with the largest sample to date of individual star tracers out to r ~ 150 kpc. We find that the radial velocity dispersion of these tracers falls rapidly at large distances and is surprisingly cold (sigma_r ~ 50-60 km/s) between 100-150 kpc. Relating the measured radial velocities to the mass of the Milky Way requires knowledge of the (unknown) tracer density profile and anisotropy at these distances. Nonetheless, by assuming the stellar halo stars between 50-150 kpc have a moderate density fall-off (with power-law slope alpha < 5) and are on radial orbits (sigma^2_t/sigma^2_r < 1), we infer that the mass within 150 kpc is less than 10^12 M_solar and suggest it probably lies in the range (5-10) x 10^11 M_solar. We discuss the implications of such a low mass for the Milky Way.

We study the evolution of the stellar and dark matter components in a galaxy cluster of $10^{15} \, \rm{M_{\odot}}$ from $z=3$ to the present epoch using the high-resolution collisionless simulations of Ruszkowski & Springel (2009). At $z=3$ the dominant progenitor halos were populated with spherical model galaxies with and without accounting for adiabatic contraction. We apply a weighting scheme which allows us to change the relative amount of dark and stellar material assigned to each simulation particle in order to produce luminous properties which agree better with abundance matching arguments and observed bulge sizes at $z=3$. This permits the study of the effect of initial compactness on the evolution of the mass-size relation. We find that for more compact initial stellar distributions the size of the final Brightest Cluster Galaxy grows with mass according to $r\propto M^{2}$, whereas for more extended initial distributions, $r\propto M$. Our results show that collisionless mergers in a cosmological context can reduce the strength of inner dark matter cusps with changes in logarithmic slope of 0.3 to 0.5 at fixed radius. Shallow cusps such as those found recently in several strong lensing clusters thus do not necessarily conflict with CDM, but may rather reflect on the initial structure of the progenitor galaxies, which was shaped at high redshift by their formation process.