Gas interior to the bar of the Milky Way has recently been shown as the closest example of a Low Ionization (Nuclear) Emission Region--LI(N)ER--in the universe. To better understand the nature of this gas, a sample of face-on galaxies with integral field spectroscopy are used to study the ionized gas conditions of 240 barred and 250 nonbarred galaxies, focusing on those that are most similar to the Milky Way. Strong optical line emission of $[NII]$ $\lambda 6584$, H$\alpha$, $[OIII]$ $\lambda 5007$, and H$\beta$ are used to diagnose the dominant ionization mechanisms of gas across galaxies and the Galaxy via Baldwin-Phillips-Terlevich (BPT) Diagrams. Barred galaxies show a strong suppression of star formation and an increase in composite and LI(N)ER like spectra in their inner regions when compared with similar nonbarred counterparts. This effect is lessened in galaxies of very low ($\log_{10}(M_\star/M_\odot) \lesssim 10.4$) or very high ($\log_{10}(M_\star/M_\odot) \gtrsim 11.1$) total stellar mass. Bar masks from Galaxy Zoo:3D show the bar's non-axisymmetric effect on the ionized gas and help predict the face-on distribution of ionized gas conditions near the bar of the Milky Way.

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

The Open Cluster Chemical Abundances and Mapping (OCCAM) survey aims to constrain key Galactic dynamical and chemical evolution parameters by the construction of a large, comprehensive, uniform, infrared-based spectroscopic data set of hundreds of open clusters. This fourth contribution from the OCCAM survey presents analysis using SDSS/APOGEE DR16 of a sample of 128 open clusters, 71 of which we designate to be "high quality" based on the appearance of their color-magnitude diagram. We find the APOGEE DR16 derived [Fe/H] abundances to be in good agreement with previous high resolution spectroscopic open cluster abundance studies. Using the high quality sample, we measure Galactic abundance gradients in 16 elements, and find evolution of some of the [X/Fe] gradients as a function of age. We find an overall Galactic [Fe/H] vs R_GC gradient of $-0.068 \pm 0.001$ dex kpc$^{-1}$ over the range of $6 <$ R_GC $< 13.9$ kpc; however, we note that this result is sensitive to the distance catalog used, varying as much as 15%. We formally derive the location a break in the [Fe/H] abundance gradient as a free parameter in the gradient fit for the first time. We also measure significant Galactic gradients in O, Mg, S, Ca, Mn, Cr, Cu, Na, Al, and K, some of which are measured for the first time. Our large sample allows us to explore four well-populated age bins to explore the time evolution of gradients for a large number of elements and comment on possible implications for Galactic chemical evolution and radial migration.

Studies of the kinematics and chemical compositions of Galactic globular clusters (GCs) enable the reconstruction of the history of star formation, chemical evolution, and mass assembly of the Galaxy. Using the latest data release (DR16) of the SDSS/APOGEE survey, we identify 3,090 stars associated with 46 GCs. Using a previously defined kinematic association, we break the sample down into eight separate groups and examine how the kinematics-based classification maps into chemical composition space, considering only $\alpha$ (mostly Si and Mg) elements and Fe. Our results show that: (i) The loci of both in situ and accreted subgroups in chemical space match those of their field counterparts; (ii) GCs from different individual accreted subgroups occupy the same locus in chemical space. This could either mean that they share a similar origin or that they are associated with distinct satellites which underwent similar chemical enrichment histories; (iii) The chemical compositions of the GCs associated with the low orbital energy subgroup defined by Massari and collaborators is broadly consistent with an in situ origin. However, at the low metallicity end, the distinction between accreted and in situ populations is blurred; (iv) Regarding the status of GCs whose origin is ambiguous, we conclude the following: the position in Si-Fe plane suggests an in situ}origin for Liller 1 and a likely accreted origin for NGC 5904 and NGC 6388. The case of NGC 288 is unclear, as its orbital properties suggest an accretion origin, its chemical composition suggests it may have formed in situ.

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

We present the first release of the MaNGA Stellar Library (MaStar), which is a large, well-calibrated, high-quality empirical library covering the wavelength range of 3,622-10,354A at a resolving power of R~1800. The spectra were obtained using the same instrument as used by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, by piggybacking on the SDSS-IV/APOGEE-2N observations. Compared to previous empirical libraries, the MaStar library will have a higher number of stars and a more comprehensive stellar-parameter coverage, especially of cool dwarfs, low-metallicity stars, and stars with different [alpha/Fe], achieved by a sophisticated target selection strategy that takes advantage of stellar-parameter catalogs from the literature. This empirical library will provide a new basis for stellar population synthesis, and is particularly well-suited for stellar-population analysis of MaNGA galaxies. The first version of the library contains 8646 high-quality per-visit spectra for 3321 unique stars. Compared to photometry, the relative flux calibration of the library is accurate to 3.9% in g-r, 2.7% in r-i, and 2.2% in i-z. The data are released as part of Sloan Digital Sky Survey Data Release 15. We expect the final release of the library to contain more than 10,000 stars.

The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.

The problem of the origin of the elements is a fundamental one in astronomy and one that has many open questions. Prominent examples include (1) the nature of Type Ia supernovae and the timescale of their contributions; (2) the observational identification of elements such as titanium and potassium with the $\alpha$-elements in conflict with core-collapse supernova predictions; (3) the number and relative importance of r-process sites; (4) the origin of carbon and nitrogen and the influence of mixing and mass loss in winds; and (5) the origin of the intermediate elements, such as Cu, Ge, As, and Se, that bridge the region between charged-particle and neutron-capture reactions. The next decade will bring to maturity many of the new tools that have recently made their mark, such as large-scale chemical cartography of the Milky Way and its satellites, the addition of astrometric and asteroseismic information, the detection and characterization of gravitational wave events, 3-D simulations of convection and model atmospheres, and improved laboratory measurements for transition probabilities and nuclear masses. All of these areas are key for continued improvement, and such improvement will benefit many areas of astrophysics.

We present elemental abundance analysis of high-resolution spectra for five giant stars, deriving Fe, Mg, Al, C, N, O, Si and Ce abundances, and spatially located within the innermost regions of the bulge globular cluster NGC 6522, based on H-band spectra taken with the multi-object APOGEE-north spectrograph from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. Of the five cluster candidates, two previously unremarked stars are confirmed to have second-generation (SG) abundance patterns, with the basic pattern of depletion in C and Mg simultaneous with enrichment in N and Al as seen in other SG globular cluster populations at similar metallicity. } In agreement with the most recent optical studies, the NGC 6522 stars analyzed exhibit (when available) only mild overabundances of the s-process element Ce, contradicting the idea of the NGC 6522 stars being formed from gas enriched by spinstars and indicating that other stellar sources such as massive AGB stars could be the primary intra-cluster medium polluters. The peculiar abundance signature of SG stars have been observed in our data, confirming the presence of multiple generations of stars in NGC 6522.

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.

The Galactic bulge has a bimodal metallicity distribution function: different kinematic, spatial, and, potentially, age distributions characterize the metal-poor and metal-rich components. Despite this observed dichotomy, which argues for different formation channels for those stars, the distribution of bulge stars in the $\alpha$-abundance versus metallicity plane has been found so far to be a rather smooth single sequence. We use data from the fourteenth data release of the APOGEE spectroscopic survey (DR14) to investigate the distribution in the Mg abundance (as tracer of the $\alpha$-elements)-versus-metallicity plane of a sample of stars selected to be in the inner region of the bulge. A clean sample has been selected from the DR14 using a set of data- and pipeline-flags to ensure the quality of their fundamental parameters and elemental abundances. An additional selection made use of computed spectro-photometric distances to select a sample of likely bulge stars as those with ${\rm R_{GC}\leq 3.5~kpc}$. We adopt magnesium abundance as an $\alpha$-abundance proxy for our clean sample as it has been proven to be the most accurate $\alpha$-element as determined by ASPCAP, the pipeline for data products from APOGEE spectra. From the distribution of our bulge sample in the [Mg/Fe]-versus-[Fe/H] plane, we found that the sequence is bimodal. This bimodality is given by the presence of a low-Mg sequence of stars parallel to the main high-Mg sequence over a range of $\sim 0.5$~dex around solar metallicity. The two sequences merge above ${\rm [Fe/H]\sim0.15~ dex}$ into a single sequence whose dispersion in [Mg/Fe] is larger than either of the two sequences visible at lower metallicity. This result is confirmed when we consider stars in our sample that are inside the bulge region according to trustworthy Gaia DR2 distances.

Galactic interstellar dust has a profound impact not only on our observations of objects throughout the Universe, but also on the morphology, star formation, and chemical evolution of the Galaxy. The advent of massive imaging and spectroscopic surveys (particularly in the infrared) places us on the threshold of being able to map the properties and dynamics of dust and the interstellar medium (ISM) in three dimensions throughout the Milky Way disk and bulge. These developments will enable a fundamentally new understanding of dust properties, including how grains respond to their local environment and how those environments affect dust attenuation of background objects of interest. Distance-resolved maps of dust motion also hold great promise for tracing the flow of interstellar material throughout the Galaxy on a variety of scales, from bar-streaming motions to the collapse and dissolution of individual molecular clouds. These advances require optical and infrared imaging of stars throughout the Galactic midplane, stretching many kiloparsecs from the Sun, matched with very dense spectroscopic coverage to probe the ISM's fine-grained structure.

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.

We present an overview of the distributions of 11 elemental abundances in the Milky Way's inner regions, as traced by APOGEE stars released as part of SDSS Data Release 14/15 (DR14/DR15), including O, Mg, Si, Ca, Cr, Mn, Co, Ni, Na, Al, and K. This sample spans ~4000 stars with R_GC<4 kpc, enabling the most comprehensive study to date of these abundances and their variations within the innermost few kiloparsecs of the Milky Way. We describe the observed abundance patterns ([X/Fe]-[Fe/H]), compare to previous literature results and to patterns in stars at the solar Galactic radius, and discuss possible trends with DR14/DR15 effective temperatures. We find that the position of the [Mg/Fe]-[Fe/H] "knee" is nearly constant with R_GC, indicating a well-mixed star-forming medium or high levels of radial migration in the early inner Galaxy. We quantify the linear correlation between pairs of elements in different subsamples of stars and find that these relationships vary; some abundance correlations are very similar between the alpha-rich and alpha-poor stars, but others differ significantly, suggesting variations in the metallicity dependencies of certain supernova yields. These empirical trends will form the basis for more detailed future explorations and for the refinement of model comparison metrics. That the inner Milky Way abundances appear dominated by a single chemical evolutionary track and that they extend to such high metallicities underscore the unique importance of this part of the Galaxy for constraining the ingredients of chemical evolution modeling and for improving our understanding of the evolution of the Galaxy as a whole.

In this work, we derive a spatially resolved map of the line-of-sight velocity of the interstellar medium and use it, along with a second map of line-of-sight velocity from Paper I of this series, to determine the nature of gaseous spiral structure in the Milky Way. This map is derived from measurements of the 1.527 $\mu$m diffuse interstellar band (DIB) in stellar spectra from the APOGEE survey and covers the nearest 4-5 kpc of the Northern Galactic plane. We cross-check this new DIB-based line-of-sight velocity map with the map derived in Paper I and find that they agree. We then compare these maps with line-of-sight velocity maps derived from simulations of quasi-stationary density wave spiral structure and dynamic, or material, spiral structure in a Milky Way-like galaxy. While none of the maps derived from these simulations is an exact match to the measured velocity field of the Milky Way, the measurements are more consistent with simulations of dynamic spiral structure. In the dynamic spiral structure simulation that best matches the measurements, the Perseus spiral arm is being disrupted.

The Open Cluster Chemical Abundances and Mapping (OCCAM) survey aims to produce a comprehensive, uniform, infrared-based spectroscopic dataset for hundreds of open clusters, and to constrain key Galactic dynamical and chemical parameters from this sample. This second contribution from the OCCAM survey presents analysis of 259 member stars with [Fe/H] determinations in 19 open clusters, using Sloan Digital Sky Survey Data Release 14 (SDSS/DR14) data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and ESA Gaia. This analysis, which includes clusters with R_{GC} ranging from 7 to 13 kpc, measures an [Fe/H] gradient of -0.061 \pm 0.004 dex/kpc. We also confirm evidence of a significant positive gradient in the \alpha-elements ([O/Fe], [Mg/Fe], and [Si/Fe]) and present evidence for a significant negative gradient in iron-peak elements ([Mn/Fe] and [Ni/Fe]).

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.

IC 166 is an intermediate-age open cluster ($\sim 1$ Gyr) which lies in the transition zone of the metallicity gradient in the outer disc. Its location, combined with our very limited knowledge of its salient features, make it an interesting object of study. We present the first high-resolution spectroscopic and precise kinematical analysis of IC 166, which lies in the outer disc with $R_{GC} \sim 12.7$ kpc. High resolution \textit{H}-band spectra were analyzed using observations from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We made use of the Brussels Automatic Stellar Parameter (BACCHUS) code to provide chemical abundances based on a line-by-line approach for up to eight chemical elements (Mg, Si, Ca, Ti, Al, K, Mn and Fe). The $\alpha-$element (Mg, Si, Ca and whenever available Ti) abundances, and their trends with Fe abundances have been analysed for a total of 13 high-likelihood cluster members. No significant abundance scatter was found in any of the chemical species studied. Combining the positional, heliocentric distance, and kinematic information we derive, for the first time, the probable orbit of IC 166 within a Galactic model including a rotating boxy bar, and found that it is likely that IC 166 formed in the Galactic disc, supporting its nature as an unremarkable Galactic open cluster with an orbit bound to the Galactic plane.

In this work we apply and expand on a recently introduced outlier detection algorithm that is based on an unsupervised random forest. We use the algorithm to calculate a similarity measure for stellar spectra from the Apache Point Observatory Galactic Evolution Experiment (APOGEE). We show that the similarity measure traces non-trivial physical properties and contains information about complex structures in the data. We use it for visualization and clustering of the dataset, and discuss its ability to find groups of highly similar objects, including spectroscopic twins. Using the similarity matrix to search the dataset for objects allows us to find objects that are impossible to find using their best fitting model parameters. This includes extreme objects for which the models fail, and rare objects that are outside the scope of the model. We use the similarity measure to detect outliers in the dataset, and find a number of previously unknown Be-type stars, spectroscopic binaries, carbon rich stars, young stars, and a few that we cannot interpret. Our work further demonstrates the potential for scientific discovery when combining machine learning methods with modern survey data.

Bars are common galactic structures in the local universe that play an important role in the secular evolution of galaxies, including the Milky Way. In particular, the velocity distribution of individual stars in our galaxy is useful to shed light on stellar dynamics, and provides information complementary to that inferred from the integrated light of external galaxies. However, since a wide variety of models reproduce the distribution of velocity and the velocity dispersion observed in the Milky Way, we look for signatures of the bar on higher-order moments of the line-of-sight velocity ($V_{los}$) distribution. We make use of two different numerical simulations --one that has developed a bar and one that remains nearly axisymmetric-- to compare them with observations in the latest APOGEE data release (SDSS DR14). This comparison reveals three interesting structures that support the notion that the Milky Way is a barred galaxy. A high skewness region found at positive longitudes constrains the orientation angle of the bar, and is incompatible with the orientation of the bar at $\ell=0^\circ$ proposed in previous studies. We also analyse the $V_{los}$ distributions in three regions, and introduce the Hellinger distance to quantify the differences among them. Our results show a strong non-Gaussian distribution both in the data and in the barred model, confirming the qualitative conclusions drawn from the velocity maps. In contrast to earlier work, we conclude it is possible to infer the presence of the bar from the kurtosis distribution.

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.

The Orion Star Forming Complex (OSFC) is a central target for the APOGEE-2 Young Cluster Survey. Existing membership catalogs span limited portions of the OSFC, reflecting the difficulty of selecting targets homogeneously across this extended, highly structured region. We have used data from wide field photometric surveys to produce a less biased parent sample of young stellar objects (YSOs) with infrared (IR) excesses indicative of warm circumstellar material or photometric variability at optical wavelengths across the full 420 square degrees extent of the OSFC. When restricted to YSO candidates with H < 12.4, to ensure S/N ~100 for a six visit source, this uniformly selected sample includes 1307 IR excess sources selected using criteria vetted by Koenig & Liesawitz and 990 optical variables identified in the Pan-STARRS1 3$\pi$ survey: 319 sources exhibit both optical variability and evidence of circumstellar disks through IR excess. Objects from this uniformly selected sample received the highest priority for targeting, but required fewer than half of the fibers on each APOGEE-2 plate. We fill the remaining fibers with previously confirmed and new color-magnitude selected candidate OSFC members. Radial velocity measurements from APOGEE-1 and new APOGEE-2 observations taken in the survey's first year indicate that ~90% of the uniformly selected targets have radial velocities consistent with Orion membership.The APOGEE-2 Orion survey will include >1100 bona fide YSOs whose uniform selection function will provide a robust sample for comparative analyses of the stellar populations and properties across all sub-regions of Orion.

Detailed chemical abundance distributions for fourteen elements are derived for eight high-probability stellar members of the solar metallicity old open cluster M67 with an age of $\sim$4 Gyr. The eight stars consist of four pairs, with each pair occupying a distinct phase of stellar evolution: two G-dwarfs, two turnoff stars, two G-subgiants, and two red clump K-giants. The abundance analysis uses near-IR high-resolution spectra ($\lambda$1.5 -- 1.7$\mu$m) from the APOGEE survey and derives abundances for C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, and Fe. Our derived stellar parameters and metallicity for 2M08510076+113115 suggest that this star is a solar-twin, exhibiting abundance differences relative to the Sun of $\leq$ 0.04 dex for all elements. Chemical homogeneity is found within each class of stars ($\sim$0.02 dex), while significant abundance variations ($\sim$0.05 -- 0.20 dex) are found across the different evolutionary phases; the turnoff stars typically have the lowest abundances, while the red clump tend to have the largest. Non-LTE corrections to the LTE-derived abundances are unlikely to explain the differences. A detailed comparison of the derived Fe, Mg, Si, and Ca abundances with recently published surface abundances from stellar models that include chemical diffusion, provides a good match between the observed and predicted abundances as a function of stellar mass. Such agreement would indicate the detection of chemical diffusion processes in the stellar members of M67.

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.

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.

A cold high-velocity (HV, $\sim$ 200 km/s) peak was first reported in several Galactic bulge fields based on the APOGEE commissioning observations. Both the existence and the nature of the high-velocity peak are still under debate. Here we revisit this feature with the latest APOGEE DR13 data. We find that most of the low latitude bulge fields display a skewed Gaussian distribution with a HV shoulder. However, only 3 out of 53 fields show distinct high-velocity peaks around 200 km/s. The velocity distribution can be well described by Gauss-Hermite polynomials, except the three fields showing clear HV peaks. We find that the correlation between the skewness parameter ($h_{3}$) and the mean velocity ($\bar{v}$), instead of a distinctive HV peak, is a strong indicator of the bar. It was recently suggested that the HV peak is composed of preferentially young stars. We choose three fields showing clear HV peaks to test this hypothesis using the metallicity, [$\alpha$/M] and [C/N] as age proxies. We find that both young and old stars show HV features. The similarity between the chemical abundances of stars in the HV peaks and the main component indicates that they are not systematically different in terms of chemical abundance or age. In contrast, there are clear differences in chemical space between stars in the Sagittarius dwarf and the bulge stars. The strong HV peaks off-plane are still to be explained properly, and could be different in nature.

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 kinematics and dynamics of young stellar populations enable us to test theories of star formation. With this aim, we continue our analysis of the SDSS-III/APOGEE IN-SYNC survey, a high resolution near infrared spectroscopic survey of young clusters. We focus on the Orion A star-forming region, for which IN-SYNC obtained spectra of $\sim2700$ stars. In Paper IV we used these data to study the young stellar population. Here we study the kinematic properties through radial velocities ($v_r$). The young stellar population remains kinematically associated with the molecular gas, following a $\sim10\:{\rm{km\:s}}^{-1}$ gradient along filament. However, near the center of the region, the $v_r$ distribution is slightly blueshifted and asymmetric; we suggest that this population, which is older, is slightly in foreground. We find evidence for kinematic subclustering, detecting statistically significant groupings of co-located stars with coherent motions. These are mostly in the lower-density regions of the cloud, while the ONC radial velocities are smoothly distributed, consistent with it being an older, more dynamically evolved cluster. The velocity dispersion $\sigma_v$ varies along the filament. The ONC appears virialized, or just slightly supervirial, consistent with an old dynamical age. Here there is also some evidence for on-going expansion, from a $v_r$--extinction correlation. In the southern filament, $\sigma_v$ is $\sim2$--$3$ times larger than virial in the L1641N region, where we infer a superposition along the line of sight of stellar sub-populations, detached from the gas. On the contrary, $\sigma_v$ decreases towards L1641S, where the population is again in agreement with a virial state.

We report the peculiar chemical abundance patterns of eleven atypical Milky Way (MW) field red giant stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). These atypical giants exhibit strong Al and N enhancements accompanied by C and Mg depletions, strikingly similar to those observed in the so-called second-generation (SG) stars of globular clusters (GCs). Remarkably, we find low-Mg abundances ([Mg/Fe]$<$0.0) together with strong Al and N overabundances in the majority (5/7) of the metal-rich ([Fe/H]$\gtrsim - 1.0$) sample stars, which is at odds with actual observations of SG stars in Galactic CGs of similar metallicities. This chemical pattern is unique and unprecedented among MW stars, posing urgent questions about its origin. These atypical stars could be former SG stars of dissolved GCs formed with intrinsically lower abundances of Mg and enriched Al (subsequently self-polluted by massive AGB stars) or the result of exotic binary systems. We speculate that the stars Mg-deficiency as well as the orbital properties suggest that they could have an extragalactic origin. This discovery should guide future dedicated spectroscopic searches of atypical stellar chemical patterns in our Galaxy; a fundamental step forward to understand the Galactic formation and evolution.

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.

We present radial velocity measurements for 70 high confidence, and 34 potential binary systems in fields containing the Perseus Molecular Cloud, Pleiades, NGC 2264, and the Orion A star forming region. 18 of these systems have been previously identified as binaries in the literature. Candidate double-lined spectroscopic binaries (SB2s) are identified by analyzing the cross-correlation functions (CCFs) computed during the reduction of each APOGEE spectrum. We identify sources whose CCFs are well fit as the sum of two Lorentzians as likely binaries, and provide an initial characterization of the system based on the radial velocities indicated by that dual fit. For systems observed over several epochs, we present mass ratios and systemic velocities; for two systems with observations on eight or more epochs, and which meet our criteria for robust orbital coverage, we derive initial orbital parameters. The distribution of mass ratios for multi-epoch sources in our sample peaks at q=1, but with a significant tail toward lower q values. Tables reporting radial velocities, systemic velocities, and mass ratios are provided online. We discuss future improvements to the radial velocity extraction method we employ, as well as limitations imposed by the number of epochs currently available in the APOGEE database. The Appendix contains brief notes from the literature on each system in the sample, and more extensive notes for select sources of interest.

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.

We utilize elemental-abundance information for Galactic red giant stars in five open clusters (NGC 7789, NGC 6819, M67, NGC 188, and NGC 6791) from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) DR13 dataset to age-date the chemical evolution of the high- and low-$\alpha$ element sequences of the Milky Way. Key to this time-stamping is the cluster NGC 6791, whose stellar members have mean abundances that place it in the high-$\alpha$, high-[Fe/H] region of the [$\alpha$/Fe]-[Fe/H] plane. Based on the cluster's age ($\sim 8$ Gyr), Galactocentric radius, and height above the Galactic plane, as well as comparable chemistry reported for APOGEE stars in Baade's Window, we suggest that the two most likely origins for NGC 6791 are as an original part of the thick-disk, or as a former member of the Galactic bulge. Moreover, because NGC 6791 lies at the \textit{high metallicity end} ([Fe/H] $\sim 0.4$) of the high-$\alpha$ sequence, the age of NGC 6791 places a limit on the \textit{youngest age} of stars in the high-metallicity, high-$\alpha$ sequence for the cluster's parent population (i.e., either the bulge or the disk). In a similar way, we can also use the age and chemistry of NGC 188 to set a limit of $\sim 7$ Gyr on the \textit{oldest age} of the low-$\alpha$ sequence of the Milky Way. Therefore, NGC 6791 and NGC 188 are potentially a pair of star clusters that bracket both the timing and the duration of an important transition point in the chemical history of the Milky Way.

In the updated APOGEE-Kepler catalog, we have asteroseismic and spectroscopic data for over 3000 first ascent red giants. Given the size and accuracy of this sample, these data offer an unprecedented test of the accuracy of stellar models on the post-main-sequence. When we compare these data to theoretical predictions, we find a metallicity dependent temperature offset with a slope of around 100 K per dex in metallicity. We find that this effect is present in all model grids tested and that theoretical uncertainties in the models, correlated spectroscopic errors, and shifts in the asteroseismic mass scale are insufficient to explain this effect. Stellar models can be brought into agreement with the data if a metallicity dependent convective mixing length is used, with $ \Delta\alpha_{\rm ML, YREC} \sim 0.2$ per dex in metallicity, a trend inconsistent with the predictions of three dimensional stellar convection simulations. If this effect is not taken into account, isochrone ages for red giants from the Gaia data will be off by as much as a factor of 2 even at modest deviations from solar metallicity ([Fe/H]=$-$0.5).

Baade's window (BW) is one of the most observed Galactic bulge fields in terms of chemical abundances. Due to its low and homogeneous interstellar absorption it is considered as a calibration field for Galactic bulge studies. In the era of large spectroscopic surveys, calibration fields such as BW are necessary to cross calibrate the stellar parameters and individual abundances of the APOGEE survey. We use the APOGEE BW stars to derive their metallicity distribution function (MDF) and individual abundances, for $\alpha$- and iron-peak elements of the APOGEE ASPCAP pipeline (DR13), as well as the age distribution for stars in BW. We determine the MDF of APOGEE stars in BW and find a remarkable agreement with that of the Gaia-ESO survey (GES). Both exhibit a clear bimodal distribution. We also find that the Mg-metallicity planes of both surveys agree well, except for the metal-rich part ([Fe/H] >0.1), where APOGEE finds systematically higher Mg abundances with respect to the GES. The ages based on the [C/N] ratio reveal a bimodal age distribution, with a major old population at 10 Gyr, with a decreasing tail towards younger stars. A comparison between APOGEE estimates and stellar parameters, and those determined by other sources reveals detectable systematic offsets, in particular for spectroscopic surface gravity estimates. In general, we find a good agreement between individual abundances of O, Na, Mg, Al, Si, K, Ca, Cr, Mn, Co, and Ni from APOGEE with that of literature values. We have shown that in general APOGEE data show a good agreement in terms of MDF and individual chemical abundances with respect to literature works. Using the [C/N] ration we found a significant fraction of young stars in BW which is in agreement with the model of Haywood et al. (2016).

We explore to which extent stars within Galactic disk open clusters resemble each other in the high-dimensional space of their photospheric element abundances, and contrast this with pairs of field stars. Our analysis is based on abundances for 20 elements, homogeneously derived from APOGEE spectra (with carefully quantified uncertainties, with a median value of $\sim 0.03$ dex). We consider 90 red giant stars in seven open clusters and find that most stars within a cluster have abundances in most elements that are indistinguishable (in a $\chi^2$-sense) from those of the other members, as expected for stellar birth siblings. An analogous analysis among pairs of $>1000$ field stars shows that highly significant abundance differences in the 20-dimensional space can be established for the vast majority of these pairs, and that the APOGEE-based abundance measurements have high discriminating power. However, pairs of field stars whose abundances are indistinguishable even at 0.03~dex precision exist: $\sim 0.3$ percent of all field star pairs, and $\sim 1.0$ percent of field star pairs at the same (solar) metallicity [Fe/H] = $0 \pm 0.02$. Most of these pairs are presumably not birth siblings from the same cluster, but rather doppelganger. Our analysis implies that 'chemical tagging' in the strict sense, identifying birth siblings for typical disk stars through their abundance similarity alone, will not work with such data. However, our approach shows that abundances have extremely valuable information for probabilistic chemo-orbital modeling and combined with velocities, we have identified new cluster members from the field.

We report chemical abundances obtained by SDSS-III/APOGEE for giant stars in five globular clusters located within 2.2 kpc of the Galactic centre. We detect the presence of multiple stellar populations in four of those clusters (NGC 6553, NGC 6528, Terzan 5, and Palomar 6) and find strong evidence for their presence in NGC 6522. All clusters present a significant spread in the abundances of N, C, Na, and Al, with the usual correlations and anti-correlations between various abundances seen in other globular clusters. Our results provide important quantitative constraints on theoretical models for self-enrichment of globular clusters, by testing their predictions for the dependence of yields of elements such as Na, N, C, and Al on metallicity. They also confirm that, under the assumption that field N-rich stars originate from globular cluster destruction, they can be used as tracers of their parental systems in the high- metallicity regime.

We report on detection, from observations obtained with the APOGEE spectroscopic survey, of a metal-poor ([Fe/H] $= -1.3$ dex) field giant star with an extreme Mg-Al abundance ratio ([Mg/Fe] $= -0.31$ dex; [Al/Fe] $= 1.49$ dex). Such low Mg/Al ratios are seen only among the second-generation population of globular clusters, and are not present among Galactic disk field stars. The light element abundances of this star, 2M16011638-1201525, suggest that it could have been born in a globular cluster. We explore several origin scenarios, in particular studying the orbit of the star to check the probability of it being kinematically related to known globular clusters. We performed simple orbital integrations assuming the estimated distance of 2M16011638-1201525 and the available six-dimensional phase-space coordinates of 63 globular clusters, looking for close encounters in the past with a minimum distance approach within the tidal radius of each cluster. We found a very low probability that 2M16011638-1201525 was ejected from most globular clusters; however, we note that the best progenitor candidate to host this star is globular cluster $\omega$ Centauri (NGC 5139). Our dynamical investigation demonstrates that 2M16011638-1201525 reaches a distance $|Z_{max}| < 3 $ kpc from the Galactic plane and a minimum and maximum approach to the Galactic center of $R_{min}<0.62$ kpc and $R_{max}<7.26$ kpc in an eccentric ($e\sim0.53$) and retrograde orbit. Since the extreme chemical anomaly of 2M16011638-1201525 has also been observed in halo field stars, this object could also be considered a halo contaminant, likely been ejected into the Milky Way disk from the halo. We conclude that, 2M16011638-20152 is also kinematically consistent with the disk but chemically consistent with halo field stars.

Much of the inner Milky Way's (MW) global rotation and velocity dispersion patterns can be reproduced by models of secularly-evolved, bar-dominated bulges. More sophisticated constraints, including the higher moments of the line-of-sight velocity distributions (LOSVDs) and limits on the chemodynamical substructure, are critical for interpreting observations of the unresolved inner regions of extragalactic systems and for placing the MW in context with other galaxies. Here, we use SDSS-APOGEE data to develop these constraints, by presenting the first maps of the LOSVD skewness and kurtosis of metal-rich and metal-poor inner MW stars (divided at [Fe/H] = -0.4), and comparing the observed patterns to those that are seen both in N-body models and in extragalactic bars. Despite closely matching the mean velocity and dispersion, the models do not reproduce the observed LOSVD skewness patterns in different ways, which demonstrates that our understanding of the detailed orbital structure of the inner MW remains an important regime for improvement. We find evidence in the MW of the skewness-velocity correlation that is used as a diagnostic of extragalactic bar/bulges. This correlation appears in metal-rich stars only, providing further evidence for different evolutionary histories of chemically differentiated populations. We connect these skewness measurements to previous work on high-velocity "peaks" in the inner Galaxy, confirming the presence of that phenomenon, and we quantify the cylindrical rotation of the inner Galaxy, finding that the latitude-independent rotation vanishes outside of lon ~ 7 deg. Finally, we evaluate the MW data in light of select extragalactic bar diagnostics and discuss progress and challenges of using the MW as a resolved analog of unresolved stellar populations.

With the advent of the space missions CoRoT and Kepler, it has become feasible to determine precise asteroseismic masses and ages for large samples of red-giant stars. In this paper, we present the CoRoGEE dataset -- obtained from CoRoT lightcurves for 606 red giant stars in two fields of the Galactic disc which have been co-observed for an ancillary project of APOGEE. We have used the Bayesian parameter estimation code PARAM to calculate distances, extinctions, masses, and ages for these stars in a homogeneous analysis, resulting in relative statistical uncertainties of $\sim2\%$ in distance, $\sim4\%$ in radius, $\sim9\%$ in mass and $\sim25\%$ in age. We also assess systematic age uncertainties due to different input physics and mass loss. We discuss the correlation between ages and chemical abundance patterns of field stars over a large radial range of the Milky Way's disc (5 kpc $<R_{\rm Gal}<$ 14 kpc), focussing on the [$\alpha$/Fe]-[Fe/H]-age plane in five radial bins of the Galactic disc. We find an overall agreement with the expectations of chemical-evolution models computed before the present data were available, especially for the outer regions. However, our data also indicate that a significant fraction of stars now observed near and beyond the Solar Neighbourhood migrated from inner regions. Mock CoRoGEE observations of a chemo-dynamical Milky Way disc model show that the number of high-metallicity stars in the outer disc is too high to be accounted for even by the strong radial mixing present in the model. The mock observations also reveal that the age distribution of the [$\alpha$/Fe]-enhanced sequence in the CoRoGEE inner-disc field is much broader than expected from a combination of radial mixing and observational errors. We suggest that a thick disc/bulge component that formed stars for more than 3 Gyr may account for these discrepancies.

Formation of globular clusters (GCs), the Galactic bulge, or galaxy bulges in general, are important unsolved problems in Galactic astronomy. Homogeneous infrared observations of large samples of stars belonging to GCs and the Galactic bulge field are one of the best ways to study these problems. We report the discovery by APOGEE of a population of field stars in the inner Galaxy with abundances of N, C, and Al that are typically found in GC stars. The newly discovered stars have high [N/Fe], which is correlated with [Al/Fe] and anti-correlated with [C/Fe]. They are homogeneously distributed across, and kinematically indistinguishable from, other field stars in the same volume. Their metallicity distribution is seemingly unimodal, peaking at [Fe/H]~-1, thus being in disagreement with that of the Galactic GC system. Our results can be understood in terms of different scenarios. N-rich stars could be former members of dissolved GCs, in which case the mass in destroyed GCs exceeds that of the surviving GC system by a factor of ~8. In that scenario, the total mass contained in so-called "first-generation" stars cannot be larger than that in "second-generation" stars by more than a factor of ~9 and was certainly smaller. Conversely, our results may imply the absence of a mandatory genetic link between "second generation" stars and GCs. Last, but not least, N-rich stars could be the oldest stars in the Galaxy, the by-products of chemical enrichment by the first stellar generations formed in the heart of the Galaxy.

We examine the distribution of the [O/Fe] abundance ratio in stars across the Galactic disk using H-band spectra from the Apache Point Galactic Evolution Experiment (APOGEE). We minimize systematic errors by considering groups of stars with similar atmospheric parameters. The APOGEE measurements in the Sloan Digital Sky Survey Data Release 12 reveal that the square root of the star-to-star cosmic variance in the oxygen-to-iron ratio at a given metallicity is about 0.03-0.04 dex in both the thin and thick disk. This is about twice as high as the spread found for solar twins in the immediate solar neighborhood and the difference is probably associated to the wider range of galactocentric distances spanned by APOGEE stars. We quantify the uncertainties by examining the spread among stars with the same parameters in clusters; these errors are a function of effective temperature and metallicity, ranging between 0.005 dex at 4000 K and solar metallicity, to about 0.03 dex at 4500 K and [Fe/H]= -0.6. We argue that measuring the spread in [O/Fe] and other abundance ratios provides strong constraints for models of Galactic chemical evolution.

Metallicity gradients provide strong constraints for understanding the chemical evolution of the Galaxy. We report on radial abundance gradients of Fe, Ni, Ca, Si, and Mg obtained from a sample of 304 red-giant members of 29 disk open clusters, mostly concentrated at galactocentric distances between ~8 - 15 kpc, but including two open clusters in the outer disk. The observations are from the APOGEE survey. The chemical abundances were derived automatically by the ASPCAP pipeline and these are part of the SDSS III Data Release 12. The gradients, obtained from least squares fits to the data, are relatively flat, with slopes ranging from -0.026 to -0.033 dex/kpc for the alpha-elements [O/H], [Ca/H], [Si/H] and [Mg/H] and -0.035 dex/kpc and -0.040 dex/kpc for [Fe/H] and [Ni/H], respectively. Our results are not at odds with the possibility that metallicity ([Fe/H]) gradients are steeper in the inner disk (R_GC ~7 - 12 kpc) and flatter towards the outer disk. The open cluster sample studied spans a significant range in age. When breaking the sample into age bins, there is some indication that the younger open cluster population in our sample (log age < 8.7) has a flatter metallicity gradient when compared with the gradients obtained from older open clusters.

We investigate interstellar extinction curve variations toward $\sim$4 deg$^{2}$ of the inner Milky Way in $VIJK_{s}$ photometry from the OGLE-III and $VVV$ surveys, with supporting evidence from diffuse interstellar bands and $F435W,F625W$ photometry. We obtain independent measurements toward $\sim$2,000 sightlines of $A_{I}$, $E(V-I)$, $E(I-J)$, and $E(J-K_{s})$, with median precision and accuracy of 2%. We find that the variations in the extinction ratios $A_{I}/E(V-I)$, $E(I-J)/E(V-I)$ and $E(J-K_{s})/E(V-I)$ are large (exceeding 20%), significant, and positively correlated, as expected. However, both the mean values and the trends in these extinction ratios are drastically shifted from the predictions of Cardelli and Fitzpatrick, regardless of how $R_{V}$ is varied. Furthermore, we demonstrate that variations in the shape of the extinction curve has at least two degrees of freedom, and not one (e.g. $R_{V}$), which we conform with a principal component analysis. We derive a median value of $<A_{V}/A_{Ks}>=13.44$, which is $\sim$60% higher than the "standard" value. We show that the Wesenheit magnitude $W_{I}=I-1.61(I-J)$ is relatively impervious to extinction curve variations. Given that these extinction curves are linchpins of observational cosmology, and that it is generally assumed that $R_{V}$ variations correctly capture variations in the extinction curve, we argue that systematic errors in the distance ladder from studies of type Ia supernovae and Cepheids may have been underestimated. Moreover, the reddening maps from the Planck experiment are shown to systematically overestimate dust extinction by $\sim$100%, and lack sensitivity to extinction curve variations.

We present the stellar kinematics across the Galactic bulge and into the disk at positive longitudes from the SDSS-III APOGEE spectroscopic survey of the Milky Way. APOGEE includes extensive coverage of the stellar populations of the bulge along the mid-plane and near-plane regions. From these data, we have produced kinematic maps of 10,000 stars across longitudes 0 deg < l < 65 deg, and primarily across latitudes of |b| < 5 deg in the bulge region. The APOGEE data reveal that the bulge is cylindrically rotating across all latitudes and is kinematically hottest at the very centre of the bulge, with the smallest gradients in both kinematic and chemical space inside the inner-most region (l,|b|) < (5,5) deg. The results from APOGEE show good agreement with data from other surveys at higher latitudes and a remarkable similarity to the rotation and dispersion maps of barred galaxies viewed edge on. The thin bar that is reported to be present in the inner disk within a narrow latitude range of |b| < 2 deg appears to have a corresponding signature in [Fe/H] and [alpha/Fe]. Stars with [Fe/H] > -0.5 have dispersion and rotation profiles that are similar to that of N-body models of boxy/peanut bulges. There is a smooth kinematic transition from the thin bar and boxy bulge (l,|b|) < (15,12) deg out into the disk for stars with [Fe/H] > -1.0, and the chemodynamics across (l,b) suggests the stars in the inner Galaxy with [Fe/H] > -1.0 have an origin in the disk.

We present the results of the SDSS APOGEE INfrared Spectroscopy of Young Nebulous Clusters program (IN-SYNC) survey of the Orion A molecular cloud. This survey obtained high resolution near infrared (NIR) spectroscopy of about 2700 young pre-main sequence stars throughout the region, acquired across five distinct fields spanning 6deg field of view (FOV). With these spectra, we have measured accurate stellar parameters (T_eff, log g, v sin i) and extinctions, and placed the sources in the Hertzsprung-Russel Diagram (HRD). We have also extracted radial velocities for the kinematic characterization of the population. We compare our measurements with literature results for a sub-sample of targets in order to assess the performances and accuracy of the survey. Source extinction shows evidence for dust grains that are larger than those in the diffuse interstellar medium (ISM): we estimate an average R_V=5.5 in the region. Importantly, we find a clear correlation between HRD inferred ages and spectroscopic surface-gravity inferred ages. This clearly indicates a real spread of stellar radii at fixed temperature, and together with additional correlations with extinction and with disk presence, strongly suggests a real spread of ages large than a few Myr. Focussing on the young population around NGC1980 iota Ori, which has previously been suggested to be a separate, foreground, older cluster, we confirm its older (5Myr) age and low A_V, but considering that its radial velocity distribution is indistinguishable from the Orion A's population, we suggest that NGC1980 is part of Orion A's star formation activity. Based on their stellar parameters and kinematic properties, we identify 383 new candidate members of Orion A, most of which are diskless sources in areas of the region poorly studied by previous works.

We present the $H$-band spectral line lists adopted by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). The APOGEE line lists comprise astrophysical, theoretical, and laboratory sources from the literature, as well as newly evaluated astrophysical oscillator strengths and damping parameters. We discuss the construction of the APOGEE line list, which is one of the critical inputs for the APOGEE Stellar Parameters and Chemical Abundances Pipeline, and present three different versions that have been used at various stages of the project. The methodology for the newly calculated astrophysical line lists is reviewed. The largest of these three line lists contains 134,457 molecular and atomic transitions. In addition to the format adopted to store the data, the line lists are available in MOOG, Synspec and Turbospectrum formats. We also present a list of $H$-band spectral features that are either poorly represented or completely missing in our line list. This list is based on the average of a large number of spectral fit residuals for APOGEE observations spanning a wide range of stellar parameters.

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 developed a code that estimates distances to stars using measured spectroscopic and photometric quantities. We employ a Bayesian approach to build the probability distribution function over stellar evolutionary models given these data, delivering estimates of model parameters for each star individually. The code was first tested on simulations, successfully recovering input distances to mock stars with <1% bias.The method-intrinsic random distance uncertainties for typical spectroscopic survey measurements amount to around 10% for dwarf stars and 20\% for giants, and are most sensitive to the quality of $\log g$ measurements. The code was validated by comparing our distance estimates to parallax measurements from the Hipparcos mission for nearby stars (< 300 pc), to asteroseismic distances of CoRoT red giant stars, and to known distances of well-studied open and globular clusters. The external comparisons confirm that our distances are subject to very small systematic biases with respect to the fundamental Hipparcos scale (+0.4 % for dwarfs, and +1.6% for giants). The typical random distance scatter is 18% for dwarfs, and 26% for giants. For the CoRoT-APOGEE sample, the typical random distance scatter is ~15%, both for the nearby and farther data. Our distances are systematically larger than the CoRoT ones by about +9%, which can mostly be attributed to the different choice of priors. The comparison to known distances of star clusters from SEGUE and APOGEE has led to significant systematic differences for many cluster stars, but with opposite signs, and with substantial scatter. Finally, we tested our distances against those previously determined for a high-quality sample of giant stars from the RAVE survey, again finding a small systematic trend of +5% and an rms scatter of 30%.

The inner Galactic Bulge has, until recently, been avoided in chemical evolution studies due to extreme extinction and stellar crowding. Large, near-IR spectroscopic surveys, such as APOGEE, allow for the first time the measurement of metallicities in the inner region of our Galaxy. We study metallicities of 33 K/M giants situated in the Galactic Center region from observations obtained with the APOGEE survey. We selected K/M giants with reliable stellar parameters from the APOGEE/ASPCAP pipeline. Distances, interstellar extinction values, and radial velocities were checked to confirm that these stars are indeed situated in the inner Galactic Bulge. We find a metal-rich population centered at [M/H] = +0.4 dex, in agreement with earlier studies of other bulge regions, but also a peak at low metallicity around $\rm [M/H] = -1.0\,dex$, suggesting the presence of a metal-poor population which has not previously been detected in the central region. Our results indicate a dominant metal-rich population with a metal-poor component that is enhanced in the $\alpha$-elements. This metal-poor population may be associated with the classical bulge and a fast formation scenario.

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.

New spectroscopic surveys offer the promise of consistent stellar parameters and abundances ('stellar labels') for hundreds of thousands of stars in the Milky Way: this poses a formidable spectral modeling challenge. In many cases, there is a sub-set of reference objects for which the stellar labels are known with high(er) fidelity. We take advantage of this with The Cannon, a new data-driven approach for determining stellar labels from spectroscopic data. The Cannon learns from the 'known' labels of reference stars how the continuum-normalized spectra depend on these labels by fitting a flexible model at each wavelength; then, The Cannon uses this model to derive labels for the remaining survey stars. We illustrate The Cannon by training the model on only 542 stars in 19 clusters as reference objects, with Teff, log g and [Fe/H] as the labels, and then applying it to the spectra of 56,000 stars from APOGEE DR10. The Cannon is very accurate. Its stellar labels compare well to the stars for which APOGEE pipeline (ASPCAP) labels are provided in DR10, with rms differences that are basically identical to the stated ASPCAP uncertainties. Beyond the reference labels, The Cannon makes no use of stellar models nor any line-list, but needs a set of reference objects that span label-space. The Cannon performs well at lower signal-to-noise, as it delivers comparably good labels even at one ninth the APOGEE observing time. We discuss the limitations of The Cannon and its future potential, particularly, to bring different spectroscopic surveys onto a consistent scale of stellar labels.

The molecular carriers of the ubiquitous absorption features called the diffuse interstellar bands (DIBs) have eluded identification for many decades, in part because of the enormous parameter space spanned by the candidates and the limited set of empirical constraints afforded by observations in the diffuse interstellar medium. Detection of these features in circumstellar regions, where the environmental properties are more easily measured, is thus a promising approach to understanding the chemical nature of the carriers themselves. Here, using high resolution spectra from the APOGEE survey, we present an analysis of the unusually asymmetric 1.5272 micron DIB feature along the sightline to the Red Square Nebula and demonstrate the likely circumstellar origin of about half of the DIB absorption in this line of sight. This interpretation is supported both by the velocities of the feature components and by the ratio of foreground to total reddening along the line of sight. The Red Square Nebula sightline offers the unique opportunity to study the behavior of DIB carriers in a constrained environment and thus to shed new light on the carriers themselves.

Most field stars will have encountered the highest stellar density and hence the largest number of interactions in their birth environment. Yet the stellar dynamics during this crucial phase are poorly understood. Here we analyze the radial velocities measured for 152 out of 380 observed stars in the 2-6 Myr old star cluster IC 348 as part of the SDSS-III APOGEE. The radial velocity distribution of these stars is fitted with one or two Gaussians, convolved with the measurement uncertainties including binary orbital motions. Including a second Gaussian improves the fit; the high-velocity outliers that are best fit by this second component may either (1) be contaminants from the nearby Perseus OB2 association, (2) be a halo of ejected or dispersing stars from IC 348, or (3) reflect that IC 348 has not relaxed to a Gaussian velocity distribution. We measure a velocity dispersion for IC 348 of $0.72 \pm 0.07$ km s$^{-1}$ (or $0.64 \pm 0.08$ km s$^{-1}$ if two Gaussians are fitted), which implies a supervirial state, unless the gas contributes more to the gravitational potential than expected. No evidence is found for a dependence of this velocity dispersion on distance from the cluster center or stellar mass. We also find that stars with lower extinction (in the front of the cloud) tend to be redshifted compared with stars with somewhat higher extinction (towards the back of the cloud). This data suggests that the stars in IC 348 are converging along the line of sight. We show that this correlation between radial velocity and extinction is unlikely to be spuriously caused by the small cluster rotation of $0.024 \pm 0.013$ km s$^{-1}$ arcmin$^{-1}$ or by correlations between the radial velocities of neighboring stars. This signature, if confirmed, will be the first detection of line-of-sight convergence in a star cluster(...)

Many astronomers now participate in large international collaborations, and it is important to examine whether these structures foster a scientific climate that is inclusive and diverse. The Committee on the Participation of Women in the Sloan Digital Sky Survey (CPWS) was formed to evaluate the demographics and gender climate within SDSS-IV, one of the largest and most geographically distributed astronomical collaborations. In April 2014, the CPWS administered a demographic survey to establish a baseline for the incipient SDSS-IV. We received responses from 250 participants (46% of the active membership). Half of the survey respondents were located in the US or Canada and 30% were based in Europe. Eleven percent of survey respondents considered themselves to be an ethnic minority at their current institution. Twenty-five percent of the SDSS-IV collaboration members are women, a fraction that is consistent with the US astronomical community, but substantially higher than the fraction of women in the IAU (16%). Approximately equal fractions of men and women report holding positions of leadership. When binned by academic age and career level, men and women also assume leadership roles at approximately equal rates, in a way that increases steadily for both genders with increasing seniority. In this sense, SDSS-IV has been successful in recruiting leaders that are representative of the collaboration. Yet, more progress needs to be made towards achieving gender balance and increasing diversity in the field of astronomy, and there is still room for improvement in the membership and leadership of SDSS-IV. For example, at the highest level of SDSS-IV leadership, women disproportionately assume roles related to education and public outreach. The CPWS plans to use these initial data to establish a baseline for tracking demographics over time as we work to assess and improve the climate of SDSS-IV.

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.

The Sloan Digital Sky Survey--III (SDSS--III) Apache Point Observatory Galactic Evolution Experiment (APOGEE) has obtained high resolution (R $\sim$ 22,500), high signal-to-noise ratio ($>$ 100) spectra in the H$-$band ($\sim$1.5$-$1.7 $\mu$m) for about 146,000 stars in the Milky Way galaxy. We have computed spectral libraries with effective temperature ($T\rm{_{eff}}$) ranging from 3500 to 8000 K for the automated chemical analy\-sis of the survey data. The libraries, used to derive stellar parameters and abundances from the APOGEE spectra in the SDSS--III data release 12 (DR12), are based on ATLAS9 model atmospheres and the ASS$\epsilon$T spectral synthesis code. We present a second set of libraries based on MARCS model atmospheres and the spectral synthesis code Turbospectrum. The ATLAS9/ASS$\epsilon$T ($T\rm{_{eff}}$ = 3500$-$8000 K) and MARCS/Turbospectrum ($T\rm{_{eff}}$ = 3500$-$5500 K) grids cover a wide range of metallicity ($-$2.5 $\leq$ [M/H] $\leq$ $+$0.5 dex), surface gravity (0 $\leq$ log $g$ $\leq$ 5 dex), microturbulence (0.5 $\leq$ $\xi$ $\leq$ 8 km~s$^{-1}$), carbon ($-$1 $\leq$ [C/M] $\leq$ $+$1 dex), nitrogen ($-$1 $\leq$ [N/M] $\leq$ $+$1 dex), and $\alpha$-element ($-$1 $\leq$ [$\alpha$/M] $\leq$ $+$1 dex) variations, having thus seven dimensions. We compare the ATLAS9/ASS$\epsilon$T and MARCS/Turbospectrum libraries and apply both of them to the analysis of the observed H$-$band spectra of the Sun and the K2 giant Arcturus, as well as to a selected sample of well-known giant stars observed at very high-resolution. The new APOGEE libraries are publicly available and can be employed for chemical studies in the H$-$band using other high-resolution spectrographs.

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.

The Kepler mission has yielded a large number of planet candidates from among the Kepler Objects of Interest (KOIs), but spectroscopic follow-up of these relatively faint stars is a serious bottleneck in confirming and characterizing these systems. We present motivation and survey design for an ongoing project with the SDSS-III multiplexed APOGEE near-infrared spectrograph to monitor hundreds of KOI host stars. We report some of our first results using representative targets from our sample, which include current planet candidates that we find to be false positives, as well as candidates listed as false positives that we do not find to be spectroscopic binaries. With this survey, KOI hosts are observed over ~20 epochs at a radial velocity precision of 100-200 m/s. These observations can easily identify a majority of false positives caused by physically-associated stellar or substellar binaries, and in many cases, fully characterize their orbits. We demonstrate that APOGEE is capable of achieving RV precision at the 100-200 m/s level over long time baselines, and that APOGEE's multiplexing capability makes it substantially more efficient at identifying false positives due to binaries than other single-object spectrographs working to confirm KOIs as planets. These APOGEE RVs enable ancillary science projects, such as studies of fundamental stellar astrophysics or intrinsically rare substellar companions. The coadded APOGEE spectra can be used to derive stellar properties (T_eff, log(g)) and chemical abundances of over a dozen elements to probe correlations of planet properties with individual elemental abundances.

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.

A Li-rich red giant star (2M19411367+4003382) recently discovered in the direction of NGC 6819 belongs to the rare subset of Li-rich stars that have not yet evolved to the luminosity bump, an evolutionary stage where models predict Li can be replenished. The currently favored model to explain Li enhancement in first-ascent red giants like 2M19411367+4003382 requires deep mixing into the stellar interior. Testing this model requires a measurement of 12C/13C, which is possible to obtain from APOGEE spectra. However, the Li-rich star also has abnormal asteroseismic properties that call into question its membership in the cluster, even though its radial velocity and location on color-magnitude diagrams are consistent with membership. To address these puzzles, we have measured a wide array of abundances in the Li-rich star and three comparison stars using spectra taken as part of the APOGEE survey to determine the degree of stellar mixing, address the question of membership, and measure the surface gravity. We confirm that the Li-rich star is a red giant with the same overall chemistry as the other cluster giants. However, its log g is significantly lower, consistent with the asteroseismology results and suggestive of a very low mass if the star is indeed a cluster member. Regardless of the cluster membership, the 12C/13C and C/N ratios of the Li-rich star are consistent with standard first dredge-up, indicating that Li dilution has already occurred, and inconsistent with internal Li enrichment scenarios that require deep mixing.

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.

The initial velocity dispersion of newborn stars is a major unconstrained aspect of star formation theory. Using near-infrared spectra obtained with the APOGEE spectrograph, we show that the velocity dispersion of young (1-2 Myr) stars in NGC 1333 is 0.92+/-0.12 km/s after correcting for measurement uncertainties and the effect of binaries. This velocity dispersion is consistent with the virial velocity of the region and the diffuse gas velocity dispersion, but significantly larger than the velocity dispersion of the dense, star-forming cores, which have a sub-virial velocity dispersion of 0.5 km/s. Since the NGC 1333 cluster is dynamically young and deeply embedded, this measurement provides a strong constraint on the initial velocity dispersion of newly-formed stars. We propose that the difference in velocity dispersion between stars and dense cores may be due to the influence of a 70 micro-Gauss magnetic field acting on the dense cores, or be the signature of a cluster with initial sub-structure undergoing global collapse.

The open cluster NGC 6791 is among the oldest, most massive and metal-rich open clusters in the Galaxy. High-resolution $H$-band spectra from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) of 11 red giants in NGC 6791 are analyzed for their chemical abundances of iron, oxygen, and sodium. The abundances of these three elements are found to be homogeneous (with abundance dispersions at the level of $\sim$ 0.05 - 0.07 dex) in these cluster red giants, which span much of the red-giant branch (T$_{\rm eff}$ $\sim$ 3500K - 4600K), and include two red-clump giants. From the infrared spectra, this cluster is confirmed to be among the most metal-rich clusters in the Galaxy ($<$[Fe/H]$>$ = 0.34 $\pm$ 0.06), and is found to have a roughly solar value of [O/Fe] and slightly enhanced [Na/Fe]. Non-LTE calculations for the studied Na I lines in the APOGEE spectral region ($\lambda$16373.86\AA\ and $\lambda$16388.85\AA) indicate only small departures from LTE ($\leq$ 0.04 dex) for the parameter range and metallicity of the studied stars. The previously reported double population of cluster members with different Na abundances is not found among the studied sample.

We map the distribution and properties of the Milky Way's interstellar medium as traced by diffuse interstellar bands (DIBs) detected in near-infrared stellar spectra from the SDSS-III/APOGEE survey. Focusing exclusively on the strongest DIB in the H-band, at ~1.527 microns, we present a projected map of the DIB absorption field in the Galactic plane, using a set of about 60,000 sightlines that reach up to 15 kpc from the Sun and probe up to 30 magnitudes of visual extinction. The strength of this DIB is linearly correlated with dust reddening over three orders of magnitude in both DIB equivalent width (W_DIB) and extinction, with a power law index of 1.01 +/- 0.01, a mean relationship of W_DIB/A_V = 0.1 Angstrom mag^-1, and a dispersion of ~0.05 Angstrom mag^-1 at extinctions characteristic of the Galactic midplane. These properties establish this DIB as a powerful, independent probe of dust extinction over a wide range of A_V values. The subset of about 14,000 robustly detected DIB features have an exponential W_DIB distribution. We empirically determine the intrinsic rest wavelength of this transition to be lambda_0 = 15,272.42 Angstrom, and then calculate absolute radial velocities of the carrier, which display the kinematical signature of the rotating Galactic disk. We probe the DIB carrier distribution in three dimensions and show that it can be characterized by an exponential disk model with a scaleheight of about 100 pc and a scalelength of about 5 kpc. Finally, we show that the DIB distribution also traces large-scale Galactic structures, including the central long bar and the warp of the outer disk.

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.

APOGEE has amassed the largest ever collection of multi-epoch, high-resolution (R~22,500), H-band spectra for B-type emission line (Be) stars. The 128/238 APOGEE Be stars for which emission had never previously been reported serve to increase the total number of known Be stars by ~6%. We focus on identification of the H-band lines and analysis of the emission peak velocity separations (v_p) and emission peak intensity ratios (V/R) of the usually double-peaked H I and non-hydrogen emission lines. H I Br11 emission is found to preferentially form in the circumstellar disks at an average distance of ~2.2 stellar radii. Increasing v_p toward the weaker Br12--Br20 lines suggests these lines are formed interior to Br11. By contrast, the observed IR Fe II emission lines present evidence of having significantly larger formation radii; distinctive phase lags between IR Fe II and H I Brackett emission lines further supports that these species arise from different radii in Be disks. Several emission lines have been identified for the first time including ~16895, a prominent feature in the spectra for almost a fifth of the sample and, as inferred from relatively large v_p compared to the Br11-Br20, a tracer of the inner regions of Be disks. Unlike the typical metallic lines observed for Be stars in the optical, the H-band metallic lines, such as Fe II 16878, never exhibit any evidence of shell absorption, even when the H I lines are clearly shell-dominated. The first known example of a quasi-triple-peaked Br11 line profile is reported for HD 253659, one of several stars exhibiting intra- and/or extra-species V/R and radial velocity variation within individual spectra. Br11 profiles are presented for all discussed stars, as are full APOGEE spectra for a portion of the sample.

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 Kepler Input Catalog (KIC) provides reddening estimates for its stars, based on the assumption of a simple exponential dusty screen. This project focuses on evaluating and improving these reddening estimates for the KIC's giant stars, for which extinction is a much more significant concern than for the nearby dwarf stars. We aim to improve the calibration (and thus consistency) amongst various photometric and spectroscopic temperatures of stars in the Kepler field by removing systematics due to incorrect extinction assumptions. The revised extinction estimates may then be used to derive improved stellar and planetary properties. We plan to eventually use the large number of KIC stars as probes into the structure and properties of the Galactic ISM.

Over two years 8,859 high-resolution H-band spectra of 3,493 young (1 - 10 Myr) stars were gathered by the multi-object spectrograph of the APOGEE project as part of the IN-SYNC ancillary program of that SDSS-III survey. Here we present the forward modeling approach used to derive effective temperatures, surface gravities, radial velocities, rotational velocities, and H-band veiling from these near-infrared spectra. We discuss in detail the statistical and systematic uncertainties in these stellar parameters. In addition we present accurate extinctions by measuring the E(J-H) of these young stars with respect to the single-star photometric locus in the Pleiades. Finally we identify an intrinsic stellar radius spread of about 25% for late-type stars in IC 348 using three (nearly) independent measures of stellar radius, namely the extinction-corrected J-band magnitude, the surface gravity and the $R \sin i$ from the rotational velocities and literature rotation periods. We exclude that this spread is caused by uncertainties in the stellar parameters by showing that the three estimators of stellar radius are correlated, so that brighter stars tend to have lower surface gravities and larger $R \sin i$ than fainter stars at the same effective temperature. Tables providing the spectral and photometric parameters for the Pleiades and IC 348 have been provided online.

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.

Galactic interstellar extinction maps are powerful and necessary tools for Milky Way structure and stellar population analyses, particularly toward the heavily-reddened bulge and in the midplane. However, due to the difficulty of obtaining reliable extinction measures and distances for a large number of stars that are independent of these maps, tests of their accuracy and systematics have been limited. Our goal is to assess a variety of photometric stellar extinction estimates, including both 2D and 3D extinction maps, using independent extinction measures based on a large spectroscopic sample of stars towards the Milky Way bulge. We employ stellar atmospheric parameters derived from high-resolution $H$-band APOGEE spectra, combined with theoretical stellar isochrones, to calculate line-of-sight extinction and distances for a sample of more than 2400 giants towards the Milky Way bulge. We compare these extinction values to those predicted by individual near-IR and near+mid-IR stellar colors, 2D bulge extinction maps and 3D extinction maps. The long baseline, near+mid-IR stellar colors are, on average, the most accurate predictors of the APOGEE extinction estimates, and the 2D and 3D extinction maps derived from different stellar populations along different sightlines show varying degrees of reliability. We present the results of all of the comparisons and discuss reasons for the observed discrepancies. We also demonstrate how the particular stellar atmospheric models adopted can have a strong impact on this type of analysis, and discuss related caveats.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) - one of the Sloan Digital Sky Survey III programs -- is using near-infrared spectra of $\sim 100,000$ red giant branch star candidates to study the structure of the Milky Way. In the course of the survey, APOGEE also acquires spectra of hot field stars to serve as telluric calibrators for the primary science targets. We report the serendipitous discovery of two rare, fast-rotating B stars of the $\sigma$ Ori E type among those blue field stars observed during the first year of APOGEE operations. Both of the discovered stars display the spectroscopic signatures of the rigidly rotating magnetospheres (RRM) common to this class of highly-magnetized ($B \sim 10$ kiloGauss) stars, increasing the number of known RRM stars by $\sim 10 \%$. One (HD 345439) is a main-sequence B star with unusually strong He absorption (similar to $\sigma$ Ori E), while the other (HD 23478) fits a "He-normal" B3IV classification. We combine the APOGEE discovery spectra with other optical and near-infrared spectra of these two stars, and of $\sigma$ Ori E itself, to show how near-infrared spectroscopy can be a uniquely powerful tool for discovering more of these rare objects, which may show little/no RRM signatures in their optical spectra. We discuss the potential for further discovery of $\sigma$ Ori E type stars, as well as the implications of our discoveries for the population of these objects and insights into their origin and evolution.

Fundamental stellar properties, such as mass, radius, and age, can be inferred using asteroseismology. Cool stars with convective envelopes have turbulent motions that can stochastically drive and damp pulsations. The properties of the oscillation frequency power spectrum can be tied to mass and radius through solar-scaled asteroseismic relations. Stellar properties derived using these scaling relations need verification over a range of metallicities. Because the age and mass of halo stars are well-constrained by astrophysical priors, they provide an independent, empirical check on asteroseismic mass estimates in the low-metallicity regime. We identify nine metal-poor red giants (including six stars that are kinematically associated with the halo) from a sample observed by both the Kepler space telescope and the Sloan Digital Sky Survey-III APOGEE spectroscopic survey. We compare masses inferred using asteroseismology to those expected for halo and thick-disk stars. Although our sample is small, standard scaling relations, combined with asteroseismic parameters from the APOKASC Catalog, produce masses that are systematically higher (<{\Delta}M>=0.17+/-0.05 Msun) than astrophysical expectations. The magnitude of the mass discrepancy is reduced by known theoretical corrections to the measured large frequency separation scaling relationship. Using alternative methods for measuring asteroseismic parameters induces systematic shifts at the 0.04 Msun level. We also compare published asteroseismic analyses with scaling relationship masses to examine the impact of using the frequency of maximum power as a constraint. Upcoming APOKASC observations will provide a larger sample of ~100 metal-poor stars, important for detailed asteroseismic characterization of Galactic stellar populations.

We investigate the chemo-kinematic properties of the Milky Way disc by exploring the first year of data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE), and compare our results to smaller optical high-resolution samples in the literature, as well as results from lower resolution surveys such as GCS, SEGUE and RAVE. We start by selecting a high-quality sample in terms of chemistry ($\sim$ 20.000 stars) and, after computing distances and orbital parameters for this sample, we employ a number of useful subsets to formulate constraints on Galactic chemical and chemodynamical evolution processes in the Solar neighbourhood and beyond (e.g., metallicity distributions -- MDFs, [$\alpha$/Fe] vs. [Fe/H] diagrams, and abundance gradients). Our red giant sample spans distances as large as 10 kpc from the Sun. We find remarkable agreement between the recently published local (d $<$ 100 pc) high-resolution high-S/N HARPS sample and our local HQ sample (d $<$ 1 kpc). The local MDF peaks slightly below solar metallicity, and exhibits an extended tail towards [Fe/H] $= -$1, whereas a sharper cut-off is seen at larger metallicities. The APOGEE data also confirm the existence of a gap in the [$\alpha$/Fe] vs. [Fe/H] abundance diagram. When expanding our sample to cover three different Galactocentric distance bins, we find the high-[$\alpha$/Fe] stars to be rare towards the outer zones, as previously suggested in the literature. For the gradients in [Fe/H] and [$\alpha$/Fe], measured over a range of 6 $ < $ R $ <$ 11 kpc in Galactocentric distance, we find a good agreement with the gradients traced by the GCS and RAVE dwarf samples. For stars with 1.5 $<$ z $<$ 3 kpc, we find a positive metallicity gradient and a negative gradient in [$\alpha$/Fe].

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

The dynamics of the core of the Sagittarius (Sgr) dwarf spheroidal (dSph) galaxy are explored using high-resolution (R~22,500), H-band, near-infrared spectra of over 1,000 giant stars in the central 3 deg^2 of the system, of which 328 are identified as Sgr members. These data, among some of the earliest observations from the SDSS-III/Apache Point Observatory Galactic Evolution Experiment (APOGEE) and the largest published sample of high resolution Sgr dSph spectra to date, reveal a distinct gradient in the velocity dispersion of Sgr from 11-14 km/s for radii >0.8 degrees from center to a dynamical cold point of 8 km/s in the Sgr center --- a trend differing from that found in previous kinematical analyses of Sgr over larger scales that suggest a more or less flat dispersion profile at these radii. Well-fitting mass models with either cored and cusped dark matter distributions can be found to match the kinematical results, although the cored profile succeeds with significantly more isotropic stellar orbits than required for a cusped profile. It is unlikely that the cold point reflects an unusual mass distribution. The dispersion gradient may arise from variations in the mixture of populations with distinct kinematics within the dSph; this explanation is suggested (e.g., by detection of a metallicity gradient across similar radii), but not confirmed, by the present data. Despite these remaining uncertainties about their interpretation, these early test data (including some from instrument commissioning) demonstrate APOGEE's usefulness for precision dynamical studies, even for fields observed at extreme airmasses.

The SDSS-III Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a three year survey that is collecting 100,000 high-resolution spectra in the near-IR across multiple Galactic populations. To derive stellar parameters and chemical compositions from this massive data set, the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP) has been developed. Here, we describe empirical calibrations of stellar parameters presented in the first SDSS-III APOGEE data release (DR10). These calibrations were enabled by observations of 559 stars in 20 globular and open clusters. The cluster observations were supplemented by observations of stars in NASA's Kepler field that have well determined surface gravities from asteroseismic analysis. We discuss the accuracy and precision of the derived stellar parameters, considering especially effective temperature, surface gravity, and metallicity; we also briefly discuss the derived results for the abundances of the alpha-elements, carbon, and nitrogen. Overall, we find that ASPCAP achieves reasonably accurate results for temperature and metallicity, but suffers from systematic errors in surface gravity. We derive calibration relations that bring the raw ASPCAP results into better agreement with independently determined stellar parameters. The internal scatter of ASPCAP parameters within clusters suggests that, metallicities are measured with a precision better than 0.1 dex, effective temperatures better than 150 K, and surface gravities better than 0.2 dex. The understanding provided by the clusters and Kepler giants on the current accuracy and precision will be invaluable for future improvements of the pipeline.

The Open Cluster Chemical Analysis and Mapping (OCCAM) Survey aims to produce a comprehensive, uniform, infrared-based dataset for hundreds of open clusters, and constrain key Galactic dynamical and chemical parameters from this sample. This first contribution from the OCCAM survey presents analysis of 141 members stars in 28 open clusters with high-resolution metallicities derived from a large uniform sample collected as part of the SDSS-III/Apache Point Observatory Galactic Evolution Experiment (APOGEE). This sample includes the first high-resolution metallicity measurements for 22 open clusters. With this largest ever uniformly observed sample of open cluster stars we investigate the Galactic disk gradients of both [M/H] and [alpha/M]. We find basically no gradient across this range in [alpha/M], but [M/H] does show a gradient for R_{GC} < 10 kpc and a significant flattening beyond R_{GC} = 10 kpc. In particular, whereas fitting a single linear trend yields an [M/H] gradient of -0.09 +/- 0.03$ dex/kpc --- similar to previously measure gradients inside 13 kpc --- by independently fitting inside and outside 10 kpc separately we find a significantly steeper gradient near the Sun (7.9 <= R_{GC} <= 10) than previously found (-0.20 +/- 0.08 dex/kpc) and a nearly flat trend beyond 10 kpc (-0.02 +/- 0.09 dex/kpc).

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution infrared spectroscopic survey spanning all Galactic environments (i.e., bulge, disk, and halo), with the principal goal of constraining dynamical and chemical evolution models of the Milky Way. APOGEE takes advantage of the reduced effects of extinction at infrared wavelengths to observe the inner Galaxy and bulge at an unprecedented level of detail. The survey's broad spatial and wavelength coverage enables users of APOGEE data to address numerous Galactic structure and stellar populations issues. In this paper we describe the APOGEE targeting scheme and document its various target classes to provide the necessary background and reference information to analyze samples of APOGEE data with awareness of the imposed selection criteria and resulting sample properties. APOGEE's primary sample consists of ~100,000 red giant stars, selected to minimize observational biases in age and metallicity. We present the methodology and considerations that drive the selection of this sample and evaluate the accuracy, efficiency, and caveats of the selection and sampling algorithms. We also describe additional target classes that contribute to the APOGEE sample, including numerous ancillary science programs, and we outline the targeting data that will be included in the public data releases.

We are carrying out a large ancillary program with the SDSS-III, using the fiber-fed multi-object NIR APOGEE spectrograph, to obtain high-resolution H-band spectra of more than 1200 M dwarfs. These observations are used to measure spectroscopic rotational velocities, radial velocities, physical stellar parameters, and variability of the target stars. Here, we describe the target selection for this survey and results from the first year of scientific observations based on spectra that is publicly available in the SDSS-III DR10 data release. As part of this paper we present RVs and vsini of over 200 M dwarfs, with a vsini precision of ~2 km/s and a measurement floor at vsini = 4 km/s. This survey significantly increases the number of M dwarfs studied for vsini and RV variability (at ~100-200 m/s), and will advance the target selection for planned RV and photometric searches for low mass exoplanets around M dwarfs, such as HPF, CARMENES, and TESS. Multiple epochs of radial velocity observations enable us to identify short period binaries, and AO imaging of a subset of stars enables the detection of possible stellar companions at larger separations. The high-resolution H-band APOGEE spectra provide the opportunity to measure physical stellar parameters such as effective temperatures and metallicities for many of these stars. At the culmination of this survey, we will have obtained multi-epoch spectra and RVs for over 1400 stars spanning spectral types of M0-L0, providing the largest set of NIR M dwarf spectra at high resolution, and more than doubling the number of known spectroscopic vsini values for M dwarfs. Furthermore, by modeling telluric lines to correct for small instrumental radial velocity shifts, we hope to achieve a relative velocity precision floor of 50 m/s for bright M dwarfs. We present preliminary results of this telluric modeling technique in this paper.

Despite its importance for understanding the nature of early stellar generations and for constraining Galactic bulge formation models, at present little is known about the metal-poor stellar content of the central Milky Way. This is a consequence of the great distances involved and intervening dust obscuration, which challenge optical studies. However, the Apache Point Observatory Galactic Evolution Experiment (APOGEE), a wide-area, multifiber, high-resolution spectroscopic survey within Sloan Digital Sky Survey III (SDSS-III), is exploring the chemistry of all Galactic stellar populations at infrared wavelengths, with particular emphasis on the disk and the bulge. An automated spectral analysis of data on 2,403 giant stars in twelve fields in the bulge obtained during APOGEE commissioning yielded five stars with low metallicity([Fe/H]$\le-1.7$), including two that are very metal-poor [Fe/H]$\sim-2.1$ by bulge standards. Luminosity-based distance estimates place the five stars within the outer bulge, where other 1,246 of the analyzed stars may reside. A manual reanalysis of the spectra verifies the low metallicities, and finds these stars to be enhanced in the $\alpha$-elements O, Mg, and Si without significant $\alpha$-pattern differences with other local halo or metal-weak thick-disk stars of similar metallicity, or even with other more metal-rich bulge stars. While neither the kinematics nor chemistry of these stars can yet definitively determine which, if any, are truly bulge members, rather than denizens of other populations co-located with the bulge, the newly-identified stars reveal that the chemistry of metal-poor stars in the central Galaxy resembles that of metal-weak thick-disk stars at similar metallicity.

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.

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.

Commissioning observations with the Apache Point Observatory Galactic Evolution Experiment (APOGEE), part of the Sloan Digital Sky Survey III, have produced radial velocities (RVs) for ~4700 K/M-giant stars in the Milky Way bulge. These high-resolution (R \sim 22,500), high-S/N (>100 per resolution element), near-infrared (1.51-1.70 um; NIR) spectra provide accurate RVs (epsilon_v~0.2 km/s) for the sample of stars in 18 Galactic bulge fields spanning -1<l<20 deg, |b|<20 deg, and dec>-32 deg. This represents the largest NIR high-resolution spectroscopic sample of giant stars ever assembled in this region of the Galaxy. A cold (sigma_v~30 km/s), high-velocity peak (V_GSR \sim +200 km/s) is found to comprise a significant fraction (~10%) of stars in many of these fields. These high RVs have not been detected in previous MW surveys and are not expected for a simple, circularly rotating disk. Preliminary distance estimates rule out an origin from the background Sagittarius tidal stream or a new stream in the MW disk. Comparison to various Galactic models suggests that these high RVs are best explained by stars in orbits of the Galactic bar potential, although some observational features remain unexplained.

We provide new, high-resolution A(Ks) extinction maps of the heavily reddened Galactic midplane based on the Rayleigh-Jeans Color Excess ("RJCE") method. RJCE determines star-by-star reddening based on a combination of near- and mid-infrared photometry. The new RJCE-generated maps have 2 x 2 arcmin pixels and span some of the most severely extinguished regions of the Galaxy -- those covered with Spitzer+IRAC imaging by the GLIMPSE-I, -II, -3D, and Vela-Carina surveys, from 256<l<65 deg and, in general, for |b| <= 1-1.5 deg (extending up to |b|<=4 deg in the bulge). Using RJCE extinction measurements, we generate dereddened color-magnitude diagrams and, in turn, create maps based on main sequence, red clump, and red giant star tracers, each probing different distances and thereby providing coarse three-dimensional information on the relative placement of dust cloud structures. The maps generated from red giant stars, which reach to ~18-20 kpc, probe beyond most of the Milky Way extinction in most directions and provide close to a "total Galactic extinction" map -- at minimum they provide high angular resolution maps of lower limits on A(Ks). Because these maps are generated directly from measurements of reddening by the very dust being mapped, rather than inferred on the basis of some less direct means, they are likely the most accurate to date for charting in detail the highly patchy differential extinction in the Galactic midplane. We provide downloadable FITS files and an IDL tool for retrieving extinction values for any line of sight within our mapped regions.

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)

The Milky Way (MW) remains a primary laboratory for understanding the structure and evolution of spiral galaxies, but typically we are denied clear views of MW stellar populations at low Galactic latitudes because of extinction by interstellar dust. However, the combination of 2MASS near-infrared (NIR) and Spitzer-IRAC mid-infrared (MIR) photometry enables a powerful method for determining the line of sight reddening to any star: the sampled wavelengths lie in the Rayleigh-Jeans part of the spectral energy distribution of most stars, where, to first order, all stars have essentially the same intrinsic color. Thus, changes in stellar NIR-MIR colors due to interstellar reddening are readily apparent, and (under an assumed extinction law) the observed colors and magnitudes of stars can be easily and accurately restored to their intrinsic values, greatly increasing their usefulness for Galactic structure studies. In this paper we explore this "Rayleigh-Jeans Color Excess" (RJCE) method and demonstrate that use of even a simple variant of the RJCE method based on a single reference color, (H-[4.5um]), can rather accurately remove dust effects from previously uninterpretable 2MASS color-magnitude diagrams of stars in fields along the heavily reddened Galactic mid-plane, with results far superior to those derived from application of other dereddening methods. We also show that "total" Galactic midplane extinction looks rather different from that predicted using 100um emission maps from the IRAS/ISSA and COBE/DIRBE instruments as presented by Schlegel et al. Instead, the Galactic mid-plane extinction strongly resembles the distribution of 13-CO (J=1->0) emission. Future papers will focus on refining the RJCE method and applying the technique to understand better not only dust and its distribution, but the distribution of stars intermixed with the dust in the low-latitude Galaxy.

The Sloan Digital Sky Survey III/Apache Point Observatory Galactic Evolution Experiment (SDSS-III/APOGEE) is a large-scale spectroscopic survey of Galactic stars and star clusters. The SDSS-III/APOGEE survey is designed to produce high-S/N, R = 27,500-31,000 spectra that cover a wavelength range of 1.51 to 1.68 microns. By utilizing APOGEE's excellent kinematics (error <= 0.5 km/s) and abundances (errors ~ 0.1 dex), we will be able to study star cluster kinematics and chemical properties in detail. Over the course of the 3-year survey beginning in 2011, APOGEE will target 25-30 key open and globular clusters. In addition, the large area coverage of the SDSS focal plane will also allow us to target stars in 100-200 additional star clusters during the main survey observations. We present the strength of APOGEE for both open and globular star cluster studies and the methods of identifying probable clusters members utilizing 2MASS and IRAC/WISE data.

We combine near-infrared (2MASS) and mid-infrared (Spitzer-IRAC) photometry to characterize the IR extinction law (1.2-8 microns) over nearly 150 degrees of contiguous Milky Way midplane longitude. The relative extinctions in 5 passbands across these wavelength and longitude ranges are derived by calculating color excess ratios for G and K giant red clump stars in contiguous midplane regions and deriving the wavelength dependence of extinction in each one. Strong, monotonic variations in the extinction law shape are found as a function of angle from the Galactic center, symmetric on either side of it. These longitudinal variations persist even when dense interstellar regions, known a priori to have a shallower extinction curve, are removed. The increasingly steep extinction curves towards the outer Galaxy indicate a steady decrease in the absolute-to-selective extinction ratio (R_V) and in the mean dust grain size at greater Galactocentric angles. We note an increasing strength of the 8 micron extinction inflection at high Galactocentric angles and, using theoretical dust models, show that this behavior is consistent with the trend in R_V. Along several lines of sight where the solution is most feasible, A_lambda/A_Ks as a function of Galactic radius is estimated and shown to have a Galactic radial dependence. Our analyses suggest that the observed relationship between extinction curve shape and Galactic longitude is due to an intrinsic dependence of the extinction law on Galactocentric radius.

We present observations of Taurus-Auriga Class I/II protostars obtained with the Spitzer InfraRed Spectrograph. Detailed spectral fits to the 6 and 15 micron features are made, using publicly-available laboratory data, to constrain the molecular composition, abundances, and levels of thermal processing along the lines of sight. We provide an inventory of the molecular environments observed, which have an average composition dominated by water ice with ~12% CO_2 (abundance relative to H_2O), >~2-9% CH_3OH, <~14% NH_3, ~4% CH_4, ~2% H_2CO, ~0.6% HCOOH, and ~0.5% SO_2. We find CO_2/H_2O ratios nearly equivalent to those observed in cold clouds and lines of sight toward the galactic center. The unidentified 6.8 micron profiles vary from source to source, and it is shown to be likely that even combinations of the most common candidates (NH_4+ and CH_3OH) are inadequate to explain the feature fully. We discuss correlations among SED spectral indices, abundance ratios, and thermally-processed ice fractions and their implications for CO_2 formation and evolution. Comparison of our spectral fits with cold molecular cloud sight-lines indicates abundant prestellar ice environments made even richer by the radiative effects of protostars. Our results add additional constraints and a finer level of detail to current full-scale models of protostellar and protoplanetary systems.