The Milky Way provides an ideal laboratory to test our understanding of galaxy evolution, owing to our ability to observe our Galaxy over fine scales. However, connecting the Galaxy to the wider galaxy population remains difficult, due to the challenges posed by our internal perspective and to the different observational techniques employed. Here, we present a sample of galaxies identified as Milky Way Analogs (MWAs) on the basis of their stellar masses and bulge-to-total ratios, observed as part of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. We analyse the galaxies in terms of their stellar kinematics and populations as well as their ionised gas contents. We find our sample to contain generally young stellar populations in their outskirts. However, we find a wide range of stellar ages in their central regions, and we detect central AGN-like or composite-like activity in roughly half of the sample galaxies, with the other half consisting of galaxies with central star-forming emission or emission consistent with old stars. We measure gradients in gas metallicity and stellar metallicity that are generally flatter in physical units than those measured for the Milky Way; however, we find far better agreement with the Milky Way when scaling gradients by galaxies' disc scale lengths. From this, we argue much of the discrepancy in metallicity gradients to be due to the relative compactness of the Milky Way, with differences in observing perspective also likely to be a factor.

On small scales there have been a number of claims of discrepancies between the standard Cold Dark Matter (CDM) model and observations. The 'missing satellites problem' infamously describes the overabundance of subhalos from CDM simulations compared to the number of satellites observed in the Milky Way. A variety of solutions to this discrepancy have been proposed; however, the impact of the specific properties of the Milky Way halo relative to the typical halo of its mass have yet to be explored. Motivated by recent studies that identified ways in which the Milky Way is atypical (e.g., Licquia et al. 2015), we investigate how the properties of dark matter halos with mass comparable to our Galaxy's --- including concentration, spin, shape, and scale factor of the last major merger --- correlate with the subhalo abundance. Using zoom-in simulations of Milky Way-like halos, we build two models of subhalo abundance as functions of host halo properties and conclude that the Milky Way should be expected to have 22%-44% fewer subhalos with low maximum rotation velocities ($V_{\rm max}^{\rm sat} \sim 10$kms$^{-1}$) at the 95% confidence level and up to 72% fewer than average subhalos with high rotation velocities ($V_{\rm max}^{\rm sat} \gtrsim 30$kms$^{-1}$, comparable to the Magellanic Clouds) than would be expected for a typical halo of the Milky Way's mass. Concentration is the most informative single parameter for predicting subhalo abundance. Our results imply that models tuned to explain the missing satellites problem assuming typical subhalo abundances for our Galaxy will be over-correcting.

It is crucial to understand how the Milky Way, the galaxy we can study in the most intimate detail, fits in amongst other galaxies. Key examples include the Tully-Fisher relation (TFR) --- i.e., the tight correlation between luminosity ($L$) and rotational velocity ($V_\textrm{rot}$) --- and the 3-dimensional luminosity-velocity-radius ($LVR$) scaling relation. Several past studies have characterized the MW as a 1--1.5$\sigma$ outlier to the TFR. This study reexamines such comparisons using new estimates of MW properties that are robust to many of the systematic uncertainties that have been a problem in the past and are based on assumptions consistent with those used for other spiral galaxies. Comparing to scaling relations derived from modern extragalactic data, we find that our Galaxy's properties are in excellent agreement with TFRs defined using any SDSS-filter absolute magnitude, stellar mass, or baryonic mass as the $L$ proxy. We next utilize disk scale length ($R_\textrm{d}$) measurements to extend this investigation to the $LVR$ relation. Here we find that our Galaxy lies farther from the relation than $\sim$90\% of other spiral galaxies, yielding $\sim$9.5$\sigma$ evidence that it is unusually compact for its $L$ and $V_\textrm{rot}$ (based on MW errors alone), a result that holds for all of the $L$ proxies considered. The expected $R_\textrm{d}$ for the MW from the $LVR$ relation is $\sim$5 kpc, nearly twice as large as the observed value, with error estimates placing the two in tension at the $\sim$1.4$\sigma$ level. The compact scale length of the Galactic disk could be related to other ways in which the MW has been found to be anomalous.

The exponential scale length ($L_d$) of the Milky Way's (MW's) disk is a critical parameter for describing the global physical size of our Galaxy, important both for interpreting other Galactic measurements and helping us to understand how our Galaxy fits into extragalactic contexts. Unfortunately, current estimates span a wide range of values and often are statistically incompatible with one another. Here, we perform a Bayesian meta-analysis to determine an improved, aggregate estimate for $L_d$, utilizing a mixture-model approach to account for the possibility that any one measurement has not properly accounted for all statistical or systematic errors. Within this machinery we explore a variety of ways of modeling the nature of problematic measurements, and then employ a Bayesian model averaging technique to derive net posterior distributions that incorporate any model-selection uncertainty. Our meta-analysis combines 29 different (15 visible and 14 infrared) photometric measurements of $L_d$ available in the literature; these involve a broad assortment of observational datasets, MW models and assumptions, and methodologies, all tabulated herein. Analyzing the visible and infrared measurements separately yields estimates for $L_d$ of $2.71^{+0.22}_{-0.20}$ kpc and $2.51^{+0.15}_{-0.13}$ kpc, respectively, whereas considering them all combined yields $2.64\pm0.13$ kpc. The ratio between the visible and infrared scale lengths determined here is very similar to that measured in external spiral galaxies. We use these results to update the model of the Galactic disk from our previous work, constraining its stellar mass to be $4.8^{+1.5}_{-1.1}\times10^{10}$ $\textrm{M}_\odot$, and the MW's total stellar mass to be $5.7^{+1.5}_{-1.1}\times10^{10}$ $\textrm{M}_\odot$.

We present improved estimates of several global properties of the Milky Way, including its current star formation rate (SFR), the stellar mass contained in its disk and bulge+bar components, as well as its total stellar mass. We do so by combining previous measurements from the literature using a hierarchical Bayesian (HB) statistical method that allows us to account for the possibility that any value may be incorrect or have underestimated errors. We show that this method is robust to a wide variety of assumptions about the nature of problems in individual measurements or error estimates. Ultimately, our analysis yields a SFR for the Galaxy of $\dot{\mathrm{M}}_\star=1.65\pm0.19$ $\textrm{M}_\odot \textrm{yr}^{-1}$, assuming a Kroupa initial mass function (IMF). By combining HB methods with Monte Carlo simulations that incorporate the latest estimates of the Galactocentric radius of the Sun, $R_0$, the exponential scale length of the disk, $L_d$, and the local surface density of stellar mass, $\Sigma_\star(R_0)$, we show that the mass of the Galactic bulge+bar is $\textrm{M}_\star^B=0.91\pm0.07\times10^{10}$ $\textrm{M}_\odot$, the disk mass is $\textrm{M}_\star^D=5.17\pm1.11\times10^{10}$ $\textrm{M}_\odot$, and their combination yields a total stellar mass of $\textrm{M}_\star=6.08\pm1.14\times10^{10}$ $\textrm{M}_\odot$ (assuming a Kroupa IMF and an exponential disk profile). This analysis is based upon a new compilation of literature bulge mass estimates, normalized to common assumptions about the stellar initial mass function and Galactic disk properties, presented herein. We additionally find a bulge-to-total mass ratio for the Milky Way of $B/T=0.150^{+0.028}_{-0.019}$ and a specific star formation rate of $\dot{\mathrm{M}}_\star/\textrm{M}_\star=2.71\pm0.59\times10^{-11} \textrm{yr}^{-1}$.

We demonstrate a new statistical method of determining the global photometric properties of the Milky Way (MW) to an unprecedented degree of accuracy, allowing our Galaxy to be compared directly to objects measured in extragalactic surveys. Capitalizing on the high-quality imaging and spectroscopy dataset from the Sloan Digital Sky Survey (SDSS), we exploit the inherent dependence of galaxies' luminosities and colors on their total stellar mass, $\mathrm{M}_\star$, and star formation rate (SFR), $\mathrm{\dot{M}}_\star$, by selecting a sample of $Milky$ $Way$ $analog$ $galaxies$ designed to reproduce the best Galactic $\mathrm{M}_\star$ and $\mathrm{\dot{M}}_\star$ measurements, including all measurement uncertainties. Making the Copernican assumption that the MW is not extraordinary amongst galaxies of similar stellar mass and SFR, we then analyze the photometric properties of this matched sample, constraining the characteristics of our Galaxy without suffering interference from interstellar dust. We explore a variety of potential systematic errors that could affect this method, and find that they are subdominant to random uncertainties. We present both SDSS $ugriz$ absolute magnitudes and colors in both rest-frame $z$=0 and $z$=0.1 passbands for the MW, which are in agreement with previous estimates but can have up to $\sim$3$\times$ lower errors. We find the MW to have absolute magnitude $^0\!M_r-5\log h=-21.00_{-0.37}^{+0.38}$ and integrated color $^0(g-r)=0.682_{-0.056}^{+0.066}$, indicating that it may belong to the green-valley region in color-magnitude space and ranking it amongst the brightest and reddest of spiral galaxies. We also present new estimates of global stellar mass-to-light ratios for our Galaxy. This work will help relate our in-depth understanding of the Galaxy to studies of more distant objects.

We describe the algorithm used to select the Luminous Red Galaxy (LRG) sample for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV) using photometric data from both the SDSS and the Wide-Field Infrared Survey Explorer (WISE). LRG targets are required to meet a set of color selection criteria and have z-band and i-band MODEL magnitudes z < 19.95 and 19.9 < i < 21.8, respectively. Our algorithm selects roughly 50 LRG targets per square degree, the great majority of which lie in the redshift range 0.6 < z < 1.0 (median redshift 0.71). We demonstrate that our methods are highly effective at eliminating stellar contamination and lower-redshift galaxies. We perform a number of tests using spectroscopic data from SDSS-III/BOSS to determine the redshift reliability of our target selection and its ability to meet the science requirements of eBOSS. The SDSS spectra are of high enough signal-to-noise ratio that at least 89% of the target sample yields secure redshift measurements. We also present tests of the uniformity and homogeneity of the sample, demonstrating that it should be clean enough for studies of the large-scale structure of the universe at higher redshifts than SDSS-III/BOSS LRGs reached.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We describe a new method of combining optical and infrared photometry to select Luminous Red Galaxies (LRGs) at redshifts $z > 0.6$. We explore this technique using a combination of optical photometry from CFHTLS and HST, infrared photometry from the WISE satellite, and spectroscopic or photometric redshifts from the DEEP2 Galaxy Redshift Survey or COSMOS. We present a variety of methods for testing the success of our selection, and present methods for optimization given a set of rest-frame color and redshift requirements. We have tested this selection in two different regions of the sky, the COSMOS and Extended Groth Strip (EGS) fields, to reduce the effect of cosmic/sample variance. We have used these methods to assemble large samples of LRGs for two different ancillary programs as a part of the SDSS-III/ BOSS spectroscopic survey. This technique is now being used to select $\sim$600,000 LRG targets for SDSS-IV/eBOSS, which began observations in Fall 2014, and will be adapted for the proposed DESI survey. We have found these methods can select high-redshift LRGs efficiently with minimal stellar contamination; this is extremely difficult to achieve with selections that rely on optical photometry alone.