Redshifts used in current cosmological supernova samples are measured using two primary techniques, one based on well-measured host galaxy spectral lines and the other based on supernova-dominated spectra. Here, we construct an updated Pantheon catalog with revised redshifts, redshift sources and estimated uncertainties for the entire sample to investigate whether these two techniques yield consistent results. The best-fit cosmological parameters using these two measurement techniques disagree, and we explore several possible sources of bias which could result from using the lower-precision supernova-dominated redshifts. In a pilot study, we show that using a host redshift-only subsample of the Pantheon catalog will result in lower $\Omega_m$ and matter density $\Omega_m h^2$ and slightly higher $H_0$ than previous analysis which, among other possibilities, could result in supernova and CMB measurements agreeing on $\Omega_m h^2$ despite tension in $H_0$. To obtain rigorous results, though, the Pantheon catalog should be improved by obtaining host spectra for supernova that have faded and future surveys should be designed to use host galaxy redshifts rather than lower-precision methods.

A major goal in the field of galaxy formation is to understand the formation of the Milky Way's disk. The first step toward doing this is to empirically describe its present state. We use the new high-dimensional dataset of 19 abundances from 27,135 red clump APOGEE stars to examine the distribution of clusters defined using abundances. We explore different dimensionality reduction techniques and implement a non-parametric agglomerate hierarchical clustering method. We see that groups defined using abundances are spatially separated, as a function of age. Furthermore, the abundance groups represent different distributions in the [Fe/H]-age plane. Ordering our clusters by age reveals patterns suggestive of the sequence of chemical enrichment in the disk over time. Our results indicate that a promising avenue to trace the details of the disk's assembly is via a full interpretation of the empirical connections we report.

Tidal debris structures formed from disrupted satellites contain important clues about the assembly histories of galaxies. To date, studies of these structures have been hampered by reliance on by-eye identification and morphological classification which leaves their interpretation significantly uncertain. In this work we present a new machine-vision technique based on the Subspace-Constrained Mean Shift (SCMS) algorithm which can perform these tasks automatically. SCMS finds the location of the high-density `ridges' that define substructure morphology. After identification, the coefficients of an orthogonal series density estimator are used to classify points on the ridges as part of a continuum between shell-like or stream-like debris, from which a global morphological classification can be determined. We dub this procedure Subspace Constrained Unsupervised Detection of Structure (SCUDS). By applying this tool to controlled N--body simulations of minor mergers we demonstrate that the extracted classifications correspond to the well--understood underlying physics of phase mixing. The application of SCUDS to resolved stellar population data from near-future surveys will inform our understanding of the buildup of galaxies stellar halos.

Observational studies of halo stars during the last two decades have placed some limits on the quantity and nature of accreted dwarf galaxy contributions to the Milky Way stellar halo by typically utilizing stellar phase-space information to identify the most recent halo accretion events. In this study we tested the prospects of using 2-D chemical abundance ratio distributions (CARDs) found in stars of the stellar halo to determine its formation history. First, we used simulated data from eleven "MW-like" halos to generate satellite template sets of 2-D CARDs of accreted dwarf satellites which are comprised of accreted dwarfs from various mass regimes and epochs of accretion. Next, we randomly drew samples of $\sim10^{3-4}$ mock observations of stellar chemical abundance ratios ([$\alpha$/Fe], [Fe/H]) from those eleven halos to generate samples of the underlying densities for our CARDs to be compared to our templates in our analysis. Finally, we used the expectation-maximization algorithm to derive accretion histories in relation to the satellite template set (STS) used and the sample size. For certain STS used we typically can identify the relative mass contributions of all accreted satellites to within a factor of 2. We also find that this method is particularly sensitive to older accretion events involving low-luminous dwarfs e.g. ultra-faint dwarfs - precisely those events that are too ancient to be seen by phase-space studies of stars and too faint to be seen by high-z studies of the early Universe. Since our results only exploit two chemical dimensions and near-future surveys promise to provide $\sim6-9$ dimensions, we conclude that these new high-resolution spectroscopic surveys of the stellar halo will allow us to recover its accretion history - and the luminosity function of infalling dwarf galaxies - across cosmic time.

One way to constrain the nature of the high-redshift progenitors of the Milky Way is to look at the low-metallicity stellar populations of the different Galactic components today. For example, high-resolution spectroscopy of very metal poor (VMP) stars demonstrates remarkable agreement between the distribution of [Ti/Fe] in the stellar populations of the Milky Way halo (MW) and ultra-faint dwarf (UFD) galaxies. In contrast, for the neutron capture (nc) abundance ratio distributions [(Sr,Ba)/Fe], the peak of the small UFD sample (6 stars) exhibits a signicant under-abundance relative to the VMP stars in the larger MW halo sample (~ 300 stars). We present a simple scenario that can simultaneously explain these similarities and differences by assuming: (i) that the MW VMP stars were predominately enriched by a prior generation of stars which possessed a higher total mass than the prior generation of stars that enriched the UFD VMP stars; and (ii) a much stronger mass-dependent yield (MDY) for nc-elements than for the (known) MDY for Ti. Simple statistical tests demonstrate that conditions (i) and (ii) are consistent with the observed abundance distributions, albeit without strong constraints on model parameters. A comparison of the broad constraints for these nc-MDY with those derived in the literature seems to rule out Ba production from low-mass SNs and affirms models that primarily generate yields from high-mass SN. Our scenario can be confirmed by a relatively modest (factor of ~ 3-4) increase in the number of high-resolution spectra of VMP stars in UFDs.

Extensive photometric studies of the globular clusters located towards the center of the Milky Way have been historically neglected. The presence of patchy differential reddening in front of these clusters has proven to be a significant obstacle to their detailed study. We present here a well-defined and reasonably homogeneous photometric database for 25 of the brightest Galactic globular clusters located in the direction of the inner Galaxy. These data were obtained in the B, V, and I bands using the Magellan 6.5m telescope and the Hubble Space Telescope. A new technique is extensively used in this paper to map the differential reddening in the individual cluster fields, and to produce cleaner, dereddened color-magnitude diagrams for all the clusters in the database. Subsequent papers will detail the astrophysical analysis of the cluster populations, and the properties of the obscuring material along the clusters' lines of sight.

A serious limitation in the study of many globular clusters -- especially those located near the Galactic Center -- has been the existence of large and differential extinction by foreground dust. In a series of papers we intend to map the differential extinction and remove its effects, using a new dereddening technique, in a sample of clusters in the direction of the inner Galaxy, observed using the Magellan 6.5m telescope and the Hubble Space Telescope. These observations and their analysis will let us produce high quality color-magnitude diagrams of these poorly studied clusters that will allow us to determine these clusters' relative ages, distances and chemistry and to address important questions about the formation and the evolution of the inner Galaxy. We also intend to use the maps of the differential extinction to sample and characterize the interstellar medium along the numerous low latitude lines of sight where the clusters in our sample lie. In this first paper we describe in detail our dereddening method along with the powerful statistics tools that allow us to apply it, and we show the kind of results that we can expect, applying the method to M62, one of the clusters in our sample. The width of the main sequence and lower red giant branch narrows by a factor of 2 after applying our dereddening technique, which will significantly help to constrain the age, distance, and metallicity of the cluster.

(abridged) We develop an algorithm for estimating parameters of a distribution sampled with contamination, employing a statistical technique known as ``expectation maximization'' (EM). Given models for both member and contaminant populations, the EM algorithm iteratively evaluates the membership probability of each discrete data point, then uses those probabilities to update parameter estimates for member and contaminant distributions. The EM approach has wide applicability to the analysis of astronomical data. Here we tailor an EM algorithm to operate on spectroscopic samples obtained with the Michigan-MIKE Fiber System (MMFS) as part of our Magellan survey of stellar radial velocities in nearby dwarf spheroidal (dSph) galaxies. These samples are presented in a companion paper and contain discrete measurements of line-of-sight velocity, projected position, and Mg index for ~1000 - 2500 stars per dSph, including some fraction of contamination by foreground Milky Way stars. The EM algorithm quantifies both dSph and contaminant distributions, returning maximum-likelihood estimates of the means and variances, as well as the probability that each star is a dSph member. Applied to our MMFS data, the EM algorithm identifies more than 5000 probable dSph members. We test the performance of the EM algorithm on simulated data sets that represent a range of sample size, level of contamination, and amount of overlap between dSph and contaminant velocity distributions. The simulations establish that for samples ranging from large (N ~3000) to small (N~30), the EM algorithm distinguishes members from contaminants and returns accurate parameter estimates much more reliably than conventional methods of contaminant removal (e.g., sigma clipping).

We present stellar velocity dispersion profiles for seven Milky Way dwarf spheroidal (dSph) satellite galaxies. We have measured 8394 line-of-sight velocities (+/- 2.5 km/s) for 6804 stars from high-resolution spectra obtained at the Magellan and MMT telescopes. We combine these new data with previously published velocities to obtain the largest available kinematic samples, which include more than 5500 dSph members. All the measured dSphs have stellar velocity dispersion of order 10 km/s that remains approximately constant with distance from the dSph center, out to and in some cases beyond the radius at which the mean surface brightness falls to the background level. Assuming dSphs reside within dark matter halos characterized by the NFW density profile, we obtain reasonable fits to the empirical velocity dispersion profiles. These fits imply that, among the seven dSphs, M_vir ~ 10^[8-9] M_sun. The mass enclosed at a radius of 600 pc, the region common to all data sets, ranges from (2-7) x 10^7 M_sun .

We introduce a stellar velocity survey of dwarf spheroidal galaxies, undertaken using the Michigan/MIKE Fiber System (MMFS) at the Magellan/Clay 6.5 m telescope at Las Campanas Observatory. As of 2006 November we have used MMFS to collect 6415 high-resolution (R= 20000-25000) spectra from 5180 stars in four dwarf spheroidal galaxies: Carina, Fornax, Sculptor and Sextans. Spectra sample the range 5140-5180 Angstroms, which includes the prominent magnesium triplet absorption feature. We measure radial velocity (RV) to a median precision of 2.0 km/s for stars as faint as V ~ 20.5. From the spectra we also are able to measure the strength of iron and magnesium absorption features using spectral indices that correlate with effective temperature, surface gravity and chemical abundance. Measurement of line strength allows us to identify interloping foreground stars independently of velocity, and to examine the metallicity distribution among dSph members. Here we present detailed descriptions of MMFS, our target selection and spectroscopic observations, the data reduction procedure, and error analysis. We compare our RV results to previously published measurements for individual stars. In some cases we find evidence for a mild, velocity-dependent offset between the RVs we measure using the magnesium triplet and previously published RV measurements derived from the infrared calcium triplet. In companion papers we will present the complete data sets and kinematic analyses of these new observations.

We present multifiber echelle radial velocity results for 551 stars in the Sextans dwarf spheroidal galaxy and identify 294 stars as probable Sextans members. The projected velocity dispersion profile of the binned data remains flat to a maximum angular radius of $30^{\prime}$. We introduce a nonparametric technique for estimating the projected velocity dispersion surface, and use this to search for kinematic substructure. Our data do not confirm previous reports of a kinematically distinct stellar population at the Sextans center. Instead we detect a region near the Sextans core radius that is kinematically colder than the overall Sextans sample with 95% confidence.