We present high resolution spectroscopy taken with the Keck Echellete Spectrograph and Imager to measure stellar velocity dispersions for eight active dwarf galaxies ($M_{\ast}<3\times10^{9}~M_{\odot}$) with virial black hole masses. We double the number of systems in this stellar mass regime with measurements of both black hole mass ($M_{\rm BH}$) and stellar velocity dispersion ($\sigma_{\ast}$), and place them on the $M_{\rm BH}-\sigma_{\ast}$ relation. The tight relation between $M_{\rm BH}$ and $\sigma_{\ast}$ for higher mass galaxies is a strong piece of evidence for the co-evolution of BHs and their host galaxies, but it has been unclear whether this relation holds in the dwarf galaxy regime. Our sample is in good agreement with the extrapolation of the $M_{\rm BH}-\sigma_{\ast}$ relation to low BH/galaxy masses, suggesting that the processes which produce $M_{\rm BH}-\sigma_{\ast}$ can also operate in dwarf galaxies. These results provide important constraints for massive black hole seed formation models and models exploring the impact of AGN feedback in dwarf galaxies.

We present a detailed abundance analysis of the three brightest member stars at the top of the giant branch of the ultra-faint dwarf galaxy Grus~II. All stars exhibit a higher than expected $\mathrm{[Mg/Ca]}$ ratio compared to metal-poor stars in other ultra-faint dwarf galaxies and in the Milky Way halo. Nucleosynthesis in high mass ($\geqslant 20$M$_\odot$) core-collapse supernovae has been shown to create this signature. The abundances of this small sample (3) stars suggest the chemical enrichment of Grus~II could have occurred through substantial high-mass stellar evolution and is consistent with the framework of a top-heavy initial mass function. However, with only three stars it can not be ruled out that the abundance pattern is the result of a stochastic chemical enrichment at early times in the galaxy. The most metal-rich of the three stars also possesses a small enhancement in rapid neutron-capture ($r$-process) elements. The abundance pattern of the $r$-process elements in this star matches the scaled $r$-process pattern of the solar system and $r$-process enhanced stars in other dwarf galaxies and in the Milky Way halo, hinting at a common origin for these elements across a range of environments. All current proposed astrophysical sites of $r$-process element production are associated with high-mass stars, thus the possible top-heavy initial mass function of Grus~II would increase the likelihood of any of these events occurring. The time delay between the $\alpha$ and $r$-process element enrichment of the galaxy favors a neutron star merger as the origin of the $r$-process elements in Grus~II.

We present a Bayesian method to identify multiple (chemodynamic) stellar populations in dwarf spheroidal galaxies (dSphs) using velocity, metallicity, and positional stellar data without the assumption of spherical symmetry. We apply this method to a new Keck/DEIMOS spectroscopic survey of the Ursa Minor (UMi) dSph. We identify 892 likely members, making this the largest UMi sample with line-of-sight velocity and metallicity measurements. Our Bayesian method detects two distinct chemodynamic populations with high significance ($\ln{B}\sim33$). The metal-rich ($[{\rm Fe/H}]=-2.05\pm0.03$) population is kinematically colder (radial velocity dispersion of $\sigma_v=4.9\pm0.8 \, {\rm km \, s^{-1}}$) and more centrally concentrated than the metal-poor ($[{\rm Fe/H}]=-2.29\pm0.05$) and kinematically hotter population ($\sigma_v =11.5\pm0.9\, {\rm km \, s^{-1}}$). Furthermore, we apply the same analysis to an independent MMT/Hectochelle data set and confirm the existence of two chemodynamic populations in UMi. In both data sets, the metal-rich population is significantly flattened ($\epsilon=0.75\pm0.03$) and the metal-poor population is closer to spherical ($\epsilon=0.33_{-0.09}^{+0.12}$). Despite the presence of two populations, we are unable to robustly estimate the slope of the dynamical mass profile. We found hints for prolate rotation of order $\sim 2 \, {\rm km \, s^{-1}}$ in the MMT data set, but further observations are required to verify this. The flattened metal-rich population invalidates assumptions built into simple dynamical mass estimators, so we computed new astrophysical dark matter annihilation (J) and decay profiles based on the rounder, hotter metal-poor population and inferred $\log_{10}{(J(0.5^{\circ})/{\rm GeV^{2} \, cm^{-5}})}\approx19.1$ for the Keck data set. Our results paint a more complex picture of the evolution of Ursa Minor than previously discussed.

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

We introduce the Southern Stellar Stream Spectroscopy Survey (${S}^5$), an on-going program to map the kinematics and chemistry of stellar streams in the Southern Hemisphere. The initial focus of ${S}^5$ has been spectroscopic observations of recently identified streams within the footprint of the Dark Energy Survey (DES), with the eventual goal of surveying streams across the entire southern sky. Stellar streams are composed of material that has been tidally striped from dwarf galaxies and globular clusters and hence are excellent dynamical probes of the gravitational potential of the Milky Way, as well as providing a detailed snapshot of its accretion history. Observing with the 3.9-m Anglo-Australian Telescope's 2-degree-Field fibre positioner and AAOmega spectrograph, and combining the precise photometry of DES DR1 with the superb proper motions from $Gaia$ DR2, allows us to conduct an efficient spectroscopic survey to map these stellar streams. So far ${S}^5$ has mapped 9 DES streams and 3 streams outside of DES; the former are the first spectroscopic observations of these recently discovered streams. In addition to the stream survey, we use spare fibres to undertake a Milky Way halo survey and a low-redshift galaxy survey. This paper presents an overview of the ${S}^5$ program, describing the scientific motivation for the survey, target selection, observation strategy, data reduction and survey validation. Finally, we describe early science results on stellar streams and Milky Way halo stars drawn from the survey. Updates on ${S}^5$, including future public data release, can be found at \url{http://s5collab.github.io}.

We present an analysis of the long-term optical variability for $\sim50,000$ nearby (z<0.055) galaxies from the NASA-Sloan Atlas, $35,000$ of which are low-mass ($M_{\ast}<10^{10}~M_{\odot}$). We use difference imaging of Palomar Transient Factory (PTF) R-band observations to construct light curves with typical baselines of several years. We then search for subtle variations in the nuclear light output. We determine whether detected variability is AGN-like by assessing the fit quality to a damped random walk model. We identify 424 variability-selected AGN, including 244 with stellar masses between $10^{7}$ and $10^{10}~M_{\odot}$. 75% of low-mass galaxies with AGN-like variability have narrow emission lines dominated by star formation. After controlling for nucleus magnitude, the fraction of variable AGN is constant down to $M_{\ast}=10^{9}~M_{\odot}$, suggesting no drastic decline in the BH occupation fraction down to this stellar mass regime. Combining our NASA-Sloan Atlas sample with samples of nearby galaxies with broad H$\alpha$ emission, we find no dependence of variability properties with black hole mass. However, we caution that the variable AGN fraction is strongly dependent on baseline. For baselines less than two years, the variable fraction for the full sample is 0.25%, compared to 1.0% for baselines longer than two years. Finally, comparing Stripe 82 light curves (Baldassare et al. 2018) to PTF light curves, we find populations of changing-look AGN: 8 galaxies that are variable in Stripe 82, but quiescent in PTF, and 15 galaxies where the reverse is true. Our PTF work demonstrates the promise of long-term optical variability searches in low-mass galaxies for finding AGNs missed by other selection techniques.

Stars that escape globular clusters form tidal tails that are predominantly shaped by the global distribution of mass in the Galaxy, but also preserve a historical record of small-scale perturbations. Using deep $grz$ photometry from DECaLS, we present highly probable members of the tidal tails associated with the disrupting globular cluster Palomar 5. These data yield the cleanest view of a stellar stream beyond $\sim20\,\rm kpc$ and reveal: (1) a wide, low surface-brightness extension of the leading tail; (2) significant density variations along the stream; and (3) sharp changes in the direction of both the leading and the trailing tail. In the fiducial Milky Way model, a rotating bar perturbs the Palomar 5 tails and can produce streams with similar width and density profiles to those observed. However, the deviations of the stream track in this simple model do not match those observed in the Palomar 5 trailing tail, indicating the need for an additional source of perturbation. These discoveries open up the possibility of measuring the population of perturbers in the Milky Way, including dark-matter subhalos, with an ensemble of stellar streams and deep photometry alone.

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.

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

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.

Nearby dwarf galaxies are local analogues of high-redshift and metal-poor stellar populations. Most of these systems ceased star formation long ago, but they retain signatures of their past that can be unraveled by detailed study of their resolved stars. Archaeological examination of dwarf galaxies with resolved stellar spectroscopy provides key insights into the first stars and galaxies, galaxy formation in the smallest dark matter halos, stellar populations in the metal-free and metal-poor universe, the nature of the first stellar explosions, and the origin of the elements. Extremely large telescopes with multi-object R=5,000-30,000 spectroscopy are needed to enable such studies for galaxies of different luminosities throughout the Local Group.

For nearly 40 years, dark matter has been widely assumed to be cold and collisionless. Cold dark matter models make fundamental predictions for the behavior of dark matter on small (<10 kpc) scales. These predictions include cuspy density profiles at the centers of dark matter halos and a halo mass function that increases as dN/dM ~ M^-1.9 down to very small masses. We suggest two observational programs relying on extremely large telescopes to critically test these predictions, and thus shed new light on the nature of dark matter. (1) Combining adaptive optics-enabled imaging with deep spectroscopy to measure the three-dimensional motions of stars within a sample of Local Group dwarf galaxies that are the cleanest dark matter laboratories known in the nearby universe. From these observations the inner slope of the dark matter density profile can be determined with an accuracy of 0.20 dex, enabling a central cusp to be distinguished from a core at 5 sigma significance. (2) Diffraction-limited AO imaging and integral field spectroscopy of gravitationally lensed galaxies and quasars to quantify the abundance of dark substructures in the halos of the lens galaxies and along the line of sight. Observations of 50 lensed arcs and 50 multiply-imaged quasars will be sufficient to measure the halo mass function over the range 10^7 < M < 10^10 Msun at cosmological scales, independent of the baryonic and stellar composition of those structures. These two observational probes provide complementary information about the small scale structure, with a joint self-consistent analysis mitigating limitations of either probe. This program will produce the strongest existing constraints on the properties of dark matter on small scales, allowing conclusive tests of alternative warm, fuzzy, and self-interacting dark matter models.

The 2020s are poised to continue the past two decades of significant advances based on observations of dwarf galaxies in the nearby universe. Upcoming wide-field photometric surveys will probe substantially deeper than previous data sets, pushing the discovery frontier for new dwarf galaxies to fainter magnitudes, lower surface brightnesses, and larger distances. These dwarfs will be compelling targets for testing models of galaxy formation and cosmology, including the properties of dark matter and possible modifications to gravity. However, most of the science that can be extracted from nearby dwarf galaxies relies on spectroscopy with large telescopes. We suggest that maximizing the scientific impact of near-future imaging surveys will require both major spectroscopic surveys on 6-10m telescopes and multiplexed spectroscopy with even larger apertures.

In this work we use structural properties of Milky Way's outer halo ($R_G > 25\,\mathrm{kpc}$) satellites (dwarf spheroidal galaxies, ultra-faint dwarf galaxies and globular clusters) derived from deep, wide-field and homogeneous data, to present evidence of a correlation in the S\'ersic index v/s effective radius plane followed by a large fraction of outer halo globular clusters and satellite dwarf galaxies. We show that this correlation can be entirely reproduced by fitting empirical relations in the central surface brightness v/s absolute magnitude and S\'ersic index v/s absolute magnitude parameter spaces, and by assuming the existence of two types of outer halo globular clusters: one of high surface brightness (HSB group), with properties similar to inner halo clusters; and another of low surface brightness (LSB group), which share characteristics with dwarf spheroidal and ultra-faint dwarf galaxies. Given the similarities of LSB clusters with dwarf spheroidal and ultra-faint dwarf galaxies, we discuss the possibility that outer halo clusters also originated inside dark matter halos and that tidal forces from different galaxy host's potentials are responsible for the different properties between HSB and LSB clusters.

We investigate the lowest-mass quiescent galaxies known to exist in isolated environments ($\mathrm{M^* = 10^{9.0-9.5} \ M_\odot}$; 1.5 Mpc from a more massive galaxy). This population may represent the lowest stellar mass galaxies in which internal feedback quenches galaxy-wide star formation. We present Keck/ESI long-slit spectroscopy for 27 isolated galaxies in this regime: 20 quiescent galaxies and 7 star-forming galaxies. We measure emission line strengths as a function of radius and place galaxies on the Baldwin Phillips Terlevich (BPT) diagram. Remarkably, 16 of 20 quiescent galaxies in our sample host central AGN-like line ratios. Only 5 of these quiescent galaxies were identified as AGN-like in SDSS due to lower spatial resolution and signal-to-noise. We find that many of the quiescent galaxies in our sample have spatially-extended emission across the non-SF regions of BPT-space. When considering only the central 1$^{\prime\prime}$, we identify a tight relationship between distance from the BPT star-forming sequence and host galaxy stellar age as traced by $\mathrm{D_n4000}$, such that older stellar ages are associated with larger distances from the star-forming locus. Our results suggest that the presence of hard ionizing radiation (AGN-like line ratios) is intrinsically tied to the quenching of what may be the lowest-mass self-quenched galaxies.

We present a new measurement of the systemic proper motion of the Small Magellanic Cloud (SMC), based on an expanded set of 30 fields containing background quasars and spanning a $\sim$3 year baseline, using the \textit{Hubble Space Telescope} (\textit{HST}) Wide Field Camera 3. Combining this data with our previous 5 \textit{HST} fields, and an additional 8 measurements from the \textit{Gaia}-Tycho Astrometric Solution Catalog, brings us to a total of 43 SMC fields. We measure a systemic motion of $\mu_{W}$ = $-0.82$ $\pm$ 0.02 (random) $\pm$ 0.10 (systematic) mas yr$^{-1}$ and $\mu_{N}$ = $-1.21$ $\pm$ 0.01 (random) $\pm$ 0.03 (systematic) mas yr$^{-1}$. After subtraction of the systemic motion, we find little evidence for rotation, but find an ordered mean motion radially away from the SMC in the outer regions of the galaxy, indicating that the SMC is in the process of tidal disruption. We model the past interactions of the Clouds with each other based on the measured present-day relative velocity between them of $103 \pm 26$ km s$^{-1}$. We find that in 97\% of our considered cases, the Clouds experienced a direct collision $147 \pm 33$ Myr ago, with a mean impact parameter of $7.5 \pm 2.5$ kpc.

A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models.

We present an analysis of the nuclear variability of $\sim28,000$ nearby ($z<0.15$) galaxies with Sloan Digital Sky Survey (SDSS) spectroscopy in Stripe 82. We construct light curves using difference imaging of SDSS g-band images, which allows us to detect subtle variations in the central light output. We select variable AGN by assessing whether detected variability is well-described by a damped random walk model. We find 135 galaxies with AGN-like nuclear variability. While most of the variability-selected AGNs have narrow emission lines consistent with the presence of an AGN, a small fraction have narrow emission lines dominated by star formation. The star-forming systems with nuclear AGN-like variability tend to be low-mass ($M_{\ast}<10^{10}~M_{\odot}$), and may be AGNs missed by other selection techniques due to star formation dilution or low-metallicities. We explore the AGN fraction as a function of stellar mass, and find that the fraction of variable AGN increases with stellar mass, even after taking into account the fact that lower mass systems are fainter. There are several possible explanations for an observed decline in the fraction of variable AGN with decreasing stellar mass, including a drop in the supermassive black hole occupation fraction, a decrease in the ratio of black hole mass to galaxy stellar mass, or a change in the variability properties of lower-mass AGNs. We demonstrate that optical photometric variability is a promising avenue for detecting AGNs in low-mass, star formation-dominated galaxies, which has implications for the upcoming Large Synoptic Survey Telescope.

We present a detailed study of the faint Milky Way satellite Draco II (Dra II) from deep CFHT/MegaCam broadband $g$ and $i$ photometry and narrow-band metallicity-sensitive CaHK observations, along with follow-up Keck II/DEIMOS multi-object spectroscopy. Forward modeling of the deep photometry allows us to refine the structural and photometric properties of Dra II: the distribution of stars in colour-magnitude space implies Dra II is old (13.5 $ \pm 0.5 $ Gyr), very metal poor, very faint ($ L_V = 180 ^{+124}_{-72} L_\odot $), and at a distance $d = 21.5 \pm 0.4$ kpc. The narrow-band, metallicity-sensitive CaHK Pristine photometry confirms this very low metallicity ([Fe/H]$ = -2.7 \pm 0.1$ dex). Even though our study benefits from a doubling of the spectroscopic sample size compared to previous investigations, the velocity dispersion of the system is still only marginally resolved ($\sigma_{vr}<5.9$ km s$^{-1}$ at the 95 per cent confidence level) and confirms that Dra II is a dynamically cold stellar system with a large recessional velocity ($\langle v_{r}\rangle = -342.5^{+1.1}_{-1.2} $ km s$^{-1}$). We further show that the spectroscopically confirmed members of Dra~II have a mean proper motion of $(\mu_\alpha^*,\mu_\delta)=(1.26 \pm 0.27,0.94 \pm 0.28) $ mas yr$^{-1}$ in the Gaia DR2 data, which translates to an orbit with a pericenter and an apocenter of $21.3 ^{+0.7}_{-1.0}$ and $153.8 ^{+56.7}_{-34.7}$ kpc, respectively. Taken altogether, these properties favour the scenario of Dra~II being a potentially disrupting dwarf galaxy. The low-significance extra-tidal features we map around the satellite tentatively support this scenario.

We use deep $HST$ WFC3/IR imaging to study the Initial Mass Function (IMF) of the ultra faint dwarf galaxy Coma Berenices (Com Ber). Our observations reach the lowest stellar mass ever probed in a resolved galaxy, with 50\% completeness at $\sim 0.17$ M$_{\odot}$. Unresolved background galaxies however limit our purity below $\sim 0.23$ M$_{\odot}$. If modeled with a single power law, we find that the IMF slope is $-1.45^{+0.29}_{-0.3}$ (68\% credible intervals), compared to a Milky Way value of $-2.3$. For a broken power law, we obtain a low-mass slope of $-1.18_{-0.33}^{+0.49}$, a high-mass slope of $-1.88_{-0.49}^{+0.43}$ and a break mass of $0.57_{-0.08}^{+0.12}$ M$_{\odot}$, compared to $-1.3$, $-2.3$ and 0.5 M$_{\odot}$ for a Kroupa IMF. For a log-normal IMF model we obtain values of $0.33_{-0.16}^{+0.15}$ M$_{\odot}$ for the location parameter and of $0.68_{-0.12}^{+0.17}$ for $\sigma$ (0.22 M$_{\odot}$ and 0.57 for the Chabrier system IMF). All three parametrizations produce similar agreement with the data. Our results agree with previous analysis of shallower optical HST data. However analysis of similar optical data of other dwarfs finds IMFs significantly more bottom-light than in the Milky Way. These results suggest two, non mutually exclusive, possibilities: that the discrepancy of the dwarf galaxies IMF with respect to the Milky Way is, at least partly, an artifact of using a single power law model, and that there is real variance in the IMF at low masses between the currently studied nearby dwarfs, with Com Ber being similar to the Milky Way, but other dwarfs differing significantly.

We present structural parameters from a wide-field homogeneous imaging survey of Milky Way satellites carried out with the MegaCam imagers on the 3.6m Canada-France-Hawaii Telescope (CFHT) and 6.5m Magellan-Clay telescope. Our survey targets an unbiased sample of "outer halo" satellites (i.e., substructures having Galactocentric distances greater than 25 kpc) and includes classical dSph galaxies, ultra-faint dwarfs, and remote globular clusters. We combine deep, panoramic $gr$ imaging for 44 satellites and archival $gr$ imaging for 14 additional objects (primarily obtained with the DECam instrument as part of the Dark Energy Survey), to measure photometric and structural parameters for 58 outer halo satellites. This is the largest and most uniform analysis of Milky Way satellites undertaken to date and represents roughly three quarters ($58/81 \simeq$72\%) of all known outer halo satellites. We use a maximum-likelihood method to fit four density laws to each object in our survey: exponential, Plummer, King and Sersic models. We examine systematically the isodensity contour maps and color magnitude diagrams for each of our program objects, present a comparison with previous results, and tabulate our best-fit photometric and structural parameters, including ellipticities, position angles, effective radii, Sersic indices, absolute magnitudes, and surface brightness measurements. We investigate the distribution of outer halo satellites in the size-magnitude diagram, and show that the current sample of outer halo substructures spans a wide range in effective radius, luminosity and surface brightness, with little evidence for a clean separation into star cluster and galaxy populations at the faintest luminosities and surface brightnesses.

We describe a deep, systematic imaging study of satellites in the outer halo of the Milky Way. Our sample consists of 58 stellar overdensities --- i.e., substructures classified as either globular clusters, classical dwarf galaxies, or ultra-faint dwarf galaxies --- that are located at Galactocentric distances of R$_{\rm GC}$ > 25 kpc (outer halo) and out to ~400 kpc. This includes 44 objects for which we have acquired deep, wide-field, $g-$ and $r-$band imaging with the MegaCam mosaic cameras on the 3.6m Canada-France-Hawaii Telescope and the 6.5m Magellan-Clay telescope. These data are supplemented by archival imaging, or published $gr$ photometry, for an additional 14 objects, most of which were discovered recently in the Dark Energy Survey (DES). We describe the scientific motivation for our survey, including sample selection, observing strategy, data reduction pipeline, calibration procedures, and the depth and precision of the photometry. The typical 5$\sigma$ point-source limiting magnitudes for our MegaCam imaging --- which collectively covers an area of ~52 deg$^{2}$ --- are $g_{\rm lim}$ ~25.6 and $r_{\rm lim}$ ~25.3 AB mag. These limits are comparable to those from the coadded DES images and are roughly a half-magnitude deeper than will be reached in a single visit with LSST. Our photometric catalog thus provides the deepest and most uniform photometric database of Milky Way satellites available for the foreseeable future. In other papers in this series, we have used these data to explore the blue straggler populations in these objects, their density distributions, star formation histories, scaling relations and possible foreground structures.

We investigate the global neutral hydrogen (HI) content of isolated galaxies selected from the SDSS spectroscopic survey with optical evidence of Active Galactic Nuclei (AGN). Our sample includes galaxies with unresolved HI observations from the ALFALFA 70% data release combined with deeper HI observations of low-mass galaxies with 7.0 < log(M*) < 9.5. We examine the HI masses of this sample using the distance from the star-forming sequence on the OIII\Hb and NII\Ha Baldwin Phillips Terlevich (BPT) diagram as a measurement of AGN activity. In agreement with previous studies, we find that, for galaxies with log(M*) > 9.5, AGN activity does not correlate with the global HI content. However, for galaxies with 9.2 < log(M*) < 9.5, we identify a set of objects at large distances from the BPT star-forming sequence and lower than expected HI masses. This gas-depleted sample is red in both g-r and NUV-r colors and compact without distinguishable signs of star formation. This is surprising because the vast majority of isolated galaxies in this stellar mass regime are both star-forming and gas-rich. These galaxies are greater than 1.5 Mpc from any massive galaxy, ruling out environmental processes as a source of the gas-depletion. We suggest that either black hole feedback or shocks from extremely bursty star formation cause the emission lines and have destroyed or otherwise consumed the cold gas.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

The shape of a galaxy's spatially unresolved, globally integrated 21-cm emission line depends on its internal gas kinematics: galaxies with rotation-supported gas disks produce double-horned profiles with steep wings, while galaxies with dispersion-supported gas produce Gaussian-like profiles with sloped wings. Using mock observations of simulated galaxies from the FIRE project, we show that one can therefore constrain a galaxy's gas kinematics from its unresolved 21-cm line profile. In particular, we find that the kurtosis of the 21-cm line increases with decreasing $V/\sigma$, and that this trend is robust across a wide range of masses, signal-to-noise ratios, and inclinations. We then quantify the shapes of 21-cm line profiles from a morphologically unbiased sample of $\sim$2000 low-redshift, HI-detected galaxies with $M_{\rm star} = 10^{7-11} M_{\odot}$ and compare to the simulated galaxies. At $M_{\rm star} \gtrsim 10^{10} M_{\odot}$, both the observed and simulated galaxies produce double-horned profiles with low kurtosis and steep wings, consistent with rotation-supported disks. Both the observed and simulated line profiles become more Gaussian-like (higher kurtosis and less-steep wings) at lower masses, indicating increased dispersion support. However, the simulated galaxies transition from rotation to dispersion support more strongly: at $M_{\rm star} = 10^{8-10}M_{\odot}$, most of the simulations produce more Gaussian-like profiles than typical observed galaxies with similar mass, indicating that gas in the low-mass simulated galaxies is, on average, overly dispersion-supported. Most of the lower-mass simulated galaxies also have somewhat lower gas fractions than the median of the observed population. The simulations nevertheless reproduce the observed line-width baryonic Tully-Fisher relation, which is insensitive to rotation vs. dispersion support.

Using deep HST/ACS observations, we demonstrate that the sub-solar stellar initial mass function (IMF) of 6 ultra-faint dwarf Milky Way Satellites (UFDs) is more bottom light than the IMF of the Milky Way disk. Our data have a lower mass limit of about 0.45 M$_{\odot}$, while the upper limit is $\sim 0.8$ M$_\odot$, set by the turn-off mass of these old, metal poor systems. If formulated as a single power law, we obtain a shallower IMF slope than the "Salpeter" value of $-2.3$, ranging from $-1.01$ for Leo IV, to $-1.87$ for Bo\"otes I. The significance of such deviations depends on the galaxy and is typically 95\% or more. When modeled as a log-normal, the IMF fit results in a larger peak mass than in the Milky Way disk, however a Milky Way disk value for the characteristic system mass ($\sim0.22$ M$_{\odot}$) is excluded only at 68\% significance, and only for some UFDs in the sample. We find that the IMF slope correlates well with the galaxy mean metallicity and, to a lesser degree, with the velocity dispersion and the total mass. The strength of the observed correlations is limited by shot noise in the number of observed stars, but future space-based missions like JWST and WFIRST will both enhance the number of dwarf Milky Way Satellites that can be studied in such detail, and the observation depth for individual galaxies.

We present the survey strategy and early results of the "Satellites Around Galactic Analogs" (SAGA) Survey. The SAGA Survey's goal is to measure the distribution of satellite galaxies around 100 systems analogous to the Milky Way down to the luminosity of the Leo I dwarf galaxy ($ M_r < -12.3 $). We define a Milky Way analog based on $K$-band luminosity and local environment. Here, we present satellite luminosity functions for 8 Milky Way analog galaxies between 20 to 40 Mpc. These systems have nearly complete spectroscopic coverage of candidate satellites within the projected host virial radius down to $ r_o < 20.75 $ using low redshift $gri$ color criteria. We have discovered a total of 25 new satellite galaxies: 14 new satellite galaxies meet our formal criteria around our complete host systems, plus 11 additional satellites in either incompletely surveyed hosts or below our formal magnitude limit. Combined with 13 previously known satellites, there are a total of 27 satellites around 8 complete Milky Way analog hosts. We find a wide distribution in the number of satellites per host, from 1 to 9, in the luminosity range for which there are five Milky Way satellites. Standard abundance matching extrapolated from higher luminosities predicts less scatter between hosts and a steeper luminosity function slope than observed. We find that the majority of satellites (26 of 27) are star-forming. These early results indicate that the Milky Way has a different satellite population than typical in our sample, potentially changing the physical interpretation of measurements based only on the Milky Way's satellite galaxies.

In low-mass galaxies, stellar feedback can drive gas outflows that generate non-equilibrium fluctuations in the gravitational potential. Using cosmological zoom-in baryonic simulations from the Feedback in Realistic Environments (FIRE) project, we investigate how these fluctuations affect stellar kinematics and the reliability of Jeans dynamical modeling in low-mass galaxies. We find that stellar velocity dispersion and anisotropy profiles fluctuate significantly over the course of galaxies' starburst cycles. We therefore predict an observable correlation between star formation rate and stellar kinematics: dwarf galaxies with higher recent star formation rates should have systemically higher stellar velocity dispersions. This prediction provides an observational test of the role of stellar feedback in regulating both stellar and dark-matter densities in dwarf galaxies. We find that Jeans modeling, which treats galaxies as virialized systems in dynamical equilibrium, overestimates a galaxy's dynamical mass during periods of post-starburst gas outflow and underestimates it during periods of net inflow. Short-timescale potential fluctuations lead to typical errors of $\sim 20\%$ in dynamical mass estimates, even if full 3-dimensional stellar kinematics -- including the orbital anisotropy -- are known exactly. When orbital anisotropy is $\textit{not}$ known a priori, typical mass errors arising from non-equilibrium fluctuations in the potential are larger than those arising from the mass-anisotropy degeneracy. However, Jeans modeling alone $\textit{cannot}$ reliably constrain the orbital anisotropy, and problematically, it often favors anisotropy models that do not reflect the true profile. If galaxies completely lose their gas and cease forming stars, fluctuations in the potential subside, and Jeans modeling becomes much more reliable.

The baryonic Tully-Fisher relation (BTFR) is a valuable observational tool and a critical test of galaxy formation theory. We explore the systematic uncertainty in the slope and the scatter of the observed line width BTFR utilizing homogeneously measured, unresolved HI observations for 930 isolated galaxies. We measure a fiducial relation of $\log_{10}{M_{\rm baryon}} = 3.24\log_{10}{V_{\rm rot}}+3.21$ with observed scatter of 0.25 dex with $10^{7.4} < M_{\rm baryon} < 10^{11.3}M_{\odot}$ where $V_{\rm rot}$ is measured from 20\% HI line widths. We vary the definitions of $M_{\rm baryon}$ and $V_{\rm rot}$, the sample selection, and the linear fitting algorithm. We obtain slopes ranging from 2.64 to 3.53 and scatter measurements ranging from 0.14 to 0.41~dex, indicating a systematic uncertainty of 0.25 in the BTFR slope derived from unresolved HI line widths. We compare our fiducial slope to literature measurements, where reported slopes range from 3.0 to 4.3 and scatter is either unmeasured, immeasurable or as large as 0.4~dex. Measurements derived from unresolved HI line widths tend to produce slopes of 3.3, while measurements derived from asymptotic rotation curves tend to produce slopes of 3.9. The largest factor affecting the BTFR slope is the definition of $V_{\rm rot}$. The sample definition, the mass range, and the linear fitting algorithm also significantly affect the measured BTFR. We find that galaxies in our sample with $V_{\rm rot} < 100~\rm{km~s^{-1}}$ are consistent with the line width BTFR of more massive galaxies, but these galaxies drive most of the observed scatter. It is critical when comparing predictions to an observed BTFR that the $V_{\rm rot}$ definition, the sample selection and the fitting algorithm are similarly defined. We recommend direct statistical comparisons between data sets with commensurable properties as opposed to comparing BTFR power-law fits.

Pegasus III (Peg III) is one of the few known ultra-faint stellar systems in the outer halo of the Milky Way. We present the results from a follow-up campaign with Magellan/IMACS and Keck/DEIMOS. Deep stellar photometry down to $r_0\approx 25$ mag at 50% completeness level has allowed accurate measurements of its photometric and structural properties. The color-magnitude diagram of Peg III confirms that the stellar system is well described by an old ($\gtrsim12$ Gyr) and metal-poor ([Fe/H]$\lesssim-2.0$ dex) stellar population at a heliocentric distance of $215\pm12$ kpc. The revised half-light radius $r_{h}=53\pm14$ pc, ellipticity $\epsilon=0.38^{+0.22}_{-0.38}$, and total luminosity $M_{V}=-3.4\pm0.4$ are in good agreement with the values quoted in our previous paper. We further report on the spectroscopic identification of seven, possibly eight member stars of Peg III. The Ca II triplet lines of the brightest member stars indicate that Peg III contains stars with metallicity as low as [Fe/H]=$-2.55\pm0.15$ dex. Peg III has a systemic velocity of $-222.9 \pm 2.6$ km s$^{-1}$ and a velocity dispersion of $5.4^{+3.0}_{-2.5}$ km s$^{-1}$. The inferred dynamical mass within the half-light radius is $1.4^{+3.0}_{-1.1} \times 10^6\rm{M_{\odot}}$ and the mass-to-light ratio $\rm{M/L}$$_{V} = 1470^{+5660}_{-1240}$ $\rm{M_{\odot}/L_{\odot}}$, providing further evidence that Peg III is a dwarf galaxy satellite. We find that Peg III and another distant dwarf satellite Pisces II lie relatively close to each other ($\Delta d_{spatial}=43\pm19$ kpc) and share similar radial velocities in the Galactic standard-of-rest frame ($\Delta v_{GSR}=12.3\pm3.7$ km s$^{-1}$). This suggests that they may share a common origin.

We use time-domain optical spectroscopy to distinguish between broad emission lines powered by accreting black holes (BHs) or stellar processes (i.e., supernovae) for 16 galaxies identified as AGN candidates by Reines \etal (2013). Our study is primarily focused on those objects with narrow emission-line ratios dominated by star formation. Based on follow-up spectra taken with the Magellan Echellette Spectrograph (MagE), the Dual Imaging Spectrograph, and the Ohio State Multi-Object Spectrograph, we find that the broad H$\alpha$ emission has faded or was ambiguous for all of the star-forming objects (14/16) over baselines ranging from 5 to 14 years. For the two objects in our follow-up sample with narrow-line AGN signatures (RGG 9 and RGG 119), we find persistent broad H$\alpha$ emission consistent with an AGN origin. Additionally, we use our MagE observations to measure stellar velocity dispersions for 15 objects in the Reines et al. (2013) sample, all with narrow-line ratios indicating the presence of an AGN. Stellar masses range from $\sim5\times10^{8}$ to $3\times10^{9}$~\msun, and we measure $\sigma_{\ast}$ ranging from $28-71~{\rm km~s^{-1}}$. These $\sigma_{\ast}$ correspond to some of the lowest-mass galaxies with optical signatures of AGN activity. We show that RGG 119, the one object which has both a measured $\sigma_{\ast}$ and persistent broad H$\alpha$ emission, falls near the extrapolation of the $\rm M_{BH}-\sigma_{\star}$ relation to the low-mass end.

We present the spatially resolved star formation history (SFH) of the Carina dwarf spheroidal galaxy, obtained from deep, wide-field g,r imaging and a metallicity distribution from the literature. Our photometry covers $\sim2$ deg$^2$, reaching up to $\sim10$ times the half-light radius of Carina with a completeness higher than $50\%$ at $g\sim24.5$, more than one magnitude fainter than the oldest turnoff. This is the first time a combination of depth and coverage of this quality has been used to derive the SFH of Carina, enabling us to trace its different populations with unprecedented accuracy. We find that Carina's SFH consists of two episodes well separated by a star formation temporal gap. These episodes occurred at old ($>10$ Gyr) and intermediate ($2$-$8$ Gyr) ages. Our measurements show that the old episode comprises the majority of the population, accounting for $54\pm5\%$ of the stellar mass within $1.3$ times the King tidal radius, while the total stellar mass derived for Carina is $1.60\pm0.09\times 10^{6} M_{\rm{\odot}}$, and the stellar mass-to-light ratio $1.8\pm0.2$. The SFH derived is consistent with no recent star formation which hints that the observed blue plume is due to blue stragglers. We conclude that the SFH of Carina evolved independently of the tidal field of the Milky Way, since the frequency and duration of its star formation events do not correlate with its orbital parameters. This result is supported by the age/metallicity relation observed in Carina, and the gradients calculated indicating that outer regions are older and more metal poor.

We present observations of the Pisces A and B galaxies with the Advanced Camera for Surveys on the Hubble Space Telescope. Photometry from these images clearly resolve a Red Giant Branch for both objects, demonstrating that they are nearby dwarf galaxies. We describe a Bayesian inferential approach to determining the distance to these galaxies using the magnitude of the tip of the RGB, and then apply this approach to these galaxies. We also provide the full probability distributions for parameters derived using this approach. This reveals the distance to these galaxies as $5.64^{+0.13}_{-0.15} \, {\rm Mpc}$ and $8.89^{+0.75}_{-0.85} \, {\rm Mpc}$ for Pisces A and B, respectively, placing both within the Local Volume but not the Local Group. We estimate the star formation histories of these galaxies, which suggests that they have recently undergone an increase in their star formation rates. Together these yield luminosities for Pisces A and B of $M_V=-11.57^{+0.06}_{-0.05}$ and $-12.9 \pm 0.2$, respectively, and estimated stellar masses of $\log(M_*/M_{\odot})= 7.0^{+0.4}_{-1.7}$ and $7.5^{+0.3}_{-1.8}$. We further show that these galaxies are likely at the boundary between nearby voids and higher-density filamentary structure. This suggests that they are entering a higher-density region from voids, where they would have experienced delayed evolution, consistent with their recent increased star formation rates. If this is indeed the case, they are useful for study as proxies of the galaxies that later evolved into typical Local Group satellite galaxies.

We present high-resolution Magellan/MIKE spectra of the four brightest confirmed red giant stars in the ultra-faint dwarf galaxy Bootes II (Boo II). These stars all inhabit the metal-poor tail of the Boo II metallicity distribution function. The chemical abundance pattern of all detectable elements in these stars is consistent with that of the Galactic halo. However, all four stars have undetectable amounts of neutron-capture elements Sr and Ba, with upper limits comparable to the lowest ever detected in the halo or in other dwarf galaxies. One star exhibits significant radial velocity variations over time, suggesting it to be in a binary system. Its variable velocity has likely increased past determinations of the Boo II velocity dispersion. Our four stars span a limited metallicity range, but their enhanced {\alpha}-abundances and low neutron-capture abundances are consistent with the interpretation that Boo II has been enriched by very few generations of stars. The chemical abundance pattern in Boo II confirms the emerging trend that the faintest dwarf galaxies have neutron-capture abundances distinct from the halo, suggesting the dominant source of neutron-capture elements in halo stars may be different than in ultra-faint dwarfs.

We examine the effects of stellar feedback and bursty star formation on low-mass galaxies ($M_{\rm star}=2\times10^6-5\times10^{10}{\rm M_{\odot}}$) using the FIRE (Feedback in Realistic Environments) simulations. While previous studies emphasized the impact of feedback on dark matter profiles, we investigate the impact on the stellar component: kinematics, radial migration, size evolution, and population gradients. Feedback-driven outflows/inflows drive significant radial stellar migration over both short and long timescales via two processes: (1) outflowing/infalling gas can remain star-forming, producing young stars that migrate $\sim1{\rm\,kpc}$ within their first $100 {\rm\,Myr}$, and (2) gas outflows/inflows drive strong fluctuations in the global potential, transferring energy to all stars. These processes produce several dramatic effects. First, galaxies' effective radii can fluctuate by factors of $>2$ over $\sim200 {\rm\,Myr}$, and these rapid size fluctuations can account for much of the observed scatter in radius at fixed $M_{\rm star}.$ Second, the cumulative effects of many outflow/infall episodes steadily heat stellar orbits, causing old stars to migrate outward most strongly. This age-dependent radial migration mixes---and even inverts---intrinsic age and metallicity gradients. Thus, the galactic-archaeology approach of calculating radial star-formation histories from stellar populations at $z=0$ can be severely biased. These effects are strongest at $M_{\rm star}\approx10^{7-9.6}{\rm M_{\odot}}$, the same regime where feedback most efficiently cores galaxies. Thus, detailed measurements of stellar kinematics in low-mass galaxies can strongly constrain feedback models and test baryonic solutions to small-scale problems in $\Lambda$CDM.

We present the first spectroscopic analysis of the faint and compact stellar system Draco II (Dra II, M_V=-2.9+/-0.8, r_h=19^{+8}_{-6} pc), recently discovered in the Pan-STARRS1 3{\pi} survey. The observations, conducted with DEIMOS on the Keck II telescope, establish some of its basic characteristics: the velocity data reveal a narrow peak with 9 member stars at a systemic heliocentric velocity <v_r>=-347.6^{+1.7}_{-1.8} km/s, thereby confirming Dra II is a satellite of the Milky Way; we infer a velocity dispersion with \sigma_{vr}=2.9+/-2.1 km/s (<8.4 km/s at the 95% confidence level), which implies log_{10}(M_{1/2})=5.5^{+0.4}_{-0.6} and log_{10}((M/L)_{1/2})=2.7^{+0.5}_{-0.8}, in Solar units; furthermore, very weak Calcium triplet lines in the spectra of the high signal-to-noise member stars imply [Fe/H]<-2.1, whilst variations in the line strengths of two stars with similar colours and magnitudes suggest a metallicity spread in Dra II. These new data cannot clearly discriminate whether Draco II is a star cluster or amongst the faintest, most compact, and closest dwarf galaxies. However, the sum of the three --- individually inconclusive --- pieces of evidence presented here, seems to favour the dwarf galaxy interpretation.

We present Magellan/M2FS, VLT/GIRAFFE, and Gemini South/GMOS spectroscopy of the newly discovered Milky Way satellite Reticulum II. Based on the spectra of 25 Ret II member stars selected from Dark Energy Survey imaging, we measure a mean heliocentric velocity of 62.8 +/- 0.5 km/s and a velocity dispersion of 3.3 +/- 0.7 km/s. The mass-to-light ratio of Ret II within its half-light radius is 470 +/- 210 Msun/Lsun, demonstrating that it is a strongly dark matter-dominated system. Despite its spatial proximity to the Magellanic Clouds, the radial velocity of Ret II differs from that of the LMC and SMC by 199 and 83 km/s, respectively, suggesting that it is not gravitationally bound to the Magellanic system. The likely member stars of Ret II span 1.3 dex in metallicity, with a dispersion of 0.28 +/- 0.09 dex, and we identify several extremely metal-poor stars with [Fe/H] < -3. In combination with its luminosity, size, and ellipticity, these results confirm that Ret II is an ultra-faint dwarf galaxy. With a mean metallicity of [Fe/H] = -2.65 +/- 0.07, Ret II matches Segue~1 as the most metal-poor galaxy known. Although Ret II is the third-closest dwarf galaxy to the Milky Way, the line-of-sight integral of the dark matter density squared is log J = 18.8 +/- 0.6 Gev^2/cm^5 within 0.2 degrees, indicating that the predicted gamma-ray flux from dark matter annihilation in Ret II is lower than that of several other dwarf galaxies.

We study the baryon content of low mass galaxies selected from the Sloan Digital Sky Survey (SDSS DR8), focusing on galaxies in isolated environments where the complicating physics of galaxy-galaxy interactions are minimized. We measure neutral hydrogen (HI) gas masses and line-widths for 148 isolated galaxies with stellar mass between $10^7$ and $10^{9.5} M_{\odot}$. We compare isolated low mass galaxies to more massive galaxies and galaxies in denser environments by remeasuring HI emission lines from the Arecibo Legacy Fast ALFA (ALFALFA) survey 40% data release. All isolated low mass galaxies either have large atomic gas fractions or large atomic gas fractions cannot be ruled out via their upper limits. We measure a median atomic gas fraction of $f_{\rm gas} = 0.82 \pm 0.13$ for our isolated low mass sample with no systems below 0.30. At all stellar masses, the correlations between galaxy radius, baryonic mass and velocity width are not significantly affected by environment. Finally, we estimate a median baryon to total dynamical mass fraction of $f_{\rm baryon, disk} = 0.15 \pm 0.18$. We also estimate two different median baryon to halo mass fractions using the results of semi-analytic models $(f_{\rm baryon, halo} = 0.04 \pm 0.06)$ and abundance matching $(f_{\rm baryon, halo} = 0.04 \pm 0.02)$. Baryon fractions estimated directly using HI observations appear independent of environment and maximum circular velocity, while baryon fractions estimated using abundance matching show a significant depletion of baryons at low maximum circular velocities.

We recently identified a population of low surface brightness objects in the field of the z=0.023 Coma cluster, using the Dragonfly Telephoto Array. Here we present Keck spectroscopy of one of the largest of these "ultra-diffuse galaxies" (UDGs), confirming that it is a member of the cluster. The galaxy has prominent absorption features, including the Ca II H+K lines and the G-band, and no detected emission lines. Its radial velocity of cz=6280 +- 120 km/s is within the 1 sigma velocity dispersion of the Coma cluster. The galaxy has an effective radius of 4.3 +- 0.3 kpc and a Sersic index of 0.89 +- 0.06, as measured from Keck imaging. We find no indications of tidal tails or other distortions, at least out to a radius of ~2 r_e. We show that UDGs are located in a previously sparsely populated region of the size - magnitude plane of quiescent stellar systems, as they are ~6 magnitudes fainter than normal early-type galaxies of the same size. It appears that the luminosity distribution of large quiescent galaxies is not continuous, although this could largely be due to selection effects. Dynamical measurements are needed to determine whether the dark matter halos of UDGs are similar to those of galaxies with the same luminosity or to those of galaxies with the same size.

The origin of the Galactic halo stellar structure known as the Monoceros ring is still under debate. In this work, we study that halo substructure using deep CFHT wide-field photometry obtained for the globular clusters NGC2419 and Koposov2, where the presence of Monoceros becomes significant because of their coincident projected position. Using Sloan Digital Sky Survey photometry and spectroscopy in the area surrounding these globulars and beyond, where the same Monoceros population is detected, we conclude that a second feature, not likely to be associated with Milky Way disk stars along the line-of-sight, is present as foreground population. Our analysis suggests that the Monoceros ring might be composed of an old stellar population of age t ~ 9Gyr and a new component ~ 4Gyr younger at the same heliocentric distance. Alternatively, this detection might be associated with a second wrap of Monoceros in that direction of the sky and also indicate a metallicity spread in the ring. The detection of such a low-density feature in other sections of this halo substructure will shed light on its nature.

We present the deepest optical photometry for any dwarf elliptical (dE) galaxy based on Hubble Space Telescope ACS observations of the Local Group dE galaxies NGC 147 and NGC 185. The resulting F606W and F814W color-magnitude diagrams are the first to reach below the main sequence turnoff in a dE galaxy, allowing us to determine full star formation histories in these systems. The ACS fields are located ~1.5 effective radii from the galaxy center to avoid photometric crowding. While our ACS pointings in both dEs show unambiguous evidence for old and intermediate age stars, the mean age in NGC 147 is ~ 4 Gyr younger as compared to NGC 185. In NGC 147, only 40% of stars were in place 12.5 Gyrs ago (z~5), with the bulk of the remaining stellar population forming between 5 to 7 Gyr. In contrast, 70% of stars were formed in NGC 185 field more than12.5 Gyr ago with the majority of the remaining population forming between 8 to 10 Gyr. Star formation ceased in both ACS fields at least 3 Gyr ago. Previous observations in the central regions of NGC 185 show evidence for star formation as recent as 100 Myr ago and a strong metallicity gradient with radius. We suggest that the orbit of NGC 185 has a larger pericenter as compared to NGC 147, allowing it to preserve radial gradients and maintain a small central reservoir of recycled gas. We interpret the inferred differences in star formation histories to imply an earlier infall time into the M31 environment for NGC 185 as compared to NGC 147.

High-precision astrometry throughout the Local Group is a unique capability of the Hubble Space Telescope (HST), with potential for transformative science, including constraining the nature of dark matter, probing the epoch of reionization, and understanding key physics of galaxy evolution. While Gaia will provide unparalleled astrometric precision for bright stars in the inner halo of the Milky Way, HST is the only current mission capable of measuring accurate proper motions for systems at greater distances (> 80 kpc), which represents the vast majority of galaxies in the Local Group. The next generation of proper-motion measurements will require long time baselines, spanning many years to decades and possibly multiple telescopes, combining HST with the James Webb Space Telescope (JWST) or the Wide-Field Infrared Survey Telescope (WFIRST). However, the current HST allocation process is not conducive to such multi-cycle/multi-mission science, which will bear fruit primarily over many years. We propose an HST astrometry initiative to enable long-time-baseline, multi-mission science, which we suggest could be used to provide comprehensive kinematic measurements of all dwarf galaxies and high surface-density stellar streams in the Local Group with HST's Advanced Camera for Surveys (ACS) or Wide Field Camera 3 (WFC3). Such an initiative not only would produce forefront scientific results within the next 5 years of HST's life, but also would serve as a critical anchor point for future missions to obtain unprecedented astrometric accuracy, ensuring that HST leaves a unique and lasting legacy for decades to come.

We report the discovery of 47 low surface brightness objects in deep images of a 3 x 3 degree field centered on the Coma cluster, obtained with the Dragonfly Telephoto Array. The objects have central surface brightness mu(g,0) ranging from 24 - 26 mag/arcsec^2 and effective radii r_e = 3"-10", as measured from archival Canada France Hawaii Telescope images. From their spatial distribution we infer that most or all of the objects are galaxies in the Coma cluster. This relatively large distance is surprising as it implies that the galaxies are very large: with r_e = 1.5 - 4.6 kpc their sizes are similar to those of L* galaxies even though their median stellar mass is only ~6 x 10^7 Solar masses. The galaxies are relatively red and round, with <g-i> = 0.8 and <b/a> = 0.74. One of the 47 galaxies is fortuitously covered by a deep Hubble Space Telescope ACS observation. The ACS imaging shows a large spheroidal object with a central surface brightness mu(g,0) = 25.8 mag/arcsec^2, a Sersic index n=0.6, and an effective radius of 7", corresponding to 3.4 kpc at the distance of Coma. The galaxy is not resolved into stars, consistent with expectations for a Coma cluster object. We speculate that UDGs may have lost their gas supply at early times, possibly resulting in very high dark matter fractions.

We report the discovery of two dwarf galaxies, Pisces A and B, from a blind 21 cm HI search. These were the only two galaxies found via optical imaging and spectroscopy of 22 HI clouds identified in the GALFA-HI survey as dwarf galaxy candidates. They have properties consistent with being in the Local Volume ($<10$ Mpc), and one has resolved stellar populations such that it may be on the outer edge of the Local Group ($\sim 1 \, {\rm Mpc}$ from M31). While the distance uncertainty makes interpretation ambiguous, these may be among the faintest starforming galaxies known. Additionally, rough estimates comparing these galaxies to $\Lambda$CDM dark matter simulations suggest consistency in number density, implying that dark matter halos likely to host these galaxies are primarily HI-rich. The galaxies may thus be indicative of a large population of dwarfs at the limit of detectability that are comparable to the faint satellites of the Local Group. Because they are outside the influence of a large dark matter halo to alter their evolution, these galaxies can provide critical anchors to dwarf galaxy formation models.

We present alpha element to iron abundance ratios, [$\alpha$/Fe], for four stars in the outer stellar halo of the Andromeda Galaxy (M 31). The stars were identified as high-likelihood field halo stars by Gilbert et al. (2012) and lie at projected distances between 70 and 140 kpc from M 31's center. These are the first alpha abundances measured for a halo star in a galaxy beyond the Milky Way. The stars range in metallicity between [Fe/H]= -2.2 and [Fe/H]= -1.4. The sample's average [$\alpha$/Fe] ratio is +0.20+/-0.20. The best-fit average value is elevated above solar which is consistent with rapid chemical enrichment from Type II supernovae. The mean [$\alpha$/Fe] ratio of our M31 outer halo sample agrees (within the uncertainties) with that of Milky Way inner/outer halo stars that have a comparable range of [Fe/H].

We present new constraints on the star formation histories of six ultra-faint dwarf galaxies: Bootes I, Canes Venatici II, Coma Berenices, Hercules, Leo IV, and Ursa Major I. Our analysis employs a combination of high-precision photometry obtained with the Advanced Camera for Surveys on the Hubble Space Telescope, medium-resolution spectroscopy obtained with the DEep Imaging Multi-Object Spectrograph on the W.M. Keck Observatory, and updated Victoria-Regina isochrones tailored to the abundance patterns appropriate for these galaxies. The data for five of these Milky Way satellites are best fit by a star formation history where at least 75% of the stars formed by z~10 (13.3 Gyr ago). All of the galaxies are consistent with 80% of the stars forming by z~6 (12.8 Gyr ago) and 100% of the stars forming by z~3 (11.6 Gyr ago). The similarly ancient populations of these galaxies support the hypothesis that star formation in the smallest dark matter sub-halos was suppressed by a global outside influence, such as the reionization of the universe.

We present the metallicity distribution of red giant branch (RGB) stars in M31's stellar halo, derived from photometric metallicity estimates for over 1500 spectroscopically confirmed RGB halo stars. The stellar sample comes from 38 halo fields observed with the Keck/DEIMOS spectrograph, ranging from 9 to 175 kpc in projected distance from M31's center, and includes 52 confirmed M31 halo stars beyond 100 kpc. While a wide range of metallicities is seen throughout the halo, the metal-rich peak of the metallicity distribution function becomes significantly less prominent with increasing radius. The metallicity profile of M31's stellar halo shows a continuous gradient from 9 to ~100 kpc, with a magnitude of -0.01 dex/kpc. The stellar velocity distributions in each field are used to identify stars that are likely associated with tidal debris features. The removal of tidal debris features does not significantly alter the metallicity gradient in M31's halo: a gradient is maintained in fields spanning 10 to 90 kpc. We analyze the halo metallicity profile, as well as the relative metallicities of stars associated with tidal debris features and the underlying halo population, in the context of current simulations of stellar halo formation. We argue that the large scale gradient in M31's halo implies M31 accreted at least one relatively massive progenitor in the past, while the field to field variation seen in the metallicity profile indicates that multiple smaller progenitors are likely to have contributed substantially to M31's outer halo.

We present a new method for determining the Galactic gravitational potential based on forward modeling of tidal stellar streams. We use this method to test the performance of smooth and static analytic potentials in representing realistic dark matter halos, which have substructure and are continually evolving by accretion. Our FAST-FORWARD method uses a Markov Chain Monte Carlo algorithm to compare, in 6D phase space, an "observed" stream to models created in trial analytic potentials. We analyze a large sample of streams evolved in the Via Lactea II (VL2) simulation, which represents a realistic Galactic halo potential. The recovered potential parameters are in agreement with the best fit to the global, present-day VL2 potential. However, merely assuming an analytic potential limits the dark matter halo mass measurement to an accuracy of 5 to 20%, depending on the choice of analytic parametrization. Collectively, mass estimates using streams from our sample reach this fundamental limit, but individually they can be highly biased. Individual streams can both under- and overestimate the mass, and the bias is progressively worse for those with smaller perigalacticons, motivating the search for tidal streams at galactocentric distances larger than 70 kpc. We estimate that the assumption of a static and smooth dark matter potential in modeling of the GD-1 and Pal5-like streams introduces an error of up to 50% in the Milky Way mass estimates.

We present alpha to iron abundance ratios for 226 individual red giant branch stars in nine dwarf galaxies of the Andromeda (M 31) satellite system. The abundances are measured from the combined signal of Mg, Si, Ca, and Ti lines in Keck/DEIMOS medium-resolution spectra. This constitutes the first large sample of alpha abundance ratios measured in the M 31 satellite system. The dwarf galaxies in our sample exhibit a variety of alpha abundance ratios, with the average values in each galaxy ranging from approximately solar ([alpha/Fe] ~ +0.0) to alpha-enhanced ([alpha/Fe] ~ +0.5). These variations do not show a correlation with internal kinematics, environment, or stellar density. We confirm radial gradients in the iron abundance of two galaxies out of the five with sufficient data (NGC 185 and And II). There is only tentative evidence for an alpha abundance radial gradient in NGC 185. We homogeneously compare our results to the Milky Way classical dwarf spheroidals, finding evidence for wider variation in average alpha abundance. In the absence of chemical abundances for the M 31 stellar halo, we compare to the Milky Way stellar halo. A stellar halo comprised of disrupted M31 satellites is too metal-rich and inconsistent with the Milky Way halo alpha abundance distribution even if considering only satellites with predominantly old stellar populations. The M 31 satellite population provides a second system in which to study chemical abundances of dwarf galaxies and reveals a wider variety of abundance patterns than the Milky Way.

We present global metallicity properties, metallicity distribution functions (MDFs) and radial metallicity profiles for the six most luminous M31 dwarf galaxy satellites: M32, NGC 205, NGC 185, NGC 147, Andromeda VII, and Andromeda II. The results presented are the first spectroscopic MDFs for dwarf systems surrounding a host galaxy other than the Milky Way. Our sample consists of individual metallicity measurements for 1243 red giant branch (RGB) member stars spread across these six systems. We determine metallicities based on the strength of the Ca II triplet lines using the empirical calibration of Carrera et al.(2013) which is calibrated over the metallicity range -4 < [Fe/H] <+0.5. We find that these M31 satellites lie on the same luminosity-metallicity relationship as the Milky Way dwarf satellites. We do not find a trend between the internal metallicity spread and galaxy luminosity, contrary to previous studies. The MDF widths of And II and And VII are similar to the MW dwarfs of comparable luminosity, however, our four brightest M31 dwarf are more luminous than any of the MW dwarf spheroidals and have broader MDFs. The MDFs of our six M31 dwarfs are consistent with the leaky box model of chemical evolution, although our metallicity errors allow a wide range of evolution models. We find a significant radial gradient in metallicity in only two of our six systems, NGC 185 and Andromeda II, and flat radial metallicity gradients in the rest of our sample with no observed correlation between rotational support and radial metallicity gradients. While the average properties and radial trends of the M31 dwarf galaxies agree with MW counterparts at similar luminosity, the detailed MDFs are different, particularly at the metal-rich end.

We present optical imaging and spectroscopy and HI imaging of the Virgo Cluster galaxy IC 3418, which is likely a "smoking gun" example of the transformation of a dwarf irregular into a dwarf elliptical galaxy by ram pressure stripping. IC 3418 has a spectacular 17 kpc length UV-bright tail comprised of knots, head-tail, and linear stellar features. The only H-alpha emission arises from a few HII regions in the tail, the brightest of which are at the heads of head-tail UV sources whose tails point toward the galaxy ("fireballs"). Several of the elongated tail sources have H-alpha peaks outwardly offset by 80-150 pc from the UV peaks, suggesting that gas clumps continue to accelerate through ram pressure, leaving behind streams of newly formed stars which have decoupled from the gas. Absorption line strengths, measured from Keck DEIMOS spectra, together with UV colors, show star formation stopped 300+/-100 Myr ago in the main body, and a strong starburst occurred prior to quenching. While neither H-alpha nor HI emission are detected in the main body of the galaxy, we have detected 4x10^7 M_sun of HI from the tail with the VLA. The gas consumption timescale in the tail is relatively long, implying that most of the stripped gas does not form stars but joins the ICM. The velocities of tail HII regions, measured from Keck LRIS spectra, extend only a small fraction of the way to the cluster velocity, suggesting that star formation does not happen in more distant parts of the tail. Stars in the outer tail have velocities exceeding the escape speed, but some in the inner tail should fall back into the galaxy, forming halo streams. One likely fallback stream is identified.

We present early results from a Hubble Space Telescope survey of the ultra-faint dwarf galaxies. These Milky Way satellites were discovered in the Sloan Digital Sky Survey, and appear to be an extension of the classical dwarf spheroidals to low luminosities, offering a new front in the efforts to understand the missing satellite problem. Because they are the least luminous, most dark matter dominated, and least chemically evolved galaxies known, the ultra-faint dwarfs are the best candidate fossils from the early universe. The primary goal of the survey is to measure the star-formation histories of these galaxies and discern any synchronization due to the reionization of the universe. We find that the six galaxies of our survey have very similar star-formation histories, and that each is dominated by stars older than 12 Gyr.

We present a sample of 151 dwarf galaxies (10^8.5 < M_stellar < 10^9.5 Msun) that exhibit optical spectroscopic signatures of accreting massive black holes (BHs), increasing the number of known active galaxies in this stellar mass range by more than an order of magnitude. Utilizing data from the Sloan Digital Sky Survey Data Release 8 and stellar masses from the NASA-Sloan Atlas, we have systematically searched for active BHs in ~25,000 emission-line galaxies with stellar masses comparable to the Magellanic Clouds and redshifts z<0.055. Using the narrow-line [OIII]/H-beta versus [NII]/H-alpha diagnostic diagram, we find photoionization signatures of BH accretion in 136 galaxies, a small fraction of which also exhibit broad H-alpha emission. For these broad-line AGN candidates, we estimate BH masses using standard virial techniques and find a range of 10^5 < M_BH < 10^6 Msun and a median of M_BH ~ 2 x 10^5 Msun. We also detect broad H-alpha in 15 galaxies that have narrow-line ratios consistent with star-forming galaxies. Follow-up observations are required to determine if these are true type 1 AGN or if the broad H-alpha is from stellar processes. The median absolute magnitude of the host galaxies in our active sample is Mg = -18.1 mag, which is ~1-2 magnitudes fainter than previous samples of AGN hosts with low-mass BHs. This work constrains the smallest galaxies that can form a massive BH, with implications for BH feedback in low-mass galaxies and the origin of the first supermassive BH seeds.

We present Keck/DEIMOS spectroscopy and CFHT/MegaCam photometry for the Milky Way globular cluster Palomar 13. We triple the number of spectroscopically confirmed members, including many repeat velocity measurements. Palomar 13 is the only known globular cluster with possible evidence for dark matter, based on a Keck/HIRES 21 star velocity dispersion of sigma=2.2+/-0.4 km/s. We reproduce this measurement, but demonstrate that it is inflated by unresolved binary stars. For our sample of 61 stars, the velocity dispersion is sigma=0.7(+0.6/-0.5) km/s. Combining our DEIMOS data with literature values, our final velocity dispersion is sigma=0.4(+0.4/-0.3) km/s. We determine a spectroscopic metallicity of [Fe/H]=-1.6+/-0.1 dex, placing a 1-sigma upper limit of sigma_[Fe/H]~0.2 dex on any internal metallicity spread. We determine Palomar 13's total luminosity to be M_V=-2.8+/-0.4, making it among the least luminous known globular clusters. The photometric isophotes are regular out to the half-light radius and mildly irregular outside this radius. The outer surface brightness profile slope is shallower than typical globular clusters (eta=-2.8+/-0.3). Thus at large radius, tidal debris is likely affecting the appearance of Palomar 13. Combining our luminosity with the intrinsic velocity dispersion, we find a dynamical mass of of M_1/2=1.3(+2.7/-1.3)x10^3 M_sun and a mass-to-light ratio of M/L_V=2.4(+5.0/-2.4) M_sun/L_sun. Within our measurement errors, the mass-to-light ratio agrees with the theoretical predictions for a single stellar population. We conclude that, while there is some evidence for tidal stripping at large radius, the dynamical mass of Palomar 13 is consistent with its stellar mass and neither significant dark matter, nor extreme tidal heating, is required to explain the cluster dynamics.

We present a homogeneous study of blue straggler stars across ten outer halo globular clusters, three classical dwarf spheroidal and nine ultra-faint galaxies based on deep and wide-field photometric data taken with MegaCam on the Canada-France-Hawaii Telescope. We find blue straggler stars to be ubiquitous among these Milky Way satellites. Based on these data, we can test the importance of primordial binaries or multiple systems on blue straggler star formation in low density environments. For the outer halo globular clusters we find an anti-correlation between the specific frequency of blue straggler and absolute magnitude, similar to that previously observed for inner halo clusters. When plotted against density and encounter rate, the frequency of blue stragglers are well fitted by single trends with smooth transitions between dwarf galaxies and globular clusters, which points to a common origin for their blue stragglers. The fraction of blue stragglers stays constant and high in the low encounter rate regime spanned by our dwarf galaxies, and decreases with density and encounter rate in the range spanned by our globular clusters. We find that young stars can mimic blue stragglers in dwarf galaxies only if their ages are 2.5+/-0.5 Gyr and they represent ~1-7% of the total number of stars, which we deem highly unlikely. These results point to mass-transfer or mergers of primordial binaries or multiple systems as the dominant blue straggler formation mechanism in low density systems.

Segue 2, discovered by Belokurov et al. (2009), is a galaxy with a luminosity of only 900 L_sun. We present Keck/DEIMOS spectroscopy of 25 members of Segue 2--a threefold increase in spectroscopic sample size. The velocity dispersion is too small to be measured with our data. The upper limit with 90% (95%) confidence is sigma_v < 2.2 (2.6) km/s, the most stringent limit for any galaxy. The corresponding limit on the mass within the 3-D half-light radius (46 pc) is M_1/2 < 1.5 (2.1) x 10^5 M_sun. Segue 2 is the least massive galaxy known. We identify Segue 2 as a galaxy rather than a star cluster based the wide dispersion in [Fe/H] (from -2.85 to -1.33) among the member stars. The stars' [alpha/Fe] ratios decline with increasing [Fe/H], indicating that Segue 2 retained Type Ia supernova ejecta despite its presently small mass and that star formation lasted for at least 100 Myr. The mean metallicity, <[Fe/H]> = -2.22 +/- 0.13 (about the same as the Ursa Minor galaxy, 330 times more luminous than Segue 2), is higher than expected from the luminosity-metallicity relation defined by more luminous dwarf galaxy satellites of the Milky Way. Segue 2 may be the barest remnant of a tidally stripped, Ursa Minor-sized galaxy. If so, it is the best example of an ultra-faint dwarf galaxy that came to be ultra-faint through tidal stripping. Alternatively, Segue 2 could have been born in a very low-mass dark matter subhalo (v_max < 10 km/s), below the atomic hydrogen cooling limit.

We present constraints on the stellar initial mass function (IMF) in two ultra-faint dwarf (UFD) galaxies, Hercules and Leo IV, based on deep HST/ACS imaging. The Hercules and Leo IV galaxies are extremely low luminosity (M_V = -6.2, -5.5), metal-poor (<[Fe/H]>= -2.4, -2.5) systems that have old stellar populations (> 11 Gyr). Because they have long relaxation times, we can directly measure the low-mass stellar IMF by counting stars below the main-sequence turnoff without correcting for dynamical evolution. Over the stellar mass range probed by our data, 0.52 - 0.77 Msun, the IMF is best fit by a power-law slope of alpha = 1.2^{+0.4}_{-0.5} for Hercules and alpha = 1.3 +/- 0.8 for Leo IV. For Hercules, the IMF slope is more shallow than a Salpeter IMF (alpha=2.35) at the 5.8-sigma level, and a Kroupa IMF (alpha=2.3 above 0.5 Msun) at 5.4-sigma level. We simultaneously fit for the binary fraction, finding f_binary = 0.47^{+0.16}_{-0.14} for Hercules, and 0.47^{+0.37}_{-0.17} for Leo IV. The UFD binary fractions are consistent with that inferred for Milky Way stars in the same mass range, despite very different metallicities. In contrast, the IMF slopes in the UFDs are shallower than other galactic environments. In the mass range 0.5 - 0.8 Msun, we see a trend across the handful of galaxies with directly measured IMFs such that the power-law slopes become shallower (more bottom-light) with decreasing galactic velocity dispersion and metallicity. This trend is qualitatively consistent with results in elliptical galaxies inferred via indirect methods and is direct evidence for IMF variations with galactic environment.

We present deep CMDs for two Subaru Suprime-Cam fields in the Virgo Stellar Stream(VSS)/Virgo Overdensity(VOD) and compare them to a field centred on the highest concentration of Sagittarius (Sgr) Tidal Stream stars in the leading arm, Branch A of the bifurcation. A prominent population of MS stars is detected in all three fields and can be traced as faint as g~24 mag. Using theoretical isochrone fitting we derive an age of 9.1(+1.0;-1.1)Gyr, a median abundance of [Fe/H]=-0.70 (+0.15; -0.20)dex and d_helio of 30.9+-3.0kpc for the MS of the Sgr Stream Branch A. The dominant main sequence populations in the two VSS/VOD fields (Lsun ~265 deg, Bsun ~13 deg) are located at a mean distance of 23.3+-1.6kpc and have an age ~8.2Gyr and an abundance [Fe/H]=-0.67(+0.16;-0.12)dex similar to the Sgr Stream stars. These parameters are also in good agreement with the age of the main population in the Sagittarius dwarf (8.0+-1.5Gyr). They also agree with the peak in the metallicity distribution of 2-3Gyr old M-giants, [Fe/H] ~ -0.6dex, in the Sgr leading arm north. We then compare the results from the VSS/VOD fields with the Sgr Tidal Stream model by Law & Majewski based on a triaxial Galactic halo that is calibrated with SDSS Sgr A-branch and 2MASS M-giant stars. We find that the most prominent feature in the CMDs, the MS population at 23kpc, is not explained by the model. Instead the model predicts in these directions a low density filamentary structure of Sgr debris stars at ~9kpc and a slightly higher concentration of Sgr stars spread from 42-53kpc. At best there is only marginal evidence for the presence of these populations in our data. Our findings then suggest that while there are probably some Sgr debris stars present, the dominant stellar population in the VOD originates from a different halo structure that has almost identical age and metallicity as some sections of the Sgr tidal stream.

We present Keck/DEIMOS spectroscopy of resolved stars in the M31 satellites And XXVIII & And XXIX. We show that these are likely self-bound galaxies based on 18 and 24 members in And XXVIII & And XXIX, respectively. And XXVIII has a systemic velocity of -331.1 +/- 1.8 km/s and velocity dispersion of 4.9+/-1.6 km/s, implying a mass-to-light ratio (within r_1/2) of ~ 44 +/- 41. And XXIX has a systemic velocity of -194.4 +/- 1.5 km/s and velocity dispersion of 5.7+/-1.2 km/s, implying a mass-to-light ratio (within r_1/2) of ~ 124 +/- 72. The internal kinematics and implied masses of And XXVIII & And XXIX are similar to dwarf spheroidals (dSphs) of comparable luminosities, implying that these objects are dark matter-dominated dwarf galaxies. Despite the large projected distances from their host (380 and 188 kpc), the kinematics of these dSph suggest that they are bound M31 satellites.

The Milky Way ultra-faint dwarf galaxies (UFDs) contain some of the oldest, most metal-poor stars in the Universe. We present [Mg/Fe], [Si/Fe], [Ca/Fe], [Ti/Fe], and mean [alpha/Fe], abundance ratios for 61 individual red giant branch stars across 8 UFDs. This is the largest sample of alpha abundances published to date in galaxies with absolute magnitudes M_V > -8, including the first measurements for Segue 1, Canes Venatici II, Ursa Major I, and Leo T. Abundances were determined via medium-resolution Keck/DEIMOS spectroscopy and spectral synthesis. The sample spans the metallicity range -3.4 < [Fe/H] < -1.1. With the possible exception of Segue 1 and Ursa Major II, the individual UFDs show on average lower [alpha/Fe] at higher metallicities, consistent with enrichment from Type Ia supernovae. Thus even the faintest galaxies have undergone at least a limited level of chemical self-enrichment. Together with recent photometric studies, this suggests that star formation in the UFDs was not a single burst, but instead lasted at least as much as the minimum time delay of the onset of Type Ia supernovae (~100 Myr) and less than ~2 Gyr. We further show that the combined population of UFDs has an [alpha/Fe] abundance pattern that is inconsistent with a flat, Galactic halo-like alpha abundance trend, and is also qualitatively different from that of the more luminous CVn I dSph, which does show a hint of a plateau at very low [Fe/H].

We present the surface brightness profile of M31's stellar halo out to a projected radius of 175 kpc. The surface brightness estimates are based on confirmed samples of M31 red giant branch stars derived from Keck/DEIMOS spectroscopic observations. A set of empirical spectroscopic and photometric M31 membership diagnostics is used to identify and reject foreground and background contaminants. This enables us to trace the stellar halo of M31 to larger projected distances and fainter surface brightnesses than previous photometric studies. The surface brightness profile of M31's halo follows a power-law with index -2.2 +/- 0.2 and extends to a projected distance of at least ~175 kpc (~ 2/3 of M31's virial radius), with no evidence of a downward break at large radii. The best-fit elliptical isophotes have b/a=0.94 with the major axis of the halo aligned along the minor axis of M31's disk, consistent with a prolate halo, although the data are also consistent with M31's halo having spherical symmetry. The fact that tidal debris features are kinematically cold is used to identify substructure in the spectroscopic fields out to projected radii of 90 kpc, and investigate the effect of this substructure on the surface brightness profile. The scatter in the surface brightness profile is reduced when kinematically identified tidal debris features in M31 are statistically subtracted; the remaining profile indicates a comparatively diffuse stellar component to M31's stellar halo exists to large distances. Beyond 90 kpc, kinematically cold tidal debris features can not be identified due to small number statistics; nevertheless, the significant field-to-field variation in surface brightness beyond 90 kpc suggests that the outermost region of M31's halo is also comprised to a significant degree of stars stripped from accreted objects.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

We present evidence for a new Milky Way stellar tidal stream in the direction of the Andromeda and Triangulum (M31 and M33) galaxies. Using a matched-filter technique, we search the Sloan Digital Sky Survey Data Release 8 by creating stellar density maps which probe the Milky Way halo at distances between 8 and 40 kpc. A visual search of these maps recovers all of the major known stellar streams, as well as a new stream in the direction of M31/M33 which we name the Triangulum stream. The stream spans 0.2 deg by 12 deg on the sky, or 75 pc by 5.5 kpc in physical units with a best fitting distance of 26+/-4 kpc. The width of the stream is consistent with being the tidal remnant of a globular cluster. A color magnitude diagram of the stream region shows an overdensity which, if identified as a main sequence turn-off, corresponds to an old (~12 Gyr) and metal-poor ([Fe/H] \sim -1.0 dex) stellar population. Future kinematic studies of this and similar cold streams will provide tight constraints on the shape of the Galactic gravitational potential.

We present kinematical profiles and metallicity for the M31 dwarf spheroidal (dSph) satellite galaxy Andromeda II (And II) based on Keck DEIMOS spectroscopy of 531 red giant branch stars. Our kinematical sample is among the largest for any M31 satellite and extends out to two effective radii (r_eff = 5.3' = 1.1 kpc). We find a mean systemic velocity of -192.4+-0.5 km/s and an average velocity dispersion of sigma_v = 7.8+-1.1 km/s. While the rotation velocity along the major axis of And II is nearly zero (<1 km/s), the rotation along the minor axis is significant with a maximum rotational velocity of v_max=8.6+-1.8 km/s. We find a kinematical major axis, with a maximum rotational velocity of v_max=10.9+-2.4 km/s, misaligned by 67 degrees to the isophotal major axis. And II is thus the first dwarf galaxy with evidence for nearly prolate rotation with a v_max/sigma_v = 1.1, although given its ellipticity of epsilon = 0.10, this object may be triaxial. We measured metallicities for a subsample of our data, finding a mean metallicity of [Fe/H] = -1.39+- 0.03 dex and an internal metallicity dispersion of 0.72+-0.03 dex. We find a radial metallicity gradient with metal-rich stars more centrally concentrated, but do not observe a significant difference in the dynamics of two metallicity populations. And II is the only known dwarf galaxy to show minor axis rotation making it a unique system whose existence offers important clues on the processes responsible for the formation of dSphs.

We demonstrate that dwarf galaxies (10^7 < M_stellar < 10^9 Msun) with no active star formation are extremely rare (<0.06%) in the field. Our sample is based on the NASA-Sloan Atlas which is a re-analysis of the Sloan Digital Sky Survey Data Release 8. We examine the relative number of quenched versus star forming dwarf galaxies, defining quenched galaxies as having no Halpha emission (EW_Halpha < 2 AA) and a strong 4000AA-break. The fraction of quenched dwarf galaxies decreases rapidly with increasing distance from a massive host, leveling off for distances beyond 1.5 Mpc. We define galaxies beyond 1.5 Mpc of a massive host galaxy to be in the field. We demonstrate that there is a stellar mass threshold of M_stellar < 1.0x10^9 Msun below which quenched galaxies do not exist in the field. Below this threshold, we find that none of the 2951 field dwarf galaxies are quenched; all field dwarf galaxies show evidence for recent star formation. Correcting for volume effects, this corresponds to a 1-sigma upper limit on the quenched fraction of 0.06%. In more dense environments, quenched galaxies account for 23% of the dwarf population over the same stellar mass range. The majority of quenched dwarf galaxies (often classified as dwarf elliptical galaxies) are within 2 virial radii of a massive galaxy, and only a few percent of quenched dwarf galaxies exist beyond 4 virial radii. Thus, for galaxies with stellar mass less than 1.0x10^9 Msun, ending star-formation requires the presence of a more massive neighbor, providing a stringent constraint on models of star formation feedback.

We present new constraints on the star formation histories of the ultra-faint dwarf (UFD) galaxies, using deep photometry obtained with the Hubble Space Telescope (HST). A galaxy class recently discovered in the Sloan Digital Sky Survey, the UFDs appear to be an extension of the classical dwarf spheroidals to low luminosities, offering a new front in efforts to understand the missing satellite problem. They are the least luminous, most dark-matter dominated, and least chemically-evolved galaxies known. Our HST survey of six UFDs seeks to determine if these galaxies are true fossils from the early universe. We present here the preliminary analysis of three UFD galaxies: Hercules, Leo IV, and Ursa Major I. Classical dwarf spheroidals of the Local Group exhibit extended star formation histories, but these three Milky Way satellites are at least as old as the ancient globular cluster M92, with no evidence for intermediate-age populations. Their ages also appear to be synchronized to within ~1 Gyr of each other, as might be expected if their star formation was truncated by a global event, such as reionization.

We present a resolved-star spectroscopic survey of 15 dwarf spheroidal (dSph) satellites of the Andromeda Galaxy (M31) as part of the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) project. We filter foreground contamination from Milky Way (MW) stars, noting that MW substructure is evident in this contaminant sample. We also filter M31 halo field giant stars, and identify the remainder as probable dSph members. We then use these members to determine the kinematical properties of the dSphs. For the first time, we confirm that And XVIII, XXI, and XXII show kinematics consistent with bound, dark matter-dominated galaxies. From the velocity dispersions for the full sample of dSphs we determine masses, which we combine with the size and luminosity of the galaxies to produce mass-size-luminosity scaling relations. With these scalings we determine that the M31 dSphs are fully consistent with the MW dSphs, suggesting that the well-studied MW satellite population provides a fair sample for broader conclusions. We also estimate dark matter halo masses of the satellites, and find that there is no sign that the luminosity of these galaxies depends on their dark halo mass, a result consistent with what is seen for MW dwarfs. Two of the M31 dSphs (And XV, XVI) have estimated maximum circular velocities smaller than 12 km/s (to 1sigma), which likely places them within the lowest mass dark matter halos known to host stars (along with Bootes I of the MW). Finally, we use the systemic velocities of the M31 satellites to estimate the total mass of the M31 halo: within 139 kpc, M31's mass is 8 +4.1 - 3.7 x 10^11 corresponding to a virial mass for M31's dark matter halo of 1.2 +0.9 -0.7 x 10^12, consistent with previous results.

We report the discovery of a new ultra-faint globular cluster in the constellation of Ursa Minor, based on stellar photometry from the MegaCam imager at the Canada-France-Hawaii Telescope (CFHT). We find that this cluster, Munoz 1, is located at a distance of 45 +/- 5 kpc and at a projected distance of only 45 arcmin from the center of the Ursa Minor dSph galaxy. Using a Maximum Likelihood technique we measure a half-light radius of 0.5 arcmin, or equivalently 7 pc and an ellipticity consistent with being zero. We estimate its absolute magnitude to be M_V=-0.4 +/- 0.9, which corresponds to L_V=120 (+160, -65) L_sun and we measure a heliocentric radial velocity of -137 +/- 4 km/s based on Keck/DEIMOS spectroscopy. This new satellite is separate from Ursa Minor by ~30 kpc and 110 km/s suggesting the cluster is not obviously associated with the dSph, despite the very close angular separation. Based on its photometric properties and structural parameters we conclude that Munoz 1 is a new ultra-faint stellar cluster. Along with Segue 3 this is one of the faintest stellar clusters known to date.

The combination of large size, high stellar density, high metallicity, and Sersic surface brightness profile of the spheroidal component of the Andromeda galaxy (M31) within R_proj ~ 20 kpc suggest that it is unlike any subcomponent of the Milky Way. In this work we capitalize on our proximity to and external view of M31 to probe the kinematical properties of this "inner spheroid." We employ a Markov chain Monte Carlo (MCMC) analysis of resolved stellar kinematics from Keck/DEIMOS spectra of 5651 red giant branch stars to disentangle M31's inner spheroid from its stellar disk. We measure the mean velocity and dispersion of the spheroid in each of five spatial bins after accounting for a locally cold stellar disk as well as the Giant Southern Stream and associated tidal debris. For the first time, we detect significant spheroid rotation (v_rot ~ 50 km/s) beyond R_proj ~ 5 kpc. The velocity dispersion decreases from about 140 km/s at R_proj = 7 kpc to 120 km/s at R_proj = 14 kpc, consistent to 2 sigma with existing measurements and models. We calculate the probability that a given star is a member of the spheroid and find that the spheroid has a significant presence throughout the spatial extent of our sample. Lastly, we show that the flattening of the spheroid is due to velocity anisotropy in addition to rotation. Though this suggests that the inner spheroid of M31 more closely resembles an elliptical galaxy than a typical spiral galaxy bulge, it should be cautioned that our measurements are much farther out (2 - 14 r_eff) than for the comparison samples.

As part of the SPLASH survey of the Andromeda galaxy (M31) and its neighbors, we have obtained Keck/DEIMOS spectra of the compact elliptical (cE) satellite M32. This is the first resolved-star kinematical study of any cE galaxy. In contrast to previous studies that extended out to r<30"~1Re~100pc, we measure the rotation curve and velocity dispersion profile out to r~250" and higher order Gauss-Hermite moments out to r~70". We achieve this by combining integrated-light spectroscopy at small radii (where crowding/blending are severe) with resolved stellar spectroscopy at larger radii, using spatial and kinematical information to statistically account for M31 contamination. The rotation curve and velocity dispersion profile extend well beyond the radius (r~150") where the isophotes are distorted. Unlike NGC 205, another close dwarf companion of M31, M32's kinematic are regular and symmetric and do not show obvious sharp gradients across the region of isophotal elongation and twists. We interpret M32's kinematics using three-integral axisymmetric dynamical equilibrium models constructed using Schwarzschild's orbit superposition technique. Models with a constant M/L can fit the data remarkably well. However, since such a model requires an increasing tangential anisotropy with radius, invoking the presence of an extended dark halo may be more plausible. Such an extended dark halo is definitely required to bind a half-dozen fast-moving stars observed at the largest radii, but these stars may not be an equilibrium component of M32. The observed regularity of the stellar kinematics, as well as the possible detection of an extended dark halo, are unexpected if M31 tides are significant at large radii. While these findings by themselves do not rule out tidal models for cE formation, they suggest that tidal stripping may not be as significant for shaping cE galaxies as has often been argued.

The discovery of Ultra-Faint Dwarf (UFD) galaxies in the halo of the Milky Way extends the faint end of the galaxy luminosity function to a few hundred solar luminosities. This extremely low luminosity regime poses a significant challenge for the photometric characterization of these systems. We present a suite of simulations aimed at understanding how different observational choices related to the properties of a low luminosity system impact our ability to determine its true structural parameters such as half-light radius and central surface brightness. We focus on estimating half-light radii (on which mass estimates depend linearly) and find that these numbers can have up to 100% uncertainties when relatively shallow photometric surveys, such as SDSS, are used. Our simulations suggest that to recover structural parameters within 10% or better of their true values: (a) the ratio of the field-of-view to the half-light radius of the satellite must be greater than three, (b) the total number of stars, including background objects should be larger than 1000, and (c) the central to background stellar density ratio must be higher than 20. If one or more of these criteria are not met, the accuracy of the resulting structural parameters can be significantly compromised. In the context of future surveys such as LSST, the latter condition will be closely tied to our ability to remove unresolved background galaxies. Assessing the reliability of measured structural parameters will become increasingly critical as the next generation of deep wide-field surveys detects UFDs beyond the reach of current spectroscopic limits.

We introduce a comprehensive analysis of multi-epoch stellar line-of-sight velocities to determine the intrinsic velocity dispersion of the ultrafaint satellites of the Milky Way. Our method includes a simultaneous Bayesian analysis of both membership probabilities and the contribution of binary orbital motion to the observed velocity dispersion within a 14-parameter likelihood. We apply our method to the Segue 1 dwarf galaxy and conclude that Segue 1 is a dark-matter-dominated galaxy at high probability with an intrinsic velocity dispersion of 3.7^{+1.4}_{-1.1} km/sec. The dark matter halo required to produce this dispersion must have an average density of 2.5^{+4.1}_{-1.9} solar mass/pc^3 within a sphere that encloses half the galaxy's stellar luminosity. This is the highest measured density of dark matter in the Local Group. Our results show that a significant fraction of the stars in Segue 1 may be binaries with the most probable mean period close to 10 years, but also consistent with the 180 year mean period seen in the solar vicinity at about 1 sigma. Despite this binary population, the possibility that Segue 1 is a bound star cluster with the observed velocity dispersion arising from the orbital motion of binary stars is disfavored by the multi-epoch stellar velocity data at greater than 99% C.L. Finally, our treatment yields a projected (two-dimensional) half-light radius for the stellar profile of Segue 1 of 28^{+5}_{-4} pc, in excellent agreement with photometric measurements.

We investigate the kinematic properties and stellar population of the Galactic satellite Willman 1 (Wil 1) by combining Keck/DEIMOS spectroscopy with KPNO mosaic camera imaging. Wil 1 is an ultra-low luminosity Milky Way companion. This object lies in a region of size-luminosity space (M_V ~ -2 mag, d ~ 38 kpc, r_half ~ 20 pc) also occupied by the Galactic satellites Bo\"otes II and Segue 1 and 2, but no other known old stellar system. We use kinematic and color-magnitude criteria to identify 45 stars as possible members of Wil 1. With a systemic velocity of v_helio = -12.8 +/- 1.0 km/s, Wil 1 stars have velocities similar to those of foreground Milky Way stars. Informed by Monte-Carlo simulations, we identify 5 of the 45 candidate member stars as likely foreground contaminants. We confirm a significant spread in the abundances of the likely Wil 1 red giant branch members ([Fe/H] = -1.73 +/- 0.12 and -2.65 +/- 0.12, [Ca/Fe] = -0.4 +/- 0.18 and +0.13 +/- 0.28). This spread supports the scenario that Wil 1 is an ultra-low luminosity dwarf galaxy rather than a star cluster. Wil 1's innermost stars move with radial velocities offset by 8 km/s from its outer stars and have a velocity dispersion consistent with 0 km/s, suggesting that Wil 1 may not be in dynamical equilibrium. The combination of the foreground contamination and unusual kinematic distribution make it difficult to robustly determine the dark matter mass of Wil 1. As a result, X-ray or gamma-ray observations of Wil 1 that attempt to constrain models of particle dark matter using an equilibrium mass model are strongly affected by the systematics in the observations presented here. We conclude that, despite the unusual features in the Wil 1 kinematic distribution, evidence indicates that this object is, or at least once was, a dwarf galaxy.

We investigate the kinematic and photometric properties of the Segue 3 Milky Way companion using Keck/DEIMOS spectroscopy and Magellan/IMACS g and r-band imaging. Using maximum likelihood methods to analyze the photometry, we study the structure and stellar population of Segue 3. We find the half-light radius of Segue 3 is 26" +/- 5" (2.1 +/- 0.4 pc, for a distance of 17 kpc) and the absolute magnitude is a mere M_V = 0.0 +/- 0.8 mag, making Segue 3 the least luminous old stellar system known. We find Segue 3 to be consistent with a single stellar population, with an age of 12.0 +1.5/-0.4 Gyr and an [Fe/H] of -1.7 +0.07/-0.27. Line-of-sight velocities from the spectra are combined with the photometry to determine a sample of 32 stars which are likely associated with Segue 3. The member stars within three half-light radii have a velocity dispersion of 1.2 +/- 2.6 km/s. Photometry of the members indicates the stellar population has a spread in [Fe/H] of <0.3 dex. These facts, together with the small physical size of Segue 3, imply the object is likely an old, faint stellar cluster which contains no significant dark matter. We find tentative evidence for stellar mass loss in Segue 3 through the eleven candidate member stars outside of three half-light radii, as expected from dynamical arguments. Interpretation of the data outside of three half-light radii, is complicated by the object's spatial coincidence with a previously known halo substructure, which may enhance contamination of our member sample.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

We present the results of a comprehensive Keck/DEIMOS spectroscopic survey of the ultra-faint Milky Way satellite galaxy Segue 1. We have obtained velocity measurements for 98.2% of the stars within 67 pc (10 arcmin, or 2.3 half-light radii) of the center of Segue 1 that have colors and magnitudes consistent with membership, down to a magnitude limit of r=21.7. Based on photometric, kinematic, and metallicity information, we identify 71 stars as probable Segue 1 members, including some as far out as 87 pc. After correcting for the influence of binary stars using repeated velocity measurements, we determine a velocity dispersion of 3.7^{+1.4}_{-1.1} km/s, with a corresponding mass within the half-light radius of 5.8^{+8.2}_{-3.1} x 10^5 Msun. The stellar kinematics of Segue 1 require very high mass-to-light ratios unless the system is far from dynamical equilibrium, even if the period distribution of unresolved binary stars is skewed toward implausibly short periods. With a total luminosity less than that of a single bright red giant and a V-band mass-to-light ratio of 3400 Msun/Lsun, Segue 1 is the darkest galaxy currently known. We critically re-examine recent claims that Segue 1 is a tidally disrupting star cluster and that kinematic samples are contaminated by the Sagittarius stream. The extremely low metallicities ([Fe/H] < -3) of two Segue 1 stars and the large metallicity spread among the members demonstrate conclusively that Segue 1 is a dwarf galaxy, and we find no evidence in favor of tidal effects. We also show that contamination by the Sagittarius stream has been overestimated. Segue 1 has the highest measured dark matter density of any known galaxy and will therefore be a prime testing ground for dark matter physics and galaxy formation on small scales.

We present a catalog of Fe, Mg, Si, Ca, and Ti abundances for 2961 red giant stars that are likely members of eight dwarf satellite galaxies of the Milky Way (MW): Sculptor, Fornax, Leo I, Sextans, Leo II, Canes Venatici I, Ursa Minor, and Draco. For the purposes of validating our measurements, we also observed 445 red giants in MW globular clusters and 21 field red giants in the MW halo. The measurements are based on Keck/DEIMOS medium-resolution spectroscopy combined with spectral synthesis. We estimate uncertainties in [Fe/H] by quantifying the dispersion of [Fe/H] measurements in a sample of stars in monometallic globular clusters. We estimate uncertainties in Mg, Si, Ca, and Ti abundances by comparing our medium-resolution spectroscopic measurements to high-resolution spectroscopic abundances of the same stars. For this purpose, our DEIMOS sample included 132 red giants with published high-resolution spectroscopy in globular clusters, the MW halo field, and dwarf galaxies. The standard deviations of the differences in [Fe/H] and [alpha/Fe] (the average of [Mg/Fe], [Si/Fe], [Ca/Fe], and [Ti/Fe]) between the two samples is 0.15 and 0.16, respectively. This catalog represents the largest sample of multi-element abundances in dwarf galaxies to date. The next papers in this series draw conclusions on the chemical evolution, gas dynamics, and star formation histories from the catalog presented here. The wide range of dwarf galaxy luminosity reveals the dependence of dwarf galaxy chemical evolution on galaxy stellar mass.

We isolate a sample of 43 upper RGB stars in the extreme outer halo Galactic globular cluster NGC 2419 from two Keck/DEIMOS slitmasks. The probability that there is more than one contaminating halo field star in this sample is extremely low. Analysis of moderate resolution spectra of these cluster members, as well as of our Keck/HIRES high resolution spectra of a subsample of them, demonstrates that there is a small but real spread in Ca abundance of ~ 0.2 dex within this massive metal-poor globular cluster. This provides additional support to earlier suggestions that NGC 2419 is the remnant of a dwarf galaxy accreted long ago by the Milky Way.

We derive an accurate mass estimator for dispersion-supported stellar systems and demonstrate its validity by analyzing resolved line-of-sight velocity data for globular clusters, dwarf galaxies, and elliptical galaxies. Specifically, by manipulating the spherical Jeans equation we show that the dynamical mass enclosed within the 3D deprojected half-light radius r_1/2 can be determined with only mild assumptions about the spatial variation of the stellar velocity dispersion anisotropy. We find M_1/2 = 3 \sigma_los^2 r_1/2 / G ~ 4 \sigma_los^2 R_eff / G, where \sigma_los^2 is the luminosity-weighted square of the line-of-sight velocity dispersion and R_eff is the 2D projected half-light radius. While deceptively familiar in form, this formula is not the virial theorem, which cannot be used to determine accurate masses unless the radial profile of the total mass is known a priori. We utilize this finding to show that all of the Milky Way dwarf spheroidal galaxies (MW dSphs) are consistent with having formed within a halo of mass approximately 3 x 10^9 M_sun in Lambda CDM cosmology. The faintest MW dSphs seem to have formed in dark matter halos that are at least as massive as those of the brightest MW dSphs, despite the almost five orders of magnitude spread in luminosity. We expand our analysis to the full range of observed dispersion-supported stellar systems and examine their I-band mass-to-light ratios (M/L). The M/L vs. M_1/2 relation for dispersion-supported galaxies follows a U-shape, with a broad minimum near M/L ~ 3 that spans dwarf elliptical galaxies to normal ellipticals, a steep rise to M/L ~ 3,200 for ultra-faint dSphs, and a more shallow rise to M/L ~ 800 for galaxy cluster spheroids.

We use new kinematic data from the ultra-faint Milky Way satellite Segue 1 to model its dark matter distribution and derive upper limits on the dark matter annihilation cross-section. Using gamma-ray flux upper limits from the Fermi satellite and MAGIC, we determine cross-section exclusion regions for dark matter annihilation into a variety of different particles including charged leptons. We show that these exclusion regions are beginning to probe the regions of interest for a dark matter interpretation of the electron and positron fluxes from PAMELA, Fermi, and HESS, and that future observations of Segue 1 have strong prospects for testing such an interpretation. We additionally discuss prospects for detecting annihilation with neutrinos using the IceCube detector, finding that in an optimistic scenario a few neutrino events may be detected. Finally we use the kinematic data to model the Segue 1 dark matter velocity dispersion and constrain Sommerfeld enhanced models.

We present new Keck/DEIMOS spectroscopic observations of hundreds of individual stars along the sightline to Andromeda's first three discovered dwarf spheroidal galaxies (dSphs) - And I, II, and III, and leverage recent observations by our team of three additional dSphs, And VII, X, and XIV, as a part of the SPLASH Survey. Member stars of each dSph are isolated from foreground Milky Way dwarf and M31 field contamination using a variety of photometric and spectroscopic diagnostics. Our final spectroscopic sample of member stars in each dSph, for which we measure accurate radial velocities with a median uncertainty (random plus systematic errors) of 4 - 5 km/s, includes 80 red giants in And I, 95 in And II, 43 in And III, 18 in And VII, 22 in And X, and 38 in And XIV. The sample of confirmed members in the six dSphs are used to derive each system's mean radial velocity, intrinsic central velocity dispersion, mean abundance, abundance spread, and dynamical mass. This combined data set presents us with a unique opportunity to perform the first systematic comparison of the global properties (e.g., metallicities, sizes, and dark matter masses) of one-third of Andromeda's total known dSph population with Milky Way counterparts of the same luminosity. We discuss both the luminosity-metallicity relation and the luminosity-size relation of these satellites, and find that the chemical evolution histories of each host's satellites is similar. The dynamical mass estimates of M31's dSphs are similar or smaller than Milky Way dSphs of the same luminosity despite their sizes being similar or larger, suggesting M31 dSphs are less dense than Milky Way counterparts. The implications of these results for general understanding of galaxy formation and evolution is summarized. Abridged.

We present kinematic and metallicity profiles for the M31 dwarf elliptical (dE) satellite galaxies NGC 147 and NGC 185. The profiles represent the most extensive spectroscopic radial coverage for any dE galaxy, extending to a projected distance of eight half-light radii (8 r_eff = 14'). We achieve this coverage via Keck/DEIMOS multislit spectroscopic observations of 520 and 442 member red giant branch stars in NGC 147 and NGC 185, respectively. In contrast to previous studies, we find that both dEs have significant internal rotation. We measure a maximum rotational velocity of 17+/-2 km/s for NGC 147 and 15+/-5 km/s for NGC 185. The velocity dispersions decrease gently with radius with an average dispersion of 16+/-1 km/s for NGC 147 and 24+/-1 km/s for NGC 185. Both dEs have internal metallicity dispersions of 0.5 dex, but show no evidence for a radial metallicity gradient. We construct two-integral axisymmetric dynamical models and find that the observed kinematical profiles cannot be explained without modest amounts of non-baryonic dark matter. We measure central mass-to-light ratios of ML_V = 4.2+/-0.6 and ML_V = 4.6+/-0.6 for NGC 147 and NGC 185, respectively. Both dE galaxies are consistent with being primarily flattened by their rotational motions, although some anisotropic velocity dispersion is needed to fully explain their observed shapes. The velocity profiles of all three Local Group dEs (NGC 147, NGC 185 and NGC 205) suggest that rotation is more prevalent in the dE galaxy class than previously assumed, but is often manifest only at several times the effective radius. Since all dEs outside the Local Group have been probed to only inside the effective radius, this opens the door for formation mechanisms in which dEs are transformed or stripped versions of gas-rich rotating progenitor galaxies.

We present Keck/HIRES observations of six metal-poor stars in two of the ultra-faint dwarf galaxies orbiting the Milky Way, Ursa Major II and Coma Berenices. These observations include the first high-resolution spectroscopic observations of extremely metal-poor stars ([Fe/H]<-3.0) stars not belonging to the Milky Way (MW) halo field star population. We obtain abundance measurements and upper limits for 26 elements between carbon and europium. The entire sample of stars spans a range of -3.2<[Fe/H]<-2.3, and we confirm that each galaxy contains a large intrinsic spread of Fe abundances. A comparison with MW halo stars of similar metallicities reveals substantial agreement between the abundance patterns of the ultra-faint dwarf galaxies and the MW halo for the light, alpha and iron-peak elements (C to Zn). This agreement contrasts with the results of earlier studies of more metal-rich stars (-2.5<[Fe/H]<-1.0) in more luminous dwarf spheroidal galaxies (dSphs), which found significant abundance discrepancies with respect to the MW halo data. The abundances of neutron-capture elements (Sr to Eu) in the ultra-faint dwarf galaxies are extremely low, consistent with the most metal-poor halo stars, but not with the typical halo abundance pattern at [Fe/H]>-3.0. Our results are broadly consistent with a galaxy formation model that predicts that massive dwarf galaxies are the source of the metal-rich component ([Fe/H]>-2.5) of the MW halo, but we also suggest that the faintest known dwarfs may be the primary contributors to the metal-poor end of the MW halo metallicity distribution.

We present deep CFHT/MegaCam photometry of the ultra-faint Milky Way satellite galaxies Coma Berenices (ComBer) and Ursa Major II (UMa II). These data extend to r~25, corresponding to three magnitudes below the main sequence turn-offs in these galaxies. We robustly calculate a total luminosity of M_V=-3.8 +/- 0.6 for ComBer and M_V=-3.9 +/- 0.5 for UMa II, in agreement with previous results. ComBer shows a fairly regular morphology with no signs of active tidal stripping down to a surface brightness limit of 32.4 magarcsec^-2. Using a maximum likelihood analysis, we calculate the half-light radius of ComBer to be r_half=74 +/- 4 pc (5.8 +/- 0.3 arcmin) and its ellipticity e=0.36 +/- 0.04. In contrast, UMa II shows signs of on-going disruption. We map its morphology down to mu_V=32.6 mag arcsec^-2 and found that UMa II is larger than previously determined, extending at least ~700 pc (1.2 deg on the sky) and it is also quite elongated with an ellipticity of e=0.50 +/- 0.2. However, our estimate for the half-light radius, 123 +/- 3 pc (14.1 +/- 0.3 arcmin) is similar to previous results. We discuss the implications of these findings in the context of potential indirect dark matter detections and galaxy formation. We conclude that while ComBer appears to be a stable dwarf galaxy, UMa II shows signs of on-going tidal interaction.

The giant southern stream (GSS) is the most prominent tidal debris feature in M31's stellar halo. The GSS is composed of a relatively metal-rich, high surface-brightness "core" and a lower metallicity, lower surface brightness "envelope." We present Keck/DEIMOS spectroscopy of red giant stars in six fields in the vicinity of M31's GSS and one field on Stream C, an arc-like feature on M31's SE minor axis at R=60 kpc. Several GSS-related findings and measurements are presented here. We present the innermost kinematical detection of the GSS core to date (R=17 kpc). This field also contains the continuation of a second kinematically cold component originally seen in a GSS core field at R=21 kpc. The velocity gradients of the GSS and the second component in the combined data set are parallel over a radial range of 7 kpc, suggesting a possible bifurcation in the line-of-sight velocities of GSS stars. We also present the first kinematical detection of substructure in the GSS envelope. Using kinematically identified samples, we show that the envelope debris has a ~0.7 dex lower mean photometric metallicity and possibly higher intrinsic velocity dispersion than the GSS core. The GSS is also identified in the field of the M31 dSph satellite And I; the GSS in this field has a metallicity distribution identical to that of the GSS core. We confirm the presence of two kinematically cold components in Stream C, and measure intrinsic velocity dispersions of ~10 and ~4 km/s. This compilation of the kinematical (mean velocity, intrinsic velocity dispersion) and chemical properties of stars in the GSS core and envelope, coupled with published surface brightness measurements and wide-area star-count maps, will improve constraints on the orbit and internal structure of the dwarf satellite progenitor.

Andromeda X (And X) is a newly discovered low-luminosity M31 dwarf spheroidal galaxy (dSph) found by Zucker et al. (2007) in the Sloan Digital Sky Survey (SDSS - York et al. 2000). In this paper, we present the first spectroscopic study of individual red giant branch stars in And X, as a part of the SPLASH Survey (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo). Using the Keck II telescope and multiobject DEIMOS spectrograph, we target two spectroscopic masks over the face of the galaxy and measure radial velocities for ~100 stars with a median accuracy of sigma_v ~ 3 km/s. The velocity histogram for this field confirms three populations of stars along the sight line: foreground Milky Way dwarfs at small negative velocities, M31 halo red giants over a broad range of velocities, and a very cold velocity ``spike'' consisting of 22 stars belonging to And X with v_rad = -163.8 +/- 1.2 km/s. By carefully considering both the random and systematic velocity errors of these stars (e.g., through duplicate star measurements), we derive an intrinsic velocity dispersion of just sigma_v = 3.9 +/- 1.2 km/s for And X, which for its size, implies a minimum mass-to-light ratio of M/L =37^{+26}_{-19} assuming the mass traces the light. Based on the clean sample of member stars, we measure the median metallicity of And X to be [Fe/H] = -1.93 +/- 0.11, with a slight radial metallicity gradient. The dispersion in metallicity is large, sigma([Fe/H]) = 0.48, possibly hinting that the galaxy retained much of its chemical enrichment products. We discuss the potential for better understanding the formation and evolution mechanisms for M31's system of dSphs through (current) kinematic and chemical abundance studies, especially in relation to the Milky Way sample. (abridged version)

We present measurements of Fe, Mg, Si, Ca, and Ti abundances for 388 radial velocity member stars in the Sculptor dwarf spheroidal galaxy (dSph), a satellite of the Milky Way. This is the largest sample of individual alpha element (Mg, Si, Ca, Ti) abundance measurements in any single dSph. The measurements are made from Keck/DEIMOS medium-resolution spectra (6400-9000 A, R ~ 6500). Based on comparisons to published high-resolution (R >~ 20000) spectroscopic measurements, our measurements have uncertainties of sigma([Fe/H]) = 0.14 and sigma([alpha/Fe]) = 0.13. The Sculptor [Fe/H] distribution has a mean <[Fe/H]> = -1.58 and is asymmetric with a long, metal-poor tail, indicative of a history of extended star formation. Sculptor has a larger fraction of stars with [Fe/H] < -2 than the Milky Way halo. We have discovered one star with [Fe/H] = -3.80 +/- 0.28, which is the most metal-poor star known anywhere except the Milky Way halo, but high-resolution spectroscopy is needed to measure this star's detailed abundances. As has been previously reported based on high-resolution spectroscopy, [alpha/Fe] in Sculptor falls as [Fe/H] increases. The metal-rich stars ([Fe/H] ~ -1.5) have lower [alpha/Fe] than Galactic halo field stars of comparable metallicity. This indicates that star formation proceeded more gradually in Sculptor than in the Galactic halo. We also observe radial abundance gradients of -0.030 +/- 0.003 dex per arcmin in [Fe/H] and +0.013 +/- 0.003 dex per arcmin in [alpha/Fe] out to 11 arcmin (275 pc). Together, these measurements cast Sculptor and possibly other surviving dSphs as representative of the dwarf galaxies from which the metal-poor tail of the Galactic halo formed.

Over the past five years, searches in Sloan Digital Sky Survey data have more than doubled the number of known dwarf satellite galaxies of the Milky Way, and have revealed a population of ultra-faint galaxies with luminosities smaller than typical globular clusters, L ~ 1000 Lsun. These systems are the faintest, most dark matter dominated, and most metal poor galaxies in the universe. Completeness corrections suggest that we are poised on the edge of a vast discovery space in galaxy phenomenology, with hundreds more of these extreme galaxies to be discovered as future instruments hunt for the low-luminosity threshold of galaxy formation. Dark matter dominated dwarfs of this kind probe the small-scale power-spectrum, provide the most stringent limits on the phase-space packing of dark matter, and offer a particularly useful target for dark matter indirect detection experiments. Full use of dwarfs as dark matter laboratories will require synergy between deep, large-area photometric searches; spectroscopic and astrometric follow-up with next-generation optical telescopes; and subsequent observations with gamma-ray telescopes for dark matter indirect detection.

The hierarchical theory of galaxy formation rests on the idea that smaller galactic structures merge to form the galaxies that we see today. The past decade has provided remarkable observational support for this scenario, driven in part by advances in spectroscopic instrumentation. Multi-object spectroscopy enabled the discovery of kinematically cold substructures around the Milky Way and M31 that are likely the debris of disrupting satellites. Improvements in high-resolution spectroscopy have produced key evidence that the abundance patterns of the Milky Way halo and its dwarf satellites can be explained by Galactic chemical evolution models based on hierarchical assembly. These breakthroughs have depended almost entirely on observations of nearby stars in the Milky Way and luminous red giant stars in M31 and Local Group dwarf satellites. In the next decade, extremely large telescopes will allow observations far down the luminosity function in the known dwarf galaxies, and they will enable observations of individual stars far out in the Galactic halo. The chemical abundance census now available for the Milky Way will become possible for our nearest neighbor, M31. Velocity dispersion measurements now available in M31 will become possible for systems beyond the Local Group such as Sculptor and M81 Group galaxies. Detailed studies of a greater number of individual stars in a greater number of spiral galaxies and their satellites will test hierarchical assembly in new ways because dynamical and chemical evolution models predict different outcomes for halos of different masses in different environments.

We present Hubble Space Telescope (HST) WFPC2 photometry of 13 microlensed source stars from the 5.7 year Large Magellanic Cloud (LMC) survey conducted by the MACHO Project. The microlensing source stars are identified by deriving accurate centroids in the ground-based MACHO images using difference image analysis (DIA) and then transforming the DIA coordinates to the HST frame. None of these sources is coincident with a background galaxy, which rules out the possibility that the MACHO LMC microlensing sample is contaminated with misidentified supernovae or AGN in galaxies behind the LMC. This supports the conclusion that the MACHO LMC microlensing sample has only a small amount of contamination due to non-microlensing forms of variability. We compare the WFPC2 source star magnitudes with the lensed flux predictions derived from microlensing fits to the light curve data. In most cases the source star brightness is accurately predicted. Finally, we develop a statistic which constrains the location of the Large Magellanic Cloud (LMC) microlensing source stars with respect to the distributions of stars and dust in the LMC and compare this to the predictions of various models of LMC microlensing. This test excludes at > 90% confidence level models where more than 80% of the source stars lie behind the LMC. Exotic models that attempt to explain the excess LMC microlensing optical depth seen by MACHO with a population of background sources are disfavored or excluded by this test. Models in which most of the lenses reside in a halo or spheroid distribution associated with either the Milky Way or the LMC are consistent which these data, but LMC halo or spheroid models are favored by the combined MACHO and EROS microlensing results.

We present Keck/DEIMOS spectroscopy of Segue 1, an ultra-low luminosity (M_V = -1.5) Milky Way satellite companion. While the combined size and luminosity of Segue 1 are consistent with either a globular cluster or a dwarf galaxy, we present spectroscopic evidence that this object is a dark matter-dominated dwarf galaxy. We identify 24 stars as members of Segue 1 with a mean heliocentric recession velocity of 206 +/- 1.3 kms. We measure an internal velocity dispersion of 4.3+/-1.2 kms. Under the assumption that these stars are in dynamical equilibrium, we infer a total mass of 4.5^{+4.7}_{-2.5} x 10^5 Msun in the case where mass-follow-light; using a two-component maximum likelihood model, we determine a similar mass within the stellar radius of 50 pc. This implies a mass-to-light ratio of ln(M/L_V) = 7.2^{+1.1}_{-1.2} or M/L_V = 1320^{+2680}_{-940}. The error distribution of the mass-to-light ratio is nearly log-normal, thus Segue 1 is dark matter-dominated at a high significance. Using spectral synthesis modeling, we derive a metallicity for the single red giant branch star in our sample of [Fe/H]=-3.3 +/- 0.2 dex. Finally, we discuss the prospects for detecting gamma-rays from annihilation of dark matter particles and show that Segue 1 is the most promising satellite for indirect dark matter detection. We conclude that Segue 1 is the least luminous of the ultra-faint galaxies recently discovered around the Milky Way, and is thus the least luminous known galaxy.

The Milky Way has at least twenty-three known satellite galaxies that shine with luminosities ranging from about a thousand to a billion times that of the Sun. Half of these galaxies were discovered in the past few years in the Sloan Digital Sky Survey, and they are among the least luminous galaxies in the known Universe. A determination of the mass of these galaxies provides a test of galaxy formation at the smallest scales and probes the nature of the dark matter that dominates the mass density of the Universe. Here we use new measurements of the velocities of the stars in these galaxies to show that they are consistent with them having a common mass of about 10^7 M_\odot within their central 300 parsecs. This result demonstrates that the faintest of the Milky Way satellites are the most dark matter-dominated galaxies known, and could be a hint of a new scale in galaxy formation or a characteristic scale for the clustering of dark matter.

We present new metallicity measurements for 298 individual red giant branch stars in eight of the least luminous dwarf spheroidal galaxies (dSphs) in the Milky Way (MW) system. Our technique is based on medium resolution Keck/DEIMOS spectroscopy coupled with spectral synthesis. We present the first spectroscopic metallicities at [Fe/H] < -3.0 of stars in a dwarf galaxy, with individual stellar metallicities as low as [Fe/H] = -3.3. Because our [Fe/H] measurements are not tied to empirical metallicity calibrators and are sensitive to arbitrarily low metallicities, we are able to probe this extremely metal-poor regime accurately. The metallicity distribution of stars in these dSphs is similar to the MW halo at the metal-poor end. We also demonstrate that the luminosity-metallicity relation previously seen in more luminous dSph galaxies (M_V = -13.4 to -8.8) extends smoothly down to an absolute magnitude of M_V = -3.7. The discovery of extremely metal-poor stars in dSphs lends support to the LCDM galaxy assembly paradigm wherein dwarf galaxies dissolve to form the stellar halo of the MW.

NGC 205, a close satellite of the M31 galaxy, is our nearest example of a dwarf elliptical galaxy. Photometric and kinematic observations suggest that NGC 205 is undergoing tidal distortion from its interaction with M31. Despite earlier attempts, the orbit and progenitor properties of NGC 205 are not well known. We perform an optimized search for these unknowns by combining a genetic algorithm with restricted N-body simulations of the interaction. This approach, coupled with photometric and kinematic observations as constraints, allows for an effective exploration of the parameter space. We represent NGC 205 as a static Hernquist potential with embedded massless test particles that serve as tracers of surface brightness. We explore 3 distinct, initially stable configurations of test particles: cold rotating disk, warm rotating disk, and hot, pressure-supported spheroid. Each model reproduces some, but not all, of the observed features of NGC 205, leading us to speculate that a rotating progenitor with substantial pressure support could match all of the observables. Furthermore, plausible combinations of mass and scale length for the pressure-supported spheroid progenitor model reproduce the observed velocity dispersion profile. For all 3 models, orbits that best match the observables place the satellite 11+/-9 kpc behind M31 moving at very large velocities: 300-500 km/s on primarily radial orbits. Given that the observed radial component is only 54 km/s, this implies a large tangential motion for NGC 205, moving from the NW to the SE. These results suggest NGC 205 is not associated with the stellar arc observed to the NE of NGC 205. Furthermore, NGC 205's velocity appears to be near or greater than its escape velocity, signifying that the satellite is likely on its first M31 passage.

In the course of our survey of the outer halo of the Andromeda Galaxy we have discovered a remote, possible satellite of that system at a projected 162 kpc (11.7 degrees) radius. The fairly elongated (0.31 +/- 0.09 ellipticity) dwarf can be fit with a King profile of 1.07 kpc (d/784 kpc) limiting radius, where the satellite distance, d, is estimated at ~630-850 kpc from the tip of the red giant branch. The newfound galaxy, ``Andromeda XIV'' (``AndXIV''), distinguishes itself from other Local Group galaxies by its extreme dynamics: Keck/DEIMOS spectroscopy reveals it to have a large heliocentric radial velocity (-481 km/s), or -206 km/s velocity relative to M31. Even at its projected radius AndXIV already is at the M31 escape velocity based on the latest M31 mass models. If AndXIV is bound to M31, then recent models with reduced M31 virial masses need revision upward. If not bound to M31, then AndXIV is just now falling into the Local Group for the first time and represents a dwarf galaxy that formed and spent almost its entire life in isolation.

We use kinematic data from three new, nearby, extremely low-luminosity Milky Way dwarf galaxies (Ursa Major II, Willman 1, and Coma Berenices) to constrain the properties of their dark matter halos, and from these make predictions for the gamma-ray flux from annihilation of dark matter particles in these halos. We show that these 10^3 solar luminosity dwarfs are the most dark matter dominated galaxies in the Universe, with total masses within 100 pc in excess of 10^6 solar masses. Coupled with their relative proximity, their large masses imply that they should have mean gamma-ray fluxes comparable to or greater than any other known satellite galaxy of the Milky Way. Our results are robust to both variations of the inner slope of the density profile and the effect of tidal interactions. The fluxes could be boosted by up to two orders of magnitude if we include the density enhancements caused by surviving dark matter substructure.

We present Keck/DEIMOS spectroscopy of stars in 8 of the newly discovered ultra-faint dwarf galaxies around the Milky Way. We measure the velocity dispersions of Canes Venatici I and II, Ursa Major I and II, Coma Berenices, Hercules, Leo IV and Leo T from the velocities of 18 - 214 stars in each galaxy and find dispersions ranging from 3.3 to 7.6 km/s. The 6 galaxies with absolute magnitudes M_V < -4 are highly dark matter-dominated, with mass-to-light ratios approaching 1000. The measured velocity dispersions are inversely correlated with their luminosities, indicating that a minimum mass for luminous galactic systems may not yet have been reached. We also measure the metallicities of the observed stars and find that the 6 brightest of the ultra-faint dwarfs extend the luminosity-metallicity relationship followed by brighter dwarfs by 2 orders of magnitude in luminosity; several of these objects have mean metallicities as low as [Fe/H] = -2.3 and therefore represent some of the most metal-poor known stellar systems. We detect metallicity spreads of up to 0.5 dex in several objects, suggesting multiple star formation epochs. Having established the masses of the ultra-faint dwarfs, we re-examine the missing satellite problem. After correcting for the sky coverage of the SDSS, we find that the ultra-faint dwarfs substantially alleviate the discrepancy between the predicted and observed numbers of satellites around the Milky Way, but there are still a factor of ~4 too few dwarf galaxies over a significant range of masses. We show that if galaxy formation in low-mass dark matter halos is strongly suppressed after reionization, the simulated circular velocity function of CDM subhalos can be brought into approximate agreement with the observed circular velocity function of Milky Way satellite galaxies. [slightly abridged]

The galaxy circular velocity function at small masses is related to the matter power spectrum on small scales. Although this function is well-studied for Local Group dwarfs, theoretical predictions and observational measurements are difficult for satellite galaxies, because of ram pressure and tidal stripping. By contrast, isolated dwarf galaxies are less affected by these processes, and almost always have enough 21cm emission to trace their dynamics robustly. Here, we test cold dark matter cosmology using isolated low mass dwarf galaxies from the SDSS with measured 21cm widths. We find consistency between the predicted and observed number density of isolated galaxies down to circular velocities of 50 km/s. Our technique yields a direct test of small-scale cosmology independent of the Lyman-alpha forest power spectrum, but our sample is currently statistically less powerful: warm dark matter particles heavier than 0.5 keV cannot be ruled out. Our major systematic uncertainty is the surface brightness limit of the SDSS. Blind HI surveys, such as the ALFALFA survey on Arecibo, will uncover a larger number of isolated low mass galaxies and increase the power of our constraints. With our sample, we also find that the Tully-Fisher relation for dwarf galaxies is a strong function of environment, and that the baryonic fraction is only a weak function of mass. These results suggest that for dwarf galaxies, gas loss is dominated by external, not internal, processes. [abridged]

We present the discovery of a kinematically-cold stellar population along the SE minor axis of the Andromeda galaxy (M31) that is likely the forward continuation of M31's giant southern stream. This discovery was made in the course of an on-going spectroscopic survey of red giant branch (RGB) stars in M31 using the DEIMOS instrument on the Keck II 10-m telescope. Stellar kinematics are investigated in eight fields located 9-30 kpc from M31's center (in projection). A likelihood method based on photometric and spectroscopic diagnostics is used to isolate confirmed M31 RGB stars from foreground Milky Way dwarf stars: for the first time, this is done without using radial velocity as a selection criterion, allowing an unbiased study of M31's stellar kinematics. The radial velocity distribution of the 1013 M31 RGB stars shows evidence for the presence of two components. The broad (hot) component has a velocity dispersion of 129 km/s and presumably represents M31's virialized spheroid. A significant fraction (19%) of the population is in a narrow (cold) component centered near M31's systemic velocity with a velocity dispersion that decreases with increasing radial distance, from 55.5 km/s at R_proj=12 kpc to 10.6 km/s at R_proj=18 kpc. The spatial and velocity distribution of the cold component matches that of the "Southeast shelf" predicted by the Fardal et al. (2007) orbital model of the progenitor of the giant southern stream. The metallicity distribution of the cold component matches that of the giant southern stream, but is about 0.2 dex more metal rich on average than that of the hot spheroidal component. We discuss the implications of our discovery on the interpretation of the intermediate-age spheroid population found in this region in recent ultra-deep HST imaging studies.

Using the HST ACS, we have obtained deep optical images of a southeast minor-axis field in the Andromeda Galaxy, 21 kpc from the nucleus. In both star counts and metallicity, this field represents a transition zone between the metal-rich, highly-disturbed inner spheroid that dominates within 15 kpc and the metal-poor, diffuse population that dominates beyond 30 kpc. The color-magnitude diagram reaches well below the oldest main-sequence turnoff in the population, allowing a reconstruction of the star formation history in this field. Compared to the spheroid population at 11 kpc, the population at 21 kpc is ~1.3 Gyr older and ~0.2 dex more metal-poor, on average. However, like the population at 11 kpc, the population at 21 kpc exhibits an extended star formation history; one third of the stars are younger than 10 Gyr, although only a few percent are younger than 8 Gyr. The relatively wide range of metallicity and age is inconsistent with a single, rapid star-formation episode, and instead suggests that the spheroid even at 21 kpc is dominated by the debris of earlier merging events likely occurring more than 8 Gyr ago.

We investigate the gas content and baryonic Tully-Fisher relationship for extremely low luminosity dwarf galaxies in the absolute magnitude range -13.5 > Mr > -16. The sample is selected from the Sloan Digital Sky Survey and consists of 101 galaxies for which we have obtained follow-up HI observations using the Arecibo Observatory and Green Bank Telescope. This represents the largest homogeneous sample of dwarfs at low luminosities with well-measured HI and optical properties. The sample spans a range of environments, from dense groups to truly isolated galaxies. The average neutral gas fraction is f_gas=0.6, significantly exceeding that of typical gas-rich galaxies at higher luminosities. Dwarf galaxies are therefore less efficient at turning gas into stars over their lifetimes. The strong environmental dependence of the gas fraction distribution demonstrates that while internal processes can reduce the gas fractions to roughly f_gas=0.4, external processes are required to fully remove gas from a dwarf galaxy. The average rotational velocity of our sample is vrot=50 km/s. Including more massive galaxies from the literature, we fit a baryonic Tully-Fisher slope of M_baryon \propto vrot^(3.70+/- 0.15). This slope compares well with CDM models that assume an equal baryon to dark matter ratio at all masses. While gas stripping or other processes may modify the baryon to dark matter ratio for dwarfs in the densest environments, the majority of dwarf galaxies in our sample have not preferentially lost significant baryonic mass relative to more massive galaxies.

We have analysed the globular cluster (GC) system of the nucleated dwarf elliptical galaxy VCC 1087 in the Virgo cluster, based on Keck/LRIS spectroscopy and archival HST/ACS imaging. We estimate VCC 1087 hosts a total population of 77+/-19 GCs, which corresponds to a relatively high V-band specific frequency of 5.8+/-1.4. The g-z color distribution of the GCs shows a blue (metal-poor) peak with a tail of redder (metal-rich) clusters similar in color to those seen in luminous ellipticals. Spectroscopy of a subsample of 12 GCs suggests that the GC system is old and coeval (~10 Gyr), with a fairly broad metallicity distribution (-1.8<[m/H]<-0.8). In contrast, an integrated spectrum of the underlying galaxy starlight reveals that its optical luminosity is dominated by metal-rich, intermediate-aged stars. Radial velocities of the GCs suggest rotation close to the major axis of the galaxy, and this rotation is dynamically significant with (v/sigma)^* >1. A compilation of the kinematics of the GC systems of 9 early-type galaxies shows surprising diversity in the v/sigma parameter for GC systems. In this context, the GC system of VCC 1087 exhibits the most significant rotation to velocity dispersion signature. Modeling the velocity dispersion profile of the GCs and galaxy stars suggest fairly constant mass-to-light ratios of ~3 out to 6.5 kpc. The present observations can entertain both baryonic and non-baryonic solutions, and GC velocities at larger radii would be most valuable with regard to this issue. We discuss the evolution of VCC 1087 in terms of the galaxy ``harassment'' scenario, and conclude that this galaxy may well be the remains of a faded, tidally perturbed Sc spiral [abridged].

NGC 205 is the nearest example of a dwarf elliptical (dE) galaxy and the prototype of this enigmatic galaxy class. Photometric evidence suggests that NGC 205, a close satellite of the M31 galaxy, is tidally interacting with its parent galaxy. We present stellar radial velocity measurements out to a projected radius of 20' (5 kpc) in NGC 205 based on Keck/DEIMOS multislit spectroscopic observations of 725 individual red giant branch stars. Our kinematic measurements extend from the center out to six times the effective radius of NGC 205, well past the expected tidal radius. The contamination in our kinematic sample from M31 field stars is estimated to be a few percent based on maximum likelihood fits. We measure a maximum major-axis rotation speed for the body of NGC 205 of 11+/- 5 kms; this is the first dE galaxy in which the maximum rotation velocity has been measured. Combined with the velocity dispersion, we conclude that NGC 205 is supported by a combination of rotation and random motions. At a major-axis distance of 4.5' (1 kpc), the velocity profile of NGC 205 turns over; stars beyond this radius are moving counter to the rotation of the inner part of the galaxy. The turnover radius is coincident with the onset of isophotal twisting and the estimated tidal radius, suggesting that the outer kinematics of NGC 205 are dominated by gravitational interactions with the nearby M31 galaxy. The motion of stars beyond a radius of 4.5' implies that NGC 205 is in a prograde encounter with its parent galaxy M31.

We present a measurement of the systemic proper motion of the Large Magellanic Cloud (LMC) from astrometry with the High Resolution Camera (HRC) of the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). We observed LMC fields centered on 21 background QSOs that were discovered from their optical variability in the MACHO database. The QSOs are distributed homogeneously behind the central few degrees of the LMC. With 2 epochs of HRC data and a ~2 year baseline we determine the proper motion of the LMC to better than 5% accuracy: mu_W = -2.03 +/- 0.08 mas/yr; mu_N = 0.44 +/- 0.05 mas/yr. This is the most accurate proper motion measurement for any Milky Way satellite thus far. When combined with HI data from the Magellanic Stream this should provide new constraints on both the mass distribution of the Galactic Halo and models of the Stream.

Using 7 years of MACHO survey data, we present a new determination of the optical depth to microlensing towards the Galactic bulge. We select the sample of 62 microlensing events (60 unique) on clump giant sources and perform a detailed efficiency analysis. We use only the clump giant sources because these are bright bulge stars and are not as strongly affected by blending as other events. Using a subsample of 42 clump events concentrated in just 4.5 square degrees, we find tau = 2.17^{+0.47}_{-0.38} x 10^{-6} at (l,b) = (1.50, -2.68), somewhat smaller than found in most previous MACHO studies, but in excellent agreement with recent theoretical predictions. We also present the optical depth in each of the 19 fields in which we detected events, and find limits on optical depth for fields with no events. The errors in optical depth in individual fields are dominated by Poisson noise. We measure optical depth gradients (1.06 +/- 0.71) x 10^{-6} deg^{-1} and (0.29 +/- 0.43) x 10^{-6} deg^{-1} in the galactic latitude b and longitude l directions, respectively. Finally, we discuss the possibility of anomalous duration distribution of events in the field 104 centered on (l,b) = (3.11, -3.01) as well as investigate spatial clustering of events in all fields.

We present evidence for a counter-rotating core in the low-luminosity (M_B = -18.2) elliptical galaxy NGC 770 based on internal stellar kinematic data. This counter-rotating core is unusual as NGC 770 is not the primary galaxy in the region and it lies in an environment with evidence of on-going tidal interactions. We discovered the counter-rotating core via single-slit Keck/ESI echelle spectroscopy; subsequent integral field spectroscopy was obtained with the Gemini/GMOS IFU. The counter-rotating region has a peak rotation velocity of 21 km/s as compared to the main galaxy's rotation speed of greater than 45 km/s in the opposite direction. The counter-rotating region extends to a radius of 4'' (0.6 kpc), slightly smaller than the half-light radius of the galaxy which is 5.3'' (0.8 kpc) and is confined to a disk whose scale height is less than 0.8'' (0.1 kpc). We compute an age and metallicity of the inner counter-rotating region of 3 +/- 0.5 Gyr and [Fe/H] = 0.2 +/- 0.2 dex, based on Lick absorption-line indices. The lack of other large galaxies in this region limits possible scenarios for the formation of the counter-rotating core. We discuss several scenarios and favor one in which NGC 770 accreted a small gas-rich dwarf galaxy during a very minor merging event. If this scenario is correct, it represents one of the few known examples of merging between two dwarf-sized galaxies.

We present a catalog of 450 high signal-to-noise microlensing events observed by the MACHO collaboration between 1993 and 1999. The events are distributed throughout our fields and, as expected, they show clear concentration toward the Galactic center. No optical depth is given for this sample since no blending efficiency calculation has been performed, and we find evidence for substantial blending. In a companion paper we give optical depths for the sub-sample of events on clump giant source stars, where blending is not a significant effect. Several events with sources that may belong to the Sagittarius dwarf galaxy are identified. For these events even relatively low dispersion spectra could suffice to classify these events as either consistent with Sagittarius membership or as non-Sagittarius sources. Several unusual events, such as microlensing of periodic variable source stars, binary lens events, and an event showing extended source effects are identified. We also identify a number of contaminating background events as cataclysmic variable stars.

We present results from a large spectroscopic survey of M31 red giants using the Keck telescope/DEIMOS. Photometric pre-screening, based on the 100A-wide DDO51 band centered on the Mgb/MgH feature, was used to select spectroscopic targets. Red giant candidates were targeted in a small field on M31's giant southern tidal stream at a projected distance of 31kpc from the galaxy center. We isolate a clean sample of 68 giants by removing contaminants (foreground Galactic dwarfs and background galaxies) using spectroscopic, imaging, and photometric methods. About 65% of the M31 stars are found to be members of the stream, while the rest appear to be members of the general halo population. The mean (heliocentric) radial velocity of the stream in our field is -458 km/s, or -158 km/s relative to M31's systemic velocity, in good agreement with recent measurements at other stream locations. The intrinsic velocity dispersion of the stream is constrained to be 15_{-15}^{+8} km/s (90% confidence limits). The companion paper by Font et al. (2004, astro-ph/0406146) discusses possible orbits, implications of the coldness of the stream, and progenitor satellite properties. The kinematics (and perhaps [Fe/H] distribution) of our halo sample indicate that it is different from other M31 halo samples; this may be an indication of substructure in the halo. The stream seems to have a higher mean [Fe/H] than the halo, -0.51 vs -0.74 dex, and a smaller [Fe/H] spread. The stream's high metallicity implies that its progenitor must have been a luminous dwarf galaxy. The CaII triplet strengths of the M31 giants are generally consistent with photometric estimates of their metallicity (derived by fitting RGB fiducials in the color-magnitude diagram). There is indirect evidence of intermediate-age stars in the stream.

Fourier coefficents have been derived for the $V$ and $R$ light curves of 785 overtone RR Lyrae variables in 16 MACHO fields near the bar of the LMC. The $\phi_{31}$ and $R_{21}$ coefficients have been compared with those of the first overtone RR Lyrae variables in the Galactic globular clusters NGC 6441, M107, M5, M3, M2, $\omega$ Centauri and M68. The results indicate that many of the LMC variables have properties similar to the ones in M2, M3, M5 and the Oosterhoff type I variables in $\omega$ Cen, but they are different from the Oosterhoff type II variables in $\omega$ Cen. Equations derived from hydrodynamic pulsation models have been used to calculate the luminosity and temperature for the 330 bona fide first-overtone variables. The results indicate that they have $\log L$ in the range 1.6 to $1.8\lsun$ and $\log T_{eff}$ between 3.85 and 3.87. Based on these temperatures, a mean color excess $E(V-R) =0.08$ mag, equivalent to $E(B-V)=0.14$ mag, has been estimated for these 330 stars. The 80 M5-like variables (selected according to their location in the $\phi_{31}-\log P$ plot) are used to determine a LMC distance. After correcting for the effects of extinction and crowding, a mean apparent magnitude $<V_0>=18.99 \pm 0.02$ (statistical) $\pm 0.16$ (systematic) has been estimated for these 80 stars. Combining this with a mean absolute magnitude $M_V=0.56\pm 0.06$ for M5-like stars derived from Baade-Wesselink analyses, main sequence fitting, Fourier parameters and the trigonometric parallax of RR Lyrae, we derive an LMC distance modulus $\mu=18.43\pm 0.06$ (statistical) $\pm 0.16$ (systematic) mag. The large systematic error arises from the difficulties of correcting for interstellar extinction and for crowding.

We have frequency analyzed 6391 variables classified earlier as fundamental mode RR Lyrae (RR0) stars in the MACHO database on the Large Magellanic Cloud (LMC). The overwhelming majority of these variables have been proved to be indeed RR0 stars, whereas the remaining ones have fallen in one of the following categories: single- and double-mode Cepheids, binaries, first overtone and double-mode RR Lyrae stars and non-classified variables. Special attention has been paid to the Blazhko stars. We found altogether 731 Blazhko variables showing either a doublet or an equidistant triplet pattern at the main pulsation component in their frequency spectra. The incidence rate of the Blazhko variables among the RR0 stars in the LMC is 11.9%, which is three times higher than their rate among the first overtone RR Lyrae stars. From the frequency spectra we found that variables with larger modulation amplitudes at the higher frequency side of the main pulsation component are three times more frequent than the ones showing the opposite amplitude pattern. A search made for a modulation component with the Blazhko period in the average brightness of the individual variables showed the existence of such a modulation with an overall amplitude of ~0.006 mag. On the other hand, a similar search for quadruple modulation patterns around the main pulsation component have failed to clearly detect such components at the ~0.004 mag level. This means that the amplitudes of the quadruple components should be, on the average, at least ten times smaller than those of the triplet components. This finding and the existence of Blazhko variables with highly asymmetric modulation amplitudes not only question the validity of the magnetic oblique rotator model, but also puts stringent constraints on models based on mode coupling theories.

We present spatially-resolved internal kinematics and stellar chemical abundances for a sample of dwarf elliptical (dE) galaxies in the Virgo Cluster observed with Keck/ESI. We find that 4 out of 17 dEs have major axis rotation velocities consistent with rotational flattening, while the remaining dEs have no detectable major axis rotation. Despite this difference in internal kinematics, rotating and non-rotating dEs are remarkably similar in terms of their position in the Fundamental Plane, morphological structure, stellar populations, and local environment. We present evidence for faint underlying disks and/or weak substructure in a fraction of both rotating and non-rotating dEs, but a comparable number of counter-examples exist for both types which show no evidence of such structure. Absorption-line strengths were determined based on the Lick/IDS system (Hbeta, Mgb, Fe5270, Fe5335) for the central region of each galaxy. We find no difference in the line-strength indices, and hence stellar populations, between rotating and non-rotating dE galaxies. The best-fitting mean age and metallicity for our 17 dE sample are 5 Gyr and Fe/H = -0.3 dex, respectively, with rms spreads of 3 Gyr and 0.1 dex. The majority of dEs are consistent with solar alpha/Fe abundance ratios. By contrast, the stellar populations of classical elliptical galaxies are, on average, older, more metal rich, and alpha-enhanced relative to our dE sample. The local environments of both dEs types appear to be diverse in terms of their proximity to larger galaxies in real or velocity space within the Virgo Cluster. Thus, rotating and non-rotating dEs are remarkably similar in terms of their structure, stellar content, and local environments, presenting a significant challenge to theoretical models of their formation. (abridged)

We describe a few recent microlensing results from the MACHO Collaboration. The aim of the MACHO Project was the identification and quantitative description of dark and luminous matter in the Milky Way using microlensing toward the Magellanic Clouds and Galactic bulge. We start with a discussion of the HST follow-up observations of the microlensing events toward the LMC detected in the first 5 years of the experiment. Using color-magnitude diagrams we attempt to distinguish between two possible locations of the microlensing sources: 1) in the LMC or 2) behind the LMC. We conclude that unless the extinction is extremely patchy, it is very unlikely that most of the LMC events have sources behind the LMC. During an examination of the HST images of the 13 LMC events we found a very red object next to the source star of event LMC-5. Astrometry, microlensing parallax fit, and a spectrum suggest that in this case we directly image the lens - a low-mass disk star. Then we focus on the majority of events observed by the MACHO Project, which are detected toward the Galactic bulge. We argue that the microlensing optical depth toward the bulge is best measured using events that have clump giant sources, which are almost unaffected by blending. From this sample we derive a low optical depth toward the Galactic bulge of (1.4 +/- 0.3) x 10^{-6}, in good agreement with other observational constraints and with theoretical models. The presence of many long-duration events among the bulge candidates allows us to investigate the microlensing parallax effect. Events with the strongest parallax signal are probably due to massive remnants. Events MACHO-96-BLG-5 and MACHO-98-BLG-6 might have been caused by the 6-solar-mass black holes.

Keck/ESI spectroscopy is presented for six Virgo Cluster dwarf elliptical (dE) galaxies. The mean line-of-sight velocity and velocity dispersion are resolved as a function of radius along the major axis of each galaxy, nearly doubling the total number of dEs with spatially-resolved stellar kinematics. None of the observed objects shows evidence of strong rotation: upper limits on v_rot/sigma are well below those expected for rotationally-flattened objects. Such limits place strong constraints on dE galaxy formation models. Although these galaxies continue the trend of low rotation velocities observed in Local Group dEs, they are in contrast to recent observations of large rotation velocities in slightly brighter cluster dEs. Using WFPC2 surface photometry and spherically-symmetric dynamical models, we determine global mass-to-light ratios 3 < M/L_V < 6. These ratios are comparable to those expected for an intermediate-age stellar population and are broadly consistent with the (V-I) colors of the galaxies. This implies that these dEs do not have a significant dark matter component inside an effective radius. Central black holes more massive than 10^7 M_sun can be ruled out. For the 5 nucleated dEs in our sample, we determine kinematic and photometric properties for the central nucleus separately from the underlying host dE galaxy. These nuclei are as bright or brighter than the most luminous Galactic globular clusters and lie near the region of Fundamental Plane space occupied by globular clusters. In this space, the Virgo dE galaxies lie in the same general region as Local Group and other nearby dEs, although non-rotating dEs appear to have a slightly higher mean mass and mass-to-light ratio than their rotating counterparts.

We present 47 spectroscopically-confirmed quasars discovered behind the Magellanic Clouds identified via photometric variability in the MACHO database. Thirty-eight quasars lie behind the Large Magellanic Cloud and nine behind the Small Magellanic Cloud, more than tripling the number of quasars previously known in this region. The quasars cover the redshift interval 0.2 < z < 2.8 and apparent mean magnitudes 16.6 < V < 20.1. We discuss the details of quasar candidate selection based on time variability in the MACHO database and present results of spectroscopic follow-up observations. Our follow-up detection efficiency was 20%; the primary contaminants were emission-line Be stars in Magellanic Clouds. For the 47 quasars discovered behind the Magellanic Clouds plus an additional 12 objects previously identified in this region, we present 7.5-year MACHO V- and R-band lightcurves with average sampling times of 2-10 days.

We present the photometric properties of 1279 blue variable stars within the LMC. Photometry is derived from the MACHO database. The lightcurves of the sample exhibit a variety of quasi-periodic and aperiodic outburst behavior. A characteristic feature of the photometric variation is that the objects are reddest when at maximum outburst. A subset of 102 objects were examined spectroscopically. Within this subset, 91% exhibited Balmer emission in at least one epoch, in some cases with spectacular spectral variability. The variability observed in the sample is consistent with the establishment and maintenance of the Be phenomenon.

We present an analysis of the longest timescale microlensing events discovered by the MACHO Collaboration during a 7 year survey of the Galactic bulge. We find 6 events that exhibit very strong microlensing parallax signals due, in part, to accurate photometric data from the GMAN and MPS collaborations. The microlensing parallax fit parameters are used in a likelihood analysis, which is able to estimate the distance and masses of the lens objects based upon a standard model of the Galactic velocity distribution. This analysis indicates that the most likely masses of 5 of the 6 lenses are > 1 Msun, which suggests that a substantial fraction of the Galactic lenses may be massive stellar remnants. This could explain the observed excess of long timescale microlensing events. The lenses for events MACHO-96-BLG-5 and MACHO-98-BLG-6 are the most massive, with mass estimates of M/Msun = 6 +10/-3 and M/Msun = 6 +7/-3, respectively. The observed upper limits on the absolute brightness of main sequence stars for these lenses are < 1 Lsun, so both lenses are black hole candidates. The black hole interpretation is also favored by a likelihood analysis with a Bayesian prior using a conventional model for the lens mass function. We consider the possibility that the source stars for some of these 6 events may lie in the foreground or background of the Galactic bulge, but we find that this is unlikely. Future HST observations of these events can either confirm the black hole lens hypothesis or detect the lens stars and provide a direct measurement of their masses. Future observations of similar events by SIM or the Keck or VLTI interferometers will allow direct measurements of the lens masses for stellar remnant lenses as well.

We present preliminary results of the analysis of 5 years of MACHO data on the Galactic bulge microlensing events with clump giants as sources. This class of events allows one to obtain robust conclusions because relatively bright clump stars are not strongly affected by blending. We discuss: 1) the selection of `giant' events, 2) the distribution of event durations, 3) the anomalous character of event durations and optical depth in the MACHO field 104 centered on (l,b) = (3.1, -3.0). We report the preliminary average optical depth of (2.0 +/- 0.4) x 10^{-6} (internal) at (l,b) = (3.9, -3.8), and present a map of the spatial distribution of the optical depth. When field 104 is removed from the sample, the optical depth drops to (1.4 +/- 0.3) x 10^{-6}, which is in excellent agreement with infrared-based models of the central Galactic region.

Keck spectroscopy is presented for four dwarf elliptical galaxies in the Virgo Cluster. At this distance, the mean velocity and velocity dispersion are well resolved as a function of radius between 100 to 1000 pc, allowing a clear separation between nuclear and surrounding galaxy light. We find a variety of dispersion profiles for the inner regions of these objects, and show that none of these galaxies is rotationally flattened.

We report the discovery of eight new R Coronae Borealis (RCB) stars in the Large Magellanic Cloud (LMC) using the MACHO project photometry database. The discovery of these new stars increases the number of known RCB stars in the LMC to thirteen. We have also discovered four stars similar to the Galactic variable DY Per. These stars decline much more slowly and are cooler than the RCB stars. The absolute luminosities of the Galactic RCB stars are unknown since there is no direct measurement of the distance to any Galactic RCB star. Hence, the importance of the LMC RCB stars. We find a much larger range of absolute magnitudes (M(V) = -2.5 to -5 mag) than inferred from the small pre-MACHO sample of LMC RCB stars. It is likely that there is a temperature - M(V)relationship with the cooler stars being intrinsically fainter. Cool (~5000 K) RCB stars are much more common than previously thought based on the Galactic RCB star sample. Using the fairly complete sample of RCB stars discovered in the MACHO fields, we have estimated the likely number of RCB stars in the Galaxy to be ~3,200. The SMC MACHO fields were also searched for RCB stars but none were found.

By cross-correlating the results of two recent large-scale surveys, the general properties of a well defined sample of semi-regular variable stars have been determined. ISOGAL mid-infrared photometry and MACHO lightcurves are assembled for approximately 300 stars in the Baade's Windows of low extinction towards the Galactic bulge. These stars are mainly giants of late M spectral type, evolving along the asymptotic giant branch (AGB). They are found to possess a wide and continuous distribution of pulsation periods and to obey an approximate log~period -- bolometric magnitude relation or set of such relations. Approximate mass-loss rates in the range of 1e-8 to 5e-7 M_sun per year are derived from ISOGAL mid-infrared photometry and models of stellar spectra adjusted for the presence of optically-thin circumstellar silicate dust. Mass-loss rates depend on luminosity and pulsation period. Some stars lose mass as rapidly as short-period Miras but do not show Mira-like amplitudes. A period of 70 days or longer is a necessary but not a sufficient condition for mass loss to occur. For AGB stars in the mass-loss ranges that we observe, the functional dependence of mass-loss rate on temperature and luminosity is found to be in agreement with recent theoretical predictions. If we include our mass-loss rates with a sample of extreme mass-losing AGB stars in the Large Magellanic Cloud, we get the general result for AGB stars that mass-loss rate is proportional to luminosity^{2.7}, valid for AGB stars with 10^{-8} to 10^{-4} M_sun per year (Abridged).

We attempt to determine whether the MACHO microlensing source stars are drawn from the average population of the LMC or from a population behind the LMC by examining the HST color-magnitude diagram (CMD) of microlensing source stars. We present WFPC2 HST photometry of eight MACHO microlensing source stars and the surrounding fields in the LMC. The microlensing source stars are identified by deriving accurate centroids in the ground-based MACHO images using difference image analysis (DIA) and then transforming the DIA coordinates to the HST frame. We consider in detail a model for the background population of source stars based on that presented by Zhao, Graff & Guhathakurta. In this model, the source stars have an additional reddening <E(B-V)> = 0.13 mag and a slightly larger distance modulus <Delta mu> ~ 0.3 mag than the average LMC population. We also investigate a series of source star models, varying the relative fraction of source stars drawn from the average and background populations and the displacement of the background population from the LMC. Due to the small number of analyzed events the distribution of probabilities of different models is rather flat. A shallow maximum occurs at a fraction s_LMC ~ 0.8 of the source stars in the LMC. This is consistent with the interpretation that a significant fraction of observed microlensing events are due to lenses in the Milky Way halo, but does not definitively exclude other models.

We present photometry and analysis of the microlensing alert MACHO 96-LMC-2. The ~3% photometry provided by the Global Microlensing Alert Network follow--up effort reveals a periodic modulation in the lightcurve. We attribute this to binarity of the lensed source. Microlensing fits to a rotating binary source magnified by a single lens converge on two minima, separated by delta chi^2 ~ 1. The most significant fit X1 predicts a primary which contributes ~100% of the light, a dark secondary, and an orbital period (T) of 9.2 days. The second fit X2 yields a binary source with two stars of roughly equal mass and luminosity, and T = 21.2 days. The lensed object appears to lie on the upper LMC main sequence. We estimate the mass of the primary component of the binary system, M ~2 M_sun. For the preferred model X1, we explore the range of dark companions by assuming 0.1 M_sun and 1.4 M_sun objects in models X1a and X1b, respectively. We find lens velocities projected to the LMC in these models of v^hat_X1a = 18.3 +/- 3.1 km/s and v^hat_X1b = 188 +/- 32 k/ms. In both these cases, a likelihood analysis suggests an LMC lens is preferred over a Galactic halo lens, although only marginally so in model X1b. We also find v^hat_X2 = 39.6 +/- 6.1 k/ms, where the likelihood for the lens location is strongly dominated by the LMC disk. In all cases, the lens mass is consistent with that of an M-dwarf. The LMC self-lensing rate contributed by 96-LMC-2 is consistent with model self-lensing rates. (Abridged)

We report on a search for long duration microlensing events towards the Large Magellanic Cloud. We find none, and therefore put limits on the contribution of high mass objects to the Galactic dark matter. At 95% confidence level we exclude objects in the mass range 0.3 solar masses to 30.0 solar masses from contributing more than 4 times 10^11 solar masses to the Galactic halo. Combined with earlier results, this means that objects with masses under 30 solar masses cannot make up the entire dark matter halo if the halo is of typical size. For a typical dark halo, objects with masses under 10 solar masses contribute less than 40% of the dark matter.

We present a ~ 5-yr optical light curve of the recurrent Be/X-ray transient A0538-66 obtained as a by-product of the MACHO Project. These data reveal both a long-term modulation at P = 420.8 +/- 0.8 d and a short-term modulation at 16.6510 +/- 0.0022 d which, within errors, confirms the previously found orbital period. Furthermore, the orbital activity is only seen at certain phases of the 421 d cycle suggesting that the long-term modulation is related to variations in the Be star envelope.

We present preliminary results of the analysis of 5 years of MACHO data on the Galactic bulge microlensing events with clump giants as sources. In particular, we discuss: 1) the selection of `giant' events, 2) distribution of impact parameters, 3) distribution of event durations, 4) the concentration of long duration events in MACHO field 104 centered on (l,b) = (3.1,-3.0). We report the preliminary average optical depth of (2.0 +/- 0.4) x 10^{-6} (internal) at (l,b) = (3.9,-3.8). We discuss future work and prospects for building a coherent model of the Galaxy.

Using the exceptional long-term monitoring capabilities of the MACHO project, we present here the optical history of LMC X-2 for a continuous 6-yr period. These data were used to investigate the previously claimed periodicities for this source of 8.15 h and 12.54 d : we find upper amplitude limits of 0.10 mag and 0.09 mag, respectively.

More than 1300 variables classified provisionally as first overtone RR Lyrae pulsators in the MACHO variable star database of the Large Magellanic Cloud (LMC) have been subjected to standard frequency analysis. Based on the remnant power in the prewhitened spectra, we found 70% of the total population to be monoperiodic. The remaining 30% (411 stars) are classified as one of 9 types according to their frequency spectra. Several types of RR Lyrae pulsational behavior are clearly identified here for the first time. Together with the earlier discovered double-mode (fundamental & first overtone) variables this study increased the number of the known double-mode stars in the LMC to 181. During the total 6.5yr time span of the data, 10% of the stars show strong period changes. We also discovered two additional types of multifrequency pulsators with low occurrence rates of 2% for each. In the first type there remains one closely spaced component after prewhitening by the main pulsation frequency. In the second type the number of remnant components is two, they are also closely spaced, and, in addition, they are symmetric in their frequency spacing relative to the central component. This latter type of variables is associated with their relatives among the fundamental pulsators, known as Blazhko variables. Their high frequency (~20%) among the fundamental mode variables versus the low occurrence rate of their first overtone counterparts makes it more difficult to explain Blazhko phenomenon by any theory depending mainly on the role of aspect angle or magnetic field. Current theoretical models invoke nonradial pulsation components in these stars.

We present the microlensing optical depth towards the Galactic bulge based on the detection of 99 events found in our Difference Image Analysis (DIA) survey. This analysis encompasses three years of data, covering ~ 17 million stars in ~ 4 deg^2, to a source star baseline magnitude limit of V = 23. The DIA technique improves the quality of photometry in crowded fields, and allows us to detect more microlensing events with faint source stars. We find this method increases the number of detection events by 85% compared with the standard analysis technique. DIA light curves of the events are presented and the microlensing fit parameters are given. The total microlensing optical depth is estimated to be tau_(total)= 2.43^(+0.39/-0.38) x 10^(-6) averaged over 8 fields centered at l=2.68 and b=-3.35. For the bulge component we find tau_(bulge)=3.23^(+0.52/-0.50) x 10^(-6) assuming a 25% stellar contribution from disk sources. These optical depths are in good agreement with the past determinations of the MACHO Alcock et al. (1997) and OGLE Udalski et al. (1994) groups, and are higher than predicted by contemporary Galactic models. We show that our observed event timescale distribution is consistent with the distribution expected from normal mass stars, if we adopt the stellar mass function of Scalo (1986) as our lens mass function. However, we note that as there is still disagreement about the exact form of the stellar mass function, there is uncertainty in this conclusion. Based on our event timescale distribution we find no evidence for the existence of a large population of brown dwarfs in the direction of the Galactic bulge.

The MACHO project is a search for dark matter in the form of massive compact halo objects (MACHOs). The project has photometrically monitored tens of millions of stars in the Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and Galactic bulge in search of rare gravitational microlensing events caused by these otherwise invisible objects. In 5.7 years of observations toward the LMC some 13-17 microlensing events have been observed by the MACHO survey, allowing powerful statements to be made about the nature of the dark population in the halo of our Galaxy. A critical component of these statements is an accurate determination of the survey's detection efficiency. The detection efficiency is a complicated function of temporal sampling, stellar crowding (the luminosity function), image quality, photometry, time-series analysis, and criteria used to select the microlensing candidates. Such a complex interdependence is most naturally solved using a Monte Carlo approach. Here we describe the details of the Monte Carlo used to calculate the efficiency presented in the MACHO 5.7-year LMC results. Here we correct several shortcomings of past determinations, including (1) adding fainter source stars (2.5 magnitudes below our faintest detected "stars"), (2) an up-to-date luminosity function for the LMC, (3) better sampling of real images in both stellar density and observing conditions, (4) an improved scheme for adding artificial microlensing onto a random sample of real lightcurves, and many other improvements. [Abridged]

We present a 9 million star color-magnitude diagram (9M CMD) of the LMC bar. The 9M CMD reveals a complex superposition of different age and metallicity stellar populations, with important stellar evolutionary phases occurring over 3 orders of magnitude in number density. First, we count the non-variable supergiants, the associated Cepheids, and measure the effective temperatures defining the instability strip. Lifetime predictions of stellar evolution theory are tested, with implications for the origin of low-luminosity Cepheids. The highly-evolved AGB stars have a bimodal distribution in brightness, which we interpret as discrete old populations (>1 Gyr). The faint AGB may be metal-poor and very old. We identify the clusters NGC 411 and M3 as templates for the admixture of old stellar populations. However, there are indications that the old and metal-poor field population has a red HB morphology: the RR Lyraes lie on the red edge of the instability strip, the AGB-bump is very red, and the ratio of AGB-bump stars to RR Lyraes is quite large. If the HB second parameter is age, the old and metal-poor field population likely formed after the oldest clusters. Lifetime predictions of stellar evolution theory lead us to associate a significant fraction of the red HB clump giants with the same old and metal-poor population producing the RR Lyraes and the AGB-bump. In this case, compared to the age-dependent luminosity predictions of stellar evolution theory, the red HB clump is too bright relative to the RR Lyraes and AGB-bump. Last, the surface density profile of RR Lyraes is fit by an exponential, favoring a disk-like rather than spheroidal distribution. We conclude that the age of the LMC disk is probably similar to the age of the Galactic disk. (ABRIDGED)

We report on our search for microlensing towards the Large Magellanic Cloud (LMC). Analysis of 5.7 years of photometry on 11.9 million stars in the LMC reveals 13 - 17 microlensing events. This is significantly more than the $\sim$ 2 to 4 events expected from lensing by known stellar populations. The timescales ($\that$) of the events range from 34 to 230 days. We estimate the microlensing optical depth towards the LMC from events with $2 < \that < 400$ days to be 1.2 ^{+0.4}_ {-0.3} \ten{-7}$, with an additional 20% to 30% of systematic error. The spatial distribution of events is mildly inconsistent with LMC/LMC disk self-lensing, but is consistent with an extended lens distribution such as a Milky Way or LMC halo. Interpreted in the context of a Galactic dark matter halo, consisting partially of compact objects, a maximum likelihood analysis gives a MACHO halo fraction of 20% for a typical halo model with a 95% confidence interval of 8% to 50%. A 100% MACHO halo is ruled out at the 95% C.L. for all except our most extreme halo model. Interpreted as a Galactic halo population, the most likely MACHO mass is between $ 0.15 \msun$ and $ 0.9 \msun$, depending on the halo model, and the total mass in MACHOs out to 50 kpc is found to be 9+4-3 10^{10} msun, independent of the halo model. These results are marginally consistent with our previous results, but are lower by about a factor of two. Besides a larger data set, this work also includes an improved efficiency determination, improved likelihood analysis, and more thorough testing of systematic errors, especially with respect to the treatment of potential backgrounds to microlensing, such as supernovae in galaxies behind the LMC. [Abridged]

We present HST photometry for three fields in the outer disk of the LMC extending approximately four magnitudes below the faintest main sequence turnoff. We cannot detect any strongly significant differences in the stellar populations of the three fields based on the morphologies of the color-magnitude diagrams, the luminosity functions, and the relative numbers of stars in different evolutionary stages. Our observations therefore suggest similar star formation histories in these regions, although some variations are certainly allowed. The fields are located in two regions of the LMC: one is in the north-east field and two are located in the north-west. Under the assumption of a common star formation history, we combine the three fields with ground-based data at the same location as one of the fields to improve statistics for the brightest stars. We compare this stellar population with those predicted from several simple star formation histories suggested in the literature, using a combination of the R-method of Bertelli et al (1992) and comparisons with the observed luminosity function. The only model which we consider that is not rejected by the observations is one in which the star formation rate is roughly constant for most of the LMC's history and then increases by a factor of three about 2 Gyr ago. Such a model has roughly equal numbers of stars older and younger than 4 Gyr, and thus is not dominated by young stars. This star formation history, combined with a closed box chemical evolution model, is consistent with observations that the metallicity of the LMC has doubled in the past 2 Gyr.