We present the first asymmetric drift (AD) measurements for unresolved stellar populations of different characteristic ages above and below 1.5 Gyr. These measurements sample the age-velocity relation (AVR) in galaxy disks. In this first paper we develop two efficient algorithms to extract AD on a spaxel-by-spaxel basis from optical integral-field spectroscopic (IFS) data-cubes. The algorithms apply different spectral templates, one using simple stellar populations and the other a stellar library; their comparison allows us to assess systematic errors in derived multi-component velocities, such as template-mismatch. We test algorithm reliability using mock spectra and Monte Carlo Markov Chains on real data from the MaNGA survey in SDSS-IV. We quantify random and systematic errors in AD as a function of signal-to-noise and stellar population properties with the aim of applying this technique to large subsets of the MaNGA galaxy sample. As a demonstration of our methods, we apply them to an initial sample of seven galaxies with comparable stellar mass and color to the Milky Way. We find a wide range of distinct AD radial profiles for young and old stellar populations.

PypeIt is a Python package for semi-automated reduction of astronomical, spectroscopic data. Its algorithms build on decades-long development of previous data reduction pipelines by the developers (Bernstein, Burles, & Prochaska, 2015; Bochanski et al., 2009). The reduction procedure -- including a complete list of the input parameters and available functionality -- is provided as online documentation hosted by Read the Docs, which is regularly updated. (https://pypeit.readthedocs.io/en/latest/). Release v1.0.3 serves the following spectrographs: Gemini/GNIRS, Gemini/GMOS, Gemini/FLAMINGOS 2, Lick/Kast, Magellan/MagE, Magellan/Fire, MDM/OSMOS, Keck/DEIMOS (600ZD, 830G, 1200G), Keck/LRIS, Keck/MOSFIRE (J and Y gratings tested), Keck/NIRES, Keck/NIRSPEC (low-dispersion), LBT/Luci-I, Luci-II, LBT/MODS (beta), NOT/ALFOSC (grism4), VLT/X-Shooter (VIS, NIR), VLT/FORS2 (300I, 300V), WHT/ISIS.

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.

SDSS-IV MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is the largest integral-field spectroscopy survey to date, aiming to observe a statistically representative sample of 10,000 low-redshift galaxies. In this paper we study the reliability of the emission-line fluxes and kinematic properties derived by the MaNGA Data Analysis Pipeline (DAP). We describe the algorithmic choices made in the DAP with regards to measuring emission-line properties, and the effect of our adopted strategy of simultaneously fitting the continuum and line emission. The effect of random errors are quantified by studying various fit-quality metrics, idealized recovery simulations and repeat observations. This analysis demonstrates that the emission lines are well-fit in the vast majority of the MaNGA dataset and the derived fluxes and errors are statistically robust. The systematic uncertainty on emission-line properties introduced by the choice of continuum templates is also discussed. In particular, we test the effect of using different stellar libraries and simple stellar-population models on the derived emission-line fluxes and the effect of introducing different tying prescriptions for the emission-line kinematics. We show that these effects can generate large ($>$ 0.2 dex) discrepancies at low signal-to-noise and for lines with low equivalent width (EW); however, the combined effect is noticeable even for H$\alpha$ EW $>$ 6~\AA. We provide suggestions for optimal use of the data provided by SDSS data release 15 and propose refinements on the \DAP\ for future MaNGA data releases.

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.

High-multiplex and deep spectroscopic follow-up of upcoming panoramic deep-imaging surveys like LSST, Euclid, and WFIRST is a widely recognized and increasingly urgent necessity. No current or planned facility at a U.S. observatory meets the sensitivity, multiplex, and rapid-response time needed to exploit these future datasets. FOBOS, the Fiber-Optic Broadband Optical Spectrograph, is a near-term fiber-based facility that addresses these spectroscopic needs by optimizing depth over area and exploiting the aperture advantage of the existing 10m Keck II Telescope. The result is an instrument with a uniquely blue-sensitive wavelength range (0.31-1.0 um) at R~3500, high-multiplex (1800 fibers), and a factor 1.7 greater survey speed and order-of-magnitude greater sampling density than Subaru's Prime Focus Spectrograph (PFS). In the era of panoramic deep imaging, FOBOS will excel at building the deep, spectroscopic reference data sets needed to interpret vast imaging data. At the same time, its flexible focal plane, including a mode with 25 deployable integral-field units (IFUs) across a 20 arcmin diameter field, enables an expansive range of scientific investigations. Its key programmatic areas include (1) nested stellar-parameter training sets that enable studies of the Milky Way and M31 halo sub-structure, as well as local group dwarf galaxies, (2) a comprehensive picture of galaxy formation thanks to detailed mapping of the baryonic environment at z~2 and statistical linking of evolving populations to the present day, and (3) dramatic enhancements in cosmological constraints via precise photometric redshifts and determined redshift distributions. In combination with Keck I instrumentation, FOBOS also provides instant access to medium-resolution spectroscopy for transient sources with full coverage from the UV to the K-band.

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.

The exact nature of the arms of spiral galaxies is still an open question. It has been widely assumed that spiral arms in galaxies with two distinct symmetrical arms are the products of density waves that propagate around the disk, with the spiral arms being visibly enhanced by the star formation that is triggered as the passing wave compresses gas in the galaxy disk. Such a persistent wave would propagate with an approximately constant angular speed, its pattern speed Op. The quasi-stationary density wave theory can be tested by measuring this quantity and showing that it does not vary with radius in the galaxy. Unfortunately, this measurement is difficult because Op is only indirectly connected to observables such as the stellar rotation speed. Here, we use the detailed information on stellar populations of the grand-design spiral galaxy UGC 3825, extracted from spectral mapping, to measure the offset between young stars of a known age and the spiral arm in which they formed, allowing the first direct measure of Op at a range of radii. The offset in this galaxy is found to be as expected for a pattern speed that varies little with radius, indicating consistency with a quasi-stationary density wave, and lending credence to this new method.

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.

A new class of quiescent galaxies harboring possible AGN-driven winds has been discovered using spatially resolved optical spectroscopy from the ongoing SDSS-IV MaNGA survey. These galaxies, termed "red geysers", constitute $5-10\%$ of the local quiescent population and are characterized by narrow bisymmetric patterns in ionized gas emission features. Cheung et al. argued that these galaxies host large-scale AGN-driven winds that may play a role in suppressing star formation at late times. In this work, we test the hypothesis that AGN activity is ultimately responsible for the red geyser phenomenon. We compare the nuclear radio activity of the red geysers to a matched control sample with similar stellar mass, redshift, rest frame $NUV-r$ color, axis ratio and presence of ionized gas. We have used the 1.4 GHz radio continuum data from VLA FIRST survey to stack the radio flux from the red geyser and control samples. In addition to a 3 times higher FIRST detection rate, we find that red geysers have a 5$\sigma$ higher level of average radio flux than control galaxies. After restricting to rest-frame $NUV-r$ color $>$ 5 and checking mid-IR WISE photometry, we rule out star formation contamination and conclude that red geysers are associated with more active AGN. Red geysers and a possibly-related class with disturbed H$\alpha$ emission account for 40\% of all radio-detected red galaxies with $\rm log~(M_\star/M_\odot) < 11$. Our results support a picture in which episodic AGN activity drives large-scale-relatively weak ionized winds that may provide a feedback mechanism for many early-type galaxies.

We show how the Radon transform (defined as a series of line integrals through an image at different orientations and offsets from the origin) can be used as a simple, non-parametric tool to characterize galaxy velocity fields, specifically their global kinematic position angles (PA_k) and any radial variation or asymmetry in PA_k. This method is fast and easily automated, making it particularly beneficial in an era where IFU and interferometric surveys are yielding samples of thousands of galaxies. We demonstrate the Radon transform by applying it to gas and stellar velocity fields from the first ~2800 galaxies of the SDSS-IV MaNGA IFU survey. We separately classify gas and stellar velocity fields into five categories based on the shape of their radial PA_k profiles. At least half of stellar velocity fields and two-thirds of gas velocity fields are found to show detectable deviations from uniform coplanar circular motion, although most of these variations are symmetric about the center of the galaxy. The behavior of gas and stellar velocity fields is largely independent, even when PA_k profiles for both components are measured over the same radii. We present evidence that one class of symmetric PA_k variations is likely associated with bars and/or oval distortions, while another class is more consistent with warped disks. This analysis sets the stage for more in-depth future studies which explore the origin of diverse kinematic behavior in the galaxy population.

Ground-layer adaptive optics (GLAO) systems offer the possibility of improving the "seeing" of large ground-based telescopes and increasing the efficiency and sensitivity of observations over a wide field-of-view. We explore the utility and feasibility of deploying a GLAO system at the W. M. Keck Observatory in order to feed existing and future multi-object spectrographs and wide-field imagers. We also briefly summarize science cases spanning exoplanets to high-redshift galaxy evolution that would benefit from a Keck GLAO system. Initial simulations indicate that a Keck GLAO system would deliver a 1.5x and 2x improvement in FWHM at optical (500 nm) and infrared (1.5 micron), respectively. The infrared instrument, MOSFIRE, is ideally suited for a Keck GLAO feed in that it has excellent image quality and is on the telescope's optical axis. However, it lacks an atmospheric dispersion compensator, which would limit the minimum usable slit size for long-exposure science cases. Similarly, while LRIS and DEIMOS may be able to accept a GLAO feed based on their internal image quality, they lack either an atmospheric dispersion compensator (DEIMOS) or flexure compensation (LRIS) to utilize narrower slits matched to the GLAO image quality. However, some science cases needing shorter exposures may still benefit from Keck GLAO and we will investigate the possibility of installing an ADC.

Einstein's theory of gravity, General Relativity, has been precisely tested on Solar System scales, but the long-range nature of gravity is still poorly constrained. The nearby strong gravitational lens, ESO 325-G004, provides a laboratory to probe the weak-field regime of gravity and measure the spatial curvature generated per unit mass, $\gamma$. By reconstructing the observed light profile of the lensed arcs and the observed spatially resolved stellar kinematics with a single self-consistent model, we conclude that $\gamma = 0.97 \pm 0.09$ at 68% confidence. Our result is consistent with the prediction of 1 from General Relativity and provides a strong extragalactic constraint on the weak-field metric of gravity.

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 measure $\lambda_{R_e}$, a proxy for galaxy specific stellar angular momentum within one effective radius, and the ellipticity, $\epsilon$, for about 2300 galaxies of all morphological types observed with integral field spectroscopy as part of the MaNGA survey, the largest such sample to date. We use the $(\lambda_{R_e}, \epsilon)$ diagram to separate early-type galaxies into fast and slow rotators. We also visually classify each galaxy according to its optical morphology and two-dimensional stellar velocity field. Comparing these classifications to quantitative $\lambda_{R_e}$ measurements reveals tight relationships between angular momentum and galaxy structure. In order to account for atmospheric seeing, we use realistic models of galaxy kinematics to derive a general approximate analytic correction for $\lambda_{R_e}$. Thanks to the size of the sample and the large number of massive galaxies, we unambiguously detect a clear bimodality in the $(\lambda_{R_e}, \epsilon)$ diagram which may result from fundamental differences in galaxy assembly history. There is a sharp secondary density peak inside the region of the diagram with low $\lambda_{R_e}$ and $\epsilon < 0.4$, previously suggested as the definition for slow rotators. Most of these galaxies are visually classified as non-regular rotators and have high velocity dispersion. The intrinsic bimodality must be stronger, as it tends to be smoothed by noise and inclination. The large sample of slow rotators allows us for the first time to unveil a secondary peak at +/-90 degrees in their distribution of the misalignments between the photometric and kinematic position angles. We confirm that genuine slow rotators start appearing above a stellar mass of 2\times10^{11} M_{\odot}$ where a significant number of high-mass fast rotators also exist.

We study the spatially resolved star formation of 1494 galaxies in the SDSSIV-MaNGA Survey. SFRs are calculated using a two-step process, using $H_\alpha$ in star forming regions and $D_n4000$ in regions identified as AGN/LI(N)ER or lineless. The roles of secular and environmental quenching processes are investigated by studying the dependence of the radial profiles of specific star formation rate on stellar mass, galaxy structure and environment. We report on the existence of `Centrally Suppressed' galaxies, which have suppressed SSFR in their cores compared to their disks. The profiles of centrally suppressed and unsuppressed galaxies are distibuted in a bimodal way. Galaxies with high stellar mass and core velocity dispersion are found to be much more likely to be centrally suppressed than low mass galaxies, and we show that this is related to morphology and the presence of AGN/LI(N)ER like emission. Centrally suppressed galaxies also display lower star formation at all radii compared to unsuppressed galaxies. The profiles of central and satellite galaxies are also compared, and we find that satellite galaxies experience lower specific star formation rates at all radii than central galaxies. This uniform suppression could be a signal of the stripping of hot halo gas in the process known as strangulation. We find that satellites are not more likely to be suppressed in their cores than centrals, indicating that the core suppression is an entirely internal process. We find no correlation between the local environment density and the profiles of star formation rate surface density.

Regions of disc galaxies with widespread star formation tend to be both gravitationally unstable and self-shielded against ionizing radiation, whereas extended outer discs with little or no star formation tend to be stable and unshielded on average. We explore what drives the transition between these two regimes, specifically whether discs first meet the conditions for self-shielding (parameterized by dust optical depth, $\tau$) or gravitational instability (parameterized by a modified version of Toomre's instability parameters, $Q_{\rm thermal}$, which quantifies the stability of a gas disc that is thermally supported at $T=10^4$ K). We first introduce a new metric formed by the product of these quantities, $Q_{\rm thermal}\tau$, which indicates whether the conditions for disk instability or self-shielding are easier to meet in a given region of a galaxy, and we discuss how $Q_{\rm thermal}\tau$ can be constrained even in the absence of direct gas information. We then analyse a sample of 13 galaxies with resolved gas measurements and find that on average galaxies will reach the threshold for disk instabilities ($Q_{\rm thermal}<1$) before reaching the threshold for self-shielding ($\tau>1$). Using integral field spectroscopic observations of a sample of 236 galaxies from the MaNGA survey, we find that the value of $Q_{\rm thermal}\tau$ in star-forming discs is consistent with similar behavior. These results support a scenario where disc fragmentation and collapse occurs before self-shielding, suggesting that gravitational instabilities are the primary condition for widespread star formation in galaxy discs. Our results support similar conclusions based on recent galaxy simulations.

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.

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 study the gas phase metallicity (O/H) and nitrogen abundance gradients traced by star forming regions in a representative sample of 550 nearby galaxies in the stellar mass range $\rm 10^9-10^{11.5} M_\odot$ with resolved spectroscopic data from the SDSS-IV MaNGA survey. Using strong-line ratio diagnostics (R23 and O3N2 for metallicity and N2O2 for N/O) and referencing to the effective (half-light) radius ($\rm R_e$), we find that the metallicity gradient steepens with stellar mass, lying roughly flat among galaxies with $\rm log(M_\star/M_\odot) = 9.0$ but exhibiting slopes as steep as -0.14 dex $\rm R_e^{-1}$ at $\rm log(M_\star/M_\odot) = 10.5$ (using R23, but equivalent results are obtained using O3N2). At higher masses, these slopes remain typical in the outer regions of our sample ($\rm R > 1.5 ~R_e$), but a flattening is observed in the central regions ($\rm R < 1~ R_e$). In the outer regions ($\rm R > 2.0 ~R_e$) we detect a mild flattening of the metallicity gradient in stacked profiles, although with low significance. The N/O ratio gradient provides complementary constraints on the average chemical enrichment history. Unlike the oxygen abundance, the average N/O profiles do not flatten out in the central regions of massive galaxies. The metallicity and N/O profiles both depart significantly from an exponential form, suggesting a disconnect between chemical enrichment and stellar mass surface density on local scales. In the context of inside-out growth of discs, our findings suggest that central regions of massive galaxies today have evolved to an equilibrium metallicity, while the nitrogen abundance continues to increase as a consequence of delayed secondary nucleosynthetic production.

Using spatially resolved spectroscopy from SDSS-IV MaNGA we have demonstrated that low ionisation emission line regions (LIERs) in local galaxies result from photoionisation by hot evolved stars, not active galactic nuclei. LIERs are ubiquitous in both quiescent galaxies and in the central regions of galaxies where star formation takes place at larger radii. We refer to these two classes of galaxies as extended LIER (eLIER) and central LIER (cLIER) galaxies respectively. cLIERs are late type galaxies located around the green valley, in the transition region between the star formation main sequence and quiescent galaxies. These galaxies display regular disc rotation in both stars and gas, although featuring a higher central stellar velocity dispersion than star forming galaxies of the same mass. cLIERs are consistent with being slowly quenched inside-out; the transformation is associated with massive bulges, pointing towards the importance of bulge growth via secular evolution. eLIERs are morphologically early types and are indistinguishable from passive galaxies devoid of line emission in terms of their stellar populations, morphology and central stellar velocity dispersion. Ionised gas in eLIERs shows both disturbed and disc-like kinematics. When a large-scale flow/rotation is observed in the gas, it is often misaligned relative to the stellar component. These features indicate that eLIERs are passive galaxies harbouring a residual cold gas component, acquired mostly via external accretion. Importantly, quiescent galaxies devoid of line emission reside in denser environments and have significantly higher satellite fraction than eLIERs. Environmental effects thus represent the likely cause for the existence of line-less galaxies on the red sequence.

Spiral structure is the most distinctive feature of disk galaxies and yet debate persists about which theory of spiral structure is the correct one. Many versions of the density wave theory demand that the pitch angle is uniquely determined by the distribution of mass in the bulge and disk of the galaxy. We present evidence that the tangent of the pitch angle of logarithmic spiral arms in disk galaxies correlates strongly with the density of neutral atomic hydrogen in the disk and with the central stellar bulge mass of the galaxy. These three quantities, when plotted against each other, form a planar relationship which, we argue should be fundamental to our understanding of spiral structure in disk galaxies. We further argue that any successful theory of spiral structure must be able to explain this relationship.

We study the internal radial gradients of stellar population properties within $1.5\;R_{\rm e}$ and analyse the impact of galaxy environment. We use a representative sample of 721 galaxies with masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV survey MaNGA. We split this sample by morphology into early-type and late-type galaxies. Using the full spectral fitting code FIREFLY, we derive the light and mass-weighted stellar population properties age and metallicity, and calculate the gradients of these properties. We use three independent methods to quantify galaxy environment, namely the $N^{th}$ nearest neighbour, the tidal strength parameter $Q$ and distinguish between central and satellite galaxies. In our analysis, we find that early-type galaxies generally exhibit shallow light-weighted age gradients in agreement with the literature and mass-weighted median age gradients tend to be slightly positive. Late-type galaxies, instead, have negative light-weighted age gradients. We detect negative metallicity gradients in both early and late-type galaxies that correlate with galaxy mass, with the gradients being steeper and the correlation with mass being stronger in late-types. We find, however, that stellar population gradients, for both morphological classifications, have no significant correlation with galaxy environment for all three characterisations of environment. Our results suggest that galaxy mass is the main driver of stellar population gradients in both early and late-type galaxies, and any environmental dependence, if present at all, must be very subtle.

We study the internal gradients of stellar population properties within $1.5\;R_{\rm e}$ for a representative sample of 721 galaxies with stellar masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV MaNGA IFU survey. Through the use of our full spectral fitting code FIREFLY, we derive light and mass-weighted stellar population properties and their radial gradients, as well as full star formation and metal enrichment histories. We also quanfify the impact that different stellar population models and full spectral fitting routines have on the derived stellar population properties, and the radial gradient measurements. In our analysis, we find that age gradients tend to be shallow for both early-type and late-type galaxies. {\em Mass-weighted} age gradients of early-types are positive ($\sim 0.09\; {\rm dex}/R_{\rm e}$) pointing to "outside-in" progression of star formation, while late-type galaxies have negative {\em light-weighted} age gradients ($\sim -0.11\; {\rm dex}/R_{\rm e}$), suggesting an "inside-out" formation of discs. We detect negative metallicity gradients in both early and late-type galaxies, but these are significantly steeper in late-types, suggesting that radial dependence of chemical enrichment processes and the effect of gas inflow and metal transport are far more pronounced in discs. Metallicity gradients of both morphological classes correlate with galaxy mass, with negative metallicity gradients becoming steeper with increasing galaxy mass. The correlation with mass is stronger for late-type galaxies, with a slope of $d(\nabla [Z/H])/d(\log M)\sim -0.2\pm 0.05\;$, compared to $d(\nabla [Z/H])/d(\log M)\sim -0.05\pm 0.05\;$ for early-types. This result suggests that the merger history plays a relatively small role in shaping metallicity gradients of galaxies.

Galaxies grow through both internal and external processes. In about 10% of nearby red galaxies with little star formation, gas and stars are counter-rotating, demonstrating the importance of external gas acquisition in these galaxies. However, systematic studies of such phenomena in blue, star-forming galaxies are rare, leaving uncertain the role of external gas acquisition in driving evolution of blue galaxies. Based on new measurements with integral field spectroscopy of a large representative galaxy sample, we find an appreciable fraction of counter-rotators among blue galaxies (9 out of 489 galaxies). The central regions of blue counter-rotators show younger stellar populations and more intense, ongoing star formation than their outer parts, indicating ongoing growth of the central regions. The result offers observational evidence that the acquisition of external gas in blue galaxies is possible; the interaction with pre-existing gas funnels the gas into nuclear regions (< 1 kpc) to form new stars.

We study the properties of 66 galaxies with kinematically misaligned gas and stars from MaNGA survey. The fraction of kinematically misaligned galaxies varies with galaxy physical parameters, i.e. M*, SFR and sSFR. According to their sSFR, we further classify these 66 galaxies into three categories, 10 star-forming, 26 "Green Valley" and 30 quiescent ones. The properties of different types of kinematically misaligned galaxies are different in that the star-forming ones have positive gradient in D4000 and higher gas-phase metallicity, while the green valley/quiescent ones have negative D4000 gradients and lower gas-phase metallicity on average. There is evidence that all types of the kinematically misaligned galaxies tend to live in more isolated environment. Based on all these observational results, we propose a scenario for the formation of star forming galaxies with kinematically misaligned gas and stars - the progenitor accretes misaligned gas from a gas-rich dwarf or cosmic web, the cancellation of angular momentum from gas-gas collisions between the pre-existing gas and the accreted gas largely accelerates gas inflow, leading to fast centrally-concentrated star-formation. The higher metallicity is due to enrichment from this star formation. For the kinematically misaligned green valley and quiescent galaxies, they might be formed through gas-poor progenitors accreting kinematically misaligned gas from satellites which are smaller in mass.

Using kinematic maps from the Sloan Digital Sky Survey (SDSS) Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we reveal that the majority of low-mass quenched galaxies exhibit coherent rotation in their stellar kinematics. Our sample includes all 39 quenched low-mass galaxies observed in the first year of MaNGA. The galaxies are selected with $M_{r} > -19.1$, stellar masses $10^{9}$ M$_{\odot} < M_{\star} < 5\times10^{9}$ M$_{\odot}$, EW$_{H\alpha} <2$ \AA, and all have red colours $(u-r)>1.9$. They lie on the size-magnitude and $\sigma$-luminosity relations for previously studied dwarf galaxies. Just six ($15\pm5.7$ per cent) are found to have rotation speeds $v_{e,rot} < 15$ km s$^{-1}$ at $\sim1$ $R_{e}$, and may be dominated by pressure support at all radii. Two galaxies in our sample have kinematically distinct cores in their stellar component, likely the result of accretion. Six contain ionised gas despite not hosting ongoing star formation, and this gas is typically kinematically misaligned from their stellar component. This is the first large-scale Integral Field Unit (IFU) study of low mass galaxies selected without bias against low-density environments. Nevertheless, we find the majority of these galaxies are within $\sim1.5$ Mpc of a bright neighbour ($M_{K} < -23$; or M$_{\star} > 5\times10^{10}$ M$_{\odot}$), supporting the hypothesis that galaxy-galaxy or galaxy-group interactions quench star formation in low-mass galaxies. The local bright galaxy density for our sample is $\rho_{proj} = 8.2\pm2.0$ Mpc$^{-2}$, compared to $\rho_{proj} = 2.1\pm0.4$ Mpc$^{-2}$ for a star forming comparison sample, confirming that the quenched low mass galaxies are preferentially found in higher density environments.

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.

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.

Quiescent galaxies with little or no ongoing star formation dominate the galaxy population above $M_{*}\sim 2 \times 10^{10}~M_{\odot}$, where their numbers have increased by a factor of $\sim25$ since $z\sim2$. Once star formation is initially shut down, perhaps during the quasar phase of rapid accretion onto a supermassive black hole, an unknown mechanism must remove or heat subsequently accreted gas from stellar mass loss or mergers that would otherwise cool to form stars. Energy output from a black hole accreting at a low rate has been proposed, but observational evidence for this in the form of expanding hot gas shells is indirect and limited to radio galaxies at the centers of clusters, which are too rare to explain the vast majority of the quiescent population. Here we report bisymmetric emission features co-aligned with strong ionized gas velocity gradients from which we infer the presence of centrally-driven winds in typical quiescent galaxies that host low-luminosity active nuclei. These galaxies are surprisingly common, accounting for as much as $10\%$ of the population at $M_* \sim 2 \times 10^{10}~ M_{\odot}$. In a prototypical example, we calculate that the energy input from the galaxy's low-level active nucleus is capable of driving the observed wind, which contains sufficient mechanical energy to heat ambient, cooler gas (also detected) and thereby suppress star formation.

We study the spatially resolved excitation properties of the ionised gas in a sample of 646 galaxies using integral field spectroscopy data from SDSS-IV MaNGA. Making use of Baldwin-Philips-Terlevich diagnostic diagrams we demonstrate the ubiquitous presence of extended (kpc scale) low ionisation emission-line regions (LIERs) in both star forming and quiescent galaxies. In star forming galaxies LIER emission can be associated with diffuse ionised gas, most evident as extra-planar emission in edge-on systems. In addition, we identify two main classes of galaxies displaying LIER emission: `central LIER' (cLIER) galaxies, where central LIER emission is spatially extended, but accompanied by star formation at larger galactocentric distances, and `extended LIER' (eLIER) galaxies, where LIER emission is extended throughout the whole galaxy. In eLIER and cLIER galaxies, LIER emission is associated with radially flat, low H$\alpha$ equivalent width of line emission ($<$ 3 \AA) and stellar population indices demonstrating the lack of young stellar populations, implying that line emission follows tightly the continuum due to the underlying old stellar population. The H$\alpha$ surface brightness radial profiles are always shallower than $\rm 1/r^{2}$ and the line ratio [OIII]$\lambda$5007/[OII]$\lambda$3727,29 (a tracer of the ionisation parameter of the gas) shows a flat gradient. This combined evidence strongly supports the scenario in which LIER emission is not due to a central point source but to diffuse stellar sources, the most likely candidates being hot, evolved (post-asymptotic giant branch) stars. Shocks are observed to play a significant role in the ionisation of the gas only in rare merging and interacting systems.

We present results from 21 cm radio synthesis imaging of 28 spiral galaxies from the DiskMass Survey obtained with the VLA, WSRT, and GMRT facilities. We detail the observations and data reduction procedures and present a brief analysis of the radio data. We construct 21 cm continuum images, global HI emission-line profiles, column-density maps, velocity fields, and position-velocity diagrams. From these we determine star formation rates (SFRs), HI line widths, total HI masses, rotation curves, and azimuthally-averaged radial HI column-density profiles. All galaxies have an HI disk that extends beyond the readily observable stellar disk, with an average ratio and scatter of R_{HI}/R_{25}=1.35+/-0.22, and a majority of the galaxies appear to have a warped HI disk. A tight correlation exists between total HI mass and HI diameter, with the largest disks having a slightly lower average column density. Galaxies with relatively large HI disks tend to exhibit an enhanced stellar velocity dispersion at larger radii, suggesting the influence of the gas disk on the stellar dynamics in the outer regions of disk galaxies. We find a striking similarity among the radial HI surface density profiles, where the average, normalized radial profile of the late-type spirals is described surprisingly well with a Gaussian profile. These results can be used to estimate HI surface density profiles in galaxies that only have a total HI flux measurement. We compare our 21 cm radio continuum luminosities with 60 micron luminosities from IRAS observations for a subsample of 15 galaxies and find that these follow a tight radio-infrared relation, with a hint of a deviation from this relation at low luminosities. We also find a strong correlation between the average SFR surface density and the K-band surface brightness of the stellar disk.

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

We present the correlation between the extrapolated central disk surface brightness (mu) and extrapolated central surface mass density (Sigma) for galaxies in the DiskMass sample. This mu-Sigma-relation has a small scatter of 30% at the high-surface-brightness (HSB) end. At the low surface brightness (LSB) end, galaxies fall above the mu-Sigma-relation, which we attribute to their higher dark matter content. After correcting for the dark matter, as well as for the contribution of gas and the effects of radial gradients in the disk, the LSB end falls back on the linear mu-Sigma-relation. The resulting scatter about the corrected mu-Sigma-relation is 25% at the HSB end, and about 50% at the LSB end. The intrinsic scatter in the mu-Sigma-relation is estimated to be 10% to 20%. Thus, if the surface brightness is known, the stellar surface mass density is known to within 10-20% (random error). Assuming disks have an exponential vertical distribution of mass, the average (M_L)_K is 0.24 Msun/Lsun, with an intrinsic scatter around the mean of at most 0.05 Msun/Lsun. This value for (M/L)_K is 20% smaller than we found in Martinsson et al., mainly due to the correction for dark matter applied here. This small scatter means that among the galaxies in our sample variations in scale height, vertical density profile shape, and/or the ratio of vertical over radial velocity dispersion must be small.

We study the relationship between the stability level of late-type galaxy disks and their star-formation activity using integral-field gaseous and stellar kinematic data. Specifically, we compare the two-component (gas+stars) stability parameter from Romeo & Wiegert (Q_RW), incorporating stellar kinematic data for the first time, and the star-formation rate estimated from 21cm continuum emission. We determine the stability level of each disk probabilistically using a Bayesian analysis of our data and a simple dynamical model. Our method incorporates the shape of the stellar velocity ellipsoid (SVE) and yields robust SVE measurements for over 90% of our sample. Averaging over this subsample, we find a meridional shape of sigma_z/sigma_R = 0.51^{+0.36}_{-0.25} for the SVE and, at 1.5 disk scale lengths, a stability parameter of Q_RW = 2.0 +/- 0.9. We also find that the disk-averaged star-formation-rate surface density (Sigma-dot_e,*) is correlated with the disk-averaged gas and stellar mass surface densities (Sigma_e,g and Sigma_e,*) and anti-correlated with Q_RW. We show that an anti-correlation between Sigma-dot_e,* and Q_RW can be predicted using empirical scaling relations, such that this outcome is consistent with well-established statistical properties of star-forming galaxies. Interestingly, Sigma-dot_e,* is not correlated with the gas-only or star-only Toomre parameters, demonstrating the merit of calculating a multi-component stability parameter when comparing to star-formation activity. Finally, our results are consistent with the Ostriker et al. model of self-regulated star-formation, which predicts Sigma-dot_e,*/Sigma_e,g/sqrt(Sigma_e,*). Based on this and other theoretical expectations, we discuss the possibility of a physical link between disk stability level and star-formation rate in light of our empirical results.

We calculate the stellar surface mass density (Sigma_*) and two-component (gas+stars) disk stability (Q_RW) for 25 late-type galaxies from the DiskMass Survey. These calculations are based on fits of a dynamical model to our ionized-gas and stellar kinematic data performed using a Markov Chain Monte Carlo sampling of the Bayesian posterior. Marginalizing over all galaxies, we find a median value of Q_RW=2.0+/-0.9 at 1.5 scale lengths. We also find that Q_RW is anti-correlated with the star-formation rate surface density (Sigma_SFR), which can be predicted using a closed set of empirical scaling relations. Finally, we find that the star-formation efficiency (Sigma_SFR/Sigma_g) is correlated with Sigma_* and weakly anti-correlated with Q_RW. The former is consistent with an equilibrium prediction of Sigma_SFR/Sigma_g propto Sigma_*^{1/2}. Despite its order-of-magnitude range, we find no correlation of Sigma_SFR/Sigma_g/Sigma_*^{1/2} with any other physical quantity derived by our study.

We present dynamically-determined rotation-curve mass decompositions of 30 spiral galaxies, which were carried out to test the maximum-disk hypothesis and to quantify properties of their dark-matter (DM) halos. We used measured vertical velocity dispersions of the disk stars to calculate dynamical mass surface densities. Together with our atomic and molecular gas mass surface densities, we derived the stellar mass surface densities, and thus have absolute measurements of all dominant baryonic components. Using K-band surface brightness profiles, we calculated the K-band mass-to-light ratio of the stellar disks (M/L). Our result is consistent with all galaxies in the sample having equal M/L, with a sample average and scatter of <M/L>=0.31+/-0.07. Rotation-curves of the baryonic components were calculated from their mass surface densities, and used with circular-speed measurements to derive the structural parameters of the DM halos, modeled as either a pseudo-isothermal sphere (pISO) or an NFW halo. All galaxies in our sample are submaximal, such that at 2.2 disk scale lengths (hR) the ratios between the baryonic and total rotation-curves (Fb^{2.2hR}) are less than 0.75. We find this ratio to be nearly constant between 1-6 hR within individual galaxies. We find a sample average and scatter of <Fb^{2.2hR}>=0.57+/-0.07, with trends of larger Fb^{2.2hR} for more luminous and higher-surface-brightness galaxies. To enforce these being maximal, we need to scale M/L by a factor 3.6 on average. The DM rotation curves are marginally better fit by a pISO than by an NFW halo. For the nominal-M/L (submaximal) case, the derived NFW-halo parameters have values consistent with LCDM N-body simulations, suggesting that the baryonic matter has only had a minor effect on the DM distribution. In contrast, maximum-M/L decompositions yield halo concentrations that are too low compared to the LCDM simulations.

We present ionized-gas (OIII) and stellar kinematics (velocities and velocity dispersions) for 30 nearly face-on spiral galaxies out to as much as three disk scale lengths (h_R). These data have been derived from PPak IFU spectroscopy (4980-5370A), observed at a mean resolution of R=7700 (sigma_inst=17km/s). These data are a fundamental product of our survey and will be used in companion papers to, e.g., derive the detailed (baryonic+dark) mass budget of each galaxy in our sample. Our presentation provides a comprehensive description of the observing strategy, data reduction, and analysis. Along with a clear presentation of the data, we demonstrate: (1) The OIII and stellar rotation curves exhibit a clear signature of asymmetric drift with a rotation difference that is 11% of the maximum rotation speed of the galaxy disk, comparable to measurements in the solar neighborhood in the Milky Way. (2) The e-folding length of the stellar velocity dispersion is two times h_R on average, as expected for a disk with a constant scale height and mass-to-light ratio, with a scatter that is notably smaller for massive, high-surface-brightness disks in the most luminous galaxies. (3) At radii larger than 1.5 h_R, the stellar velocity dispersion tends to decline slower than the best-fitting exponential function, which may be due to an increase in the disk mass-to-light ratio, disk flaring, or disk heating by the dark-matter halo. (4) A strong correlation exists between the central vertical stellar velocity dispersion of the disks and their circular rotational speed at 2.2 h_R, with a zero point indicating that galaxy disks are submaximal. Moreover, weak but consistent correlations exist such that disks with a fainter central surface brightness in bluer and less luminous galaxies of later morphological types are kinematically colder with respect to their rotational velocities.

We measure the contribution of galaxy disks to the overall gravitational potential of 30 nearly face-on intermediate-to-late-type spirals from the DiskMass Survey. The central vertical velocity dispersion of the disk stars, sigma(z,R=0), is related to the maximum rotation speed (Vmax) as sigma(z,R=0) ~ 0.26 Vmax, consistent with previous measurements for edge-on disk galaxies and a mean stellar velocity ellipsoid axial ratio sigma(z) / sigma(R) = 0.6. For reasonable values of disk oblateness, this relation implies these galaxy disks are submaximal. We find disks in our sample contribute only 15% to 30% of the dynamical mass within 2.2 disk scale-lengths (hR), with percentages increasing systematically with luminosity, rotation speed and redder color. These trends indicate the mass ratio of disk-to-total matter remains at or below 50% at 2.2 hR even for the most extreme, fast-rotating disks (Vmax > 300 km/s), of the reddest rest-frame, face-on color (B-K ~ 4 mag), and highest luminosity (M(K)<-26.5 mag). Therefore, spiral disks in general should be submaximal. Our results imply that the stellar mass-to-light ratio and hence the accounting of baryons in stars should be lowered by at least a factor of 3.

We present a detailed and unique mass budget for the high-surface-brightness galaxy UGC 463, showing it is dominated by dark matter (DM) at radii beyond one scale length (h_R) and has a baryonic-to-DM mass ratio of approximately 1:3 within 4.2 h_R. Assuming a constant scale height (h_z, calculated via an empirical oblateness relation), we calculate dynamical disk mass surface densities from stellar kinematics, which provide vertical velocity dispersions after correcting for the shape of the stellar velocity ellipsoid (measured to have sigma_theta/sigma_R=1.04 +/- 0.22 and sigma_z/sigma_R=0.48 +/- 0.09). We isolate the stellar mass surface density by accounting for all gas mass components and find an average K-band mass-to-light ratio of 0.22 +/- 0.09 (ran) ^{+0.16}_{-0.15} (sys) M_{sun}/L_{sun}^{K}; Zibetti et al. and Bell et al. predict, respectively, 0.56 and 3.6 times our dynamical value based on stellar-population-synthesis modeling. The baryonic matter is submaximal by a factor of ~3 in mass and the baryonic-to-total circular-speed ratio is 0.61^{+0.07}_{-0.09} (ran) ^{+0.12}_{-0.18} (sys) at 2.2 h_R; however, the disk is globally stable with a multi-component stability that decreases asymptotically with radius to Q~2. We directly calculate the circular speed of the DM halo by subtracting the baryonic contribution to the total circular speed; the result is equally well described by either a Navarro-Frenk-White halo or a pseudo-isothermal sphere. The volume density is dominated by DM at heights of |z|>1.6 h_z for radii of R > h_R. As is shown in follow-up papers, UGC 463 is just one example among nearly all galaxies we have observed that contradict the hypothesis that high-surface-brightness spiral galaxies have maximal disks.

Wide-field multi-object spectroscopy is a high priority for European astronomy over the next decade. Most 8-10m telescopes have a small field of view, making 4-m class telescopes a particularly attractive option for wide-field instruments. We present a science case and design drivers for a wide-field multi-object spectrograph (MOS) with integral field units for the 4.2-m William Herschel Telescope (WHT) on La Palma. The instrument intends to take advantage of a future prime-focus corrector and atmospheric-dispersion corrector that will deliver a field of view 2 deg in diameter, with good throughput from 370 to 1,000 nm. The science programs cluster into three groups needing three different resolving powers R: (1) high-precision radial-velocities for Gaia-related Milky Way dynamics, cosmological redshift surveys, and galaxy evolution studies (R = 5,000), (2) galaxy disk velocity dispersions (R = 10,000) and (3) high-precision stellar element abundances for Milky Way archaeology (R = 20,000). The multiplex requirements of the different science cases range from a few hundred to a few thousand, and a range of fibre-positioner technologies are considered. Several options for the spectrograph are discussed, building in part on published design studies for E-ELT spectrographs. Indeed, a WHT MOS will not only efficiently deliver data for exploitation of important imaging surveys planned for the coming decade, but will also serve as a test-bed to optimize the design of MOS instruments for the future E-ELT.

We present a performance analysis of the DiskMass Survey. The survey uses collisionless tracers in the form of disk stars to measure the surface-density of spiral disks, to provide an absolute calibration of the stellar mass-to-light ratio, and to yield robust estimates of the dark-matter halo density profile in the inner regions of galaxies. We find a disk inclination range of 25-35 degrees is optimal for our measurements, consistent with our survey design to select nearly face-on galaxies. Uncertainties in disk scale-heights are significant, but can be estimated from radial scale-lengths to 25% now, and more precisely in the future. We detail the spectroscopic analysis used to derive line-of-sight velocity dispersions, precise at low surface-brightness, and accurate in the presence of composite stellar populations. Our methods take full advantage of large-grasp integral-field spectroscopy and an extensive library of observed stars. We show that the baryon-to-total mass fraction (F_b) is not a well-defined observational quantity because it is coupled to the halo mass model. This remains true even when the disk mass is known and spatially-extended rotation curves are available. In contrast, the fraction of the rotation speed supplied by the disk at 2.2 scale lengths (disk maximality) is a robust observational indicator of the baryonic disk contribution to the potential. We construct the error-budget for the key quantities: dynamical disk mass surface-density, disk stellar mass-to-light ratio, and disk maximality (V_disk / V_circular). Random and systematic errors in these quantities for individual galaxies will be ~25%, while survey precision for sample quartiles are reduced to 10%, largely devoid of systematic errors outside of distance uncertainties.

We present a survey of the mass surface-density of spiral disks, motivated by outstanding uncertainties in rotation-curve decompositions. Our method exploits integral-field spectroscopy to measure stellar and gas kinematics in nearly face-on galaxies sampled at 515, 660, and 860 nm, using the custom-built SparsePak and PPak instruments. A two-tiered sample, selected from the UGC, includes 146 nearly face-on galaxies, with B<14.7 and disk scale-lengths between 10 and 20 arcsec, for which we have obtained H-alpha velocity-fields; and a representative 46-galaxy subset for which we have obtained stellar velocities and velocity dispersions. Based on re-calibration of extant photometric and spectroscopic data, we show these galaxies span factors of 100 in L(K) (0.03 < L/L(K)* < 3), 8 in L(B)/L(K), 10 in R-band disk central surface-brightness, with distances between 15 and 200 Mpc. The survey is augmented by 4-70 micron Spitzer IRAC and MIPS photometry, ground-based UBVRIJHK photometry, and HI aperture-synthesis imaging. We outline the spectroscopic analysis protocol for deriving precise and accurate line-of-sight stellar velocity dispersions. Our key measurement is the dynamical disk-mass surface-density. Star-formation rates and kinematic and photometric regularity of galaxy disks are also central products of the study. The survey is designed to yield random and systematic errors small enough (i) to confirm or disprove the maximum-disk hypothesis for intermediate-type disk galaxies, (ii) to provide an absolute calibration of the stellar mass-to-light ratio well below uncertainties in present-day stellar-population synthesis models, and (iii) to make significant progress in defining the shape of dark halos in the inner regions of disk galaxies.

We analyze integral-field ionized gas and stellar line-of-sight kinematics in the context of determining the stellar velocity ellipsoid for spiral galaxies observed by the Disk-Mass Survey. Our new methodology enables us to measure, for the first time, a radial gradient in the ellipsoid ratio sigma_z / sigma_R. Random errors in this decomposition are 15% at two disk scale-lengths.

We report the discovery of optical ghosts generated when using Volume Phase Holographic (VPH) gratings in spectrographs employing the Littrow configuration. The ghost is caused by light reflected off the detector surface, recollimated by the camera, recombined by, and reflected from, the grating and reimaged by the camera onto the detector. This recombination can occur in two different ways. We observe this ghost in two spectrographs being developed by the University of Wisconsin - Madison: the Robert Stobie Spectrograph for the Southern African Large Telescope and the Bench Spectrograph for the WIYN 3.5m telescope. The typical ratio of the brightness of the ghost relative to the integrated flux of the spectrum is of order 10^-4, implying a recombination efficiency of the VPH gratings of order 10^-3 or higher, consistent with the output of rigorous coupled wave analysis. Any spectrograph employing VPH gratings, including grisms, in Littrow configuration will suffer from this ghost, though the general effect is not intrinsic to VPH gratings themselves and has been observed in systems with conventional gratings in non-Littrow configurations. We explain the geometric configurations that can result in the ghost as well as a more general prescription for predicting its position and brightness on the detector. We make recommendations for mitigating the ghost effects for spectrographs and gratings currently built. We further suggest design modifications for future VPH gratings to eliminate the problem entirely, including tilted fringes and/or prismatic substrates. We discuss the resultant implications on the spectrograph performance metrics.

NGC 3310 is a local galaxy with an intense, ongoing starburst thought to result from a merger with a companion galaxy. It has several known tidal features in the northwest and southern regions around the main galactic disc, as well as a closed, tidal loop emerging from the eastern side of the disc and rejoining in the north. This loop appears to be distinct from the rest of the shells surrounding NGC3310 and is the first of its kind to be detected in a starburst galaxy. In this work, we present UBVR photometry to faint surface brightness levels of this debris network, and we explore various strategies for modelling NGC 3310's disc and subtracting its contribution from each region of debris. We then compare these photometric results with the GALEV spectral synthesis models, and find possible material from the intruder galaxy, suggesting that the recent accretion of several small galaxies is driving the evolution of NGC 3310.

Little is known about the content and distribution of dark matter in spiral galaxies. To break the degeneracy in galaxy rotation curve decompositions, which allows a wide range of dark matter halo density profiles, an independent measure of the mass surface density of stellar disks is needed. Here, we present our ongoing Disk Mass project, using two custom-built Integral Field Units, to measure the vertical velocity dispersion of stars in ~40 spiral galaxies. This will provide a kinematic measurement of the stellar disk mass required to break the degeneracy, enabling us to determine the dark matter properties in spiral galaxies with unprecedented accuracy. Here we present preliminary results for three galaxies with different central disk surface brightness levels.

We present the decomposition of the stellar velocity ellipsoid using stellar velocity dispersions within a 40 deg wedge about the major-axis (sigma_maj), the epicycle approximation, and the asymmetric drift equation. Thus, we employ no fitted forms for sigma_maj and escape interpolation errors resulting from comparisons of the major and minor axes. We apply the theoretical construction of the method to integral field data taken for NGC 3949 and NGC 3982. We derive the vertical-to-radial velocity dispersion ratio (sigma_z / sigma_R) and find (1) our decomposition method is accurate and reasonable, (2) NGC 3982 appears to be rather typical of an Sb type galaxy with sigma_z / sigma_R = 0.73 (+0.13/-0.11) despite its high surface brightness and small size, and (3) NGC 3949 has a hot disk with sigma_z / sigma_R = 1.18 (+0.36/-0.28).

(Abridged) We explore the spatial distribution of stars in the Sculptor dwarf spheroidal (dSph) galaxy over an area of 7.82 deg^2. We identify red giant branch (RGB) starts via Washington M, T_2+DDO51 photometry and a blue horizontal branch (BHB) population to map the spatial structure of the dSph. A spectroscopically observed subset of Sculptor candidate stars yield a systemic heliocentric velocity for the system of v_{hel}=110.43 km/s, in good agreement with previous studies, and a global velocity dispersion of sigma_v=8.8 km/s, which may rise slightly past 0.4r_{lim}. To a limit of M~19, we find 94% of the photometrically-selected Sculptor giant star candidates with spectrocopic measurements are kinematically associated with Sculptor, and four of ten stars not selected photometrically are selected kinematically; our candidate samples are likely to be very pure. We take considerable care in assessing the contaminating background level in our photometric sample to ensure accurate density profiles. These assessments verify that we detect a considerable stellar density of Sculptor stars to the limits of our survey area in both the RGB and BHB samples. We find the Sculptor density profile is well-fit by a King profile of limiting radius r_{lim} = 79.6 within ~60 arcmin, beyond which a "break" to a power law profile occurs. This break population must be either a bound group of "halo stars" around the Sculptor dSph or unbound tidal debris. The latter is supported by 2D distribution analyses and, if true, implies a fractional mass-loss rate of ~0.042 Gyr^{-1} for Sculptor. Finally, likely more metal-poor RGB stars (as selected by color and magnitude) are significantly less centrally concentrated and, therefore, constitute the primary contributor to the likely tidally-stripped parts of the dSph.

We present a performance analysis of SparsePak and the WIYN Bench Spectrograph for precision studies of stellar and ionized gas kinematics of external galaxies. We focus on spectrograph configurations with echelle and low-order gratings yielding spectral resolutions of ~10000 between 500-900nm. These configurations are of general relevance to the spectrograph performance. Benchmarks include spectral resolution, sampling, vignetting, scattered light, and an estimate of the system absolute throughput. Comparisons are made to other, existing, fiber feeds on the WIYN Bench Spectrograph. Vignetting and relative throughput are found to agree with a geometric model of the optical system. An aperture-correction protocol for spectrophotometric standard-star calibrations has been established using independent WIYN imaging data and the unique capabilities of the SparsePak fiber array. The WIYN point-spread-function is well-fit by a Moffat profile with a constant power-law outer slope of index -4.4. We use SparsePak commissioning data to debunk a long-standing myth concerning sky-subtraction with fibers: By properly treating the multi-fiber data as a ``long-slit'' it is possible to achieve precision sky subtraction with a signal-to-noise performance as good or better than conventional long-slit spectroscopy. No beam-switching is required, and hence the method is efficient. Finally, we give several examples of science measurements which SparsePak now makes routine. These include H$\alpha$ velocity fields of low surface-brightness disks, gas and stellar velocity-fields of nearly face-on disks, and stellar absorption-line profiles of galaxy disks at spectral resolutions of ~24,000.

The combination of dynamical and photometric properties of galaxies offers a largely un-tapped source of information on how galaxies assembled and where stars formed. Bi-dimensional kinematic measurements have been the stumbling block. The light-gathering power of SALT coupled with the high-throughput performance of the Prime Focus Imaging Spectrograph (PFIS) yield a superb facility for measuring velocity-ellipsoids of stars and gas in galaxies out to gigaparsec distances. From these data dynamical asymmetries arising from lopsided or elliptical halos may be probed; disk-mass and mass-decompositions may be uniquely determined; mechanisms for disk heating constrained; and a zeropoint for the mass-to-light ratios of stellar populations set. A number of groups within the SALT consortium are interested in making these measurements using a variety of different, but complementary approaches. The scientific potential from their synthesis is very promising. We describe some unusual observational modes in which PFIS may be used to probe the shape of dark-matter halos and the content of galaxy disks.

We present our Disk Mass project as the main science case for building a new fiber IFU-module for the PMAS spectrograph, currently mounted at the Cassegrain focus of the 3.5m telescope on Calar Alto. Compared to traditional long-slit observations, the large light collecting power of 2-dimensional Integral Field Units dramatically improves the prospects for performing spectroscopy on extended low surface brightness objects with high spectral resolution. This enables us to measure stellar velocity dispersions in the outer disk of normal spiral galaxies. We describe some results from a PMAS pilot study using the existing lenslet array, and provide a basic description of the new fiber IFU-module for PMAS.

We explore the metallicity distribution function (MDF) of red giant stars in Omega Centauri from a catalogue of Washington M, T_2 and DDO51 photometry covering over 1.1 deg^2 outside the cluster core. Using updated calibrations of giant branch isometallicity loci in this filter system, photometric metallicities, guided by previously published spectroscopic abundances, are derived. Several methods are employed to correct the MDF for contamination by Galactic stars, including: (1) use of the surface gravity sensitivity of the (M-DDO51) color index to eliminate foreground dwarf stars, (2) radial velocities, and (3) membership probabilities from proper motions. The contamination-corrected MDF for Omega Cen shows a range of enrichment levels spanning nearly 2 dex in [Fe/H], and with peaks at [Fe/H]=-1.6, -1.2, and -0.9.

We are exploring the extended stellar distributions of Galactic satellite galaxies and globular clusters. For seven objects studied thus far, the observed profile departs from a King function at large r, revealing a ``break population'' of stars. In our sample, the relative density of the ``break'' correlates to the inferred M/L of these objects. We discuss opposing hypotheses for this trend: (1) Higher M/L objects harbor more extended dark matter halos that support secondary, bound, stellar ``halos''. (2) The extended populations around dwarf spheroidals (and some clusters) consist of unbound, extratidal debris from their parent objects, which are undergoing various degrees of tidal disruption. In this scenario, higher M/L ratios reflect higher degrees of virial non-equilibrium in the parent objects, thus invalidating a precept underlying the use of core radial velocities to obtain masses.