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.

$Gaia$ DR2 has provided an unprecedented wealth of information about the positions and motions of stars in our Galaxy, and has highlighted the degree of disequilibria in the disc. As we collect data over a wider area of the disc it becomes increasingly appealing to start analysing stellar actions and angles, which specifically label orbit space, instead of their current phase space location. Conceptually, while $\bar{x}$ and $\bar{v}$ tell us about the potential and local interactions, grouping in action puts together stars that have similar frequencies and hence similar responses to dynamical effects occurring over several orbits. Grouping in actions and angles refines this further to isolate stars which are travelling together through space and hence have shared histories. Mixing these coordinate systems can confuse the interpretation. For example, it has been suggested that by moving stars to their guiding radius, the Milky Way spiral structure is visible as ridge-like overdensities in the $Gaia$ data \citep{Khoperskov+19b}. However, in this work, we show that these features are in fact the known kinematic moving groups, both in the $L_z-\phi$ and the $v_{\mathrm{R}}-v_{\phi}$ planes. Using simulations we show how this distinction will become even more important as we move to a global view of the Milky Way. As an example, we show that the radial velocity wave seen in the Galactic disc in $Gaia$ and APOGEE should become stronger in the action-angle frame, and that it can be reproduced by transient spiral structure.

We measured [Fe/H] and [$\alpha$/Fe] using spectral synthesis of low-resolution stellar spectroscopy for 70 individual red giant branch stars across four fields spanning the outer disk, Giant Stellar Stream (GSS), and inner halo of M31. Fields at M31-centric projected distances of 23 kpc in the halo, 12 kpc in the halo, 22 kpc in the GSS, and 26 kpc in the outer disk are $\alpha$-enhanced, with $\langle$[$\alpha$/Fe]$\rangle$ = 0.43, 0.50, 0.41, and 0.58, respectively. The 23 kpc and 12 kpc halo fields are relatively metal-poor, with $\langle$[Fe/H]$\rangle$ = $-$1.54 and $-$1.30, whereas the 22 kpc GSS and 26 kpc outer disk fields are relatively metal-rich with $\langle$[Fe/H]$\rangle$ = $-$0.84 and $-$0.92, respectively. For fields with substructure, we separated the stellar populations into kinematically hot stellar halo components and kinematically cold components. We did not find any evidence of an [$\alpha$/Fe] gradient along the high surface brightness core of the GSS between $\sim$17$-$22 kpc. However, we found tentative suggestions of a negative [$\alpha$/Fe] gradient in the stellar halo, which may indicate that different progenitor(s) or formation mechanisms contributed to the build up of the inner versus outer halo. Additionally, the [$\alpha$/Fe] distribution of the metal-rich ([Fe/H] $>$ $-$1.5), smooth inner stellar halo (r$_{\rm{proj}}$ $\lesssim$ 26 kpc) is inconsistent with having formed from the disruption of progenitor(s) similar to present-day M31 satellite galaxies. The 26 kpc outer disk is most likely associated with the extended disk of M31, where its high $\alpha$-enhancement provides support for an episode of rapid star formation in M31's disk, possibly induced by a major merger.

Measurements of [Fe/H] and [$\alpha$/Fe] can probe the minor merging history of a galaxy, providing a direct way to test the hierarchical assembly paradigm. While measurements of [$\alpha$/Fe] have been made in the stellar halo of the Milky Way, little is known about detailed chemical abundances in the stellar halo of M31. To make progress with existing telescopes, we apply spectral synthesis to low-resolution DEIMOS spectroscopy (R $\sim$ 2500 at 7000 Angstroms) across a wide spectral range (4500 Angstroms $<$ $\lambda$ $<$ 9100 Angstroms). By applying our technique to low-resolution spectra of 170 giant stars in 5 MW globular clusters, we demonstrate that our technique reproduces previous measurements from higher resolution spectroscopy. Based on the intrinsic dispersion in [Fe/H] and [$\alpha$/Fe] of individual stars in our combined cluster sample, we estimate systematic uncertainties of $\sim$0.11 dex and $\sim$0.09 dex in [Fe/H] and [$\alpha$/Fe], respectively. We apply our method to deep, low-resolution spectra of 11 red giant branch stars in the smooth halo of M31, resulting in higher signal-to-noise per spectral resolution element compared to DEIMOS medium-resolution spectroscopy, given the same exposure time and conditions. We find $\langle$[$\alpha$/Fe]$\rangle$ = 0.49 $\pm$ 0.29 dex and $\langle$[Fe/H]$\rangle$ = 1.59 $\pm$ 0.56 dex for our sample. This implies that---much like the Milky Way---the smooth halo of M31 is likely composed of disrupted dwarf galaxies with truncated star formation histories that were accreted early in the halo's formation.

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.

The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way (MW) halo main sequence turnoff stars in the CANDELS fields. In this paper, the second in the HALO7D series, we present the proper motions for the HALO7D sample. We discuss our measurement methodology, which makes use of a Bayesian mixture modeling approach for creating the stationary reference frame of distant galaxies. Using the 3D kinematic HALO7D sample, we estimate the parameters of the halo velocity ellipsoid, $\langle v_{\phi} \rangle, \sigma_r, \sigma_{\phi}, \sigma_{\theta}$, and the velocity anisotropy $\beta$. Using the full HALO7D sample, we find $\beta=0.63 \pm 0.05$ at $\langle r \rangle =24$ kpc. We also estimate the ellipsoid parameters for our sample split into three apparent magnitude bins; the posterior medians for these estimates of $\beta$, while consistent with one another, increase as a function of mean sample distance. Finally, we estimate $\beta$ in each of the individual HALO7D fields. We find that the velocity anisotropy $\beta$ can vary from field to field, which suggests that the halo is not phase mixed at $\langle r \rangle =24$ kpc. We explore the $\beta$ variation across the skies of two stellar halos from the \textit{Latte} suite of FIRE-2 simulations, finding that both simulated galaxies show $\beta$ variation over a similar range to the variation observed across the four HALO7D fields. The accretion histories of the two simulated galaxies result in different $\beta$ variation patterns; spatially mapping $\beta$ is thus a way forward in characterizing the accretion history of the Galaxy.

The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way halo main sequence turnoff stars in the CANDELS fields. In this paper, we present the spectroscopic component of this dataset, and discuss target selection, observing strategy, and survey properties. We present a new method of measuring line-of-sight (LOS) velocities by combining multiple spectroscopic observations of a given star, utilizing Bayesian hierarchical modeling. We present the LOS velocity distributions of the four HALO7D fields, and estimate their means and dispersions. All of the LOS distributions are dominated by the "hot halo": none of our fields are dominated by substructure that is kinematically cold in the LOS velocity component. Our estimates of the LOS velocity dispersions are consistent across the different fields, and these estimates are consistent with studies using other types of tracers. To complement our observations, we perform mock HALO7D surveys using the synthetic survey software Galaxia to "observe'" the Bullock & Johnston (2005) accreted stellar halos. Based on these simulated datasets, the consistent LOS velocity distributions across the four HALO7D fields indicates that the HALO7D sample is dominated by stars from the same massive (or few relatively massive) accretion event(s).

We study winds in 12 X-ray AGN host galaxies at z ~ 1. We find, using the low-ionization Fe II 2586 absorption in the stacked spectra, that the probability distribution function (PDF) of the centroid velocity shift in AGN has a median, 16th and 84th percentiles of (-87, -251, +86) km/s respectively. The PDF of the velocity dispersion in AGN has a median, 84th and 16th percentile of (139, 253, 52) km/s respectively. The centroid velocity and the velocity dispersions are obtained from a two component (ISM+wind) absorption line model. The equivalent width PDF of the outflow in AGN has median, 84th and 16th percentiles of (0.4, 0.8, 0.1) Angstrom. There is a strong ISM component in Fe II 2586 absorption with (1.2, 1.5, 0.8) Angstrom, implying presence of substantial amount cold gas in the host galaxies. For comparison, star-forming and X-ray undetected galaxies at a similar redshift, matched roughly in stellar mass and galaxy inclination, have a centroid velocity PDF with percentiles of (-74, -258, +90) km/s, and a velocity dispersion PDF percentiles of (150, 259, 57) km/s. Thus, winds in the AGN are similar to star-formation-driven winds, and are too weak to escape and expel substantial cool gas from galaxies. Our sample doubles the previous sample of AGN studied at z ~ 0.5 and extends the analysis to z ~ 1. A joint reanalysis of the z ~ 0.5 AGN sample and our sample yields consistent results to the measurements above.

The dwarf carbon (dC) star SDSS J112801.67+004034.6 has an unusually high radial velocity, 531$\pm 4$ km s$^{-1}$. We present proper motion and new spectroscopic observations which imply a large Galactic rest frame velocity, 425$\pm 9$ km s$^{-1}$. Several other SDSS dC stars are also inferred to have very high galactocentric velocities, again each based on both high heliocentric radial velocity and also confidently detected proper motions. Extreme velocities and the presence of $C_2$ bands in the spectra of dwarf stars are both rare. Passage near the Galactic center can accelerate stars to such extreme velocities, but the large orbital angular momentum of SDSS J1128 precludes this explanation. Ejection from a supernova in a binary system or disruption of a binary by other stars are possibilities, particularly as dC stars are thought to obtain their photospheric $C_2$ via mass transfer from an evolved companion.

We present the line-of-sight (LOS) velocities for 13 distant main sequence Milky Way halo stars with published proper motions. The proper motions were measured using long baseline (5-7 years) multi-epoch HST/ACS photometry, and the LOS velocities were extracted from deep (5-6 hour integrations) Keck II/DEIMOS spectra. We estimate the parameters of the velocity ellipsoid of the stellar halo using a Markov chain Monte Carlo ensembler sampler method. The velocity second moments in the directions of the Galactic $(l,b,$ LOS) coordinate system are $\langle v^2_l \rangle^{1/2} = 138^{+43}_{-26}$ km/s, $\langle v^2_b \rangle^{1/2} = 88^{+28}_{-17}$ km/s, and $\langle v^2_{\rm{LOS}} \rangle^{1/2} = 91^{+27}_{-14}$ km/s. We use these ellipsoid parameters to constrain the velocity anisotropy of the stellar halo. Ours is the first measurement of the anisotropy parameter $\beta$ using 3D kinematics outside of the solar neighborhood. We find $\beta=-0.3^{+0.4}_{-0.9}$, consistent with isotropy and lower than solar neighborhood $\beta$ measurements by 2$\sigma$ ($\beta_{SN} \sim 0.5-0.7$). We identify two stars in our sample that are likely members of the known TriAnd substructure, and excluding these objects from our sample increases our estimate of the anisotropy to $\beta=0.1^{+0.4}_{-1.0}$, which is still lower than solar neighborhood measurements by $1\sigma$. The potential decrease in $\beta$ with Galactocentric radius is inconsistent with theoretical predictions, though consistent with recent observational studies, and may indicate the presence of large, shell-type structure (or structures) at $r \sim 25$ kpc. The methods described in this paper will be applied to a much larger sample of stars with 3D kinematics observed through the ongoing HALO7D program.

We report on the properties of the most massive ultra-compact dwarf galaxy (UCD) in the nearby Virgo Cluster of galaxies using imaging from the Next Generation Virgo Cluster Survey (NGVS) and spectroscopy from Keck/DEIMOS. This object (M59-UCD3) appears to be associated with the massive Virgo galaxy M59 (NGC 4621), has an integrated velocity dispersion of 78 km/s, a dynamical mass of $3.7\times10^8 M_\odot$, and an effective radius ($R_e$) of 25 pc. With an effective surface mass density of $9.4\times10^{10} M_\odot/kpc^2$, it is the densest galaxy in the local Universe discovered to date, surpassing the density of the luminous Virgo UCD, M60-UCD1. M59-UCD3 has a total luminosity of $M_{g'}=-14.2$ mag, and a spectral energy distribution consistent with an old (14 Gyr) stellar population with [Fe/H]=0.0 and [$\alpha$/Fe]=+0.2. We also examine deep imaging around M59 and find a broad low surface brightness stream pointing towards M59-UCD3, which may represent a tidal remnant of the UCD progenitor. This UCD, along with similar objects like M60-UCD1 and M59cO, likely represents an extreme population of tidally stripped galaxies more akin to larger and more massive compact early-type galaxies than to nuclear star clusters in present-day dwarf galaxies.

We explore the Triangulum-Andromeda (TriAnd) overdensity in the SPLASH (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo) and SEGUE (the Sloan Extension for Galactic Understanding and Exploration) spectroscopic surveys. Milky Way main sequence turn-off stars in the SPLASH survey reveal that the TriAnd overdensity and the recently discovered PAndAS stream (Martin et al. 2014) share a common heliocentric distance (D ~ 20 kpc), position on the sky, and line-of-sight velocity (V_GSR ~ 50 km/s). Similarly, A-type, giant, and main sequence turn-off stars selected from the SEGUE survey in the vicinity of the Segue 2 satellite show that TriAnd is prevalent in these fields, with a velocity and distance similar to Segue 2. The coincidence of the PAndAS stream and Segue 2 satellite in positional and velocity space to TriAnd suggests that these substructures are all associated, and may be a fossil record of group-infall onto the Milky Way halo. In this scenario, the Segue 2 satellite and PAndAS stream are "satellites of satellites", and the large, metal-rich TriAnd overdensity is the remains of the group central.

We present an extensive search for RR Lyrae stars in and around the ultra-faint Milky Way companions Segue 2 and Segue 3. The former (M_V = -2.5, Belokurov et al. 2009) appears to be an extremely faint dwarf galaxy companion of the Milky Way. The latter (M_V = 0.0, Fadely et al. 2011) is among the faintest star clusters known. We use B and V band time-series imaging obtained at the WIYN 0.9 meter telescope at Kitt Peak National Observatory to search for RR Lyrae in these objects. In our Segue 2 observations, we present a previously unknown fundamental mode (RRab) RR Lyrae star with a period of P_ab = 0.748 days. With this measurement, we revisit the inverse correlation between <P_ab> and <[Fe/H]> established in the literature for Milky Way dwarf galaxies and their RR Lyrae. In this context, the long period of Segue 2's RRab star as well as the known significant spread in metallicity in this dwarf galaxy are consistent with the observed trend in <P_ab> and <[Fe/H]>. We derive the first robust distance to Segue 2, using both its RRab star and spectroscopically confirmed blue horizontal branch stars. Using [Fe/H] = -2.16 and -2.44 dex, we find d_RRL = 36.6 +2.5/-2.4 and 37.7 +/- 2.7 kpc; assuming [Fe/H] = -2.257 dex, we find d_BHB = 34.4 +/- 2.6 kpc. Although no RR Lyrae were present in the Segue 3 field, we found a candidate eclipsing binary star system.