We present [Fe/H] and [$\alpha$/Fe] abundances, derived using spectral synthesis techniques, for stars in M31's outer stellar halo. The 21 [Fe/H] measurements and 7 [$\alpha$/Fe] measurements are drawn from fields ranging from 43 to 165 kpc in projected distance from M31. We combine our measurements with existing literature measurements, and compare the resulting sample of 23 stars with [Fe/H] and 9 stars with [$\alpha$/Fe] measurements in M31's outer halo with [$\alpha$/Fe] and [Fe/H] measurements, also derived from spectral synthesis, in M31's inner stellar halo ($r < $26 kpc) and dSph galaxies. The stars in M31's outer halo have [$\alpha$/Fe] patterns that are consistent with the largest of M31's dSph satellites (And I and And VII). These abundances provide tentative evidence that the [$\alpha$/Fe] abundances of stars in M31's outer halo are more similar to the abundances of Milky Way halo stars than to the abundances of stars in M31's inner halo. We also compare the spectral synthesis-based [Fe/H] measurements of stars in M31's halo with previous photometric [Fe/H] estimates, as a function of projected distance from M31. The spectral synthesis-based [Fe/H] measurements are consistent with a large-scale metallicity gradient previously observed in M31's stellar halo to projected distances as large as 100 kpc.

We present chemical abundances of red giant branch (RGB) stars in the dwarf spheroidal (dSph) satellite system of Andromeda (M31), using spectral synthesis of medium resolution (R $\sim 6000$) spectra obtained with the Keck II telescope and DEIMOS spectrograph via the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. We coadd stars according to their similarity in photometric metallicity or effective temperature to obtain a signal-to-noise ratio (S/N) high enough to measure average [Fe/H] and [$\alpha$/Fe] abundances. We validate our method using high S/N spectra of RGB stars in Milky Way globular clusters as well as deep observations for a subset of the M31 dSphs in our sample. For this set of validation coadds, we compare the weighted average abundance of the individual stars with the abundance determined from the coadd. We present individual and coadded measurements of [Fe/H] and [$\alpha$/Fe] for stars in ten M31 dSphs, including the first [$\alpha$/Fe] measurements for And IX, XIV, XV, and XVIII. These fainter, less massive dSphs show declining [$\alpha$/Fe] relative to [Fe/H], implying an extended star formation history. In addition, these dSphs also follow the same mass-metallicity relation found in other Local Group satellites. The conclusions we infer from coadded spectra agree with those from previous measurements in brighter M31 dSphs with individual abundance measurements, as well as conclusions from photometric studies. These abundances greatly increase the number of spectroscopic measurements of the chemical composition of M31's less massive dwarf satellites, which are crucial to understanding their star formation history and interaction with the M31 system.

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.

We present deep spectroscopy from Keck/DEIMOS of Andromeda I, III, V, VII, and X, all of which are dwarf spheroidal satellites of M31. The sample includes 256 spectroscopic members across all five dSphs. We confirm previous measurements of the velocity dispersions and dynamical masses, and we provide upper limits on bulk rotation. Our measurements confirm that M31 satellites obey the same relation between stellar mass and stellar metallicity as Milky Way (MW) satellites and other dwarf galaxies in the Local Group. The metallicity distributions show similar trends with stellar mass as MW satellites, including evidence in massive satellites for external influence, like pre-enrichment or gas accretion. We present the first measurements of individual element ratios, like [Si/Fe], in the M31 system, as well as measurements of the average [alpha/Fe] ratio. The trends of [alpha/Fe] with [Fe/H] also follow the same galaxy mass-dependent patterns as MW satellites. Less massive galaxies have more steeply declining slopes of [alpha/Fe] that begin at lower [Fe/H]. Finally, we compare the chemical evolution of M31 satellites to M31's Giant Stellar Stream and smooth halo. The properties of the M31 system support the theoretical prediction that the inner halo is composed primarily of massive galaxies that were accreted early. As a result, the inner halo exhibits higher [Fe/H] and [alpha/Fe] than surviving satellite galaxies.

We use cosmological hydrodynamical simulations of Milky-Way-mass galaxies from the FIRE project to evaluate various strategies for estimating the mass of a galaxy's stellar halo from deep, integrated-light images. We find good agreement with integrated-light observations if we mimic observational methods to measure the mass of the stellar halo by selecting regions of an image via projected radius relative to the disk scale length or by their surface density in stellar mass . However, these observational methods systematically underestimate the accreted stellar component, defined in our (and most) simulations as the mass of stars formed outside of the host galaxy, by up to a factor of ten, since the accreted component is centrally concentrated and therefore substantially obscured by the galactic disk. Furthermore, these observational methods introduce spurious dependencies of the estimated accreted stellar component on the stellar mass and size of galaxies that can obscure the trends in accreted stellar mass predicted by cosmological simulations, since we find that in our simulations the size and shape of the central galaxy is not strongly correlated with the assembly history of the accreted stellar halo. This effect persists whether galaxies are viewed edge-on or face-on. We show that metallicity or color information may provide a way to more cleanly delineate in observations the regions dominated by accreted stars. Absent additional data, we caution that estimates of the mass of the accreted stellar component from single-band images alone should be taken as lower limits.

We present the first measurements of [Fe/H] and [$\alpha$/Fe] abundances, obtained using spectral synthesis modeling, for red giant branch stars in M31's giant stellar stream. The spectroscopic observations, obtained at a projected distance of 17 kpc from M31's center, yielded 61 stars with [Fe/H] measurements, including 21 stars with [$\alpha$/Fe] measurements, from 112 targets identified as M31 stars. The [Fe/H] measurements confirm the expectation from photometric metallicity estimates that stars in this region of M31's halo are relatively metal-rich compared to stars in the MW's inner halo: more than half the stars in the field, including those not associated with kinematically identified substructure, have [Fe/H] abundances $> -1.0$. The stars in this field are $\alpha$-enhanced at lower metallicities, while [$\alpha$/Fe] decreases with increasing [Fe/H] above metallicities of [Fe/H] $\gtrsim -0.9$. Three kinematical components have been previously identified in this field: the giant stellar stream, a second kinematically cold feature of unknown origin, and M31's kinematically hot halo. We compare probabilistic [Fe/H] and [$\alpha$/Fe] distribution functions for each of the components. The giant stellar stream and the second kinematically cold feature have very similar abundance distributions, while the halo component is more metal-poor. Although the current sample sizes are small, a comparison of the abundances of stars in the giant stellar stream field with abundances of M31 halo and dSph stars from the literature indicate that the progenitor of the stream was likely more massive, and experienced a higher efficiency of star formation, than M31's existing dSphs or the dEs NGC147 and NGC185.

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.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

We investigate stellar metallicity distribution functions (MDFs), including Fe and ${\alpha}$-element abundances, in dwarf galaxies from the Feedback in Realistic Environments (FIRE) project. We examine both isolated dwarf galaxies and those that are satellites of a Milky Way-mass galaxy. In particular, we study the effects of including a sub-grid turbulent model for the diffusion of metals in gas. Simulations that include diffusion have narrower MDFs and abundance ratio distributions, because diffusion drives individual gas and star particles toward the average metallicity. This effect provides significantly better agreement with observed abundance distributions of dwarf galaxies in the Local Group, including the small intrinsic scatter in [${\alpha}$/Fe] vs. [Fe/H] (less than 0.1 dex). This small intrinsic scatter arises in our simulations because the interstellar medium (ISM) in dwarf galaxies is well-mixed at nearly all cosmic times, such that stars that form at a given time have similar abundances to within 0.1 dex. Thus, most of the scatter in abundances at z = 0 arises from redshift evolution and not from instantaneous scatter in the ISM. We find similar MDF widths and intrinsic scatter for satellite and isolated dwarf galaxies, which suggests that environmental effects play a minor role compared with internal chemical evolution in our simulations. Overall, with the inclusion of metal diffusion, our simulations reproduce abundance distribution widths of observed low-mass galaxies, enabling detailed studies of chemical evolution in galaxy formation.

The 32 Orionis group is a co-moving group of roughly 20 young (24 Myr) M3-B5 stars 100 pc from the Sun. Here we report the discovery of its first substellar member, WISE J052857.69+090104.2. This source was previously reported to be an M giant star based on its unusual near-infrared spectrum and lack of measurable proper motion. We re-analyze previous data and new moderate-resolution spectroscopy from Magellan/FIRE to demonstrate that this source is a young near-infrared L1 brown dwarf with very low surface gravity features. Spectral model fits indicate T$_{eff}$ = 1880$^{+150}_{-70}$ K and $\log{g}$ = 3.8$^{+0.2}_{-0.2}$ (cgs), consistent with a 15-22 Myr object with a mass near the deuterium-burning limit. Its sky position, estimated distance, kinematics (both proper motion and radial velocity), and spectral characteristics are all consistent with membership in 32 Orionis, and its temperature and age imply a mass (M = 14$^{+4}_{-3}$ M$_{Jup}$) that straddles the brown dwarf/planetary-mass object boundary. The source has a somewhat red $J-W2$ color compared to other L1 dwarfs, but this is likely a low-gravity-related temperature offset; we find no evidence of significant excess reddening from a disk or cool companion in the 3-5 $\mu$m waveband.

We present spectroscopic analysis of the low mass binary star system GJ 660.1AB, a pair of nearby M dwarfs for which we have obtained separated near-infrared spectra (0.9-2.5 $\mu$m) with the SpeX spectrograph. The spectrum of GJ 660.1B is distinctly peculiar, with a triangular-shaped 1.7 $\mu$m peak that initially suggests it to be a low surface gravity, young brown dwarf. However, we rule out this hypothesis and determine instead that this companion is a mild subdwarf (d/sdM7) based on the subsolar metallicity of the primary, [Fe/H] = $-$0.63$\pm$0.06. Comparison of the near-infrared spectrum of GJ 660.1B to two sets of spectral models yields conflicting results, with a common effective temperature T$_{eff}$ = 2550-2650 K, but alternately low surface gravity ($\log{g}$ = 4.4$\pm$0.5 cgs) and very low metallicity ([M/H] = $-$0.96$^{+0.19}_{-0.24}$), or high surface gravity ($\log{g}$ = 5.0-5.5) and slightly subsolar metallicity ([M/H] =$-$0.20$^{+0.13}_{-0.19}$). We conjecture that insufficient condensate opacity and excessive collision induced H$_2$ absorption in the models bias them toward low surface gravities and a metallicity inconsistent with the primary, and points toward improvements needed in the spectral modeling of metal-poor, very-low mass dwarfs. The peculiar spectral characteristics of GJ 660.1B emphasize that care is needed when interpreting surface gravity features in the spectra of ultracool dwarfs.

We have observed 6 late-L and T dwarfs with the Karl G. Jansky Very Large Array (VLA) to investigate the presence of highly circularly polarized radio emission, associated with large-scale auroral currents. Previous surveys encompassing ~60 L6 or later targets in this spectral range have yielded only one detection. Our sample includes the previously detected T6.5 dwarf 2MASS 10475385+2124234 as well as 5 new targets selected for the presence of H-alpha emission or optical/infrared photometric variability, which are possible manifestations of auroral activity. We detect 2MASS 10475385+2124234, as well as 4 of the 5 targets in our biased sample, including the strong IR variable SIMP J01365662+0933473 and bright H-alpha emitter 2MASS 12373919+6526148, reinforcing the possibility that activity at these disparate wavelengths is related. The radio emission frequency corresponds to a precise determination of the lower-bound magnetic field strength near the surface of each dwarf and this new sample provides robust constraints on dynamo theory in the low mass brown dwarf regime. Magnetic fields >2.5 kG are confirmed for 5/6 targets. Our results provide tentative evidence that the dynamo operating in this mass regime may be inconsistent with predicted values from a recently proposed model. Further observations at higher radio frequencies are essential for verifying this assertion.