The spectral analysis and data products in Data Release 16 (DR16; December 2019) from the high-resolution near-infrared APOGEE-2/SDSS-IV survey are described. Compared to the previous APOGEE data release (DR14; July 2017), APOGEE DR16 includes about 200000 new stellar spectra, of which 100000 are from a new southern APOGEE instrument mounted on the 2.5 m du Pont telescope at Las Campanas Observatory in Chile. DR16 includes all data taken up to August 2018, including data released in previous data releases. All of the data have been re-reduced and re-analyzed using the latest pipelines, resulting in a total of 473307 spectra of 437445 stars. Changes to the analysis methods for this release include, but are not limited to, the use of MARCS model atmospheres for calculation of the entire main grid of synthetic spectra used in the analysis, a new method for filling "holes" in the grids due to unconverged model atmospheres, and a new scheme for continuum normalization. Abundances of the neutron capture element Ce are included for the first time. A new scheme for estimating uncertainties of the derived quantities using stars with multiple observations has been applied, and calibrated values of surface gravities for dwarf stars are now supplied. Compared to DR14, the radial velocities derived for this release more closely match those in the Gaia DR2 data base, and a clear improvement in the spectral analysis of the coolest giants can be seen. The reduced spectra as well as the result of the analysis can be downloaded using links provided in the SDSS DR16 web page.

Chemically tagging groups of stars born in the same birth cluster is a major goal of spectroscopic surveys. To investigate the feasibility of such strong chemical tagging, we perform a blind chemical tagging experiment on abundances measured from APOGEE survey spectra. We apply a density-based clustering algorithm to the eight dimensional chemical space defined by [Mg/Fe], [Al/Fe], [Si/Fe], [K/Fe], [Ti/Fe], [Mn/Fe], [Fe/H], and [Ni/Fe], abundances ratios which together span multiple nucleosynthetic channels. In a high quality sample of 182,538 giant stars, we detect twenty-one candidate clusters with more than fifteen members. Our candidate clusters are more chemically homogeneous than a population of non-member stars with similar [Mg/Fe] and [Fe/H], even in abundances not used for tagging. Group members are consistent with having the same age and fall along a single stellar-population track in logg vs. Teff space. Each group's members are distributed over multiple kpc, and the spread in their radial and azimuthal actions increases with age. We qualitatively reproduce this increase using N-body simulations of cluster dissolution in Galactic potentials that include transient winding spiral arms. Observing our candidate birth clusters with high-resolution spectroscopy in other wavebands to investigate their chemical homogeneity in other nucleosynthetic groups will be essential to confirming the efficacy of strong chemical tagging. Our initially spatially-compact but now widely dispersed candidate clusters will provide novel limits on chemical evolution and orbital diffusion in the Galactic disc, and constraints on star formation in loosely-bound groups.

We compare the results of high-resolution simulations of individual dark matter subhalos evolving in external tidal fields with and without baryonic bulge and disk components, where the average dark matter particle mass is three orders of magnitude smaller than cosmological zoom-in simulations of galaxy formation. The Via Lactea II simulation is used to setup our initial conditions and provides a basis for our simulations of subhalos in a dark matter-only tidal field, while an observationally motivated model for the Milky Way is used for the tidal field that is comprised of a dark matter halo, a stellar disk, and a stellar bulge. Our simulations indicate that including stellar components in the tidal field results in the number of subhalos in Milky Way-like galaxies being only $65\%$ of what is predicted by $\Lambda$ Cold Dark Matter ($\Lambda$CDM). For subhalos with small pericentres $(r_p \lesssim 25$ kpc), the subhalo abundance is reduced further to $40\%$, with the surviving subhalos being less dense and having a tangentially-anisotropic orbital distribution. Conversely, subhalos with larger pericentres are minimally affected by the inclusion of a stellar component in the tidal field, with the total number of outer subhalos $\approx 75\%$ of the $\Lambda$CDM prediction. The densities of outer subhalos are comparable to predictions from $\Lambda$CDM, with the subhalos having an isotropic distribution of orbits. These ratios are higher than those found in previous studies that include the effects baryonic matter, which are affected by spurious disruption caused by low resolution.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We make use of APOGEE and $Gaia$ data to identify stars that are consistent with being born in the same association or star cluster as the Sun. We limit our analysis to stars that match solar abundances within their uncertainties, as they could have formed from the same Giant Molecular Cloud (GMC) as the Sun. We constrain the range of orbital actions that solar siblings can have with a suite of simulations of solar birth clusters evolved in static and time-dependent tidal fields. The static components of each galaxy model are the bulge, disk, and halo, while the various time-dependent components include a bar, spiral arms, and GMCs. In galaxy models without GMCs, simulated solar siblings all have $J_R < 122$ km $\rm s^{-1}$ kpc, $990 < L_z < 1986$ km $\rm s^{-1}$ kpc, and $0.15 < J_z < 0.58$ km $\rm s^{-1}$ kpc. Given the actions of stars in APOGEE and $Gaia$, we find 104 stars that fall within this range. One candidate in particular, Solar Sibling 1, has both chemistry and actions similar enough to the solar values that strong interactions with the bar or spiral arms are not required for it to be dynamically associated with the Sun. Adding GMCs to the potential can eject solar siblings out of the plane of the disk and increase their $J_z$, resulting in a final candidate list of 296 stars. The entire suite of simulations indicate that solar siblings should have $J_R < 122$ km $\rm s^{-1}$ kpc, $353 < L_z < 2110$ km $\rm s^{-1}$ kpc, and $J_z < 0.8$ km $\rm s^{-1}$ kpc. Given these criteria, it is most likely that the association or cluster that the Sun was born in has reached dissolution and is not the commonly cited open cluster M67.

The systematic offset of Gaia parallaxes has been widely reported with Gaia's second data release, and it is expected to persist in future Gaia data. In order to use Gaia parallaxes to infer distances to high precision, we develop a hierarchical probabilistic model to determine the Gaia parallax zero point offset along with the calibration of an empirical model for luminosity of red clump stars by combining astrometric and photometric measurements. Using a cross-matched sample of red clump stars from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and Gaia Data Release 2 (DR2), we report the parallax zero point offset in DR2 to be $ \varpi_0 = -48 \pm 1~\mu\text{as} $. We infer the red clump absolute magnitude to be $ -1.622 \pm 0.004 $ in $ K_s $, $ 0.517 \pm 0.004 $ in $ G $, $ -1.019 \pm 0.004 $ in $ J $, and $ -1.516 \pm 0.004 $ in $ H $. The intrinsic scatter of the red clump is $ \sim 0.09 $ mag in $ J $, $ H $ and $ K_s $, or $ \sim 0.17 $ mag in $ G $. We tailor our models to accommodate more complex analyses such as investigating the variations of the parallax zero point with each source's observed magnitude, observed colour, and sky position. In particular, we find fluctuations of the zero point across the sky to be of order or less than a few 10s of $ \mu\text{as} $.

We summarise recent results from analysis of APOGEE/Gaia data for stellar populations in the Galactic halo, disk, and bulge, leading to constraints on the contribution of dwarf galaxies and globular clusters to the stellar content of the Milky Way halo. Interpretation of the extant data in light of cosmological numerical simulations suggests that the Milky Way has been subject to an unusually intense accretion history at z >~ 1.5.

The stellar mass in the halo of the Milky Way is notoriously difficult to determine, owing to the paucity of its stars in the solar neighbourhood. With tentative evidence from $\it{Gaia}$ that the nearby stellar halo is dominated by a massive accretion event -- referred to as $\it{Gaia-Enceladus}$ or Sausage -- these constraints are now increasingly urgent. We measure the mass in kinematically selected mono-abundance populations (MAPs) of the stellar halo between $-3 <$ [Fe/H] $< -1$ and $0.0 <$ [Mg/Fe] $< 0.4$ using red giant star counts from APOGEE DR14. We find that MAPs are well fit by single power laws on triaxial ellipsoidal surfaces, and we show that that the power law slope $\alpha$ changes such that high $\mathrm{[Mg/Fe]}$ populations have $\alpha \sim 4$, whereas low [Mg/Fe] MAPs are more extended with shallow slopes, $\alpha \sim 2$. We estimate the total stellar mass to be $M_{*,\mathrm{tot}} = 1.3^{+0.3}_{-0.2}\times10^{9}\ \mathrm{M_{\odot}}$, of which we estimate $\sim 0.9^{+0.2}_{-0.1} \times 10^{9} \mathrm{M_{\odot}}$ to be accreted. We estimate that the mass of accreted stars with $e > 0.7$ is $M_{*,\mathrm{accreted},e>0.7} = 3 \pm 1\ \mathrm{(stat.)}\pm 1\ \mathrm{(syst.)}\times10^{8}\ \mathrm{M_\odot}$, or $\sim 30-50\%$ of the accreted halo mass. If the majority of these stars $\it{are}$ the progeny of a massive accreted dwarf, this places an upper limit on its stellar mass, and implies a halo mass for the progenitor of $\sim 10^{10.2\pm0.2}\ \mathrm{M_\odot}$. This constraint not only shows that the $\it{Gaia-Enceladus}$/Sausage progenitor may not be as massive as originally suggested, but that the majority of the Milky Way stellar halo was accreted. These measurements are an important step towards fully reconstructing the assembly history of the Milky Way.

New data from the $\textit{Gaia}$ satellite, when combined with accurate photometry from the Pan-STARRS survey, allow us to accurately estimate the properties of the GD-1 stream. Here, we analyze the stellar density perturbations in the GD-1 stream and show that they cannot be due to known baryonic structures like giant molecular clouds, globular clusters, or the Milky Way's bar or spiral arms. A joint analysis of the GD-1 and Pal 5 streams instead requires a population of dark substructures with masses $\approx 10^{7}$ to $10^9 \ M_{\rm{\odot}}$. We infer a total abundance of dark subhalos normalised to standard cold dark matter $n_{\rm sub}/n_{\rm sub, CDM} = 0.4 ^{+0.3}_{-0.2}$ ($68 \%$), which corresponds to a mass fraction contained in the subhalos $f_{\rm{sub}} = 0.14 ^{+0.11}_{-0.07} \%$, compatible with the predictions of hydrodynamical simulation of cold dark matter with baryons.

We analyze the distribution of stars along the GD-1 stream with a combination of data from the ${\it Gaia}$ satellite and the Pan-STARRS survey, and we show that the population of subhalos predicted by the cold dark matter paradigm are necessary and sufficient to explain the perturbations observed in the linear density of stars. This allows us to set novel constraints on alternative dark matter scenarios that predict a suppression of the subhalo mass function on scales smaller than the mass of dwarf galaxies. A combined analysis of the density perturbations in the GD-1 and Pal 5 streams leads to a $95\%$ lower limit on the mass of warm dark matter thermal relics $m_{\rm WDM}>4.6$ keV; adding dwarf satellite counts strengthens this to $m_{\rm WDM}>6.3$ keV.

Studies of the ages, abundances, and motions of individual stars in the Milky Way provide one of the best ways to study the evolution of disk galaxies over cosmic time. The formation of the Milky Way's barred inner region in particular is a crucial piece of the puzzle of disk galaxy evolution. Using data from APOGEE and Gaia, we present maps of the kinematics, elemental abundances, and age of the Milky Way bulge and disk that show the barred structure of the inner Milky Way in unprecedented detail. The kinematic maps allow a direct, purely kinematic determination of the bar's pattern speed of 41+/-3 km/s/kpc and of its shape and radial profile. We find the bar's age, metallicity, and abundance ratios to be the same as those of the oldest stars in the disk that are formed in its turbulent beginnings, while stars in the bulge outside of the bar are younger and more metal-rich. This implies that the bar likely formed ~8 Gyr ago, when the decrease in turbulence in the gas disk allowed a thin disk to form that quickly became bar-unstable. The bar's formation therefore stands as a crucial epoch in the evolution of the Milky Way, a picture that is in line with the evolutionary path that emerges from observations of the gas kinematics in external disk galaxies over the last ~10 Gyr.

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

(Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.

The application of machine learning (ML) methods to the analysis of astrophysical datasets is on the rise, particularly as the computing power and complex algorithms become more powerful and accessible. As the field of ML enjoys a continuous stream of breakthroughs, its applications demonstrate the great potential of ML, ranging from achieving tens of millions of times increase in analysis speed (e.g., modeling of gravitational lenses or analysing spectroscopic surveys) to solutions of previously unsolved problems (e.g., foreground subtraction or efficient telescope operations). The number of astronomical publications that include ML has been steadily increasing since 2010. With the advent of extremely large datasets from a new generation of surveys in the 2020s, ML methods will become an indispensable tool in astrophysics. Canada is an unambiguous world leader in the development of the field of machine learning, attracting large investments and skilled researchers to its prestigious AI Research Institutions. This provides a unique opportunity for Canada to also be a world leader in the application of machine learning in the field of astrophysics, and foster the training of a new generation of highly skilled researchers.

The era of data-intensive astronomy is being ushered in with the increasing size and complexity of observational data across wavelength and time domains, the development of algorithms to extract information from this complexity, and the computational power to apply these algorithms to the growing repositories of data. Data-intensive approaches are pushing the boundaries of nearly all fields of astronomy, from exoplanet science to cosmology, and they are becoming a critical modality for how we understand the universe. The success of these approaches range from the discovery of rare or unexpected phenomena, to characterizing processes that are now accessible with precision astrophysics and a deep statistical understanding of the datasets, to developing algorithms that maximize the science that can be extracted from any set of observations. In this white paper, we propose a number of initiatives to maximize Canada's ability to compete in this data-intensive era. We propose joining international collaborations and leveraging Canadian facilities for legacy data potential. We propose continuing to build a more agile computing infrastructure that's responsive to the needs of tackling larger and more complex data, as well as enabling quick prototyping and scaling of algorithms. We recognize that developing the fundamental skills of the field will be critical for Canadian astronomers, and discuss avenues through with the appropriate computational and statistical training could occur. Finally, we note that the transition to data-intensive techniques is not limited to astronomy, and we should coordinate with other disciplines to develop and make use of best practises in methods, infrastructure, and education.

We present the results of a detailed search for members of the Pal 5 tidal tail system in Gaia Data Release 2 (DR2). Tidal tails provide a sensitive method for measuring the current and past gravitational potential of their host galaxy as well as for testing predictions for the abundance of dark matter subhalos. The Pal 5 globular cluster and its associated tails are an excellent candidate for such analysis; however, only ~23 degrees of arc are currently known, with in particular the leading tail much shorter than the trailing. Using Gaia DR2 and its precise astrometry, we extend the known extent of the Pal 5 tail to ~30 degrees, 7 degrees of which are newly detected along the leading arm. The detected leading and trailing arms are symmetric in length and remain near constant width. This detection constrains proposed models in which the Galactic bar truncates Pal 5's leading arm. Follow-up spectroscopic observations are necessary to verify the candidate stream stars are consistent with the known tidal tails. If confirmed, this Pal 5 stream extension opens up new possibilities to constrain the Galactic potential.

We present direct $N$-body simulations of tidally filling 30,000 ${\rm M}_\odot$ star clusters orbiting between 10 kpc and 100 kpc in galaxies with a range of dark matter substructure properties. The time-dependent tidal force is determined based on the combined tidal tensor of the galaxy's smooth and clumpy dark matter components, the latter of which causes fluctuations in the tidal field that can heat clusters. The strength and duration of these fluctuations are sensitive to the local dark matter density, substructure fraction, sub-halo mass function, and the sub-halo mass-size relation. Based on the cold dark matter framework, we initially assume sub-halos are Hernquist spheres following a power-law mass function between $10^5$ and $10^{11} {\rm M}_\odot$ and find that tidal fluctuations are too weak and too short to affect star cluster evolution. Treating sub-halos as point masses, to explore how denser sub-halos affect clusters, we find that only sub-halos with masses greater than $10^{6} {\rm M}_\odot$ will cause cluster dissolution times to decrease. These interactions can also decrease the size of a cluster while increasing the velocity dispersion and tangential anisotropy in the outer regions via tidal heating. Hence increased fluctuations in the tidal tensor, especially fluctuations that are due to low-mass halos, do not necessarily translate into mass loss. We further conclude that the tidal approximation can be used to model cluster evolution in the tidal fields of cosmological simulations with a minimum cold dark matter sub-halo mass of $10^{6} {\rm M}_\odot$, as the effect of lower-mass sub-halos on star clusters is negligible.

The kinematics of the Milky Way disc as a function of age are well measured at the solar radius, but have not been studied over a wider range of Galactocentric radii. Here, we measure the kinematics of mono-age, mono-$\mathrm{[Fe/H]}$ populations in the low and high $\mathrm{[\alpha/Fe]}$ discs between $4 \lesssim R \lesssim 13$ kpc and $|z| \lesssim 2$ kpc using 65,719 stars in common between APOGEE DR14 and $\it{Gaia}$ DR2 for which we estimate ages using a Bayesian neural network model trained on asteroseismic ages. We determine the vertical and radial velocity dispersions, finding that the low and high $\mathrm{[\alpha/Fe]}$ discs display markedly different age--velocity-dispersion relations (AVRs) and shapes $\sigma_z/\sigma_R$. The high $\mathrm{[\alpha/Fe]}$ disc has roughly flat AVRs and constant $\sigma_z/\sigma_R = 0.64\pm 0.04$, whereas the low $\mathrm{[\alpha/Fe]}$ disc has large variations in this ratio which positively correlate with the mean orbital radius of the population at fixed age. The high $\mathrm{[\alpha/Fe]}$ disc component's flat AVRs and constant $\sigma_z/\sigma_R$ clearly indicates an entirely different heating history. Outer disc populations also have flatter radial AVRs than those in the inner disc, likely due to the waning effect of spiral arms. Our detailed measurements of AVRs and $\sigma_z/\sigma_R$ across the disc indicate that low $\mathrm{[\alpha/Fe]}$, inner disc ($R \lesssim 10\,\mathrm{kpc}$) stellar populations are likely dynamically heated by both giant molecular clouds and spiral arms, while the observed trends for outer disc populations require a significant contribution from another heating mechanism such as satellite perturbations. We also find that outer disc populations have slightly positive mean vertical and radial velocities, likely because they are part of the warped disc.

Chemical tagging has great promise as a technique to unveil our Galaxy's history. Grouping stars based on their similar chemistry can establish details of the star formation and merger history of the Milky Way. With precise measurements of stellar chemistry, chemical tagging may be able to group together stars born from the same gas cloud, regardless of their current positions and kinematics. Successfully tagging these birth clusters requires high quality chemical space information and a good cluster-finding algorithm. To test the feasibility of chemical tagging on data from current and upcoming spectroscopic surveys, we construct a realistic set of synthetic clusters, creating both observed spectra and derived chemical abundances for each star. We use Density-Based Spatial Clustering of Applications with Noise (DBSCAN) to group stars based on their spectra or abundances; these groups are matched to input clusters and are found to be highly homogeneous and complete. The percentage of clusters with more than 10 members recovered is 40% when tagging on abundances with uncertainties achievable with current techniques. Based on our fiducial model for the Milky Way, we predict recovering over 600 clusters with at least 10 observed members and 70% membership homogeneity in a sample similar to the APOGEE survey. Tagging larger surveys like the GALAH survey and the future Milky Way Mapper in SDSS V could recover tens of thousands of clusters at high homogeneity. Access to so many unique co-eval clusters will transform how we understand the star formation history and chemical evolution of our Galaxy.

The Ophiuchus stream is a short arc-like stellar feature of uncertain origin located $\sim 5$ kpc North of the Galactic centre. New proper motions from the second $Gaia$ data release reconcile the direction of motion of stream members with the stream arc, resolving a puzzling mismatch reported in earlier work. We use N-body simulations to show that the stream is likely only on its second pericentric passage, and thus was formed recently. The simulations suggest that the entire disrupted progenitor is visible in the observed stream today, and that little further tidal debris lies beyond the ends of the stream. The luminosity, length, width, and velocity dispersion of the stream suggest a globular cluster (GC) progenitor substantially fainter and of lower surface brightness than estimated in previous work, and unlike any other known globulars in the Galaxy. This result suggests the existence of clusters that would extend the known GC population to fainter and more weakly bound systems than hitherto known. How such a weakly-bound cluster of old stars survived until it was disrupted so recently, however, remains a mystery. Integrating backwards in time, we find that the orbits of Sagittarius and Ophiuchus passed within $\sim 5$ kpc of each other about $\sim 100$ Myrs ago, an interaction that might help resolve this puzzle.

Gaia DR2 has provided an unprecedented wealth of information about the kinematics of stars in the Solar neighbourhood, and has highlighted the degree of features in the Galactic disc. We confront the data with a range of bar and spiral models in both action-angle space, and the $R_{\mathrm{G}}-v_{\phi}$ plane. We find that the phase mixing induced by transient spiral structure creates ridges and arches in the local kinematics which are consistent with the Gaia data. We are able to produce a qualitatively good match to the data when combined with a bar with a variety of pattern speeds, and show that it is non trivial to decouple the effects of the bar and the spiral structure.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

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

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.

Gaia Data Release 2 provides a wealth of data to study the internal structure of nearby globular clusters. We use this data to investigate the internal kinematics of 11 nearby globular clusters, with a particular focus on their poorly-studied outer regions. We apply a strict set of selection criteria to remove contaminating sources and create pure cluster-member samples over a significant fraction of the radial range of each cluster. We confirm previous measurements of rotation (or a lack thereof) in the inner regions of several clusters, while extending the detection of rotation well beyond where it was previously measured and finding a steady decrease in rotation with radius. We also determine the orbital anisotropy profile and determine that clusters have isotropic cores, are radially anisotropic out to $\approx$ 4 half-light radii or $35\%$ of their limiting radii, and are then isotropic out to the limits of our datasets. We detect for the first time the presence of radial anisotropy in M 22, while confirming previous detections of radial anisotropy in 47 Tuc, M 3, M 13, M 15, and $\omega$ Cen's innermost regions. The implications of these measurements are that clusters can be separated into two categories: 1) clusters with observed radial anisotropy that likely formed tidally under-filling or are dynamically young, and 2) clusters that are primarily isotropic that likely formed tidally filling or are dynamically old.

We perform a large suite of direct N-body simulations aimed at revealing the location of the progenitor, or its remnant, of the GD-1 stream. Data from \gaia\ DR2 reveals the GD-1 stream extends over $\approx 100^\circ$, allowing us to determine the stream's leading and trailing ends. Our models suggest the length of the stream is consistent with a dynamical age of between 2-3 Gyr and the exact length, width and location of the GD-1 stream correspond to the stream's progenitor being located between $-30^\circ < \phi_{1,\mathrm{pro}} < -45^{\circ}$ in the standard GD-1 coordinate system. The model stream density profiles reveal that intact progenitors leave a strong over-density, recently-dissolved progenitors appear as gaps in the stream as escaped stars continue to move away from the remnant progenitor's location, and long-dissolved progenitors leave no observational signature on the remaining stream. Comparing our models to the GD-1 stream yields two possible scenarios for its progenitor's history: a) the progenitor reached dissolution approximately 500 Myr ago during the cluster's previous perigalactic pass and is both located at and responsible for the observed gap at $\phi_1=-40^{\circ}$ or b) the progenitor reached dissolution over 2.5 Gyr ago, the fully-dissolved remnant is at $-30^\circ < \phi_1 < -45^{\circ}$, and an observational signature of its location no longer exists. That the dissolved progenitor is in the range $-30^\circ < \phi_1 < -45^{\circ}$ implies that density fluctuations outside of this range, e.g., a deep gap at $\phi_1 \approx -20^\circ$, are likely produced by compact baryonic or dark-matter perturbers.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Chemical abundances for 15 elements (C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Ni) are presented for 83 stellar members of the 4 Gyr old solar-metallicity open cluster M67. The sample contains stars spanning a wide range of evolutionary phases, from G dwarfs to red clump stars. The abundances were derived from near-IR ($\lambda$1.5 -- 1.7$\mu$m) high-resolution spectra ($R$ = 22,500) from the SDSS-IV/APOGEE survey. A 1-D LTE abundance analysis was carried out using the APOGEE synthetic spectral libraries, via chi-square minimization of the synthetic and observed spectra with the qASPCAP code. We found significant abundance differences ($\sim$0.05 -- 0.30 dex) between the M67 member stars as a function of the stellar mass (or position on the HR diagram), where the abundance patterns exhibit a general depletion (in [X/H]) in stars at the main-sequence turnoff. The amount of the depletion is different for different elements. We find that atomic diffusion models provide, in general, good agreement with the abundance trends for most chemical species, supporting recent studies indicating that measurable atomic diffusion operates in M67 stars.

Gaia measures the five astrometric parameters for stars in the Milky Way, but only four of them (positions and proper motion, but not parallax) are well measured beyond a few kpc from the Sun. Modern spectroscopic surveys such as APOGEE cover a large area of the Milky Way disk and we can use the relation between spectra and luminosity to determine distances to stars beyond Gaia's parallax reach. Here, we design a deep neural network trained on stars in common between Gaia and APOGEE that determines spectro-photometric distances to APOGEE stars, while including a flexible model to calibrate parallax zero-point biases in Gaia DR2. We determine the zero-point offset to be $-52.3 \pm 2.0uas$ when modeling it as a global constant, but also train a multivariate zero-point offset model that depends on $G$, $G_{BP} - G_{RP}$ color, and $T_\mathrm{eff}$ and that can be applied to all 139 million stars in Gaia DR2 within APOGEE's color--magnitude range. Our spectro-photometric distances are more precise than Gaia at distances $\approx 2kpc$ from the Sun. We release a catalog of spectro-photometric distances for the entire APOGEE DR14 data set which covers Galactocentric radii $2kpc\lesssim R \lesssim19kpc$; $\approx 150,000$ stars have <10% uncertainty, making this a powerful sample to study the chemo-dynamical structure of the disk. We use this sample to map the mean [Fe/H] and 15 abundance ratios [X/Fe] from the Galactic center to the edge of the disk. Among many interesting trends, we find that the bulge and bar region at $R \lesssim 5kpc$ clearly stands out in [Fe/H] and most abundance ratios.

Gravitational encounters between small-scale dark matter substructure and cold stellar streams in the Milky Way halo lead to density perturbations in the latter, making streams an effective probe for detecting dark matter substructure. The Pal 5 stream is one such system for which we have some of the best data. However, Pal 5 orbits close to the center of the Milky Way and has passed through the Galactic disk many times, where its structure can be perturbed by baryonic structures such as the Galactic bar and giant molecular clouds (GMCs). In order to understand how these baryonic structures affect Pal 5's density, we present a detailed study of the effects of the Galactic bar, spiral structure, GMCs, and globular clusters on the Pal 5 stream. We estimate the effect of each perturber on the stream density by computing its power spectrum and comparing it to the power induced by a CDM-like population of dark matter subhalos. We find that the bar and GMCs can each individually create power that is comparable to the observed power on large scales, leaving little room for dark matter substructure, while spirals are subdominant on all scales. On degree scales, the power induced by the bar is small, but GMCs create small-scale density variations that are similar in amplitude to the dark-matter induced variations but otherwise indistinguishable from it. These results demonstrate that Pal 5 is a poor system for constraining the dark matter substructure fraction and that observing streams further out in the halo will be necessary to confidently detect dark matter subhalos.

Deep learning with artificial neural networks is increasingly gaining attention, because of its potential for data-driven astronomy. However, this methodology usually does not provide uncertainties and does not deal with incompleteness and noise in the training data. In this work, we design a neural network for high-resolution spectroscopic analysis using APOGEE data that mimics the methodology of standard spectroscopic analyses: stellar parameters are determined using the full wavelength range, but individual element abundances use censored portions of the spectrum. We train this network with a customized objective function that deals with incomplete and noisy training data and apply dropout variational inference to derive uncertainties on our predictions. We determine parameters and abundances for 18 individual elements at the $\approx 0.03$ dex level, even at low signal-to-noise ratio. We demonstrate that the uncertainties returned by our method are a realistic estimate of the precision and they automatically blow up when inputs or outputs outside of the training set are encountered, thus shielding users from unwanted extrapolation. By using standard deep-learning tools for GPU acceleration, our method is extremely fast, allowing analysis of the entire APOGEE data set of $\approx250,000$ spectra in ten minutes on a single, low-cost GPU. We release the stellar parameters and 18 individual-element abundances with associated uncertainty for the entire APOGEE DR14 dataset. Simultaneously, we release astroNN, a well-tested, open-source python package developed for this work, but that is also designed to be a general package for deep learning in astronomy. astroNN is available at https://github.com/henrysky/astroNN with extensive documentation at http://astroNN.readthedocs.io.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

The vertical structure and dynamics of stars in our local Galactic neighbourhood contains much information about the local distribution of visible and dark matter and of perturbations to the Milky Way disc. We use data on the positions and velocities of stars in the solar neighbourhood from \gaia\ DR2 and large spectroscopic surveys to investigate the vertical number counts and mean-velocity trend as a function of distance from the local Galactic mid-plane. We perform a detailed measurement of the wave-like North-South asymmetry in the vertical number counts, which reveals a number of deficits at heights $\approx 0.4\,\mathrm{kpc}$, $\approx 0.9\,\mathrm{kpc}$, and $\approx 1.5\,\mathrm{kpc}$, and peaks at $\approx 0.2\,\mathrm{kpc}$, $\approx 0.7\,\mathrm{kpc}$, and $\approx 1.1\,\mathrm{kpc}$. We find that the asymmetry pattern is independent of colour. The mean vertical velocity is almost constant to $<1\,\mathrm{km\,s}^{-1}$ within a few 100 pc from the mid-plane and then displays a North-South symmetric dip at $\approx0.5\,\mathrm{kpc}$ with an amplitude of $\approx 2\,\mathrm{km\,s}^{-1}$ that is a plausible velocity counterpart to the main number-count dip at a similar height. Thus, with \gaia\ DR2 we confirm at high fidelity that the local Galactic disc is undergoing a wave-like oscillation and a dynamically-consistent observational picture of the perturbed local vertical structure emerges for the first time. We also present the most precise and accurate determination of the Sun's height above the local Galactic mid-plane, correcting for any asymmetry in the vertical density: $z_\odot = 20.8 \pm 0.3\,\mathrm{pc}$.

Recent work indicates that the nearby Galactic halo is dominated by the debris from a major accretion event. We confirm that result from an analysis of APOGEE-DR14 element abundances and $\textit{Gaia}$-DR2 kinematics of halo stars. We show that $\sim$2/3 of nearby halo stars have high orbital eccentricities ($e \gtrsim 0.8$), and abundance patterns typical of massive Milky Way dwarf galaxy satellites today, characterised by relatively low [Fe/H], [Mg/Fe], [Al/Fe], and [Ni/Fe]. The trend followed by high $e$ stars in the [Mg/Fe]-[Fe/H] plane shows a change of slope at [Fe/H]$\sim-1.3$, which is also typical of stellar populations from relatively massive dwarf galaxies. Low $e$ stars exhibit no such change of slope within the observed [Fe/H] range and show slightly higher abundances of Mg, Al and Ni. Unlike their low $e$ counterparts, high $e$ stars show slightly retrograde motion, make higher vertical excursions and reach larger apocentre radii. By comparing the position in [Mg/Fe]-[Fe/H] space of high $e$ stars with those of accreted galaxies from the EAGLE suite of cosmological simulations we constrain the mass of the accreted satellite to be in the range $10^{8.5}\lesssim M_*\lesssim 10^{9}\mathrm{M_\odot}$. We show that the median orbital eccentricities of debris are largely unchanged since merger time, implying that this accretion event likely happened at $z\lesssim1.5$. The exact nature of the low $e$ population is unclear, but we hypothesise that it is a combination of $\textit{in situ}$ star formation, high $|z|$ disc stars, lower mass accretion events, and contamination by the low $e$ tail of the high $e$ population. Finally, our results imply that the accretion history of the Milky Way was quite unusual.

We apply a simple axisymmetric disc model to 218 Galactic Cepheids whose accurate measurements of the distance and velocities are obtained by cross-matching an existing Cepheids catalogue with the Gaia DR2 data. Our model fit determines the "local centrifugal speed", $V_\mathrm{c}$ $-$ defined as the rotation speed required to balance the local radial gravitational force $-$ at the Sun's location to be $V_{c}(R_0)=236\pm 3$ km s$^{-1}$ and the Sun's azimuthal and radial peculiar motions to be $V_{\odot}=12.4\pm0.7$ km s$^{-1}$ and $U_{\odot}=7.7\pm0.9$ km s$^{-1}$, respectively. These results are obtained with strong priors on the solar radius, $R_0=8.2\pm0.1$ kpc, and Sun's angular rotation velocity, $\Omega_{\odot}=30.24\pm0.12$ km s$^{-1}$ kpc$^{-1}$. We also applied the axisymmetric model to mock data from an N-body/hydrodynamic simulation of a Milky Way-like galaxy with a bar and spiral arms. We find that our axisymmetric model fit to the young stars recovers the local centrifugal speed reasonably well, even in the face of significant non-axisymmetry. However, the local centrifugal speed determined from our Cepheid sample could suffer from systematic uncertainty as large as 6 km s$^{-1}$.

We present a search for new members of the 300 km s$^{-1}$ stream (300S) near the dwarf galaxy Segue 1 using wide-field survey data. We identify 11 previously unknown bright stream members in the APOGEE-2 and SEGUE-1 and 2 spectroscopic surveys. Based on the spatial distribution of the high-velocity stars, we confirm for the first time that this kinematic structure is associated with a 24$^{\circ}$-long stream seen in SDSS and Pan-STARRS imaging data. The 300S stars display a metallicity range of $-2.17 < {\rm [Fe/H]} < -1.24$, with an intrinsic dispersion of 0.21$_{-0.09}^{+0.12}$ dex. They also have chemical abundance patterns similar to those of Local Group dwarf galaxies, as well as that of the Milky Way halo. Using the open-source code galpy to model the orbit of the stream, we find that the progenitor of the stream passed perigalacticon about 70 Myr ago, with a closest approach to the Galactic Center of about 4.1 kpc. Using Pan-STARRS DR1 data, we obtain an integrated stream luminosity of $4 \times 10^3$ L$_{\odot}$. We conclude that the progenitor of the stream was a dwarf galaxy that is probably similar to the satellites that were accreted to build the present-day Milky Way halo.

The second data release from ESA's Gaia mission has revealed many ridge-like structures in the velocity distribution of the Milky Way. We show that these can arise naturally from winding transient spiral structure that is commonly seen in N-body simulations of disk galaxies. We construct test particle models of the winding spiral structure, and compare the resulting distribution of orbits with the observed two-dimensional velocity distribution in the extended solar neighbourhood and with the distribution of rotational velocities over 8 kpc along the Sun--Galactic-centre--Galactic anti-centre line. We show that the ridges in these observations are well reproduced by the winding spiral model. Additionally, we demonstrate that the transient winding spiral potential can create a Hercules-like feature in the kinematics of the solar neighbourhood, either alone, or in combination with a long-slow bar potential.

Orbital parameters, such as eccentricity and maximum vertical excursion, of stars in the Milky Way are an important tool for understanding its dynamics and evolution, but calculation of such parameters usually relies on computationally-expensive numerical orbit integration. We present and test a fast method for estimating these parameters using an application of the St\"ackel fudge, used previously for the estimation of action-angle variables. We show that the method is highly accurate, to a level of $<1\%$ in eccentricity, over a large range of relevant orbits and in different Milky Way-like potentials, and demonstrate its validity by estimating the eccentricity distribution of the RAVE-TGAS data set and comparing it to that from orbit integration. Using the method, the orbital characteristics of the $\sim 7$ million $\textit{Gaia}$ DR2 stars with radial velocity measurements are computed with Monte Carlo sampled errors in $\sim 116$ hours of parallelised cpu time, at a speed that we estimate to be $\sim 3$ to $4$ orders of magnitude faster than using numerical orbit integration. We demonstrate using this catalogue that $\textit{Gaia}$ DR2 samples a large range of orbits in the solar vicinity, down to those with $r_\mathrm{peri} \lesssim 2.5$ kpc, and out to $r_\mathrm{ap} \gtrsim 13$ kpc. We also show that many of the features present in orbital parameter space have a low mean $z_\mathrm{max}$, suggesting that they likely result from disk dynamical effects.

Data and analysis methodology used for the SDSS/APOGEE Data Releases 13 and 14 are described, highlighting differences from the DR12 analysis presented in Holtzman (2015). Some improvement in the handling of telluric absorption and persistence is demonstrated. The derivation and calibration of stellar parameters, chemical abundances, and respective uncertainties are described, along with the ranges over which calibration was performed. Some known issues with the public data related to the calibration of the effective temperatures (DR13), surface gravity (DR13 and DR14), and C and N abundances for dwarfs (DR13 and DR14) are highlighted. We discuss how results from a data-driven technique, The Cannon (Casey 2016), are included in DR14, and compare those with results from the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP). We describe how using The Cannon in a mode that restricts the abundance analysis of each element to regions of the spectrum with known features from that element leads to Cannon abundances can lead to significantly different results for some elements than when all regions of the spectrum are used to derive abundances.

A key prediction of the standard cosmological model -- which relies on the assumption that dark matter is cold, i.e. non-relativistic at the epoch of structure formation -- is the existence of a large number of dark matter substructures on sub-galactic scales. This assumption can be tested by studying the perturbations induced by dark matter substructures on cold stellar streams. Here, we study the prospects for discriminating cold from warm dark matter by generating mock data for upcoming astronomical surveys such as the Large Synoptic Survey Telescope (LSST), and reconstructing the properties of the dark matter particle from the perturbations induced on the stellar density profile of a stream. We discuss the statistical and systematic uncertainties, and show that the method should allow to set stringent constraints on the mass of thermal dark matter relics, and possibly to yield an actual measurement of the dark matter particle mass if it is in the $\mathcal{O}(1)$ keV range.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

In the age of high-resolution spectroscopic stellar surveys of the Milky Way, the number of stars with detailed abundances of multiple elements is rapidly increasing. These elemental abundances are directly influenced by the evolutionary history of the Galaxy, but this can be difficult to interpret without an absolute timeline of the abundance enrichment. We present age-abundance trends for [M/H], [{\alpha}/M], and 17 individual elements using a sample of 721 solar neighbourhood Hipparcos red giant stars observed by APOGEE. These age trends are determined through a Bayesian hierarchical modelling method presented by Feuillet et al. (2016). We confirm that the [{\alpha}/M]- age relation in the solar neighbourhood is steep and relatively narrow (0.20 dex age dispersion), as are the [O/M]- and [Mg/M]-age relations. The age trend of [C/N] is steep and smooth, consistent with stellar evolution. The [M/H]-age relation has a mean age dispersion of 0.28 dex and a complex overall structure. The oldest stars in our sample are those with the lowest and highest metallicities, while the youngest stars are those with solar metallicity. These results provide strong constraints on theoretical models of Galactic chemical evolution (GCE). We compare them to the predictions of one-zone GCE mod- els and multi-zone mixtures, both analytic and numerical. These comparisons support the hypothesis that the solar neighbourhood is composed of stars born at a range of Galactocentric radii, and that the most metal-rich stars likely migrated from a region with earlier and more rapid star formation such as the inner Galaxy.

There are multiple groups of comoving stars in the Solar neighbourhood, which can potentially be explained as the signatures of one of the fundamental resonances of non-axisymmetric structure such as the Galactic bar or spiral arms. One such stream, Hercules, has been proposed to result from the outer Lindblad resonance (OLR) of a short fast rotating bar as shown analytically, or the corotation resonance (CR) of a longer slower rotating bar as observed in an N-body model. We show that by including an m = 4 Fourier component in an analytical long bar model, with an amplitude that is typical for bars in N-body simulations, we can reproduce a Hercules like feature in the kinematics of the Solar neighbourhood. We then describe the expected symmetry in the velocity distribution arising from such a model, which we will soon be able to test with Gaia.

Chemical tagging of stars based on their similar compositions can offer new insights about the star formation and dynamical history of the Milky Way. We investigate the feasibility of identifying groups of stars in chemical space by forgoing the use of model derived abundances in favour of direct analysis of spectra. This facilitates the propagation of measurement uncertainties and does not presuppose knowledge of which elements are important for distinguishing stars in chemical space. We use ~16,000 red-giant and red-clump H-band spectra from the Apache Point Observatory Galactic Evolution Experiment and perform polynomial fits to remove trends not due to abundance-ratio variations. Using expectation maximized principal component analysis, we find principal components with high signal in the wavelength regions most important for distinguishing between stars. Different subsamples of red-giant and red-clump stars are all consistent with needing about 10 principal components to accurately model the spectra above the level of the measurement uncertainties. The dimensionality of stellar chemical space that can be investigated in the H-band is therefore $\lesssim 10$. For APOGEE observations with typical signal-to-noise ratios of 100, the number of chemical space cells within which stars cannot be distinguished is approximately $10^{10\pm2} \times (5\pm 2)^{n-10}$ with $n$ the number of principal components. This high dimensionality and the fine-grained sampling of chemical space are a promising first step towards chemical tagging based on spectra alone.

We use the multi-epoch radial velocities acquired by the APOGEE survey to perform a large scale statistical study of stellar multiplicity for field stars in the Milky Way, spanning the evolutionary phases between the main sequence and the red clump. We show that the distribution of maximum radial velocity shifts (\drvm) for APOGEE targets is a strong function of \logg, with main sequence stars showing \drvm\ as high as $\sim$300 \kms, and steadily dropping down to $\sim$30 \kms\ for \logg$\sim$0, as stars climb up the Red Giant Branch (RGB). Red clump stars show a distribution of \drvm\ values comparable to that of stars at the tip of the RGB, implying they have similar multiplicity characteristics. The observed attrition of high \drvm\ systems in the RGB is consistent with a lognormal period distribution in the main sequence and a multiplicity fraction of 0.35, which is truncated at an increasing period as stars become physically larger and undergo mass transfer after Roche Lobe Overflow during H shell burning. The \drvm\ distributions also show that the multiplicity characteristics of field stars are metallicity dependent, with metal-poor ([Fe/H]$\lesssim-0.5$) stars having a multiplicity fraction a factor 2-3 higher than metal-rich ([Fe/H]$\gtrsim0.0$) stars. This has profound implications for the formation rates of interacting binaries observed by astronomical transient surveys and gravitational wave detectors, as well as the habitability of circumbinary planets.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) provides spectroscopic information of regions of the inner Milky Way inaccessible to optical surveys. We present the first large study of the metallicity distribution of the innermost Galactic regions based on homogeneous measurements from the SDSS Data Release 12 for 7545 red giant stars within 4.5 kpc of the Galactic center, with the goal to shed light on the structure and origin of the Galactic bulge. Stellar metallicities are found, through multiple-Gaussian decompositions, to be distributed in several components indicative of the presence of various stellar populations such as the bar, or the thin and the thick disk. A super-solar ([Fe/H]=+0.32) and a solar ([Fe/H]=+0.00) metallicity components, tentatively associated with the thin disk and the Galactic bar, respectively, seem to be the major contributors near the midplane. The solar-metallicity component extends outwards in the midplane but is not observed in the innermost regions. The central regions (within 3 kpc of the Galactic center) reveal, on the other hand, the presence of a significant metal-poor population ([Fe/H] =-0.46), tentatively associated with the thick disk, and which becomes the dominant component far from the midplane ($|Z| > +0.75 kpc). Varying contributions from these different components produce a transition region at +0.5 kpc < |Z| < +1.0 kpc characterized by a significant vertical metallicity gradient.

Among dynamical modeling techniques, the made-to-measure (M2M) method for modeling steady-state systems is among the most flexible, allowing non-parametric distribution functions in complex gravitational potentials to be modeled efficiently using N-body particles. Here we propose and test various improvements to the standard M2M method for modeling observed data, illustrated using the simple setup of a one-dimensional harmonic oscillator. We demonstrate that nuisance parameters describing the modeled system's orientation with respect to the observer---e.g., an external galaxy's inclination or the Sun's position in the Milky Way---as well as the parameters of an external gravitational field can be optimized simultaneously with the particle weights. We develop a method for sampling from the high-dimensional uncertainty distribution of the particle weights. We combine this in a Gibbs sampler with samplers for the nuisance and potential parameters to explore the uncertainty distribution of the full set of parameters. We illustrate our M2M improvements by modeling the vertical density and kinematics of F-type stars in Gaia DR1. The novel M2M method proposed here allows full probabilistic modeling of steady-state dynamical systems, allowing uncertainties on the non-parametric distribution function and on nuisance parameters to be taken into account when constraining the dark and baryonic masses of stellar systems.

SDSS-V will be an all-sky, multi-epoch spectroscopic survey of over six million objects. It is designed to decode the history of the Milky Way, trace the emergence of the chemical elements, reveal the inner workings of stars, and investigate the origin of planets. It will also create an integral-field spectroscopic map of the gas in the Galaxy and the Local Group that is 1,000x larger than the current state of the art and at high enough spatial resolution to reveal the self-regulation mechanisms of galactic ecosystems. SDSS-V will pioneer systematic, spectroscopic monitoring across the whole sky, revealing changes on timescales from 20 minutes to 20 years. The survey will thus track the flickers, flares, and radical transformations of the most luminous persistent objects in the universe: massive black holes growing at the centers of galaxies. The scope and flexibility of SDSS-V will be unique among extant and future spectroscopic surveys: it is all-sky, with matched survey infrastructures in both hemispheres; it provides near-IR and optical multi-object fiber spectroscopy that is rapidly reconfigurable to serve high target densities, targets of opportunity, and time-domain monitoring; and it provides optical, ultra-wide-field integral field spectroscopy. SDSS-V, with its programs anticipated to start in 2020, will be well-timed to multiply the scientific output from major space missions (e.g., TESS, Gaia, eROSITA) and ground-based projects. SDSS-V builds on the 25-year heritage of SDSS's advances in data analysis, collaboration infrastructure, and product deliverables. The project is now refining its science scope, optimizing the survey strategies, and developing new hardware that builds on the SDSS-IV infrastructure. We present here an overview of the current state of these developments as we seek to build our worldwide consortium of institutional and individual members.

The Hercules stream is a group of co-moving stars in the Solar neighbourhood, which can potentially be explained as a signature of either the outer Lindblad resonance (OLR) of a fast Galactic bar or the corotation resonance of a slower bar. In either case, the feature should be present over a large area of the disc. With the recent commissioning of the APOGEE-2 Southern spectrograph we can search for the Hercules stream at $(l,b)=(270^\circ,0)$, a direction in which the Hercules stream, if caused by the bar's OLR, would be strong enough to be detected using only the line-of-sight velocities. We clearly detect a narrow, Hercules-like feature in the data that can be traced from the solar neighbourhood to a distance of about 4 kpc. The detected feature matches well the line-of-sight velocity distribution from the fast-bar (OLR) model. Confronting the data with a model where the Hercules stream is caused by the corotation resonance of a slower bar leads to a poorer match, as the corotation model does not predict clearly separated modes, possibly because the slow-bar model is too hot.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

APOGEE-2 is a high-resolution, near-infrared spectroscopic survey observing roughly 300,000 stars across the entire sky. It is the successor to APOGEE and is part of the Sloan Digital Sky Survey IV (SDSS-IV). APOGEE-2 is expanding upon APOGEE's goals of addressing critical questions of stellar astrophysics, stellar populations, and Galactic chemodynamical evolution using (1) an enhanced set of target types and (2) a second spectrograph at Las Campanas Observatory in Chile. APOGEE-2 is targeting red giant branch (RGB) and red clump (RC) stars, RR Lyrae, low-mass dwarf stars, young stellar objects, and numerous other Milky Way and Local Group sources across the entire sky from both hemispheres. In this paper, we describe the APOGEE-2 observational design, target selection catalogs and algorithms, and the targeting-related documentation included in the SDSS data releases.

The best way to trace back the history of star formation and mass assembly of the Milky Way disc is by combining chemical compositions, ages and phase-space information for a large number of disc stars. With the advent of large surveys of the stellar populations of the Galaxy, such data have become available and can be used to pose constraints on sophisticated models of galaxy formation. We use SDSS-III/APOGEE data to derive the first detailed 3D map of stellar density in the Galactic disc as a function of age, [Fe/H] and $\mathrm{[\alpha/Fe]}$. We discuss the implications of our results for the formation and evolution of the disc, presenting new constraints on the disc structural parameters, stellar radial migration and its connection with disc flaring. We also discuss how our results constrain the inside out formation of the disc, and determine the surface-mass density contributions at the solar radius for mono-age, mono-[Fe/H] populations.

The absolute number and the density profiles of different types of stars in the solar neighborhood are a fundamental anchor for studies of the initial mass function, stellar evolution, and galactic structure. Using data from the Gaia DR1 Tycho-Gaia Astrometric Solution, we reconstruct Gaia's selection function and we determine Gaia's volume completeness, the local number density, and the vertical profiles of different spectral types along the main sequence from early A stars to late K stars as well as along the giant branch. We clearly detect the expected flattening of the stellar density profile near the mid-plane for all stellar types: All vertical profiles are well represented by sech^2 profiles, with scale heights ranging from ~50 pc for A stars to ~150 pc for G and K dwarfs and giants. We determine the luminosity function along the main sequence for M_V < 7 (M >~ $0.72 M_\odot$) and along the giant branch for M_J >~ -2.5. Converting this to a mass function, we find that the high-mass (M > $1\,M_\odot$) present-day mass function along the main sequence is d n / d M = 0.016 $(M/M_\odot)^{-4.7}$ stars/pc^3/$M_\odot$. Extrapolating below M = $0.72\,M_\odot$, we find a total mid-plane stellar density of 0.040+/-0.002 $M_\odot$/pc^3. Giants contribute 0.00039+/-0.00001 stars/pc^3 or about 0.00046+/-0.00005 $M_\odot$/pc^3. The star-formation rate surface density is \Sigma(t) = 7+/-1 exp(-t/[7+/-1 Gyr]) $M_\odot$/pc^2/Gyr. Overall, we find that Gaia DR1's selection biases are manageable and allow a detailed new inventory of the solar neighborhood to be made that agrees with and extends previous studies. This bodes well for mapping the Milky Way with the full Gaia data set.

Distances to individual stars in our own Galaxy are critical in order to piece together the nature of its velocity and spatial structure. Core helium burning red clump (RC) stars have similar luminosities, are abundant throughout the Galaxy, and thus constitute good standard candles. We build a hierarchical probabilistic model to quantify the quality of RC stars as standard candles using parallax measurements from the first Gaia data release. A unique aspect of our methodology is to fully account for (and marginalize over) parallax, photometry, and dust corrections uncertainties, which leads to more robust results than standard approaches. We determine the absolute magnitude and intrinsic dispersion of the RC in 2MASS bands J, H, Ks, Gaia G band, and WISE bands W1, W2, W3, and W4. We find that the absolute magnitude of the RC is $-1.61 \pm$ 0.01 (in Ks), $+0.44 \pm$ 0.01 (in G) , $-0.93 \pm$ 0.01 (in J), $-1.46 \pm$ 0.01 (in H), $-1.68 \pm$ 0.02 (in W1), $-1.69\pm$ 0.02 (in W2), $-1.67 \pm$ 0.02 (in W3), $1.76 \pm$ 0.01 mag (in W4). The mean intrinsic dispersion is $\sim 0.17 \pm$ 0.03 mag across all bands (yielding a typical distance precision of $\sim$ 8%). Thus RC stars are reliable and precise standard candles. In addition, we have also re-calibrated the zero point of the absolute magnitude of the RC in each band, which provide a benchmark for future studies to estimate distances to RC stars. Finally, the parallax error shrinkage in the hierarchical model outlined in this work can be used to obtain more precise parallaxes than Gaia for the most distant RC stars across the Galaxy.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

The SDSS-III APOGEE DR12 is a unique resource to search for stars beyond the tidal radii of star clusters. We have examined the APOGEE DR12 database for new candidates of the young star cluster Palomar 1, a system with previously reported tidal tails (Niederste-Ostholt et al. 2010). The APOGEE ASPCAP database includes spectra and stellar parameters for two known members of Pal 1 (Stars I and II), however these do not agree with the stellar parameters determined from optical spectra by Sakari et al. (2011). We find that the APOGEE analysis of these two stars is strongly affected by the known persistence problem (Majewski et al. 2015; Nidever et al. 2015). By re-examining the individual visits, and removing the blue (and sometimes green) APOGEE detector spectra affected by persistence, then we find excellent agreement in a re-analysis of the combined spectra. These methods are applied to another five stars in the APOGEE field with similar radial velocities and metallicities as those of Pal 1. Only one of these new candidates, Star F, may be a member located in the tidal tail based on its heliocentric radial velocity, metallicity, and chemistry. The other four candidates are not well aligned with the tidal tails, and comparison to the Besancon model (Robin et al. 2003) suggests that they are more likely to be non-members, i.e. part of the Galactic halo. This APOGEE field could be re-examined for other new candidates if the persistence problem can be removed from the APOGEE spectral database.

Converting a noisy parallax measurement into a posterior belief over distance requires inference with a prior. Usually this prior represents beliefs about the stellar density distribution of the Milky Way. However, multi-band photometry exists for a large fraction of the \textsl{\small{Gaia}} \textsl{\small{TGAS}} Catalog and is incredibly informative about stellar distances. Here we use \textsl{\small{2MASS}} colors for 1.4 million \textsl{\small{TGAS}} stars to build a noise-deconvolved empirical prior distribution for stars in color--magnitude space. This model contains no knowledge of stellar astrophysics or the Milky Way, but is precise because it accurately generates a large number of noisy parallax measurements under an assumption of stationarity; that is, it is capable of combining the information from many stars. We use the Extreme Deconvolution (\textsl{\small{XD}}) algorithm---an Empirical Bayes approximation to a full hierarchical model of the true parallax and photometry of every star---to construct this prior. The prior is combined with a \textsl{\small{TGAS}} likelihood to infer a precise photometric parallax estimate and uncertainty (and full posterior) for every star. Our parallax estimates are more precise than the \textsl{\small{TGAS}} catalog entries by a median factor of 1.2 (14% are more precise by a factor >2) and are more precise than previous Bayesian distance estimates that use spatial priors. We validate our parallax inferences using members of the Milky Way star cluster M67, which is not visible as a cluster in the \textsl{\small{TGAS}} parallax estimates, but appears as a cluster in our posterior parallax estimates. Our results, including a parallax posterior pdf for each of 1.4 million \textsl{\small{TGAS}} stars, are available in companion electronic tables.

The measurement of the structure of stellar populations in the Milky Way disk places fundamental constraints on models of galaxy formation and evolution. Previously, the disk's structure has been studied in terms of populations defined geometrically and/or chemically, but a decomposition based on stellar ages provides a more direct connection to the history of the disk, and stronger constraint on theory. Here, we use positions, abundances and ages for 31,244 red giant branch stars from the SDSS-APOGEE survey, spanning $3 < R_{\mathrm{gc}} < 15$ kpc, to dissect the disk into mono-age and mono-[Fe/H] populations at low and high [$\alpha$/Fe]. For each population, with $\Delta \mathrm{age} < 2$ Gyr and $\Delta \mathrm{[Fe/H]} < 0.1$ dex, we measure the structure and surface-mass density contribution. We find that low [$\alpha$/Fe] mono-age populations are fit well by a broken exponential, which increases to a peak radius and decreases thereafter. We show that this profile becomes broader with age, interpreted here as a new signal of disk heating and radial migration. High [$\alpha$/Fe] populations are well fit as single exponentials within the radial range considered, with an average scale length of $1.9\pm 0.1$ kpc. We find that the relative contribution of high to low [$\alpha$/Fe] populations at $R_0$ is $f_\Sigma = 18\% \pm 5\%$; high [$\alpha$/Fe] contributes most of the mass at old ages, and low [$\alpha$/Fe] at young ages. The low and high [$\alpha$/Fe] populations overlap in age at intermediate [Fe/H], although both contribute mass at $R_{0}$ across the full range of [Fe/H]. The mass weighted scale height $h_Z$ distribution is a smoothly declining exponential function. High [$\alpha$/Fe] populations are thicker than low [$\alpha$/Fe], and the average $h_Z$ increases steadily with age, between 200 and 600 pc.

Red clump stars are fundamental distance indicators in astrophysics, although theoretical stellar models predict a dependence of absolute magnitudes with ages. This effect is particularly strong below 2 Gyr, but even above this limit a mild age dependence is still expected. We use seismically identified red clump stars in the Kepler field for which we have reliable distances, masses and ages from the SAGA survey to first explore this effect. By excluding red clump stars with masses larger than 1.6 Msun (corresponding to ages younger than 2 Gyr), we derive robust calibrations linking intrinsic colors to absolute magnitudes in the following photometric systems: Str\"omgren $by$, Johnson $BV$, Sloan $griz$, 2MASS $JHK_s$ and WISE $W1W2W3$. With the precision achieved we also detect a slope of absolute magnitudes 0.020(0.003) mag per Gyrin the infrared, implying that distance calibrations of clump stars can be off by up to 0.2 mag in the infrared (over the range from 2 Gyr to 12 Gyr) if their ages are unknown. Even larger uncertainties affect optical bands, because of the stronger interdependency of absolute magnitudes on colors and age. Our distance calibrations are ultimately based on asteroseismology, and we show how the distance scale can be used to test the accuracy of seismic scaling relations. Within the uncertainties our calibrations are in agreement with those built upon local red clump with Hipparcos} parallaxes, although we find a tension which if confirmed would imply that scaling relations overestimate radii of red clump stars by 2(+-20%. Data-releases post Gaia DR1 will provide an important testbed for our results.

We present a proper motion measurement for the halo globular cluster Pyxis, using HST/ACS data as the first epoch, and GeMS/GSAOI Adaptive Optics data as the second, separated by a baseline of about 5 years. This is both the first measurement of the proper motion of Pyxis and the first calibration and use of Multi-Conjugate Adaptive Optics data to measure an absolute proper motion for a faint, distant halo object. Consequently, we present our analysis of the Adaptive Optics data in detail. We obtain a proper motion of mu_alpha cos(delta)=1.09+/-0.31 mas/yr and mu_delta=0.68+/-0.29 mas/yr. From the proper motion and the line-of-sight velocity we find the orbit of Pyxis is rather eccentric with its apocenter at more than 100 kpc and its pericenter at about 30 kpc. We also investigate two literature-proposed associations for Pyxis with the recently discovered ATLAS stream and the Magellanic system. Combining our measurements with dynamical modeling and cosmological numerical simulations we find it unlikely Pyxis is associated with either system. We examine other Milky Way satellites for possible association using the orbit, eccentricity, metallicity, and age as constraints and find no likely matches in satellites down to the mass of Leo II. We propose that Pyxis probably originated in an unknown galaxy, which today is fully disrupted. Assuming that Pyxis is bound and not on a first approach, we derive a 68% lower limit on the mass of the Milky Way of 0.95*10^12 M_sun.

RoadMapping is a dynamical modeling machinery developed to constrain the Milky Way's (MW) gravitational potential by simultaneously fitting an axisymmetric parametrized potential and an action-based orbit distribution function (DF) to discrete 6D phase-space measurements of stars in the Galactic disk. In this work we demonstrate RoadMapping's robustness in the presence of spiral arms by modeling data drawn from an N-body simulation snapshot of a disk-dominated galaxy of MW mass with strong spiral arms (but no bar), exploring survey volumes with radii 500pc<=r_max<=5kpc. The potential constraints are very robust, even though we use a simple action-based DF, the quasi-isothermal DF (qDF). The best-fit RoadMapping model always recovers the correct gravitational forces where most of the stars that entered the analysis are located, even for small volumes. For data from large survey volumes, RoadMapping finds axisymmetric models that average well over the spiral arms. Unsurprisingly, the models are slightly biased by the excess of stars in the spiral arms. Gravitational potential models derived from survey volumes with at least r_max=3kpc can be reliably extrapolated to larger volumes. However, a large radial survey extent, r_max~5kpc, is needed to correctly recover the halo scale length. In general, the recovery and extrapolability of potentials inferred from data sets which were drawn from inter-arm regions appear to be better than those of data sets drawn from spiral arms. Our analysis implies that building axisymmetric models for the Galaxy with upcoming Gaia data will lead to sensible and robust approximations of the MW's potential.

The spatial variations of the velocity field of local stars provide direct evidence of Galactic differential rotation. The local divergence, shear, and vorticity of the velocity field---the traditional Oort constants---can be measured based purely on astrometric measurements and in particular depend linearly on proper motion and parallax. I use data for 304,267 main-sequence stars from the Gaia DR1 Tycho-Gaia Astrometric Solution to perform a local, precise measurement of the Oort constants at a typical heliocentric distance of 230 pc. The pattern of proper motions for these stars clearly displays the expected effects from differential rotation. I measure the Oort constants to be: A = 15.3+/-0.4 km/s/kpc, B = -11.9+/-0.4 km/s/kpc, C = -3.2+/-0.4 km/s/kpc and K = -3.3+/-0.6 km/s/kpc, with no color trend over a wide range of stellar populations. These first confident measurements of C and K clearly demonstrate the importance of non-axisymmetry for the velocity field of local stars and they provide strong constraints on non-axisymmetric models of the Milky Way.

Tidal streams in the Milky Way are sensitive probes of the population of dark-matter subhalos predicted in cold-dark-matter (CDM) simulations. We present a new calculus for computing the effect of subhalo fly-bys on cold tidal streams based on the action-angle representation of streams. The heart of this calculus is a line-of-parallel-angle approach that calculates the perturbed distribution function of a given stream segment by undoing the effect of all impacts. This approach allows one to compute the perturbed stream density and track in any coordinate system in minutes for realizations of the subhalo distribution down to 10^5 Msun, accounting for the stream's internal dispersion and overlapping impacts. We study the properties of density and track fluctuations with suites of simulations. The one-dimensional density and track power spectra along the stream trace the subhalo mass function, with higher-mass subhalos producing power only on large scales, while lower mass subhalos cause structure on smaller scales. The time-dependence of impacts and of the evolution of the stream after an impact gives rise to bispectra. We further find that tidal streams are essentially corrugated sheets in the presence of subhalo perturbations: different projections of the track all reflect the same pattern of perturbations, facilitating their observational measurement. We apply this formalism to density data for the Pal 5 stream and make a first rigorous determination of 10^{+11}_{-6} dark-matter subhalos with masses between 3x10^6 and 10^9 Msun within 20 kpc from the Galactic center (corresponding to 1.4^{+1.6}_{-0.9} times the number predicted by CDM-only simulations or to f_{sub}(r<20 kpc) ~ 0.2%). Improved data will allow measurements of the subhalo mass function down to 10^5 Msun, thus definitively testing whether dark matter clumps on the smallest scales relevant for galaxy formation.

We constrain the shape of the Milky Way's halo by dynamical modeling of the observed phase-space tracks of the Pal 5 and GD-1 tidal streams. We find that the only information about the potential gleaned from the tracks of these streams are precise measurements of the shape of the gravitational potential---the ratio of vertical to radial acceleration---at the location of the streams, with weaker constraints on the radial and vertical accelerations separately. The latter will improve significantly with precise proper-motion measurements from Gaia. We measure that the overall potential flattening is 0.95 +/- 0.04 at the location of GD-1 ([R,z] ~ [12.5,6.7] kpc) and 0.94 +/- 0.05 at the position of Pal 5 ([R,z] ~ [8.4,16.8] kpc). Combined with constraints on the force field near the Galactic disk, we determine that the axis ratio of the dark-matter halo's density distribution is 1.05 +/- 0.14 within the inner 20 kpc, with a hint that the halo becomes more flattened near the edge of this volume. The halo mass within 20 kpc is 1.1 +/- 0.1 x 10^{11} M_sun. A dark-matter halo this close to spherical is in tension with the predictions from numerical simulations of the formation of dark-matter halos.

We report on the detection in the combined $Gaia$-DR1/RAVE data of a lack of disk stars in the solar neighbourhood with velocities close to zero angular momentum. We propose that this may be caused by the scattering of stars with very low angular momentum onto chaotic, halo-type orbits when they pass through the Galactic nucleus. We model the effect in a Milky-Way like potential and fit the resulting model directly to the data, finding a likelihood ($\sim2.7\sigma$) of a dip in the distribution. Using this effect, we can make a dynamical measurement of the Solar rotation velocity around the Galactic center: $v_{\odot}=239\pm9$ km s$^{-1}$. Combined with the measured proper motion of Sgr A$^*$, this measurement gives a measurement of the distance to the Galactic centre: $R_0=7.9\pm0.3$ kpc.

We present the Data Release 12 Quasar catalog (DR12Q) from the Baryon Oscillation Spectroscopic Survey (BOSS) of the SDSS-III. This catalog includes all SDSS-III/BOSS objects that were spectroscopically targeted as quasar candidates during the full survey and that are confirmed as quasars via visual inspection of the spectra, have luminosities Mi[z=2]<-20.5 (in a LCDM cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega _L=0.7), and either display at least one emission line with a full width at half maximum (FWHM)larger than 500 km/s or, if not, have interesting/complex absorption features. The catalog also includes previously known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 297,301 quasars detected over 9,376 square degrees with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with z>2.15 is about an order of magnitude greater than the number of z>2.15 quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 29,580 broad absorption line quasars and lists their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry together with some information on the optical morphology and the selection criteria. When available, the catalog also provides information on the optical variability of quasars using SDSS and PTF multi-epoch photometry. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra, covering the wavelength region 3,600-10,500A at a spectral resolution in the range 1,300<R<2,500, can be retrieved from the SDSS Catalog Archive Server.

Ample observational capabilities exist today to detect the small density perturbations that low-mass dark matter subhaloes impart on stellar streams from disrupting Galactic satellites. In anticipation of these observations, we investigate the expected number and size of gaps by combining an analytic prescription for gap evolution on circular orbits with the flux of subhaloes near the stream. We explore the distribution of gap sizes and depths for a typical cold stream around the Milky Way and find that for a given stream age and gap depth, each subhalo mass produces a characteristic gap size. For a stream with an age of a few Gyr, orbiting at a distance of 10-20 kpc from the Galactic center, even modest subhaloes with a mass of $10^6-10^7 M_\odot$ produce gaps with sizes that are on the order of several degrees. We consider the number and distribution of gap sizes created by subhaloes with masses $10^5-10^9 M_\odot$, accounting for the expected depletion of subhaloes by the Milky Way disk, and present predictions for six cold streams around the Milky Way. For Pal 5, we forecast 0.7 gaps with a density depletion of at least 25\% and a typical gap size of $8^\circ$. Thus, there appears to be {\it no tension} between the recent non-detection of density depletions in the Pal 5 tidal tails and $\Lambda$CDM expectations. These predictions can be used to guide the scale of future gap searches.

NGC 2420 is a $\sim$2 Gyr-old well-populated open cluster that lies about 2 kpc beyond the solar circle, in the general direction of the Galactic anti-center. Most previous abundance studies have found this cluster to be mildly metal-poor, but with a large scatter in the obtained metallicities for this open cluster. Detailed chemical abundance distributions are derived for 12 red-giant members of NGC 2420 via a manual abundance analysis of high-resolution (R = 22,500) near-infrared ($\lambda$1.5 - 1.7$\mu$m) spectra obtained from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. The sample analyzed contains 6 stars that are identified as members of the first-ascent red giant branch (RGB), as well as 6 members of the red clump (RC). We find small scatter in the star-to-star abundances in NGC 2420, with a mean cluster abundance of [Fe/H] = -0.16 $\pm$ 0.04 for the 12 red giants. The internal abundance dispersion for all elements (C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Co and Ni) is also very small ($\sim$0.03 - 0.06 dex), indicating a uniform cluster abundance distribution within the uncertainties. NGC 2420 is one of the clusters used to calibrate the APOGEE Stellar Parameter and Chemical Abundance Pipeline (ASPCAP). The results from this manual analysis compare well with ASPCAP abundances for most of the elements studied, although for Na, Al and V there are more significant offsets. No evidence of extra-mixing at the RGB luminosity bump is found in the $^{12}$C and $^{14}$N abundances from the pre-luminosity-bump RGB stars in comparison to the post-He core-flash RC stars.

Inferences about the spatial density or phase-space structure of stellar populations in the Milky Way require a precise determination of the effective survey volume. The volume observed by surveys such as Gaia or near-infrared spectroscopic surveys, which have good coverage of the Galactic mid-plane region, is highly complex because of the abundant small-scale structure in the three-dimensional interstellar dust extinction. We introduce a novel framework for analyzing the importance of small-scale structure in the extinction. This formalism demonstrates that the spatially-complex effect of extinction on the selection function of a pencil-beam or contiguous sky survey is equivalent to a low-pass filtering of the extinction-affected selection function with the smooth density field. We find that the angular resolution of current 3D extinction maps is sufficient for analyzing Gaia sub-samples of millions of stars. However, the current distance resolution is inadequate and needs to be improved by an order of magnitude, especially in the inner Galaxy. We also present a practical and efficient method for properly taking the effect of extinction into account in analyses of Galactic structure through an effective selection function. We illustrate its use with the selection function of red-clump stars in APOGEE using and comparing a variety of current 3D extinction maps.

The spatial structure of stellar populations with different chemical abundances in the Milky Way contains a wealth of information on Galactic evolution over cosmic time. We use data on 14,699 red-clump stars from the APOGEE survey, covering 4 kpc <~ R <~ 15 kpc, to determine the structure of mono-abundance populations (MAPs)---stars in narrow bins in [a/Fe] and [Fe/H]---accounting for the complex effects of the APOGEE selection function and the spatially-variable dust obscuration. We determine that all MAPs with enhanced [a/Fe] are centrally concentrated and are well-described as exponentials with a scale length of 2.2+/-0.2 kpc over the whole radial range of the disk. We discover that the surface-density profiles of low-[a/Fe] MAPs are complex: they do not monotonically decrease outwards, but rather display a peak radius ranging from ~5 kpc to ~13 kpc at low [Fe/H]. The extensive radial coverage of the data allows us to measure radial trends in the thickness of each MAP. While high-[a/Fe] MAPs have constant scale heights, low-[a/Fe] MAPs flare. We confirm, now with high-precision abundances, previous results that each MAP contains only a single vertical scale height and that low-[Fe/H], low-[a/Fe] and high-[Fe/H], high-[a/Fe] MAPs have intermediate (h_Z~300 to 600 pc) scale heights that smoothly bridge the traditional thin- and thick-disk divide. That the high-[a/Fe], thick disk components do not flare is strong evidence against their thickness being caused by radial migration. The correspondence between the radial structure and chemical-enrichment age of stellar populations is clear confirmation of the inside-out growth of galactic disks. The details of these relations will constrain the variety of physical conditions under which stars form throughout the MW disk.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

We present RoadMapping, a full-likelihood dynamical modelling machinery that aims to recover the Milky Way's (MW) gravitational potential from large samples of stars in the Galactic disk. RoadMapping models the observed positions and velocities of stars with a parametrized, three-integral distribution function (DF) in a parametrized axisymmetric potential. We investigate through differential test cases with idealized mock data how the breakdown of model assumptions and data properties affect constraints on the potential and DF. Our key results are: (i) If the MW's true potential is not included in the assumed model potential family, we can - in the axisymmetric case - still find a robust estimate for the potential, with only <~ 10% difference in surface density within |z| <= 1.1 kpc inside the observed volume. (ii) Modest systematic differences between the true and model DF are inconsequential. E.g, when binning stars to define sub-populations with simple DFs, binning errors do not affect the modelling as long as the DF parameters of neighbouring bins differ by < 20%. In addition, RoadMapping ensures unbiased potential estimates for either (iii) small misjudgements of the spatial selection function (i.e., <~ 15% at the survey volume's edge), (iv) if distances are known to within 10%, or (v) if proper motion uncertainties are known within 10% or are smaller than delta mu <~ 1 mas/yr. Challenges are the rapidly increasing computational costs for large sample sizes. Overall, RoadMapping is well suited to making precise new measurements of the MW's potential with data from the upcoming Gaia releases.

We develop a formalism for modelling the impact of dark matter subhaloes on cold thin streams. Our formalism models the formation of a gap in a stream in angle-frequency space and is able to handle general stream and impact geometry. We analyse an $N$-body simulation of a cold stream formed from a progenitor on an eccentric orbit in an axisymmetric potential, which is perturbed by a direct impact from a $10^8 M_\odot$ subhalo, and produce a complete generative model of the perturbed stream that matches the simulation well at a range of times. We show how the results in angle-frequency space can be related to physical properties of the gaps and that previous results for more constrained simulations are recovered. We demonstrate how our results are dependent upon the mass of the subhalo and the location of the impact along the stream. We find that gaps formed far downstream grow more rapidly than those closer to the progenitor due to the more ordered nature of the stream members far from the progenitor. Additionally, we show that the minimum gap density plateaus in time at a value that decreases with increasing subhalo mass.

Narrow stellar streams in the Milky Way halo are uniquely sensitive to dark-matter subhalos, but many of these subhalos may be tidally disrupted. I calculate the interaction between stellar and dark-matter streams using analytical and $N$-body calculations, showing that disrupting objects can be detected as low-concentration subhalos. Through this effect, we can constrain the lumpiness of the halo as well as the orbit and present position of individual dark-matter streams. This will have profound implications for the formation of halos and for direct and indirect-detection dark-matter searches.

Determining the level of chemical homogeneity in open clusters is of fundamental importance in the study of the evolution of star-forming clouds and that of the Galactic disk. Yet limiting the initial abundance spread in clusters has been hampered by difficulties in obtaining consistent spectroscopic abundances for different stellar types. Without reference to any specific model of stellar photospheres, a model for a homogeneous cluster is that it forms a one-dimensional sequence, with any differences between members due to variations in stellar mass and observational uncertainties. I present a novel method for investigating the abundance spread in open clusters that tests this one-dimensional hypothesis at the level of observed stellar spectra, rather than constraining homogeneity using derived abundances as traditionally done. Using high-resolution APOGEE spectra for 49 giants in M67, NGC 6819, and NGC 2420 I demonstrate that these spectra form one-dimensional sequences for each cluster. With detailed forward modeling of the spectra and Approximate Bayesian Computation, I derive strong limits on the initial abundance spread of 15 elements: <0.01 (0.02) dex for C and Fe, <~0.015 (0.03) dex for N, O, Mg, Si, and Ni, <~0.02 (0.03) dex for Al, Ca, and Mn, and <~0.03 (0.05) dex for Na, S, K, Ti, and V (at 68% and 95% confidence, respectively). The strong limits on C and O imply that no pollution by massive core-collapse supernovae occurred during star formation in open clusters, which, thus, need to form within <~6 Myr. Further development of this and related techniques will bring the power of differential abundances to stars other than solar twins in large spectroscopic surveys and will help unravel the history of star formation and chemical enrichment in the Milky Way through chemical tagging.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

As part of the Sloan Digital Sky Survey IV the extended Baryon Oscillation Spectroscopic Survey (eBOSS) will improve measurements of the cosmological distance scale by applying the Baryon Acoustic Oscillation (BAO) method to quasar samples. eBOSS will adopt two approaches to target quasars over 7500 sq. deg. First, a "CORE" quasar sample will combine optical selection in ugriz using a likelihood-based routine called XDQSOz, with a mid-IR-optical color-cut. eBOSS CORE selection (to g < 22 OR r < 22) should return ~ 70 quasars per sq. deg. at redshifts 0.9 < z < 2.2 and ~7 z > 2.1 quasars per sq. deg. Second, a selection based on variability in multi-epoch imaging from the Palomar Transient Factory should recover an additional ~3-4 z > 2.1 quasars per sq. deg. to g < 22.5. A linear model of how imaging systematics affect target density recovers the angular distribution of eBOSS CORE quasars over 96.7% (76.7%) of the SDSS North (South) Galactic Cap area. The eBOSS CORE quasar sample should thus be sufficiently dense and homogeneous over 0.9 < z < 2.2 to yield the first few-percent-level BAO constraint near z~1.5. eBOSS quasars at z > 2.1 will be used to improve BAO measurements in the Lyman-alpha Forest. Beyond its key cosmological goals, eBOSS should be the next-generation quasar survey, comprising > 500,000 new quasars and > 500,000 uniformly selected spectroscopically confirmed 0.9 < z < 2.2 quasars. At the conclusion of eBOSS, the SDSS will have provided unique spectra of over 800,000 quasars.

We present a sample of local red giant stars observed using the New Mexico State University 1 m telescope with the APOGEE spectrograph, for which we estimate stellar ages and the age distribution from the high-resolution spectroscopic stellar parameters and accurate distance measurements from Hipparcos. The high-resolution (R ~ 23,000), near infrared (H-band, 1.5-1.7 micron) APOGEE spectra provide measurements of the stellar atmospheric parameters (temperature, surface gravity, [M/H], and [alpha/M]). Due to the smaller uncertainties in surface gravity possible with high-resolution spectra and accurate Hipparcos distance measurements, we are able to calculate the stellar masses to within 40%. For red giants, the relatively rapid evolution of stars up the red giant branch allows the age to be constrained based on the mass. We examine methods of estimating age using both the mass-age relation directly and a Bayesian isochrone matching of measured parameters, assuming a constant star formation history (SFH). To improve the prior on the SFH, we use a hierarchical modeling approach to constrain the parameters of a model SFH from the age probability distribution functions of the data. The results of an alpha dependent Gaussian SFH model shows a clear relation between age and [alpha/M] at all ages. Using this SFH model as the prior for an empirical Bayesian analysis, we construct a full age probability distribution function and determine ages for individual stars. The age-metallicity relation is flat, with a slight decrease in [M/H] at the oldest ages and a ~ 0.5 dex spread in metallicity. For stars with ages < 1 Gyr we find a smaller spread, consistent with radial migration having a smaller effect on these young stars than on the older stars.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) has built the largest moderately high-resolution (R=22, 500) spectroscopic map of the stars across the Milky Way, and including dust-obscured areas. The APOGEE Stellar Parameter and Chemical Abundances Pipeline (ASPCAP) is the software developed for the automated analysis of these spectra. ASPCAP determines atmospheric parameters and chemical abundances from observed spectra by comparing observed spectra to libraries of theoretical spectra, using chi-2 minimization in a multidimensional parameter space. The package consists of a fortran90 code that does the actual minimization, and a wrapper IDL code for book-keeping and data handling. This paper explains in detail the ASPCAP components and functionality, and presents results from a number of tests designed to check its performance. ASPCAP provides stellar effective temperatures, surface gravities, and metallicities precise to 2%, 0.1 dex, and 0.05 dex, respectively, for most APOGEE stars, which are predominantly giants. It also provides abundances for up to 15 chemical elements with various levels of precision, typically under 0.1 dex. The final data release (DR12) of the Sloan Digital Sky Survey III contains an APOGEE database of more than 150,000 stars. ASPCAP development continues in the SDSS-IV APOGEE-2 survey.

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

The Apache Point Observatory Galactic Evolution Experiment (APOGEE), one of the programs in the Sloan Digital Sky Survey III (SDSS-III), has now completed its systematic, homogeneous spectroscopic survey sampling all major populations of the Milky Way. After a three year observing campaign on the Sloan 2.5-m Telescope, APOGEE has collected a half million high resolution (R~22,500), high S/N (>100), infrared (1.51-1.70 microns) spectra for 146,000 stars, with time series information via repeat visits to most of these stars. This paper describes the motivations for the survey and its overall design---hardware, field placement, target selection, operations---and gives an overview of these aspects as well as the data reduction, analysis and products. An index is also given to the complement of technical papers that describe various critical survey components in detail. Finally, we discuss the achieved survey performance and illustrate the variety of potential uses of the data products by way of a number of science demonstrations, which span from time series analysis of stellar spectral variations and radial velocity variations from stellar companions, to spatial maps of kinematics, metallicity and abundance patterns across the Galaxy and as a function of age, to new views of the interstellar medium, the chemistry of star clusters, and the discovery of rare stellar species. As part of SDSS-III Data Release 12, all of the APOGEE data products are now publicly available.

We characterise the radial density, metallicity and flattening profile of the Milky Way's stellar halo, based on the large sample of 1757 spectroscopically confirmed giant stars from SDSS/SEGUE-2 after excising stars that were algorithmically attributed to apparent halo substructure (including the Sagittarius stream). Compared to BHB stars or RR Lyrae, giants are more readily understood tracers of the overall halo star population, with less bias in age or metallicity. The well-characterized selection function of the sample enables forward modelling of those data, based on ellipsoidal stellar density models, $\nu_* (R,z)$, with Einasto profiles and (broken) power laws for their radial dependence, combined with a model for the metallicity gradient and the flattening profile. Among models with constant flattening, these data are reasonably well fit by an Einasto profile of $n=3.1\pm 0.5$ with an effective radius $\rm r_{eff} = 15\pm2~$kpc and a flattening of $q=0.7\pm 0.02$; or comparably well by an equally flattened broken power-law, with radial slopes of $\alpha_{in}=2.1\pm 0.3$ and $\alpha_{out}=3.8\pm 0.1$, with a break-radius of $r_{break}=18\pm1$~kpc; this is largely consistent with earlier work. We find a modest, but significant metallicity gradient within the "outer" stellar halo, $\rm [Fe/H]$ decreasing outward. If we allow for a variable flattening $q = f(r_{GC} )$, we find the distribution of halo giants to be considerably more flattened at small radii, $q({\rm 10~kpc})\sim 0.57$, compared to $q(>30{\rm kpc})\sim 0.8$. Remarkably, the data are then very well fit by a single power-law of index $\rm \sim 4.2\pm0.1$ of the variable $r_q\equiv\sqrt{R^2+(z/q(r))^2}$. In this simple and better fitting model, there is a break in flattening at $\sim 20$~kpc, instead of a break in the radial density function.

We investigate the kinematic signatures induced by spiral and bar structure in a set of simulations of Milky Way-sized spiral disc galaxies. The set includes test particle simulations that follow a quasi-stationary density wave-like scenario with rigidly rotating spiral arms, and $N$-body simulations that host a bar and transient, co-rotating spiral arms. From a location similar to that of the Sun, we calculate the radial, tangential and line-of-sight peculiar velocity fields of a patch of the disc and quantify the fluctuations by computing the power spectrum from a two-dimensional Fourier transform. We find that the peculiar velocity power spectrum of the simulation with a bar and transient, co-rotating spiral arms fits very well to that of APOGEE red clump star data, while the quasi-stationary density wave spiral model without a bar does not. We determine that the power spectrum is sensitive to the number of spiral arms, spiral arm pitch angle and position with respect to the spiral arm. However, it is necessary to go beyond the line of sight velocity field in order to distinguish fully between the various spiral models with this method. We compute the power spectrum for different regions of the spiral discs, and discuss the application of this analysis technique to external galaxies.

The Ophiuchus stream is a recently discovered stellar tidal stream in the Milky Way. We present high-quality spectroscopic data for 14 stream member stars obtained using the Keck and MMT telescopes. We confirm the stream as a fast moving ($v_{los}\sim290$ km s$^{-1}$), kinematically cold group ($\sigma_{v_{los}}\lesssim1$ km s$^{-1}$) of $\alpha$-enhanced and metal-poor stars (${\rm [\alpha/Fe]\sim0.4}$ dex, ${\rm [Fe/H]\sim-2.0}$ dex). Using a probabilistic technique, we model the stream simultaneously in line-of-sight velocity, color-magnitude, coordinate, and proper motion space, and so determine its distribution in 6D phase-space. We find that that the stream extends in distance from 7.5 to 9 kpc from the Sun; it is 50 times longer than wide, merely appearing highly foreshortened in projection. The analysis of the stellar population contained in the stream suggests that it is $\sim12$ Gyr old, and that its initial stellar mass was $\sim2\times10^4$ $M_{\odot}$ (or at least $\gtrsim7\times10^3$ $M_{\odot}$). Assuming a fiducial Milky Way potential, we fit an orbit to the stream which matches the observed phase-space distribution, except for some tension in the proper motions: the stream has an orbital period of $\sim350$ Myr, and is on a fairly eccentric orbit ($e\sim0.66$) with a pericenter of $\sim3.5$ kpc and an apocenter of $\sim17$ kpc. The phase-space structure and stellar population of the stream show that its progenitor must have been a globular cluster that was disrupted only $\sim240$ Myr ago. We do not detect any significant overdensity of stars along the stream that would indicate the presence of a progenitor, and conclude that the stream is all that is left of the progenitor.

Extreme deconvolution (XD) of broad-band photometric data can both separate stars from quasars and generate probability density functions for quasar redshifts, while incorporating flux uncertainties and missing data. Mid-infrared photometric colors are now widely used to identify hot dust intrinsic to quasars, and the release of all-sky WISE data has led to a dramatic increase in the number of IR-selected quasars. Using forced-photometry on public WISE data at the locations of SDSS point sources, we incorporate this all-sky data into the training of the XDQSOz models originally developed to select quasars from optical photometry. The combination of WISE and SDSS information is far more powerful than SDSS alone, particularly at $z>2$. The use of SDSS$+$WISE photometry is comparable to the use of SDSS$+$ultraviolet$+$near-IR data. We release a new public catalogue of 5,537,436 (total; 3,874,639 weighted by probability) potential quasars with probability $P_{\textrm{QSO}} > 0.2$. The catalogue includes redshift probabilities for all objects. We also release an updated version of the publicly available set of codes to calculate quasar and redshift probabilities for various combinations of data. Finally, we demonstrate that this method of selecting quasars using WISE data is both more complete and efficient than simple WISE color-cuts, especially at high redshift. Our fits verify that above $z \sim 3$ WISE colors become bluer than the standard cuts applied to select quasars. Currently, the analysis is limited to quasars with optical counterparts, and thus cannot be used to find highly obscured quasars that WISE color-cuts identify in significant numbers.

We derive age constraints for 1639 red giants in the APOKASC sample for which seismic parameters from Kepler, as well as effective temperatures, metallicities and [alpha/Fe] values from APOGEE DR12 are available. We investigate the relation between age and chemical abundances for these stars, using a simple and robust approach to obtain ages. We first derive stellar masses using standard seismic scaling relations, then determine the maximum possible age for each star as function of its mass and metallicity, independently of its evolutionary stage. While the overall trend between maximum age and chemical abundances is a declining fraction of young stars with increasing [alpha/Fe], at least 14 out of 241 stars with [alpha/Fe]>0.13 are younger than 6 Gyr. Five stars with [alpha/Fe]>0.2 have ages below 4 Gyr. We examine the effect of modifications in the standard seismic scaling relations, as well as the effect of very low helium fractions, but these changes are not enough to make these stars as old as usually expected for alpha-rich stars (i.e., ages greater than 8-9 Gyr). Such unusual alpha-rich young stars have also been detected by other surveys, but defy simple explanations in a galaxy evolution context.

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 investigate the light-element behavior of red giant stars in Northern globular clusters (GCs) observed by the SDSS-III Apache Point Observatory Galactic Evolution Experiment (APOGEE). We derive abundances of nine elements (Fe, C, N, O, Mg, Al, Si, Ca, and Ti) for 428 red giant stars in 10 globular clusters. The intrinsic abundance range relative to measurement errors is examined, and the well-known C-N and Mg-Al anticorrelations are explored using an extreme-deconvolution code for the first time in a consistent way. We find that Mg and Al drive the population membership in most clusters, except in M107 and M71, the two most metal-rich clusters in our study, where the grouping is most sensitive to N. We also find a diversity in the abundance distributions, with some clusters exhibiting clear abundance bimodalities (for example M3 and M53) while others show extended distributions. The spread of Al abundances increases significantly as cluster average metallicity decreases as previously found by other works, which we take as evidence that low metallicity, intermediate mass AGB polluters were more common in the more metal poor clusters. The statistically significant correlation of [Al/Fe] with [Si/Fe] in M15 suggests that $^{28}$Si leakage has occurred in this cluster. We also present C, N and O abundances for stars cooler than 4500 K and examine the behavior of A(C+N+O) in each cluster as a function of temperature and [Al/Fe]. The scatter of A(C+N+O) is close to its estimated uncertainty in all clusters and independent on stellar temperature. A(C+N+O) exhibits small correlations and anticorrelations with [Al/Fe] in M3 and M13, but we cannot be certain about these relations given the size of our abundance uncertainties. Star-to-star variations of $\alpha-$elements (Si, Ca, Ti) abundances are comparable to our estimated errors in all clusters.

Using a sample of 69,919 red giants from the SDSS-III/APOGEE Data Release 12, we measure the distribution of stars in the [$\alpha$/Fe] vs. [Fe/H] plane and the metallicity distribution functions (MDF) across an unprecedented volume of the Milky Way disk, with radius $3<R<15$ kpc and height $|z|<2$ kpc. Stars in the inner disk ($R<5$ kpc) lie along a single track in [$\alpha$/Fe] vs. [Fe/H], starting with $\alpha$-enhanced, metal-poor stars and ending at [$\alpha$/Fe]$\sim0$ and [Fe/H]$\sim+0.4$. At larger radii we find two distinct sequences in [$\alpha$/Fe] vs. [Fe/H] space, with a roughly solar-$\alpha$ sequence that spans a decade in metallicity and a high-$\alpha$ sequence that merges with the low-$\alpha$ sequence at super-solar [Fe/H]. The location of the high-$\alpha$ sequence is nearly constant across the disk, however there are very few high-$\alpha$ stars at $R>11$ kpc. The peak of the midplane MDF shifts to lower metallicity at larger $R$, reflecting the Galactic metallicity gradient. Most strikingly, the shape of the midplane MDF changes systematically with radius, with a negatively skewed distribution at $3<R<7$ kpc, to a roughly Gaussian distribution at the solar annulus, to a positively skewed shape in the outer Galaxy. For stars with $|z|>1$ kpc or [$\alpha$/Fe]$>0.18$, the MDF shows little dependence on $R$. The positive skewness of the outer disk MDF may be a signature of radial migration; we show that blurring of stellar populations by orbital eccentricities is not enough to explain the reversal of MDF shape but a simple model of radial migration can do so.

We investigate the kinematics of stars in the mid-plane of the Milky Way on scales between 25 pc and 10 kpc with data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE), the Radial Velocity Experiment (RAVE), and the Geneva-Copenhagen Survey (GCS). Using red-clump stars in APOGEE, we determine the large-scale line-of-sight velocity field out to 5 kpc from the Sun in (0.75 kpc)^2 bins. The solar motion V_{sun-c} with respect to the circular velocity V_c is the largest contribution to the power on large scales after subtracting an axisymmetric rotation field; we determine the solar motion by minimizing the large-scale power to be V_{sun-c} = 24+/-1 (ran.)+/-2 (syst [V_c])+/-5 (syst. [large-scale]) km/s, where the systematic uncertainty is due to (a) a conservative 20 km/s uncertainty in V_c and (b) the estimated power on unobserved larger scales. Combining the APOGEE peculiar-velocity field with red-clump stars in RAVE out to 2 kpc from the Sun and with local GCS stars, we determine the power spectrum of residual velocity fluctuations in the Milky Way's disk on scales between 0.2/kpc < k < 40/kpc. Most of the power is contained in a broad peak between 0.2/kpc < k < 0.9/kpc. We investigate the expected power spectrum for various non-axisymmetric perturbations and demonstrate that the central bar with commonly used parameters but of relatively high mass can explain the bulk of velocity fluctuations in the plane of the Galactic disk near the Sun. Streaming motions ~10 km/s on >~3 kpc scales in the Milky Way are in good agreement with observations of external galaxies and directly explain why local determinations of the solar motion are inconsistent with global measurements.

The SDSS-III/APOGEE survey operated from 2011-2014 using the APOGEE spectrograph, which collects high-resolution (R~22,500), near-IR (1.51-1.70 microns) spectra with a multiplexing (300 fiber-fed objects) capability. We describe the survey data products that are publicly available, which include catalogs with radial velocity, stellar parameters, and 15 elemental abundances for over 150,000 stars, as well as the more than 500,000 spectra from which these quantities are derived. Calibration relations for the stellar parameters (Teff, log g, [M/H], [alpha/M]) and abundances (C, N, O, Na, Mg, Al, Si, S, K, Ca, Ti, V, Mn, Fe, Ni) are presented and discussed. The internal scatter of the abundances within clusters indicates that abundance precision is generally between 0.05 and 0.09 dex across a broad temperature range; within more limited ranges and at high S/N, it is smaller for some elemental abundances. We assess the accuracy of the abundances using comparison of mean cluster metallicities with literature values, APOGEE observations of the solar spectrum and of Arcturus, comparison of individual star abundances with other measurements, and consideration of the locus of derived parameters and abundances of the entire sample, and find that it is challenging to determine the absolute abundance scale; external accuracy may be good to 0.1-0.2 dex. Uncertainties may be larger at cooler temperatures (Teff<4000K). Access to the public data release and data products is described, and some guidance for using the data products is provided.

We explore the practical feasibility of AGN broad-band reverberation mapping and present first results. We lay out and apply a rigorous approach for stochastic reverberation mapping of unevenly sampled multi-broad-band flux measurements, assuming that the broad-line region (BLR) line flux is contributing up to 15 % in some bands, and is directly constrained by one spectroscopical epoch. The approach describes variations of the observed flux as the continuum, modeled as a stochastic Gaussian process, and emission line contribution, modeled as a scaled, smoothed and delayed version of the continuum. This approach is capable not only to interpolate in time between measurements, but also to determine confidence limits on continuum -- line emission delays. This approach is applied to SDSS observations in 'Stripe 82' (S82) providing flux measurements precise to 2 % at $\sim$ 60 epochs over $\sim$ 10 years. The strong annual variations in the epoch sampling prove a serious limitation in practice. Also, suitable redshift ranges must be identified, where strong broad emission line contribute to one filter, but not to another. Through generating and evaluating problem-specific mock data, we verify that S82-like data can constrain $\tau_{\mathrm{delay}}$ for a simple transfer function model. In application to real data, we estimate $\tau_{\mathrm{delay}}$ for 323 AGN with $0.225 < z < 0.846$, combining information for different objects through the ensemble-scaling relationships for BLR size and BH mass. Our analysis tentatively indicates a 1.7 times larger BLR size of H$\alpha$ and MgII compared to Kaspi2000 and Vestergaard2002, but the seasonal data sampling casts doubt on the robustness of the inference.

We report a detection of the baryon acoustic oscillation (BAO) feature in the flux-correlation function of the Ly{\alpha} forest of high-redshift quasars with a statistical significance of five standard deviations. The study uses 137,562 quasars in the redshift range $2.1\le z \le 3.5$ from the Data Release 11 (DR11) of the Baryon Oscillation Spectroscopic Survey (BOSS) of SDSS-III. This sample contains three times the number of quasars used in previous studies. The measured position of the BAO peak determines the angular distance, $D_A(z=2.34)$ and expansion rate, $H(z=2.34)$, both on a scale set by the sound horizon at the drag epoch, $r_d$. We find $D_A/r_d=11.28\pm0.65(1\sigma)^{+2.8}_{-1.2}(2\sigma)$ and $D_H/r_d=9.18\pm0.28(1\sigma)\pm0.6(2\sigma)$ where $D_H=c/H$. The optimal combination, $\sim D_H^{0.7}D_A^{0.3}/r_d$ is determined with a precision of $\sim2\%$. For the value $r_d=147.4~{\rm Mpc}$, consistent with the CMB power spectrum measured by Planck, we find $D_A(z=2.34)=1662\pm96(1\sigma)~{\rm Mpc}$ and $H(z=2.34)=222\pm7(1\sigma)~{\rm km\,s^{-1}Mpc^{-1}}$. Tests with mock catalogs and variations of our analysis procedure have revealed no systematic uncertainties comparable to our statistical errors. Our results agree with the previously reported BAO measurement at the same redshift using the quasar-Ly{\alpha} forest cross-correlation. The auto-correlation and cross-correlation approaches are complementary because of the quite different impact of redshift-space distortion on the two measurements. The combined constraints from the two correlation functions imply values of $D_A/r_d$ and $D_H/r_d$ that are, respectively, 7% low and 7% high compared to the predictions of a flat $\Lambda$CDM cosmological model with the best-fit Planck parameters. With our estimated statistical errors, the significance of this discrepancy is $\approx 2.5\sigma$.

I describe the design, implementation, and usage of galpy, a Python package for galactic-dynamics calculations. At its core, galpy consists of a general framework for representing galactic potentials both in Python and in C (for accelerated computations); galpy functions, objects, and methods can generally take arbitrary combinations of these as arguments. Numerical orbit integration is supported with a variety of Runge-Kutta-type and symplectic integrators. For planar orbits, integration of the phase-space volume is also possible. galpy supports the calculation of action-angle coordinates and orbital frequencies for a given phase-space point for general spherical potentials, using state-of-the-art numerical approximations for axisymmetric potentials, and making use of a recent general approximation for any static potential. A number of different distribution functions (DFs) are also included in the current release; currently these consist of two-dimensional axisymmetric and non-axisymmetric disk DFs, a three-dimensional disk DF, and a DF framework for tidal streams. I provide several examples to illustrate the use of the code. I present a simple model for the Milky Way's gravitational potential consistent with the latest observations. I also numerically calculate the Oort functions for different tracer populations of stars and compare it to a new analytical approximation. Additionally, I characterize the response of a kinematically-warm disk to an elliptical m=2 perturbation in detail. Overall, galpy consists of about 54,000 lines, including 23,000 lines of code in the module, 11,000 lines of test code, and about 20,000 lines of documentation. The test suite covers 99.6% of the code. galpy is available at http://github.com/jobovy/galpy with extensive documentation available at http://galpy.readthedocs.org/en/latest .

We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy & Rix (2013), we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of 20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of 20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio qDM=0.7+/-0.1; this corresponds to an axis ratio qDM=0.4+/-0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior.

I present a new framework for modeling the dynamics of tidal streams. The framework consists of simple models for the initial action-angle distribution of tidal debris, which can be straightforwardly evolved forward in time. Taking advantage of the essentially one-dimensional nature of tidal streams, the transformation to position-velocity coordinates can be linearized and interpolated near a small number of points along the stream, thus allowing for efficient computations of a stream's properties in observable quantities. I illustrate how to calculate the stream's average location (its "track"') in different coordinate systems, how to quickly estimate the dispersion around its track, and how to draw mock stream data. As a generative model, this framework allows one to compute the full probability distribution function and marginalize over or condition it on certain phase-space dimensions as well as convolve it with observational uncertainties. This will be instrumental in proper data analysis of stream data. In addition to providing a computationally-efficient practical tool for modeling the dynamics of tidal streams, the action-angle nature of the framework helps elucidate how the observed width of the stream relates to the velocity dispersion or mass of the progenitor, and how the progenitors of "orphan"' streams could be located. The practical usefulness of the proposed framework crucially depends on the ability to calculate action-angle variables for any orbit in any gravitational potential. A novel method for calculating actions, frequencies, and angles in any static potential using a single orbit integration is described in an Appendix.

We present the first APOKASC catalog of spectroscopic and asteroseismic properties of 1916 red giants observed in the Kepler fields. The spectroscopic parameters provided from the Apache Point Observatory Galactic Evolution Experiment project are complemented with asteroseismic surface gravities, masses, radii, and mean densities determined by members of the Kepler Asteroseismology Science Consortium. We assess both random and systematic sources of error and include a discussion of sample selection for giants in the Kepler fields. Total uncertainties in the main catalog properties are of order 80 K in Teff , 0.06 dex in [M/H], 0.014 dex in log g, and 12% and 5% in mass and radius, respectively; these reflect a combination of systematic and random errors. Asteroseismic surface gravities are substantially more precise and accurate than spectroscopic ones, and we find good agreement between their mean values and the calibrated spectroscopic surface gravities. There are, however, systematic underlying trends with Teff and log g. Our effective temperature scale is between 0-200 K cooler than that expected from the Infrared Flux Method, depending on the adopted extinction map, which provides evidence for a lower value on average than that inferred for the Kepler Input Catalog (KIC). We find a reasonable correspondence between the photometric KIC and spectroscopic APOKASC metallicity scales, with increased dispersion in KIC metallicities as the absolute metal abundance decreases, and offsets in Teff and log g consistent with those derived in the literature. We present mean fitting relations between APOKASC and KIC observables and discuss future prospects, strengths, and limitations of the catalog data.

We present a first determination of distances and extinctions for individual stars in the first release of the APOKASC catalogue, built from the joint efforts of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and the Kepler Asteroseismic Science Consortium (KASC). Our method takes into account the spectroscopic constraints derived from the APOGEE Stellar Parameters and Chemical Abundances Pipeline, together with the asteroseismic parameters from KASC. These parameters are then employed to estimate intrinsic stellar properties, including absolute magnitudes, using the Bayesian tool PARAM. We then find the distance and extinction that best fit the observed photometry in SDSS, 2MASS, and WISE passbands. The first 1989 giants targeted by APOKASC are found at typical distances between 0.5 and 5 kpc, with individual uncertainties of just ~1.8 per cent. Our extinction estimates are systematically smaller than provided in the Kepler Input Catalogue and by the Schlegel, Finkbeiner and Davis maps. Distances to individual stars in the NGC 6791 and NGC 6819 star clusters agree to within their credible intervals. Comparison with the APOGEE red clump and SAGA catalogues provide another useful check, exhibiting agreement with our measurements to within a few percent. Overall, present methods seem to provide excellent distance and extinction determinations for the bulk of the APOKASC sample. Approximately one third of the stars present broad or multiple-peaked probability density functions and hence increased uncertainties. Uncertainties are expected to be reduced in future releases of the catalogue, when a larger fraction of the stars will have seismically-determined evolutionary status classifications.

We employ the first two years of data from the near-infrared, high-resolution SDSS-III/APOGEE spectroscopic survey to investigate the distribution of metallicity and alpha-element abundances of stars over a large part of the Milky Way disk. Using a sample of ~10,000 kinematically-unbiased red-clump stars with ~5% distance accuracy as tracers, the [alpha/Fe] vs. [Fe/H] distribution of this sample exhibits a bimodality in [alpha/Fe] at intermediate metallicities, -0.9<[Fe/H]<-0.2, but at higher metallicities ([Fe/H]=+0.2) the two sequences smoothly merge. We investigate the effects of the APOGEE selection function and volume filling fraction and find that these have little qualitative impact on the alpha-element abundance patterns. The described abundance pattern is found throughout the range 5<R<11 kpc and 0<|Z|<2 kpc across the Galaxy. The [alpha/Fe] trend of the high-alpha sequence is surprisingly constant throughout the Galaxy, with little variation from region to region (~10%). Using simple galactic chemical evolution models we derive an average star formation efficiency (SFE) in the high-alpha sequence of ~4.5E-10 1/yr, which is quite close to the nearly-constant value found in molecular-gas-dominated regions of nearby spirals. This result suggests that the early evolution of the Milky Way disk was characterized by stars that shared a similar star formation history and were formed in a well-mixed, turbulent, and molecular-dominated ISM with a gas consumption timescale (1/SFE) of ~2 Gyr. Finally, while the two alpha-element sequences in the inner Galaxy can be explained by a single chemical evolutionary track this cannot hold in the outer Galaxy, requiring instead a mix of two or more populations with distinct enrichment histories.

The Sloan Digital Sky Survey III's Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution near-infrared spectroscopic survey covering all of the major components of the Galaxy, including the dust-obscured regions of the inner Milky Way disk and bulge. Here we present a sample of 10,341 likely red-clump stars (RC) from the first two years of APOGEE operations, selected based on their position in color-metallicity-surface-gravity-effective-temperature space using a new method calibrated using stellar-evolution models and high-quality asteroseismology data. The narrowness of the RC locus in color-metallicity-luminosity space allows us to assign distances to the stars with an accuracy of 5 to 10%. The sample extends to typical distances of about 3 kpc from the Sun, with some stars out to 8 kpc, and spans a volume of approximately 100 kpc^3 over 5 kpc <~ R <~ 14 kpc, |Z| <~ 2 kpc, and -15 deg <~ Galactocentric azimuth <~ 30 deg. The APOGEE red-clump (APOGEE-RC) catalog contains photometry from 2MASS, reddening estimates, distances, line-of-sight velocities, stellar parameters and elemental abundances determined from the high-resolution APOGEE spectra, and matches to major proper motion catalogs. We determine the survey selection function for this data set and discuss how the RC selection samples the underlying stellar populations. We use this sample to limit any azimuthal variations in the median metallicity within the ~45 degree-wide azimuthal region covered by the current sample to be <= 0.02 dex, which is more than an order of magnitude smaller than the radial metallicity gradient. This result constrains coherent non-axisymmetric flows within a few kpc from the Sun.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

We present and apply rigorous dynamical modeling with which we infer unprecedented constraints on the stellar and dark matter mass distribution within our Milky Way (MW), based on large sets of phase-space data on individual stars. Specifically, we model the dynamics of 16,269 G-type dwarfs from SEGUE, which sample 5 < R_GC/kpc < 12 and 0.3 < |Z|/kpc < 3. We independently fit a parameterized MW potential and a three-integral, action-based distribution function (DF) to the phase-space data of 43 separate abundance-selected sub-populations (MAPs), accounting for the complex selection effects affecting the data. We robustly measure the total surface density within 1.1 kpc of the mid-plane to 5% over 4.5 < R_GC/kpc < 9. Using metal-poor MAPs with small radial scale lengths as dynamical tracers probes 4.5 < R_GC/kpc < 7, while MAPs with longer radial scale lengths sample 7 < R_GC/kpc < 9. We measure the mass-weighted Galactic disk scale length to be R_d = 2.15+/-0.14 kpc, in agreement with the photometrically inferred spatial distribution of stellar mass. We thereby measure dynamically the mass of the Galactic stellar disk to unprecedented accuracy: M_* = 4.6+/-0.3+3.0x(R_0/kpc-8)x10^{10}Msun and a total local surface density of \Sigma_{R_0}(Z=1.1 kpc) = 68+/-4 Msun/pc^2 of which 38+/-4 Msun/pc^2 is contributed by stars and stellar remnants. By combining our surface density measurements with the terminal velocity curve, we find that the MW's disk is maximal in that V_{c,disk} / V_{c,total} = 0.83+/-0.04 at R=2.2 R_d. We also constrain for the first time the radial profile of the dark halo at such small Galactocentric radii, finding that \rho_{DM} (r;near R_0) \propto 1 / r^\alpha with \alpha < 1.53 at 95% confidence. Our results show that action-based distribution-function modeling of complex stellar data sets is now a feasible approach that will be fruitful for interpreting Gaia data.

We present the Data Release 10 Quasar (DR10Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the first 2.5 years of the survey and that are confirmed as quasars via visual inspection of the spectra. The catalog also includes known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 166,583 quasars (74,454 are new discoveries since SDSS-DR9) detected over 6,373 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with $z>2.15$ (117,668) is $\sim$5 times greater than the number of $z>2.15$ quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII, MgII). The catalog identifies 16,461 broad absorption line quasars and gives their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry with typical accuracy of 0.03 mag and information on the optical morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3,600-10,500\AA\ at a spectral resolution in the range 1,300$<$R$<$2,500; the spectra can be retrieved from the SDSS Catalog Archive Server. We also provide a supplemental list of an additional 2,376 quasars that have been identified among the galaxy targets of the SDSS-III/BOSS.

We present Galactic mean metallicity maps derived from the first year of the SDSS-III APOGEE experiment. Mean abundances in different zones of Galactocentric radius (0 < R < 15 kpc) at a range of heights above the plane (0 < |z| < 3 kpc), are derived from a sample of nearly 20,000 stars with unprecedented coverage, including stars in the Galactic mid-plane at large distances. We also split the sample into subsamples of stars with low and high-[{\alpha}/M] abundance ratios. We assess possible biases in deriving the mean abundances, and find they are likely to be small except in the inner regions of the Galaxy. A negative radial gradient exists over much of the Galaxy; however, the gradient appears to flatten for R < 6 kpc, in particular near the Galactic mid-plane and for low-[{\alpha}/M] stars. At R > 6 kpc, the gradient flattens as one moves off of the plane, and is flatter at all heights for high-[{\alpha}/M] stars than for low-[{\alpha}/M] stars. Alternatively, these gradients can be described as vertical gradients that flatten at larger Galactocentric radius; these vertical gradients are similar for both low and high-[{\alpha}/M] populations. Stars with higher [{\alpha}/M] appear to have a flatter radial gradient than stars with lower [{\alpha}/M]. This could suggest that the metallicity gradient has grown steeper with time or, alternatively, that gradients are washed out over time by migration of stars.

We analyse the structure and chemical enrichment of a Milky Way-like galaxy with a stellar mass of 2 10^{10} M_sun, formed in a cosmological hydrodynamical simulation. It is disk-dominated with a flat rotation curve, and has a disk scale length similar to the Milky Way's, but a velocity dispersion that is ~50% higher. Examining stars in narrow [Fe/H] and [\alpha/Fe] abundance ranges, we find remarkable qualitative agreement between this simulation and observations: a) The old stars lie in a thickened distribution with a short scale length, while the young stars form a thinner disk, with scale lengths decreasing, as [Fe/H] increases. b) Consequently, there is a distinct outward metallicity gradient. c) Mono-abundance populations exist with a continuous distribution of scale heights (from thin to thick). However, the simulated galaxy has a distinct and substantive very thick disk (h_z~1.5 kpc), not seen in the Milky Way. The broad agreement between simulations and observations allows us to test the validity of observational proxies used in the literature: we find in the simulation that mono-abundance populations are good proxies for single age populations (<1 Gyr) for most abundances.

With close pairs of quasars at different redshifts, a background quasar sightline can be used to study a foreground quasar's environment in absorption. We use a sample of 650 projected quasar pairs to study the HI Lya absorption transverse to luminous, z~2 quasars at proper separations of 30kpc < R < 1Mpc. In contrast to measurements along the line-of-sight, regions transverse to quasars exhibit enhanced HI Lya absorption and a larger variance than the ambient intergalactic medium, with increasing absorption and variance toward smaller scales. Analysis of composite spectra reveals excess absorption characterized by a Lya equivalent width profile W = 2.3A (R/100kpc)^-0.46. We also observe a high (~60%) covering factor of strong, optically thick HI absorbers (HI column log NHI > 17.3) at separations R<200kpc, which decreases to ~20% at R~1Mpc, but still represents a significant excess over the cosmic average. This excess of optically thick absorption can be described by a quasar-absorber cross-correlation function xi_QA(r) = (r/r_0)^gamma with a large correlation length r_0 = 12.5+2.7-1.4 Mpc/h (comoving) and gamma = 1.68+0.14-0.30. The HI absorption measured around quasars exceeds that of any previously studied population, consistent with quasars being hosted by massive dark matter halos Mhalo~10^12.5 Msun at z~2.5. The environments of these massive halos are highly biased towards producing optically thick gas, and may even dominate the cosmic abundance of Lyman limit systems and hence the intergalactic opacity to ionizing photons at z~2.5. The anisotropic absorption around quasars implies the transverse direction is much less likely to be illuminated by ionizing radiation than the line-of-sight, which we interpret in terms of the same obscuration effects frequently invoked in unified models of active galactic nuclei.

We present a rigorous and practical way of constraining the Galactic potential based on the phase-space information for many individual stars. Such an approach is needed to dynamically model the data from ongoing spectroscopic surveys of the Galaxy and in the future Gaia. This approach describes the orbit distribution of stars by a family of parametrized distribution function (DF) proposed by McMillan and Binney, which are based on actions. We find that these parametrized DFs are flexible enough to capture well the observed phase-space distributions of individual abundance-selected Galactic subpopulations of stars (`mono-abundance populations') for a disc-like gravitational potential, which enables independent dynamical constraints from each of the Galactic mono-abundance populations. We lay out a statistically rigorous way to constrain the Galactic potential parameters by constructing the joint likelihood of potential and DF parameters, and subsequently marginalizing over the DF parameters. This approach explicitly incorporates the spatial selection function inherent to all Galactic surveys, and can account for the uncertainties of the individual position--velocity observations. On that basis, we study the precision of the parameters of the Galactic potential that can be reached with various sample sizes and realistic spatial selection functions. By creating mock samples from the DF, we show that, even under a restrictive and realistic spatial selection function, given a two-parameter gravitational potential, one can recover the true potential parameters to a few per cent with sample sizes of a few thousands. The assumptions of axisymmetry, of DFs that are smooth in the actions and of no time variation remain important limitations in our current study.

To constrain the Galactic gravitational potential near the Sun ($\sim$1.5 kpc), we derive and model the spatial and velocity distribution for a sample of 9000 K-dwarfs that have spectra from SDSS/SEGUE, which yield radial velocities and abundances ([Fe/H] & [$\alpha$/Fe]). We first derive the spatial density distribution for stars of three abundance-selected sub-populations by accounting for the survey's selection function. The vertical profile of these sub-populations are simple exponentials and their vertical dispersion profile is nearly isothermal. To model these data, we apply the `vertical' Jeans Equation, which relates the observable tracer number density and vertical velocity dispersion to the gravitational potential or vertical force. We explore a number of functional forms for the vertical force law, and fit the dispersion and density profiles of all abundance selected sub-populations simultaneously in the same potential, and explore all parameter co-variances using MCMC. Our fits constrain a disk {\it mass} scale height $\lesssim$ 300 pc and the total surface mass density to be $67 \pm 6 M_{\odot} {\rm pc^{-2}}$ at $|z| = 1.0$ kpc of which the contribution from all stars is $42 \pm 5 M_{\odot} {\rm pc^{-2}}$ (presuming a contribution from cold gas of $13 M_{\odot} {\rm pc^{-2}}$). We find significant constraints on the local dark matter density of $0.0065\pm0.0023 M_{\odot} {\rm pc^{-3}}$ ($0.25\pm0.09 {\rm GeV cm^{-3}} $). Together with recent experiments this firms up the best estimate of $0.0075\pm0.0021 M_{\odot} {\rm pc^{-3}}$ ($0.28\pm0.08 {\rm GeV cm^{-3}} $), consistent with global fits of approximately round dark matter halos to kinematic data in the outskirts of the Galaxy.

We use the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) to detect and measure the position of the Baryonic Acoustic Oscillation (BAO) feature in the three-dimensional correlation function in the Lyman-alpha forest flux fluctuations at a redshift z=2.4. The feature is clearly detected at significance between 3 and 5 sigma (depending on the broadband model and method of error covariance matrix estimation) and is consistent with predictions of the standard LCDM model. We assess the biases in our method, stability of the error covariance matrix and possible systematic effects. We fit the resulting correlation function with several models that decouple the broadband and acoustic scale information. For an isotropic dilation factor, we measure 100x(alpha_iso-1) = -1.6 ^{+2.0+4.3+7.4}_{-2.0-4.1-6.8} (stat.) +/- 1.0 (syst.) (multiple statistical errors denote 1,2 and 3 sigma confidence limits) with respect to the acoustic scale in the fiducial cosmological model (flat LCDM with Omega_m=0.27, h=0.7). When fitting separately for the radial and transversal dilation factors we find marginalised constraints 100x(alpha_par-1) = -1.3 ^{+3.5+7.6 +12.3}_{-3.3-6.7-10.2} (stat.) +/- 2.0 (syst.) and 100x(alpha_perp-1) = -2.2 ^{+7.4+17}_{-7.1-15} +/- 3.0 (syst.). The dilation factor measurements are significantly correlated with cross-correlation coefficient of ~ -0.55. Errors become significantly non-Gaussian for deviations over 3 standard deviations from best fit value. Because of the data cuts and analysis method, these measurements give tighter constraints than a previous BAO analysis of the BOSS DR9 Lyman-alpha forest sample, providing an important consistency test of the standard cosmological model in a new redshift regime.

We report a detection of the baryon acoustic oscillation (BAO) feature in the three-dimensional correlation function of the transmitted flux fraction in the \Lya forest of high-redshift quasars. The study uses 48,640 quasars in the redshift range $2.1\le z \le 3.5$ from the Baryon Oscillation Spectroscopic Survey (BOSS) of the third generation of the Sloan Digital Sky Survey (SDSS-III). At a mean redshift $z=2.3$, we measure the monopole and quadrupole components of the correlation function for separations in the range $20\hMpc<r<200\hMpc$. A peak in the correlation function is seen at a separation equal to $(1.01\pm0.03)$ times the distance expected for the BAO peak within a concordance $\Lambda$CDM cosmology. This first detection of the BAO peak at high redshift, when the universe was strongly matter dominated, results in constraints on the angular diameter distance $\da$ and the expansion rate $H$ at $z=2.3$ that, combined with priors on $H_0$ and the baryon density, require the existence of dark energy. Combined with constraints derived from Cosmic Microwave Background (CMB) observations, this result implies $H(z=2.3)=(224\pm8){\rm km\,s^{-1}Mpc^{-1}}$, indicating that the time derivative of the cosmological scale parameter $\dot{a}=H(z=2.3)/(1+z)$ is significantly greater than that measured with BAO at $z\sim0.5$. This demonstrates that the expansion was decelerating in the range $0.7<z<2.3$, as expected from the matter domination during this epoch. Combined with measurements of $H_0$, one sees the pattern of deceleration followed by acceleration characteristic of a dark-energy dominated universe.

A suite of vast stellar surveys mapping the Milky Way, culminating in the Gaia mission, is revolutionizing the empirical information about the distribution and properties of stars in the Galactic stellar disk. We review and lay out what analysis and modeling machinery needs to be in place to test mechanisms of disk galaxy evolution and to stringently constrain the Galactic gravitational potential, using such Galactic star-by-star measurements. We stress the crucial role of stellar survey selection functions in any such modeling; and we advocate the utility of viewing the Galactic stellar disk as made up from `mono-abundance populations' (MAPs), both for dynamical modeling and for constraining the Milky Way's evolutionary processes. We review recent work on the spatial and kinematical distribution of MAPs, and lay out how further study of MAPs in the Gaia era should lead to a decisively clearer picture of the Milky Way's dark matter distribution and formation history.

We explore the stellar metallicity distribution function of the Galactic halo based on SDSS ugriz photometry. A set of stellar isochrones is calibrated using observations of several star clusters and validated by comparisons with medium-resolution spectroscopic values over a wide range of metal abundance. We estimate distances and metallicities for individual main-sequence stars in the multiply scanned SDSS Stripe 82, at heliocentric distances in the range 5 - 8 kpc and |b| > 35 deg, and find that the in situ photometric metallicity distribution has a shape that matches that of the kinematically-selected local halo stars from Ryan & Norris. We also examine independent kinematic information from proper-motion measurements for high Galactic latitude stars in our sample. We find that stars with retrograde rotation in the rest frame of the Galaxy are generally more metal poor than those exhibiting prograde rotation, which is consistent with earlier arguments by Carollo et al. that the halo system comprises at least two spatially overlapping components with differing metallicity, kinematics, and spatial distributions. The observed photometric metallicity distribution and that of Ryan & Norris can be described by a simple chemical evolution model by Hartwick (or by a single Gaussian distribution); however, the suggestive metallicity-kinematic correlation contradicts the basic assumption in this model that the Milky Way halo consists primarily of a single stellar population. When the observed metallicity distribution is deconvolved using two Gaussian components with peaks at [Fe/H] ~ -1.7 and -2.3, the metal-poor component accounts for ~20% - 35% of the entire halo population in this distance range.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $\Lambda$CDM cosmology with H0 = 70 km/s/Mpc, $\Omega_{\rm M}$ = 0.3, and $\Omega_{\Lambda}$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

We measure the Milky Way's rotation curve over the Galactocentric range 4 kpc <~ R <~ 14 kpc from the first year of data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE). We model the line-of-sight velocities of 3,365 stars in fourteen fields with b = 0 deg between 30 deg < l < 210 deg out to distances of 10 kpc using an axisymmetric kinematical model that includes a correction for the asymmetric drift of the warm tracer population (\sigma_R ~ 35 km/s). We determine the local value of the circular velocity to be V_c(R_0) = 218 +/- 6 km/s and find that the rotation curve is approximately flat with a local derivative between -3.0 km/s/kpc and 0.4 km/s/kpc. We also measure the Sun's position and velocity in the Galactocentric rest frame, finding the distance to the Galactic center to be 8 kpc < R_0 < 9 kpc, radial velocity V_{R,sun} = -10 +/- 1 km/s, and rotational velocity V_{\phi,sun} = 242^{+10}_{-3} km/s, in good agreement with local measurements of the Sun's radial velocity and with the observed proper motion of Sgr A*. We investigate various systematic uncertainties and find that these are limited to offsets at the percent level, ~2 km/s in V_c. Marginalizing over all the systematics that we consider, we find that V_c(R_0) < 235 km/s at >99% confidence. We find an offset between the Sun's rotational velocity and the local circular velocity of 26 +/- 3 km/s, which is larger than the locally-measured solar motion of 12 km/s. This larger offset reconciles our value for V_c with recent claims that V_c >~ 240 km/s. Combining our results with other data, we find that the Milky Way's dark-halo mass within the virial radius is ~8x10^{11} M_sun.

We present the vertical kinematics of stars in the Milky Way's stellar disk inferred from SDSS/SEGUE G-dwarf data, deriving the vertical velocity dispersion, \sigma_z, as a function of vertical height |z| and Galactocentric radius R for a set of 'mono-abundance' sub-populations of stars with very similar elemental abundances [\alpha/Fe] and [Fe/H]. We find that all components exhibit nearly isothermal kinematics in |z|, and a slow outward decrease of the vertical velocity dispersion: \sigma_z (z,R|[\alpha/Fe],[Fe/H]) ~ \sigma_z ([\alpha/Fe],[Fe/H]) x \exp (-(R-R_0)/7 kpc}). The characteristic velocity dispersions of these components vary from ~ 15 km/s for chemically young, metal-rich stars, to >~ 50 km/s for metal poor stars. The mean \sigma_z gradient away from the mid plane is only 0.3 +/- 0.2 km/s/kpc. We find a continuum of vertical kinetic temperatures (~\sigma^2_z) as function of ([\alpha/Fe],[Fe/H]), which contribute to the stellar surface mass density as \Sigma_{R_0}(\sigma^2_z) ~ \exp(-\sigma^2_z). The existence of isothermal mono-abundance populations with intermediate dispersions reject the notion of a thin-thick disk dichotomy. This continuum of disks argues against models where the thicker disk portions arise from massive satellite infall or heating; scenarios where either the oldest disk portion was born hot, or where internal evolution plays a major role, seem the most viable. The wide range of \sigma_z ([\alpha/Fe],[Fe/H]) combined with a constant \sigma_z(z) for each abundance bin provides an independent check on the precision of the SEGUE abundances: \delta_[\alpha/Fe] ~ 0.07 dex and \delta_[Fe/H] ~ 0.15 dex. The radial decline of the vertical dispersion presumably reflects the decrease in disk surface-mass density. This measurement constitutes a first step toward a purely dynamical estimate of the mass profile the disk in our Galaxy. [abridged]

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

An analysis of the kinematics of 412 stars at 1-4 kpc from the Galactic mid-plane by Moni Bidin et al. (2012) has claimed to derive a local density of dark matter that is an order of magnitude below standard expectations. We show that this result is incorrect and that it arises from the assumption that the mean azimuthal velocity of the stellar tracers is independent of Galactocentric radius at all heights. We substitute the assumption, supported by data, that the circular speed is independent of radius in the mid-plane. We demonstrate that the assumption of constant mean azimuthal velocity is implausible by showing that it requires the circular velocity to drop more steeply than allowed by any plausible mass model, with or without dark matter, at large heights above the mid-plane. Using the approximation that the circular velocity curve is flat in the mid-plane, we find that the data imply a local dark-matter density of 0.008 +/- 0.003 Msun/pc^3 = 0.3 +/- 0.1 GeV/cm3, fully consistent with standard estimates of this quantity. This is the most robust direct measurement of the local dark-matter density to date.

We measure the quasar two-point correlation function over the redshift range 2.2<z<2.8 using data from the Baryon Oscillation Spectroscopic Survey. We use a homogeneous subset of the data consisting of 27,129 quasars with spectroscopic redshifts---by far the largest such sample used for clustering measurements at these redshifts to date. The sample covers 3,600 square degrees, corresponding to a comoving volume of 9.7(Gpc/h)^3 assuming a fiducial LambdaCDM cosmology, and it has a median absolute i-band magnitude of -26, k-corrected to z=2. After accounting for redshift errors we find that the redshift space correlation function is fit well by a power-law of slope -2 and amplitude s_0=(9.7\pm 0.5)Mpc/h over the range 3<s<25Mpc/h. The projected correlation function, which integrates out the effects of peculiar velocities and redshift errors, is fit well by a power-law of slope -1 and r_0=(8.4\pm 0.6)Mpc/h over the range 4<R<16Mpc/h. There is no evidence for strong luminosity or redshift dependence to the clustering amplitude, in part because of the limited dynamic range in our sample. Our results are consistent with, but more precise than, previous measurements at similar redshifts. Our measurement of the quasar clustering amplitude implies a bias factor of b~3.5 for our quasar sample. We compare the data to models to constrain the manner in which quasars occupy dark matter halos at z~2.4 and infer that such quasars inhabit halos with a characteristic mass of <M>~10^{12}Msun/h with a duty cycle for the quasar activity of 1 per cent.

The spatial, kinematic, and elemental-abundance structure of the Milky Way's stellar disk is complex, and has been difficult to dissect with local spectroscopic or global photometric data. Here, we develop and apply a rigorous density modeling approach for Galactic spectroscopic surveys that enables investigation of the global spatial structure of stellar sub-populations in narrow bins of [\alpha/Fe] and [Fe/H], using 23,767 G-type dwarfs from SDSS/SEGUE. We fit models for the number density of each such mono-abundance component, properly accounting for the complex spectroscopic SEGUE sampling of the underlying stellar population. We find that each mono-abundance sub-population has a simple spatial structure that can be described by a single exponential in both the vertical and radial direction, with continuously increasing scale heights (~200 pc to 1 kpc) and decreasing scale lengths (>4.5 kpc to 2 kpc) for increasingly older sub-populations, as indicated by their lower metallicities and [\alpha/Fe] enhancements. That the abundance-selected sub-components with the largest scale heights have the shortest scale lengths is in sharp contrast with purely geometric `thick--thin disk' decompositions. To the extent that [\alpha/Fe] is an adequate proxy for age, our results directly show that older disk sub-populations are more centrally concentrated, which implies inside-out formation of galactic disks. The fact that the largest scale-height sub-components are most centrally concentrated in the Milky Way is an almost inevitable consequence of explaining the vertical structure of the disk through internal evolution. Whether the simple spatial structure of the mono-abundance sub-components, and the striking correlations between age, scale length, and scale height can be plausibly explained by satellite accretion or other external heating remains to be seen.

We study the effect of the supersonic baryon--CDM flow, which has recently been shown to have a large effect on structure formation during the dark ages 10 <~ z <~ 1000, on the abundance of luminous, low-mass satellite galaxies around galaxies like the Milky Way. As the supersonic baryon--CDM flow significantly suppresses both the number of halos formed and the amount of baryons accreted onto such halos of masses 10^6 < M_{halo} / M_solar < 10^8 at z >~ 10, a large effect results on the stellar luminosity function before reionization. As halos of these masses are believed to have very little star formation after reionization due to the effects of photo-heating by the ultraviolet background, this effect persists to the present day. We calculate that the number of low-mass 10^6 < M_{halo} / M_solar < 10^8 halos that host luminous satellite galaxies today is typically suppressed by 50 percent, with values ranging up to 90 percent in regions where the initial supersonic velocity is high. We show that this previously-ignored cosmological effect resolves most of the tension between the observed and predicted number of low-mass satellites in the Milky Way, obviating the need for any other mass-dependent star-formation suppression before reionization.

Different stellar sub-populations of the Milky Way's stellar disk are known to have different vertical scale heights, their thickness increasing with age. Using SEGUE spectroscopic survey data, we have recently shown that mono-abundance sub-populations, defined in the [\alpha/Fe]-[Fe/H] space, are well described by single exponential spatial-density profiles in both the radial and the vertical direction; therefore any star of a given abundance is clearly associated with a sub-population of scale height h_z. Here, we work out how to determine the stellar surface-mass density contributions at the solar radius R_0 of each such sub-population, accounting for the survey selection function, and for the fraction of the stellar population mass that is reflected in the spectroscopic target stars given populations of different abundances and their presumed age distributions. Taken together, this enables us to derive \Sigma_{R_0}(h_z), the surface-mass contributions of stellar populations with scale height h_z. Surprisingly, we find no hint of a thin-thick disk bi-modality in this mass-weighted scale-height distribution, but a smoothly decreasing function, approximately \Sigma_{R_0}(h_z)\propto \exp(-h_z), from h_z ~ 200 pc to h_z ~ 1 kpc. As h_z is ultimately the structurally defining property of a thin or thick disk, this shows clearly that the Milky Way has a continuous and monotonic distribution of disk thicknesses: there is no 'thick disk' sensibly characterized as a distinct component. We discuss how our result is consistent with evidence for seeming bi-modality in purely geometric disk decompositions, or chemical abundances analyses. We constrain the total visible stellar surface-mass density at the Solar radius to be \Sigma^*_{R_0} = 30 +/- 1 M_\odot pc^{-2}.

We describe a technique for simultaneously classifying and estimating the redshift of quasars. It can separate quasars from stars in arbitrary redshift ranges, estimate full posterior distribution functions for the redshift, and naturally incorporate flux uncertainties, missing data, and multi-wavelength photometry. We build models of quasars in flux-redshift space by applying the extreme deconvolution technique to estimate the underlying density. By integrating this density over redshift one can obtain quasar flux-densities in different redshift ranges. This approach allows for efficient, consistent, and fast classification and photometric redshift estimation. This is achieved by combining the speed obtained by choosing simple analytical forms as the basis of our density model with the flexibility of non-parametric models through the use of many simple components with many parameters. We show that this technique is competitive with the best photometric quasar classification techniques---which are limited to fixed, broad redshift ranges and high signal-to-noise ratio data---and with the best photometric redshift techniques when applied to broadband optical data. We demonstrate that the inclusion of UV and NIR data significantly improves photometric quasar--star separation and essentially resolves all of the redshift degeneracies for quasars inherent to the ugriz filter system, even when included data have a low signal-to-noise ratio. For quasars spectroscopically confirmed by the SDSS 84 and 97 percent of the objects with GALEX UV and UKIDSS NIR data have photometric redshifts within 0.1 and 0.3, respectively, of the spectroscopic redshift; this amounts to about a factor of three improvement over ugriz-only photometric redshifts. Our code to calculate quasar probabilities and redshift probability distributions is publicly available.

We quantify quasar color-variability using an unprecedented variability database - ugriz photometry of 9093 quasars from SDSS Stripe 82, observed over 8 years at ~60 epochs each. We confirm previous reports that quasars become bluer when brightening. We find a redshift dependence of this blueing in a given set of bands (e.g. g and r), but show that it is the result of the flux contribution from less-variable or delayed emission lines in the different SDSS bands at different redshifts. After correcting for this effect, quasar color-variability is remarkably uniform, and independent not only of redshift, but also of quasar luminosity and black hole mass. The color variations of individual quasars, as they vary in brightness on year timescales, are much more pronounced than the ranges in color seen in samples of quasars across many orders of magnitude in luminosity. This indicates distinct physical mechanisms behind quasar variability and the observed range of quasar luminosities at a given black hole mass - quasar variations cannot be explained by changes in the mean accretion rate. We do find some dependence of the color variability on the characteristics of the flux variations themselves, with fast, low-amplitude, brightness variations producing more color variability. The observed behavior could arise if quasar variability results from flares or ephemeral hot spots in an accretion disc.

(Abridged) Carbon-enhanced metal-poor (CEMP) stars in the halo components of the Milky Way are explored, based on accurate determinations of the carbon-to-iron ([C/Fe]) abundance ratios and kinematic quantities for over 30000 calibration stars from the Sloan Digital Sky Survey (SDSS). Using our present criterion that low-metallicity stars exhibiting [C/Fe] ratios ("carbonicity") in excess of [C/Fe]$ = +0.7$ are considered CEMP stars, the global frequency of CEMP stars in the halo system for \feh\ $< -1.5$ is 8%; for \feh\ $< -2.0$ it is 12%; for \feh\ $<-2.5$ it is 20%. We also confirm a significant increase in the level of carbon enrichment with declining metallicity, growing from $<$[C/Fe]$>$ $\sim +1.0$ at \feh\ $= -1.5$ to $<$[C/Fe]$>$ $\sim +1.7$ at \feh\ $= -2.7$. The nature of the carbonicity distribution function (CarDF) changes dramatically with increasing distance above the Galactic plane, $|$Z$|$. For $|$Z$|$ $< 5$ kpc, relatively few CEMP stars are identified. For distances $|$Z$|$ $> 5$ kpc, the CarDF exhibits a strong tail towards high values, up to [C/Fe] $>$ +3.0. We also find a clear increase in the CEMP frequency with $|$Z$|$. For stars with $-2.0 <$ [Fe/H] $< -$1.5, the frequency grows from 5% at $|$Z$|$ $\sim 2$ kpc to 10% at $|$Z$|$ $\sim 10$ kpc. For stars with [Fe/H] $< -$2.0, the frequency grows from 8% at $|$Z$|$ $\sim 2$ kpc to 25% at $|$Z$|$ $\sim 10$ kpc. For stars with $-2.0 <$ [Fe/H] $< -$1.5, the mean carbonicity is $<$[C/Fe]$>$ $\sim +1.0$ for 0 kpc $<$ $|$Z$|$ $<$ 10 kpc, with little dependence on $|$Z$|$; for [Fe/H] $< -$2.0, $<$[C/Fe]$>$ $\sim +1.5$, again roughly independent of $|$Z$|$.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

We argue that the decomposition of gamma-ray maps in spherical harmonics is a sensitive tool to study dark matter (DM) annihilation or decay in the main Galactic halo of the Milky Way. Using the spherical harmonic decomposition in a window excluding the Galactic plane, we show for one year of Fermi data that adding a spherical template (such as a line-of-sight DM annihilation profile) to an astrophysical background significantly reduces chi^2 of the fit to the data. In some energy bins the significance of this DM-like fraction is above three sigma. This can be viewed as a hint of DM annihilation signal, although astrophysical sources cannot be ruled out at this moment. We use the derived DM fraction as a conservative upper limit on DM annihilation signal. In the case of bb-bar annihilation channel the limits are about a factor of two less constraining than the limits from dwarf galaxies. The uncertainty of our method is dominated by systematics related to modeling the astrophysical background. We show that with one year of Fermi data the statistical sensitivity would be sufficient to detect DM annihilation with thermal freeze out cross section for masses below 100 GeV.

We generalize the well-known mixtures of Gaussians approach to density estimation and the accompanying Expectation--Maximization technique for finding the maximum likelihood parameters of the mixture to the case where each data point carries an individual $d$-dimensional uncertainty covariance and has unique missing data properties. This algorithm reconstructs the error-deconvolved or "underlying" distribution function common to all samples, even when the individual data points are samples from different distributions, obtained by convolving the underlying distribution with the heteroskedastic uncertainty distribution of the data point and projecting out the missing data directions. We show how this basic algorithm can be extended with conjugate priors on all of the model parameters and a "split-and-merge" procedure designed to avoid local maxima of the likelihood. We demonstrate the full method by applying it to the problem of inferring the three-dimensional velocity distribution of stars near the Sun from noisy two-dimensional, transverse velocity measurements from the Hipparcos satellite.

The SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), a five-year spectroscopic survey of 10,000 deg^2, achieved first light in late 2009. One of the key goals of BOSS is to measure the signature of baryon acoustic oscillations in the distribution of Ly-alpha absorption from the spectra of a sample of ~150,000 z>2.2 quasars. Along with measuring the angular diameter distance at z\approx2.5, BOSS will provide the first direct measurement of the expansion rate of the Universe at z > 2. One of the biggest challenges in achieving this goal is an efficient target selection algorithm for quasars over 2.2 < z < 3.5, where their colors overlap those of stars. During the first year of the BOSS survey, quasar target selection methods were developed and tested to meet the requirement of delivering at least 15 quasars deg^-2 in this redshift range, out of 40 targets deg^-2. To achieve these surface densities, the magnitude limit of the quasar targets was set at g <= 22.0 or r<=21.85. While detection of the BAO signature in the Ly-alpha absorption in quasar spectra does not require a uniform target selection, many other astrophysical studies do. We therefore defined a uniformly-selected subsample of 20 targets deg^-2, for which the selection efficiency is just over 50%. This "CORE" subsample will be fixed for Years Two through Five of the survey. In this paper we describe the evolution and implementation of the BOSS quasar target selection algorithms during the first two years of BOSS operations. We analyze the spectra obtained during the first year. 11,263 new z>2.2 quasars were spectroscopically confirmed by BOSS. Our current algorithms select an average of 15 z > 2.2 quasars deg^-2 from 40 targets deg^-2 using single-epoch SDSS imaging. Multi-epoch optical data and data at other wavelengths can further improve the efficiency and completeness of BOSS quasar target selection. [Abridged]

We present the SDSS-XDQSO quasar targeting catalog for efficient flux-based quasar target selection down to the faint limit of the Sloan Digital Sky Survey (SDSS) catalog, even at medium redshifts (2.5 <~ z <~ 3) where the stellar contamination is significant. We build models of the distributions of stars and quasars in flux space down to the flux limit by applying the extreme-deconvolution method to estimate the underlying density. We convolve this density with the flux uncertainties when evaluating the probability that an object is a quasar. This approach results in a targeting algorithm that is more principled, more efficient, and faster than other similar methods. We apply the algorithm to derive low-redshift (z < 2.2), medium-redshift (2.2 <= z <= 3.5), and high-redshift (z > 3.5) quasar probabilities for all 160,904,060 point sources with dereddened i-band magnitude between 17.75 and 22.45 mag in the 14,555 deg^2 of imaging from SDSS Data Release 8. The catalog can be used to define a uniformly selected and efficient low- or medium-redshift quasar survey, such as that needed for the SDSS-III's Baryon Oscillation Spectroscopic Survey project. We show that the XDQSO technique performs as well as the current best photometric quasar-selection technique at low redshift, and outperforms all other flux-based methods for selecting the medium-redshift quasars of our primary interest. We make code to reproduce the XDQSO quasar target selection publicly available.

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

It has been proposed that the Hercules stream, a group of co-moving stars in the Solar neighborhood offset from the bulk of the velocity distribution, is the result of resonant interactions between stars in the outer disk and the Galactic bar. So far it has only been seen in the immediate Solar neighborhood, but the resonance model makes a prediction over a large fraction of the Galactic disk. I predict the distribution of stellar velocities and the changing Hercules feature in this distribution as a function of location in the Galactic disk in a simple model for the Galaxy and the bar that produces the observed Hercules stream. The Hercules feature is expected to be strong enough to be unambiguously detected in the distribution of line-of-sight velocities in selected directions. I identify quantitatively the most promising lines of sight for detection in line-of-sight velocities using the Kullback-Leibler divergence between the predictions of the resonance model and an axisymmetric model; these directions are at 250 deg <~ l <~ 290 deg. The predictions presented here are only weakly affected by distance uncertainties, assumptions about the distribution function in the stellar disk, and the details of the Galactic potential including the effect of spiral structure. Gaia and future spectroscopic surveys of the Galactic disk such as APOGEE and HERMES will be able to robustly test the origin of the Hercules stream and constrain the properties of the Galactic bar.

Standard maximum-likelihood estimators for binary-star and exoplanet eccentricities are biased high, in the sense that the estimated eccentricity tends to be larger than the true eccentricity. As with most non-trivial observables, a simple histogram of estimated eccentricities is not a good estimate of the true eccentricity distribution. Here we develop and test a hierarchical probabilistic method for performing the relevant meta-analysis, that is, inferring the true eccentricity distribution, taking as input the likelihood functions for the individual-star eccentricities, or samplings of the posterior probability distributions for the eccentricities (under a given, uninformative prior). The method is a simple implementation of a hierarchical Bayesian model; it can also be seen as a kind of heteroscedastic deconvolution. It can be applied to any quantity measured with finite precision--other orbital parameters, or indeed any astronomical measurements of any kind, including magnitudes, parallaxes, or photometric redshifts--so long as the measurements have been communicated as a likelihood function or a posterior sampling.

We go through the many considerations involved in fitting a model to data, using as an example the fit of a straight line to a set of points in a two-dimensional plane. Standard weighted least-squares fitting is only appropriate when there is a dimension along which the data points have negligible uncertainties, and another along which all the uncertainties can be described by Gaussians of known variance; these conditions are rarely met in practice. We consider cases of general, heterogeneous, and arbitrarily covariant two-dimensional uncertainties, and situations in which there are bad data (large outliers), unknown uncertainties, and unknown but expected intrinsic scatter in the linear relationship being fit. Above all we emphasize the importance of having a "generative model" for the data, even an approximate one. Once there is a generative model, the subsequent fitting is non-arbitrary because the model permits direct computation of the likelihood of the parameters or the posterior probability distribution. Construction of a posterior probability distribution is indispensible if there are "nuisance parameters" to marginalize away.

The velocity distribution of nearby stars contains many "moving groups" that are inconsistent with the standard assumption of an axisymmetric, time-independent, and steady-state Galaxy. We study the age and metallicity properties of the low-velocity moving groups based on the reconstruction of the local velocity distribution in Paper I of this series. We perform stringent, conservative hypothesis testing to establish for each of these moving groups whether it could conceivably consist of a coeval population of stars. We conclude that they do not: the moving groups are not trivially associated with their eponymous open clusters nor with any other inhomogeneous star formation event. Concerning a possible dynamical origin of the moving groups, we test whether any of the moving groups has a higher or lower metallicity than the background population of thin disk stars, as would generically be the case if the moving groups are associated with resonances of the bar or spiral structure. We find clear evidence that the Hyades moving group has higher than average metallicity and weak evidence that the Sirius moving group has lower than average metallicity, which could indicate that these two groups are related to the inner Lindblad resonance of the spiral structure. Further we find weak evidence that the Hercules moving group has higher than average metallicity, as would be the case if it is associated with the bar's outer Lindblad resonance. The Pleiades moving group shows no clear metallicity anomaly, arguing against a common dynamical origin for the Hyades and Pleiades groups. Overall, however, the moving groups are barely distinguishable from the background population of stars, raising the likelihood that the moving groups are associated with transient perturbations. [abridged]

If a dynamical system is long-lived and non-resonant (that is, if there is a set of tracers that have evolved independently through many orbital times), and if the system is observed at any non-special time, it is possible to infer the dynamical properties of the system (such as the gravitational force or acceleration law) from a snapshot of the positions and velocities of the tracer population at a single moment in time. In this paper we describe a general inference technique that solves this problem while allowing (1) the unknown distribution function of the tracer population to be simultaneously inferred and marginalized over, and (2) prior information about the gravitational field and distribution function to be taken into account. As an example, we consider the simplest problem of this kind: We infer the force law in the Solar System using only an instantaneous kinematic snapshot (valid at 2009 April 1.0) for the eight major planets. We consider purely radial acceleration laws of the form a_r = -A [r/r_0]^{-\alpha}, where r is the distance from the Sun. Using a probabilistic inference technique, we infer 1.989 < \alpha < 2.052 (95 percent interval), largely independent of any assumptions about the distribution of energies and eccentricities in the system beyond the assumption that the system is phase-mixed. Generalizations of the methods used here will permit, among other things, inference of Milky Way dynamics from Gaia-like observations.

Masers found in massive star-forming regions can be located precisely in six-dimensional phase space and therefore serve as a tool for studying Milky Way dynamics. The non-random orbital phases at which the masers are found and the sparseness of current samples require modeling. Here we model the phase-space distribution function of 18 precisely measured Galactic masers, permitting a mean velocity offset and a general velocity dispersion tensor relative to their local standards of rest, and accounting for different pieces of prior information. With priors only on the Sun's distance from the Galactic Center and on its motion with respect to the local standard of rest, the maser data provide a weak constraint on the circular velocity at the Sun of V_c = 246 +/- 30 km/s. Including prior information on the proper motion of Sgr A* leads to V_c = 244 +/- 13 km/s. We do not confirm the value of V_c \approx 254 km/s found in more restrictive models. This analysis shows that there is no conflict between recent determinations of V_c from Galactic Center analyses, orbital fitting of the GD-1 stellar stream, and the kinematics of Galactic masers; a combined estimate is V_c = 236 +/- 11 km/s. Apart from the dynamical parameters, we find that masers tend to occur at post-apocenter, circular-velocity-lagging phases of their orbits.

We present a three-dimensional reconstruction of the velocity distribution of nearby stars (<~ 100 pc) using a maximum likelihood density estimation technique applied to the two-dimensional tangential velocities of stars. The underlying distribution is modeled as a mixture of Gaussian components. The algorithm reconstructs the error-deconvolved distribution function, even when the individual stars have unique error and missing-data properties. We apply this technique to the tangential velocity measurements from a kinematically unbiased sample of 11,865 main sequence stars observed by the Hipparcos satellite. We explore various methods for validating the complexity of the resulting velocity distribution function, including criteria based on Bayesian model selection and how accurately our reconstruction predicts the radial velocities of a sample of stars from the Geneva-Copenhagen survey (GCS). Using this very conservative external validation test based on the GCS, we find that there is little evidence for structure in the distribution function beyond the moving groups established prior to the Hipparcos mission. This is in sharp contrast with internal tests performed here and in previous analyses, which point consistently to maximal structure in the velocity distribution. We quantify the information content of the radial velocity measurements and find that the mean amount of new information gained from a radial velocity measurement of a single star is significant. This argues for complementary radial velocity surveys to upcoming astrometric surveys.

Conservation of the phase-space density of photons plus Lorentz invariance requires that the cosmological luminosity distance be larger than the angular diameter distance by a factor of $(1+z)^2$, where $z$ is the redshift. Because this is a fundamental symmetry, this prediction--known sometimes as the "Etherington relation" or the "Tolman test"--is independent of world model, or even the assumptions of homogeneity and isotropy. It depends, however, on Lorentz invariance and transparency. Transparency can be affected by intergalactic dust or interactions between photons and the dark sector. Baryon acoustic feature and type Ia supernovae measures of the expansion history are differently sensitive to the angular diameter and luminosity distances and can therefore be used in conjunction to limit cosmic transparency. At the present day, the comparison only limits the change $\Delta\tau$ in the optical depth from redshift 0.20 to 0.35 at visible wavelengths to $\Delta\tau < 0.13$ at 95% confidence. In a model with a constant comoving number density $n$ of scatterers of constant proper cross-section $\sigma$, this limit implies $n \sigma< 2\times10^{-4} h \Mpc^{-1}$. These limits depend weakly on the cosmological world model. Within the next few years, the limits could extend to redshifts $z\approx2.5$ and improve to $n \sigma<1.1 \times10^{-5} h \Mpc^{-1}$. Cosmic variance will eventually limit the sensitivity of any test using the BAF at the $n \sigma\sim 4\times10^{-7} h \Mpc^{-1}$ level. Comparison with other measures of the transparency is provided; no other measure in the visible is as free of astrophysical assumptions.

The recently introduced Sommerfeld enhancement of the dark matter annihilation cross section has important implications for the detection of dark matter annihilation in subhalos in the Galactic halo. In addition to the boost to the dark matter annihilation cross section from the high densities of these subhalos with respect to the main halo, an additional boost caused by the Sommerfeld enhancement results from the fact that they are kinematically colder than the Galactic halo. If we further believe the generic prediction of CDM that in each subhalo there is an abundance of substructure which is approximately self-similar to that of the Galactic halo, then I show that additional boosts coming from the density enhancements of these small substructures and their small velocity dispersions enhance the dark matter annihilation cross section even further. I find that very large boost factors ($10^5$ to $10^9$) are obtained in a large class of models. The implications of these boost factors for the detection of dark matter annihilation from dwarf Spheroidal galaxies in the Galactic halo are such that, generically, they outshine the background gamma-ray flux and are detectable by the Fermi Gamma-ray Space Telescope.

If there is a fifth force in the dark sector and dark sector particles interact non-gravitationally with ordinary matter, quantum corrections generically lead to a fifth force in the visible sector. We show how the strong experimental limits on fifth forces in the visible sector constrain the direct detection cross section, and the strength of the fifth force in the dark sector. If the latter is comparable to gravity, the spin-independent direct detection cross section must typically be <~ 10^{-55} cm^2. The anomalous acceleration of ordinary matter falling towards dark matter is also constrained: \eta_{OM-DM} <~ 10^{-8}.

If galaxy clusters contain intracluster dust, the spectra of galaxies lying behind clusters should show attenuation by dust absorption. We compare the optical (3500 - 7200 \AA) spectra of 60,267 luminous, early-type galaxies selected from the Sloan Digital Sky Survey to search for the signatures of intracluster dust in z ~ 0.05 clusters. We select massive, quiescent (i.e., non-star-forming) galaxies using an EW(Halpha) <= 2 \AA cut and consider galaxies in three bins of velocity dispersion, ranging from 150 to 300 km s^{-1}. The uniformity of early-type galaxy spectra in the optical allows us to construct inverse-variance-weighted composite spectra with high signal-to-noise ratio (ranging from 10^2-10^3). We compare the composite spectra of galaxies that lie behind and adjacent to galaxy clusters and find no convincing evidence of dust attenuation on scales ~ 0.15-2 Mpc; we derive a generic limit of E(B-V) < 3 x 10^{-3} mag on scales ~ 1-2 Mpc at the 99% confidence level, using conservative jackknife error bars, corresponding to a dust mass <~ 10^8 $M_{\odot}$. On scales smaller than 1 Mpc this limit is slightly weaker, E(B-V) < 8 x 10^{-3} mag.