Bars inhabit the majority of local-Universe disk galaxies and may be important drivers of galaxy evolution through the redistribution of gas and angular momentum within disks. We investigate the star formation and gas properties of bars in galaxies spanning a wide range of masses, environments, and star formation rates using the MaNGA galaxy survey. Using a robustly-defined sample of 684 barred galaxies, we find that fractional (or scaled) bar length correlates with the host's offset from the star-formation main sequence. Considering the morphology of the H$\alpha$ emission we separate barred galaxies into different categories, including barred, ringed, and central configurations, together with H$\alpha$ detected at the ends of a bar. We find that only low-mass galaxies host star formation along their bars, and that this is located predominantly at the leading edge of the bar itself. Our results are supported by recent simulations of massive galaxies, which show that the position of star formation within a bar is regulated by a combination of shear forces, turbulence and gas flows. We conclude that the physical properties of a bar are mostly governed by the existing stellar mass of the host galaxy, but that they also play an important role in the galaxy's ongoing star formation.

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

Dust attenuation in star-forming spiral galaxies affects stars and gas in different ways due to local variations in dust geometry. We present spatially resolved measurements of dust attenuation for a sample of 232 such star-forming spiral galaxies, derived from spectra acquired by the SDSS-IV MaNGA survey. The dust attenuation affecting the stellar populations of these galaxies (obtained using full spectrum stellar population fitting methods) is compared with the dust attenuation in the gas (derived from the Balmer decrement). Both of these attenuation measures increase for local regions of galaxies with higher star formation rates; the dust attenuation affecting the stellar populations increases more so than the dust attenuation in the gas, causing the ratio of the dust attenuation affecting the stellar populations to the dust attenuation in the gas to decrease for local regions of galaxies with higher star formation rate densities. No systematic difference is discernible in any of these dust attenuation quantities between the spiral arm and inter-arm regions of the galaxies. While both the dust attenuation in the gas and the dust attenuation affecting the stellar populations decrease with galactocentric radius, the ratio of the two quantities does not vary with radius. This ratio does, however, decrease systematically as the stellar mass of the galaxy increases. Analysis of the radial profiles of the two dust attenuation measures suggests that there is a disproportionately high concentration of birth clouds (incorporating gas, young stars and clumpy dust) nearer to the centres of star-forming spiral galaxies.

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.

In this paper we use strong line calibrations of N/O and O/H in MaNGA spaxel data to explore the systematics introduced by variations in N/O on various strong-line metallicity diagnostics. We find radial variations in N/O at fixed O/H which correlate with total galaxy stellar-mass; and which can induce $\sim 40 \%$ systematic uncertainties in oxygen abundance gradients when nitrogen-dependent abundance calibrations are used. Empirically, we find that these differences are associated with variation in the local star formation efficiency, as predicted by recent chemical evolution models for galaxies, but we cannot rule out other processes such as radial migration and the accretion of passive dwarf galaxies also playing a role.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is acquiring integral-field spectroscopy for the largest sample of galaxies to date. By 2020, the MaNGA Survey --- one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) --- will have observed a statistically representative sample of 10$^4$ galaxies in the local Universe ($z\lesssim0.15$). In addition to a robust data-reduction pipeline (DRP), MaNGA has developed a data-analysis pipeline (DAP) that provides higher-level data products. To accompany the first public release of its code base and data products, we provide an overview of the MaNGA DAP, including its software design, workflow, measurement procedures and algorithms, performance, and output data model. In conjunction with our companion paper Belfiore et al., we also assess the DAP output provided for 4718 observations of 4648 unique galaxies in the recent SDSS Data Release 15 (DR15). These analysis products focus on measurements that are close to the data and require minimal model-based assumptions. Namely, we provide stellar kinematics (velocity and velocity dispersion), emission-line properties (kinematics, fluxes, and equivalent widths), and spectral indices (e.g., D4000 and the Lick indices). We find that the DAP provides robust measurements and errors for the vast majority ($>$99%) of analyzed spectra. We summarize assessments of the precision and accuracy of our measurements as a function of signal-to-noise, and provide specific guidance to users regarding the limitations of the data. The MaNGA DAP software is publicly available and we encourage community involvement in its development.

Classification of intermediate redshift ($z$ = 0.3--0.8) emission line galaxies as star-forming galaxies, composite galaxies, active galactic nuclei (AGN), or low-ionization nuclear emission regions (LINERs) using optical spectra alone was impossible because the lines used for standard optical diagnostic diagrams: [NII], H$\alpha$, and [SII] are redshifted out of the observed wavelength range. In this work, we address this problem using four supervised machine learning classification algorithms: $k$-nearest neighbors (KNN), support vector classifier (SVC), random forest (RF), and a multi-layer perceptron (MLP) neural network. For input features, we use properties that can be measured from optical galaxy spectra out to $z < 0.8$---[OIII]/H$\beta$, [OII]/H$\beta$, [OIII] line width, and stellar velocity dispersion---and four colors ($u-g$, $g-r$, $r-i$, and $i-z$) corrected to $z=0.1$. The labels for the low redshift emission line galaxy training set are determined using standard optical diagnostic diagrams. RF has the best area under curve (AUC) score for classifying all four galaxy types, meaning highest distinguishing power. Both the AUC scores and accuracies of the other algorithms are ordered as MLP$>$SVC$>$KNN. The classification accuracies with all eight features (and the four spectroscopically-determined features only) are 93.4% (92.3%) for star-forming galaxies, 69.4% (63.7%) for composite galaxies, 71.8% (67.3%) for AGNs, and 65.7% (60.8%) for LINERs. The stacked spectrum of galaxies of the same type as determined by optical diagnostic diagrams at low redshift and RF at intermediate redshift are broadly consistent. Our publicly available code (https://github.com/zkdtc/MLC_ELGs) and trained models will be instrumental for classifying emission line galaxies in upcoming wide-field spectroscopic surveys.

A galaxy's stellar mass is one of its most fundamental properties, but it remains challenging to measure reliably. With the advent of very large optical spectroscopic surveys, efficient methods that can make use of low signal-to-noise spectra are needed. With this in mind, we created a new software package for estimating effective stellar mass-to-light ratios $\log \Upsilon^*$ that uses principal component analysis(PCA) basis set to optimize the comparison between observed spectra and a large library of stellar population synthesis models. In Paper I, we showed that a with a set of six PCA basis vectors we could faithfully represent most optical spectra from the Mapping Nearby Galaxies at APO (MaNGA) survey;and we tested the accuracy of our M/L estimates using synthetic spectra. Here, we explore sources of systematic error in our mass measurements by comparing our new measurements to data from the literature. We compare our stellar mass surface density estimates to kinematics-derived dynamical mass surface density measurements from the DiskMass Survey and find some tension between the two which could be resolved if the disk scale-heights used in the kinematic analysis were overestimated by a factor of $\sim 1.5$. We formulate an aperture-corrected stellar mass catalog for the MaNGA survey, and compare to previous stellar mass estimates based on multi-band optical photometry, finding typical discrepancies of 0.1 dex. Using the spatially resolved MaNGA data, we evaluate the impact of estimating total stellar masses from spatially unresolved spectra, and we explore how the biases that result from unresolved spectra depend upon the galaxy's dust extinction and star formation rate. Finally, we describe a SDSS Value-Added Catalog which will include both spatially resolved and total (aperture-corrected) stellar masses for MaNGA galaxies.

We present a method of fitting optical spectra of galaxies using a basis set of six vectors obtained from principal component analysis (PCA) of a library of synthetic spectra of 40000 star formation histories (SFHs). Using this library, we provide estimates of resolved effective stellar mass-to-light ratio ($\log \Upsilon^*$) for thousands of galaxies from the SDSS-IV/MaNGA integral-field spectroscopic survey. Using a testing framework built on additional synthetic SFHs, we show that the estimates of stellar mass-to-light ratio are reliable (as are their uncertainties) at a variety of signal-to-noise ratios, stellar metallicities, and dust attenuation conditions. Finally, we describe the future release of the resolved stellar mass-to-light ratios as a SDSS-IV/MaNGA Value-Added Catalog (VAC) and provide a link to the software used to conduct this analysis.

The late assembly of massive galaxies is thought to be dominated by stellar accretion in their outskirts (beyond 2 effective radii Re) due to dry, minor galaxy mergers. We use observations of 1010 passive early-type galaxies (ETGs) within z<0.15 from SDSS IV MaNGA to search for evidence of this accretion. The outputs from the stellar population fitting codes FIREFLY, pPXF, and Prospector are compared to control for systematic errors in stellar metallicity (Z) estimation. We find that the average radial logZ/Zsun profiles of ETGs in various stellar mass (M) bins are not linear. As a result, these profiles are poorly characterized by a single gradient value, explaining why weak trends reported in previous work can be difficult to interpret. Instead, we examine the full radial extent of stellar metallicity profiles and find them to flatten in the outskirts of M>10^{11}Msun ETGs. This is a signature of stellar accretion. Based on a toy model for stellar metallicity profiles, we infer the ex-situ stellar mass fraction in ETGs as a function of M and galactocentric radius. We find that ex-situ stars at 2Re make up 20% of the projected stellar mass of M<10^{10.5}Msun ETGs, rising up to 80% for M>10^{11.5}Msun ETGs.

Bars in galaxies are thought to stimulate both inflow of material and radial mixing along them. Observational evidence for this mixing has been inconclusive so far however, limiting the evaluation of the impact of bars on galaxy evolution. We now use results from the MaNGA integral field spectroscopic survey to characterise radial stellar age and metallicity gradients along the bar and outside the bar in 128 strongly barred galaxies. We find that age and metallicity gradients are flatter in the barred regions of almost all barred galaxies when compared to corresponding disk regions at the same radii. Our results re-emphasize the key fact that by azimuthally averaging integral field spectroscopic data one loses important information from non-axisymmetric galaxy components such as bars and spiral arms. We interpret our results as observational evidence that bars are radially mixing material in galaxies of all stellar masses, and for all bar morphologies and evolutionary stages.

We investigate the resolved relation between local extinction and star formation surface density within nearby star-forming galaxies selected from the MaNGA survey. Balmer decrement measurements imply an extinction of the H{\alpha} line emission which scales approximately linearly with the logarithm of the star formation surface density: $ A_{H{\alpha}} = 0.46 \log(\Sigma_{SFR}) + 1.53$. Secondary dependencies are observed such that, at a given $\Sigma_{SFR}$, regions of lower metallicity and/or enhanced H{\alpha} equivalent width (EW) suffer less obscuration than regions of higher metallicity and/or lower H{\alpha} EW. Spaxels lying above the mean relation also tend to belong to galaxies that are more massive, larger and viewed under higher inclination than average. We present a simple model in which the observed trends can be accounted for by a metallicity-dependent scaling between $\Sigma_{SFR}$ and $\Sigma_{dust}$ via a super-linear Kennicutt-Schmidt relation ($n_{KS} \sim 1.47$) and a dust-to-gas ratio which scales linearly with metallicity (DGR($Z_{\odot}$) = 0.01). The relation between the resulting total dust column and observed effective extinction towards nebular regions requires a geometry for the relative distribution of H{\alpha} emitting regions and dust that deviates from a uniform foreground screen and also from an entirely homogeneous mixture of dust and emitting sources. The best-fit model features an H{\alpha} EW and galactocentric distance dependent fraction of the dust mass in a clumpy foreground screen in front of a homogeneous mixture.

We use data from 1222 late-type star-forming galaxies in the SDSS IV MaNGA survey to identify regions in which the gas-phase metallicity is anomalously-low compared to expectations from the tight empirical relation between metallicity and stellar surface mass-density at a given stellar mass. We find anomalously low metallicity (ALM) gas in 10% of the star-forming spaxels, and in 25% of the galaxies in the sample. The incidence rate of ALM gas increases strongly with both global and local measures of the specific star-formation rate, and is higher in lower mass galaxies and in the outer regions of galaxies. The incidence rate is also significantly higher in morphologically disturbed galaxies. We estimate that the lifetimes of the ALM regions are a few hundred Myr. We argue that the ALM gas has been delivered to its present location by a combination of interactions, mergers, and accretion from the halo, and that this infusion of gas stimulates star-formation. Given the estimated lifetime and duty cycle of such events, we estimate that the time-averaged accretion rate of ALM gas is similar to the star-formation rate in late type galaxies over the mass-range M$_* \sim10^9$ to 10$^{10}$ M$_{\odot}$.

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

The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, one of three core programs of the fourth-generation Sloan Digital Sky Survey (SDSS-IV), is producing a massive, high-dimensional integral field spectroscopic data set. However, leveraging the MaNGA data set to address key questions about galaxy formation presents serious data-related challenges due to the combination of its spatially inter-connected measurements and sheer volume. For each galaxy, the MaNGA pipelines produce relatively large data files to preserve the spatial correlations of the spectra and measurements, but this comes at the expense of storing the data set in a coarsely-chunked manner. The coarse chunking and total volume of the data make it time-consuming to download and curate locally-stored data. Thus, accessing, querying, visually exploring, and performing statistical analyses across the whole data set at a fine-grained scale is extremely challenging using just FITS files. To overcome these challenges, we have developed \marvin: a toolkit consisting of a Python package, Application Programming Interface (API), and web application utilizing a remote database. \marvin's robust and sustainable design minimizes maintenance, while facilitating user-contributed extensions such as high level analysis code. Finally, we are in the process of abstracting out \marvin's core functionality into a separate product so that it can serve as a foundation for others to develop \marvin-like systems for new science applications.

Chemical enrichment signatures strongly constrain galaxy formation and evolution, and a detailed understanding of abundance patterns provides clues regarding the nucleosynthetic production pathways of elements. Using the SDSS-IV MaNGA IFU survey, we study radial gradients of chemical element abundances in detail. We use stacked spectra out to 1 Re of 366 early-type galaxies with masses 9.9 - 10.8 log $M/M_{\odot}$ to probe the abundances of the elements C, N, Na, Mg, Ca, and Ti, relative to the abundance of Fe, by fitting stellar population models to a combination of Lick absorption indices. We find that C, Mg, and Ti trace each other both as a function of galaxy radius and galaxy mass. These similar C and Mg abundances within and across galaxies set a lower limit for star-formation timescales. Conversely, N and Ca are generally offset to lower abundances. The under-abundance of Ca compared to Mg implies delayed enrichment of Ca through Type Ia supernovae, whereas the correlated behaviour of Ti and the lighter $\alpha$ elements, C and Mg, suggest contributions to Ti from Type II supernovae. We obtain shallow radial gradients in [Mg/Fe], [C/Fe], and [Ti/Fe], meaning that these inferences are independent of radius. However, we measure strong negative radial gradients for [N/Fe] and [Na/Fe], of up to $-0.25\pm0.05$ and $-0.29\pm0.02$ dex/Re respectively. These gradients become shallower with decreasing galaxy mass. We find that N and Na abundances increase more steeply with velocity dispersion within galaxies than globally, while the other elements show the same relation locally and globally. This implies that the high Na and N abundances found in massive early type galaxies are generated by internal processes within galaxies. These are strongly correlated with the total metallicity, suggesting metallicity-dependent Na enrichment, and secondary N production in massive early-type galaxies.

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.

A key task of observational extragalactic astronomy is to determine where -- within galaxies of diverse masses and morphologies -- stellar mass growth occurs, how it depends on galaxy properties and what processes regulate star formation. Using spectroscopic indices derived from the stellar continuum at $\sim 4000$\AA, we determine the spatially resolved star-formation histories of 980000 spaxels in 2404 galaxies in the SDSS-IV MaNGA IFU survey. We examine the spatial distribution of star-forming, quiescent, green valley, starburst and post-starburst spaxels as a function of stellar mass and morphology to see where and in what types of galaxy star formation is occurring. The spatial distribution of star-formation is dependent primarily on stellar mass, with a noticeable change in the distribution at \mstar$>10^{10}$\msun. Galaxies above this mass have an increasing fraction of regions that are forming stars with increasing radius, whereas lower mass galaxies have a constant fraction of star forming regions with radius. Our findings support a picture of inside-out growth and quenching at high masses. We find that morphology (measured via concentration) correlates with the fraction of star-forming spaxels, but not with their radial distribution. We find (post-)starburst regions are more common outside of the galaxy centre, are preferentially found in asymmetric galaxies, and have lower gas-phase metallicity than other regions, consistent with interactions triggering starbursts and driving low metallicity gas into regions at $<1.5R_e$.

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.

MaNGA provides the opportunity to make precise spatially resolved measurements of the IMF slope in galaxies owing to its unique combination of spatial resolution, wavelength coverage and sample size. We derive radial gradients in age, element abundances and IMF slope analysing optical and near-infrared absorption features from stacked spectra out to the half-light radius of 366 early-type galaxies with masses $9.9 - 10.8\;\log M/M_{\odot}$. We find flat gradients in age and [$\alpha$/Fe] ratio, as well as negative gradients in metallicity, consistent with the literature. We further derive significant negative gradients in the [Na/Fe] ratio with galaxy centres being well enhanced in Na abundance by up to 0.5 dex. Finally, we find a gradient in IMF slope with a bottom-heavy IMF in the centre (typical mass excess factor of 1.5) and a Milky Way-type IMF at the half-light radius. This pattern is mass-dependent with the lowest mass galaxies in our sample featuring only a shallow gradient around a Milky Way IMF. Our results imply the local IMF-$\sigma$ relation within galaxies to be even steeper than the global relation and hint towards the local metallicity being the dominating factor behind the IMF variations. We also employ different stellar population models in our analysis and show that a radial IMF gradient is found independently of the stellar population model used. A similar analysis of the Wing-Ford band provides inconsistent results and further evidence of the difficulty in measuring and modelling this particular feature.

We present a catalogue of 38 spectroscopically detected strong galaxy-galaxy gravitational lens candidates identified in the Sloan Digital Sky Survey IV (SDSS-IV). We were able to simulate narrow-band images for 8 of them demonstrating evidence of multiple images. Two of our systems are compound lens candidates, each with 2 background source-planes. One of these compound systems shows clear lensing features in the narrow-band image. Our sample is based on 2812 galaxies observed by the Mapping Nearby Galaxies at APO (MaNGA) integral field unit (IFU). This Spectroscopic Identification of Lensing Objects (SILO) survey extends the methodology of the Sloan Lens ACS Survey (SLACS) and BOSS Emission-Line Survey (BELLS) to lower redshift and multiple IFU spectra. We searched ~ 1.5 million spectra, of which 3065 contained multiple high signal-to-noise background emission-lines or a resolved [OII] doublet, that are included in this catalogue. Upon manual inspection, we discovered regions with multiple spectra containing background emission-lines at the same redshift, providing evidence of a common source-plane geometry which was not possible in previous SLACS and BELLS discovery programs. We estimate more than half of our candidates have an Einstein radius > 1.7", which is significantly greater than seen in SLACS and BELLS. These larger Einstein radii produce more extended images of the background galaxy increasing the probability that a background emission-line will enter one of the IFU spectroscopic fibres, making detection more likely.

We present new evidence for AGN feedback in a subset of 69 quenched low-mass galaxies ($M_{\star} \lesssim 5\times10^{9}$ M$_{\odot}$, $M_{\rm{r}} > -19$) selected from the first two years of the SDSS-IV MaNGA survey. The majority (85 per cent) of these quenched galaxies appear to reside in a group environment. We find 6 galaxies in our sample that appear to have an active AGN that is preventing on-going star-formation; this is the first time such a feedback mechanism has been observed in this mass range. Interestingly, five of these six galaxies have an ionised gas component that is kinematically offset from their stellar component, suggesting the gas is either recently accreted or outflowing. We hypothesise these six galaxies are low-mass equivalents to the "red geysers" observed in more massive galaxies. Of the other 63 galaxies in the sample, we find 8 do appear for have some low-level, residual star formation, or emission from hot, evolved stars. The remaining galaxies in our sample have no detectable ionised gas emission throughout their structures, consistent with them being quenched. This work shows the potential for understanding the detailed physical properties of dwarf galaxies through spatially resolved spectroscopy.

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.

In this paper, we investigate 2727 galaxies observed by MaNGA as of June 2016 to develop spatially resolved techniques for identifying signatures of active galactic nuclei (AGN). We identify 303 AGN candidates. The additional spatial dimension imposes challenges in identifying AGN due to contamination from diffuse ionized gas, extra-planar gas and photoionization by hot stars. We show that the combination of spatially-resolved line diagnostic diagrams and additional cuts on H$\alpha$ surface brighness and H$\alpha$ equivalent width can distinguish between AGN-like signatures and high-metallicity galaxies with LINER-like spectra. Low mass galaxies with high specific star formation rates are particularly difficult to diagnose and routinely show diagnostic line ratios outside of the standard star-formation locus. We develop a new diagnostic -- the distance from the standard diagnostic line in the line-ratios space -- to evaluate the significance of the deviation from the star-formation locus. We find 173 galaxies that would not have been selected as AGN candidates based on single-fibre spectral measurements but exhibit photoionization signatures suggestive of AGN activity in the MaNGA resolved observations, underscoring the power of large integral field unit (IFU) surveys. A complete census of these new AGN candidates is necessary to understand their nature and probe the complex co-evolution of supermassive black holes and their hosts.

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.

We analyze the optical continuum of star-forming galaxies in SDSS by fitting stacked spectra with stellar population synthesis models to investigate the relation between stellar mass, stellar metallicity, dust attenuation and star formation rate. We fit models calculated with star formation and chemical evolution histories that are derived empirically from multi-epoch observations of the stellar mass---star formation rate and the stellar mass---gas-phase metallicity relations, respectively. We also fit linear combinations of single burst models with a range of metallicities and ages. Star formation and chemical evolution histories are unconstrained for these models. The stellar mass---stellar metallicity relations obtained from the two methods agree with the relation measured from individual supergiant stars in nearby galaxies. These relations are also consistent with the relation obtained from emission line analysis of gas-phase metallicity after accounting for systematic offsets in the gas-phase-metallicity. We measure dust attenuation of the stellar continuum and show that its dependence on stellar mass and star formation rate is consistent with previously reported results derived from nebular emission lines. However, stellar continuum attenuation is smaller than nebular emission line attenuation. The continuum-to-nebular attenuation ratio depends on stellar mass and is smaller in more massive galaxies. Our consistent analysis of stellar continuum and nebular emission lines paves the way for a comprehensive investigation of stellar metallicities of star-forming and quiescent galaxies.

The "kinematic" morphology-density relation for early-type galaxies posits that those galaxies with low angular momentum are preferentially found in the highest-density regions of the universe. We use a large sample of galaxy groups with halo masses 10^12.5 < M_halo < 10^14.5 M_sun/h observed with the Mapping Nearby Galaxies at APO (MaNGA) survey to examine whether there is a correlation between local environment and rotational support that is independent of stellar mass. We find no compelling evidence for a relationship between the angular momentum content of early-type galaxies and either local overdensity or radial position within the group at fixed stellar mass.

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

The physical process whereby a carbon--oxygen white dwarf explodes as a Type Ia supernova (SN Ia) remains highly uncertain. The degree of neutronization in SN Ia ejecta holds clues to this process because it depends on the mass and the metallicity of the stellar progenitor, and on the thermodynamic history prior to the explosion. We report on a new method to determine ejecta neutronization using Ca and S lines in the X-ray spectra of Type Ia supernova remnants (SNRs). Applying this method to \textit{Suzaku} data of Tycho, Kepler, 3C 397 and G337.2$-$0.7 in the Milky Way, and N103B in the Large Magellanic Cloud, we find that the neutronization of the ejecta in N103B is comparable to that of Tycho and Kepler, which suggests that progenitor metallicity is not the only source of neutronization in SNe Ia. We then use a grid of SN Ia explosion models to infer the metallicities of the stellar progenitors of our SNRs. The implied metallicities of 3C 397, G337.2$-$0.7, and N103B are major outliers compared to the local stellar metallicity distribution functions, indicating that progenitor metallicity can be ruled out as the origin of neutronization for these SNRs. Although the relationship between ejecta neutronization and equivalent progenitor metallicity is subject to uncertainties stemming from the $^{12}$C$\,$+$^{16}$O reaction rate, which affects the Ca/S mass ratio, our main results are not sensitive to these details.

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.

Diffuse Ionized Gas (DIG) is prevalent in star-forming galaxies. Using a sample of 365 nearly face-on star-forming galaxies observed by MaNGA, we demonstrate how DIG in star-forming galaxies impacts the measurements of emission line ratios, hence the interpretation of diagnostic diagrams and gas-phase metallicity measurements. At fixed metallicity, DIG-dominated low H\alpha\ surface brightness regions display enhanced [SII]/H\alpha, [NII]/H\alpha, [OII]/H\beta, and [OI]/H\alpha. The gradients in these line ratios are determined by metallicity gradients and H\alpha\ surface brightness. In line ratio diagnostic diagrams, contamination by DIG moves HII regions towards composite or LI(N)ER-like regions. A harder ionizing spectrum is needed to explain DIG line ratios. Leaky HII region models can only shift line ratios slightly relative to HII region models, and thus fail to explain the composite/LI(N)ER line ratios displayed by DIG. Our result favors ionization by evolved stars as a major ionization source for DIG with LI(N)ER-like emission. DIG can significantly bias the measurement of gas metallicity and metallicity gradients derived using strong-line methods. Metallicities derived using N2O2 are optimal because they exhibit the smallest bias and error. Using O3N2, R23, N2=[NII]/H\alpha, and N2S2H\alpha\ (Dopita et al. 2016) to derive metallicities introduces bias in the derived metallicity gradients as large as the gradient itself. The strong-line method of Blanc et al. (2015; IZI hereafter) cannot be applied to DIG to get an accurate metallicity because it currently contains only HII region models which fail to describe the DIG.

Chemical evolution models are powerful tools for interpreting stellar abundance surveys and understanding galaxy evolution. However, their predictions depend heavily on the treatment of inflow, outflow, star formation efficiency (SFE), the stellar initial mass function, the Type Ia supernova delay time distribution, stellar yields, and stellar population mixing. Using flexCE, a flexible one-zone chemical evolution code, we investigate the effects of and trade-offs between parameters. Two critical parameters are SFE and the outflow mass-loading parameter, which shift the knee in [O/Fe]-[Fe/H] and the equilibrium abundances that the simulations asymptotically approach, respectively. One-zone models with simple star formation histories follow narrow tracks in [O/Fe]-[Fe/H] unlike the observed bimodality (separate high-alpha and low-alpha sequences) in this plane. A mix of one-zone models with inflow timescale and outflow mass-loading parameter variations, motivated by the inside-out galaxy formation scenario with radial mixing, reproduces the two sequences better than a one-zone model with two infall epochs. We present [X/Fe]-[Fe/H] tracks for 20 elements assuming three different supernova yield models and find some significant discrepancies with solar neighborhood observations, especially for elements with strongly metallicity-dependent yields. We apply principal component abundance analysis (PCAA) to the simulations and existing data to reveal the main correlations amongst abundances and quantify their contributions to variation in abundance space. For the stellar population mixing scenario, the abundances of alpha-elements and elements with metallicity-dependent yields dominate the first and second principal components, respectively, and collectively explain 99% of the variance in the model. flexCE is a python package available at https://github.com/bretthandrews/flexCE.

We have conducted a study of extra-planar diffuse ionized gas using the first year data from the MaNGA IFU survey. We have stacked spectra from 49 edge-on, late-type galaxies as a function of distance from the midplane of the galaxy. With this technique we can detect the bright emission lines Halpha, Hbeta, [OII]3726, 3729, [OIII]5007, [NII]6549, 6584, and [SII]6717, 6731 out to about 4 kpc above the midplane. With 16 galaxies we can extend this analysis out to about 9 kpc, i.e. a distance of ~2R_e, vertically from the midplane. In the halo, the surface brightnesses of the [OII] and Halpha emission lines are comparable, unlike in the disk where Halpha dominates. When we split the sample by specific star formation rate, concentration index, and stellar mass, each subsample's emission line surface brightness profiles and ratios differ, indicating that extra-planar gas properties can vary. The emission line surface brightnesses of the gas around high specific star formation rate galaxies are higher at all distances, and the line ratios are closer to ratios characteristic of HII regions compared with low specific star formation rate galaxies. The less concentrated and lower stellar mass samples exhibit line ratios that are more like HII regions at larger distances than their more concentrated and higher stellar mass counterparts. The largest difference between different subsamples occurs when the galaxies are split by stellar mass. We additionally infer that gas far from the midplane in more massive galaxies has the highest temperatures and steepest radial temperature gradients based on their [NII]/Halpha and [OII]/Halpha ratios between the disk and the halo.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is an optical fiber-bundle integral-field unit (IFU) spectroscopic survey that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). With a spectral coverage of 3622 - 10,354 Angstroms and an average footprint of ~ 500 arcsec^2 per IFU the scientific data products derived from MaNGA will permit exploration of the internal structure of a statistically large sample of 10,000 low redshift galaxies in unprecedented detail. Comprising 174 individually pluggable science and calibration IFUs with a near-constant data stream, MaNGA is expected to obtain ~ 100 million raw-frame spectra and ~ 10 million reduced galaxy spectra over the six-year lifetime of the survey. In this contribution, we describe the MaNGA Data Reduction Pipeline (DRP) algorithms and centralized metadata framework that produces sky-subtracted, spectrophotometrically calibrated spectra and rectified 3-D data cubes that combine individual dithered observations. For the 1390 galaxy data cubes released in Summer 2016 as part of SDSS-IV Data Release 13 (DR13), we demonstrate that the MaNGA data have nearly Poisson-limited sky subtraction shortward of ~ 8500 Angstroms and reach a typical 10-sigma limiting continuum surface brightness mu = 23.5 AB/arcsec^2 in a five arcsec diameter aperture in the g band. The wavelength calibration of the MaNGA data is accurate to 5 km/s rms, with a median spatial resolution of 2.54 arcsec FWHM (1.8 kpc at the median redshift of 0.037) and a median spectral resolution of sigma = 72 km/s.

We present new analytic solutions for one-zone (fully mixed) chemical evolution models and explore their implications. In contrast to existing analytic models, we incorporate a realistic delay time distribution for Type Ia supernovae (SNIa) and can therefore track the separate evolution of $\alpha$-elements produced by core collapse supernovae (CCSNe) and iron peak elements synthesized in both CCSNe and SNIa. In generic cases, $\alpha$ and iron abundances evolve to an equilibrium at which element production is balanced by metal consumption and gas dilution, instead of continuing to increase over time. The equilibrium absolute abundances depend principally on supernova yields and the outflow mass loading parameter $\eta$, while the equilibrium abundance ratio [$\alpha$/Fe] depends mainly on yields and secondarily on star formation history. A stellar population can be metal-poor either because it has not yet evolved to equilibrium or because high outflow efficiency makes the equilibrium abundance itself low. Systems with ongoing gas accretion develop metallicity distribution functions (MDFs) that are sharply peaked, while "gas starved" systems with rapidly declining star formation have broadly peaked MDFs. A burst of star formation that consumes a significant fraction of a system's available gas can temporarily boost [$\alpha$/Fe] by 0.1-0.3 dex, a possible origin for rare, $\alpha$-enhanced stars with intermediate age or high metallicity. Other sudden transitions in system properties can produce surprising behavior, including backward evolution of a stellar population from high metallicity to low metallicity. An Appendix provides a user's guide for calculating enrichment histories, [$\alpha$/Fe] tracks, and MDFs for a wide variety of scenarios, including flexible forms of star formation history.

We use direct method oxygen abundances in combination with strong optical emission lines, stellar masses ($M_{\star}$), and star formation rates (SFRs) to recalibrate the N2, O3N2, and N2O2 oxygen abundance diagnostics. We stack spectra of $\sim$200,000 star-forming galaxies from the Sloan Digital Sky Survey in bins of $M_{\star}$ and SFR offset from the star forming main sequence to measure the weak emission lines needed to apply the direct method. All three new calibrations are reliable to within $\pm 0.10$ dex from $\log(M_{\star}/M_{\odot}) \sim 7.5 - 10.5$ and up to at least $200~M_{\odot}$ yr$^{-1}$ in SFR. The N2O2 diagnostic is the least subject to systematic biases. We apply the diagnostics to galaxies in the local universe and investigate the $M_{\star}$-$Z$-${\rm SFR}$ relation. The N2 and O3N2 diagnostics suggest the SFR dependence of the $M_{\star}$-$Z$-${\rm SFR}$ relation varies with both $M_{\star}$ and $\Delta \log(SSFR)$, whereas the N2O2 diagnostic suggests a nearly constant dependence on SFR. We apply our calibrations to a sample of high redshift galaxies from the literature, and find them to be metal poor relative to local galaxies with similar $M_{\star}$ and SFR. The calibrations do reproduce direct method abundances of the local analogs. We conclude that the $M_{\star}$-$Z$-${\rm SFR}$ relation evolves with redshift.

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.

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

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

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

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

We present Keck/MOSFIRE K-band spectroscopy of the first mass-selected sample of galaxies at $z\sim2.3$. Targets are selected from the 3D-HST Treasury survey. The six detected galaxies have a mean [NII]$\lambda$6584/H$\alpha$ ratio of $0.27\pm0.01$, with a small standard deviation of 0.05. This mean value is similar to that of UV-selected galaxies of the same mass. The mean gas-phase oxygen abundance inferred from the [NII]/H$\alpha$ ratios depends on the calibration method, and ranges from 12+log(O/H)$_{gas}=8.57$ for the {Pettini} & {Pagel} (2004) calibration to 12+log(O/H)$_{gas}= 8.87$ for the {Maiolino} {et~al.} (2008) calibration. Measurements of the stellar oxygen abundance in nearby quiescent galaxies with the same number density indicate 12+log(O/H)$_{stars}= 8.95$, similar to the gas-phase abundances of the $z\sim2.3$ galaxies if the {Maiolino} {et~al.} (2008) calibration is used. This suggests that these high-redshift star forming galaxies may be progenitors of today's massive early-type galaxies. The main uncertainties are the absolute calibration of the gas-phase oxygen abundance and the incompleteness of the $z\sim2.3$ sample: the galaxies with detected H$\alpha$ tend to be larger and have higher star formation rates than the galaxies without detected H$\alpha$, and we may still be missing the most dust-obscured progenitors.

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.

The relation between galaxy stellar mass and gas-phase metallicity is a sensitive diagnostic of the main processes that drive galaxy evolution, namely cosmological gas inflow, metal production in stars, and gas outflow via galactic winds. We employed the direct method to measure the metallicities of ~200,000 star-forming galaxies from the SDSS that were stacked in bins of (1) stellar mass and (2) both stellar mass and star formation rate (SFR) to significantly enhance the signal-to-noise ratio of the weak [O III] 4363 and [O II] 7320, 7330 auroral lines required to apply the direct method. These metallicity measurements span three decades in stellar mass from log(Mstar/Msun) = 7.4-10.5, which allows the direct method mass-metallicity relation to simultaneously capture the high-mass turnover and extend a full decade lower in mass than previous studies that employed more uncertain strong line methods. The direct method mass-metallicity relation rises steeply at low mass (O/H ~ Mstar^{1/2}) until it turns over at log(Mstar/Msun) = 8.9 and asymptotes to 12 + log(O/H) = 8.8 at high mass. The direct method mass-metallicity relation has a steeper slope, a lower turnover mass, and a factor of two to three greater dependence on SFR than strong line mass-metallicity relations. Furthermore, the SFR-dependence appears monotonic with stellar mass, unlike strong line mass-metallicity relations. We also measure the N/O abundance ratio, an important tracer of star formation history, and find the clear signature of primary and secondary nitrogen enrichment. N/O correlates tightly with oxygen abundance, and even more so with stellar mass.

Elemental abundance patterns can provide vital clues to the formation and enrichment history of a stellar population. Here we present an investigation of the Galactic bulge, where we apply principal component abundance analysis (PCAA)---a principal component decomposition of relative abundances [X/Fe]---to a sample of 35 microlensed bulge dwarf and subgiant stars, characterizing their distribution in the 12-dimensional space defined by their measured elemental abundances. The first principal component PC1, which suffices to describe the abundance patterns of most stars in the sample, shows a strong contribution from alpha-elements, reflecting the relative contributions of Type II and Type Ia supernovae. The second principal component PC2 is characterized by a Na--Ni correlation, the likely product of metallicity-dependent Type II supernova yields. The distribution in PC1 is bimodal, showing that the bimodality previously found in the [Fe/H] values of these stars is robustly and independently recovered by looking at only their relative abundance patterns. The two metal-rich stars that are alpha-enhanced have outlier values of PC2 and PC3, respectively, further evidence that they have distinctive enrichment histories. Applying PCAA to a sample of local thin and thick disk dwarfs yields a nearly identical PC1; in PC1, the metal-rich and metal-poor bulge dwarfs track kinematically selected thin and thick disk dwarfs, respectively, suggesting broadly similar alpha-enrichment histories. However, the disk PC2 is dominated by a Y--Ba correlation, likely indicating a greater contribution of s-process enrichment from long-lived asymptotic giant branch stars that is absent from the bulge PC2 because of its rapid formation.

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.

We evaluate radiation pressure from starlight on dust as a feedback mechanism in star-forming galaxies by comparing the luminosity and flux of star-forming systems to the dust Eddington limit. The linear LFIR--L'HCN correlation provides evidence that galaxies may be regulated by radiation pressure feedback. We show that star-forming galaxies approach but do not dramatically exceed Eddington, but many systems are significantly below Eddington, perhaps due to the "intermittency" of star formation. Better constraints on the dust-to-gas ratio and the CO- and HCN-to-H2 conversion factors are needed to make a definitive assessment of radiation pressure as a feedback mechanism.

Radiation pressure from the absorption and scattering of starlight by dust grains may be an important feedback mechanism in regulating star-forming galaxies. We compile data from the literature on star clusters, star-forming subregions, normal star-forming galaxies, and starbursts to assess the importance of radiation pressure on dust as a feedback mechanism, by comparing the luminosity and flux of these systems to their dust Eddington limit. This exercise motivates a novel interpretation of the Schmidt Law, the LIR-L'CO correlation, and the LIR-L'HCN correlation. In particular, the linear LIR-L'HCN correlation is a natural prediction of radiation pressure regulated star formation. Overall, we find that the Eddington limit sets a hard upper bound to the luminosity of any star-forming region. Importantly, however, many normal star-forming galaxies have luminosities significantly below the Eddington limit. We explore several explanations for this discrepancy, especially the role of "intermittency" in normal spirals - the tendency for only a small number of subregions within a galaxy to be actively forming stars at any moment because of the time-dependence of the feedback process and the luminosity evolution of the stellar population. If radiation pressure regulates star formation in dense gas, then the gas depletion timescale is 6 Myr, in good agreement with observations of the densest starbursts. Finally, we highlight the importance of observational uncertainties - namely, the dust-to-gas ratio and the CO-H2 and HCN-H2 conversion factors - that must be understood before a definitive assessment of radiation pressure as a feedback mechanism in star-forming galaxies.