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

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

Integral field spectroscopy can map astronomical objects spatially and spectroscopically. Due to instrumental and atmospheric effects, it is common for integral field instruments to yield a sampling of the sky image that is both irregular and wavelength-dependent. Most subsequent analysis procedures require a regular, wavelength independent sampling (for example a fixed rectangular grid), and thus an initial step of fundamental importance is to resample the data onto a new grid. The best possible resampling would produce a well-sampled image, with a resolution equal to that imposed by the intrinsic spatial resolution of the instrument, telescope, and atmosphere, and with no statistical correlations between neighboring pixels. A standard method in the field to produce a regular set of samples from an irregular set of samples is Shepard's method, but Shepard's method typically yields images with a degraded resolution and large statistical correlations between pixels. Here we introduce a new method, which improves on Shepard's method in both these respects. We apply this method to data from the Mapping Nearby Galaxies at Apache Point Observatory survey, part of Sloan Digital Sky Survey IV, demonstrating a full-width half maximum close to that of the intrinsic resolution (and ~ 16% better than Shepard's method) and low statistical correlations between pixels. These results nearly achieve the ideal resampling. This method can have broader applications to other integral field data sets and to other astronomical data sets (such as dithered images) with irregular sampling.

We describe the design and performance of the near-infrared (1.51--1.70 micron), fiber-fed, multi-object (300 fibers), high resolution (R = lambda/delta lambda ~ 22,500) spectrograph built for the Apache Point Observatory Galactic Evolution Experiment (APOGEE). APOGEE is a survey of ~ 10^5 red giant stars that systematically sampled all Milky Way populations (bulge, disk, and halo) to study the Galaxy's chemical and kinematical history. It was part of the Sloan Digital Sky Survey III (SDSS-III) from 2011 -- 2014 using the 2.5 m Sloan Foundation Telescope at Apache Point Observatory, New Mexico. The APOGEE-2 survey is now using the spectrograph as part of SDSS-IV, as well as a second spectrograph, a close copy of the first, operating at the 2.5 m du Pont Telescope at Las Campanas Observatory in Chile. Although several fiber-fed, multi-object, high resolution spectrographs have been built for visual wavelength spectroscopy, the APOGEE spectrograph is one of the first such instruments built for observations in the near-infrared. The instrument's successful development was enabled by several key innovations, including a "gang connector" to allow simultaneous connections of 300 fibers; hermetically sealed feedthroughs to allow fibers to pass through the cryostat wall continuously; the first cryogenically deployed mosaic volume phase holographic grating; and a large refractive camera that includes mono-crystalline silicon and fused silica elements with diameters as large as ~ 400 mm. This paper contains a comprehensive description of all aspects of the instrument including the fiber system, optics and opto-mechanics, detector arrays, mechanics and cryogenics, instrument control, calibration system, optical performance and stability, lessons learned, and design changes for the second instrument.

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.

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.

We present the Data Release 14 Quasar catalog (DR14Q) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). This catalog includes all SDSS-IV/eBOSS objects that were spectroscopically targeted as quasar candidates and that are confirmed as quasars via a new automated procedure combined with a partial visual inspection of spectra, have luminosities $M_{\rm i} \left[ z=2 \right] < -20.5$ (in a $\Lambda$CDM cosmology with $H_0 = 70 \ {\rm km \ s^{-1} \ Mpc ^{-1}}$, $\Omega_{\rm M} = 0.3$, and $\Omega_{\rm \Lambda} = 0.7$), and either display at least one emission line with a full width at half maximum (FWHM) larger than $500 \ {\rm km \ s^{-1}}$ or, if not, have interesting/complex absorption features. The catalog also includes previously spectroscopically-confirmed quasars from SDSS-I, II and III. The catalog contains 526,356 quasars 144,046 are new discoveries since the beginning of SDSS-IV) detected over 9,376 deg$^2$ (2,044 deg$^2$ having new spectroscopic data available) with robust identification and redshift measured by a combination of principal component eigenspectra. The catalog is estimated to have about 0.5% contamination. The catalog identifies 21,877 broad absorption line quasars and lists their characteristics. For each object, the catalog presents SDSS five-band CCD-based photometry with typical accuracy of 0.03 mag. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

As astronomers increasingly exploit the information available in the time domain, spectroscopic variability in particular opens broad new channels of investigation. Here we describe the selection algorithms for all targets intended for repeat spectroscopy in the Time Domain Spectroscopic Survey (TDSS), part of the extended Baryon Oscillation Spectroscopic Survey within the Sloan Digital Sky Survey-IV. Also discussed are the scientific rationale and technical constraints leading to these target selections. The TDSS includes a large "Repeat Quasar Spectroscopy" (RQS) program delivering ~13,000 repeat spectra of confirmed SDSS quasars, and several smaller "Few-Epoch Spectroscopy" (FES) programs targeting specific classes of quasars as well as stars. The RQS program aims to provide a large and diverse quasar data set for studying variations in quasar spectra on timescales of years, a comparison sample for the FES quasar programs, and opportunity for discovering rare, serendipitous events. The FES programs cover a wide variety of phenomena in both quasars and stars. Quasar FES programs target broad absorption line quasars, high signal-to-noise ratio normal broad line quasars, quasars with double-peaked or very asymmetric broad emission line profiles, binary supermassive black hole candidates, and the most photometrically variable quasars. Strongly variable stars are also targeted for repeat spectroscopy, encompassing many types of eclipsing binary systems, and classical pulsators like RR Lyrae. Other stellar FES programs allow spectroscopic variability studies of active ultracool dwarf stars, dwarf carbon stars, and white dwarf/M dwarf spectroscopic binaries. We present example TDSS spectra and describe anticipated sample sizes and results.

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

We investigate the possibility of detecting in redshift surveys a hemispherical power asymmetry similar to that first reported in CMB observations. We assume the hemispherical asymmetry arises from a linear gradient in comoving coordinates in the perturbation amplitude. We predict the resulting clustering of galaxy or galaxy cluster tracers using an excursion set approach; doing so accounts for the variation of both the underlying clustering and the tracer bias. Based on the predicted variation of the clustering of tracers, we perform a Fisher matrix forecast of the galaxy clustering amplitude and calculate the statistical significance for ideal surveys and planned surveys. The results indicate that the DESI galaxy survey would be able to detect this signal with higher than $3\sigma$ significance if the asymmetry does exist. We also investigate the amplitude and scale dependence of the above result. The DESI galaxy survey can probe the dipole amplitude higher than 0.04, which correspond to a $\pm4\%$ difference of the temperature fluctuation along and opposite the dipole direction, at least at the $2\sigma$ level. Additionally, we investigate a modulation of the power spectrum that exhibits asymmetry only for large scales. This modulation is potentially detectable. For Milky Way galaxy mass tracers, the scale-dependent modulation yields a larger change in the large scale power spectrum than does a scale-independent modulation, because the former does not alter the bias.

We compare a large suite of theoretical cosmological models to observational data from the cosmic microwave background, baryon acoustic oscillation measurements of expansion, Type Ia SNe measurements of expansion, redshift space distortion measurements of the growth of structure, and the local Hubble constant. Our theoretical models include parametrizations of dark energy as well as physical models of dark energy and modified gravity. We determine the constraints on the model parameters, incorporating the redshift space distortion data directly in the analysis. To determine whether models can be ruled out, we evaluate the $p$ value (the probability under the model of obtaining data as bad or worse than the observed data). In our comparison, we find the well known tension of H$_0$ with the other data; no model resolves this tension successfully. Among the models we consider, the large scale growth of structure data does not affect the modified gravity models as a category particularly differently than dark energy models; it matters for some modified gravity models but not others, and the same is true for dark energy models. We compute predicted observables for each model under current observational constraints, and identify models for which future observational constraints will be particularly informative.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

We present the first scientific results from the luminous red galaxy sample (LRG) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We measure the small and intermediate scale clustering from a sample of more than 61,000 galaxies in the redshift range $0.6 < z < 0.9$. We interpret these measurements in the framework of the Halo Occupation Distribution. The bias of eBOSS LRGs is $2.30 \pm 0.03$, with a satellite fraction of $13\pm3$\% and a mean halo mass of $2.5\times10^{13}h^{-1}M_{\odot}$. These results are consistent with expectations, demonstrating that eBOSS galaxies will be reliable tracers of large scale structure at $z\sim 0.7$. The eBOSS galaxy bias implies a scatter of luminosity at fixed halo mass, $\sigma_{\log L}$, of 0.19 dex. Using the clustering of massive galaxies from BOSS-CMASS, BOSS-LOWZ, and SDSS, we find that $\sigma_{\log L}=0.19$ is consistent with observations over the full redshift range that these samples cover. The addition of eBOSS to previous surveys allows investigation of the evolution of massive galaxies over the past $\sim 7$ Gyr.

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 constrain the spatial gradient of star formation history within galaxies using the colour gradients in NUV-u and u-i for a local spatially-resolved galaxy sample. By splitting each galaxy into an inner and an outer part, we find that most galaxies show negative gradients in these two colours. We first rule out dust extinction gradient and metallicity gradient as the dominant source for the colour gradient. Then using stellar population models, we explore variations in star formation history to explain the colour gradients. As shown by our earlier work, a two-phase SFH consisting of an early secular evolution (growth) phase and a subsequent rapid evolution (quenching) phase is necessary to explain the observed colour distributions among galaxies. We explore two different inside-out growth models and two different inside-out quenching models by varying parameters of the SFH between inner and outer regions of galaxies. Two of the models can explain the observed range of colour gradients in NUV-u and u-i colours. We further distinguish them using an additional constraint provided by the u-i colour gradient distribution, under the assumption of constant galaxy formation rate and a common SFH followed by most galaxies. We find the best model is an inside-out growth model in which the inner region has a shorter e-folding time scale in the growth phase than the outer region. More spatially resolved ultraviolet (UV) observations are needed to improve the significance of the result.

We describe the sample design for the SDSS-IV MaNGA survey and present the final properties of the main samples along with important considerations for using these samples for science. Our target selection criteria were developed while simultaneously optimizing the size distribution of the MaNGA integral field units (IFUs), the IFU allocation strategy, and the target density to produce a survey defined in terms of maximizing S/N, spatial resolution, and sample size. Our selection strategy makes use of redshift limits that only depend on i-band absolute magnitude ($M_i$), or, for a small subset of our sample, $M_i$ and color (NUV-i). Such a strategy ensures that all galaxies span the same range in angular size irrespective of luminosity and are therefore covered evenly by the adopted range of IFU sizes. We define three samples: the Primary and Secondary samples are selected to have a flat number density with respect to $M_i$ and are targeted to have spectroscopic coverage to 1.5 and 2.5 effective radii (Re), respectively. The Color-Enhanced supplement increases the number of galaxies in the low-density regions of color-magnitude space by extending the redshift limits of the Primary sample in the appropriate color bins. The samples cover the stellar mass range $5\times10^8 \leq M_* \leq 3\times10^{11} M_{\odot}$ and are sampled at median physical resolutions of 1.37 kpc and 2.5 kpc for the Primary and Secondary samples respectively. We provide weights that will statistically correct for our luminosity and color-dependent selection function and IFU allocation strategy, thus correcting the observed sample to a volume limited sample.

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 Sloan Digitial Sky Surveys (SDSS) have been collecting imaging and spectoscopic data since 1998. These data as well as their derived data products are made publicly available through regular data releases, of which the 13th took place summer 2016. Although public data releases can be challenging to manage, they signficantly increase the impact of a survey, both scientifically and educationally.

Fiber-fed multi-object spectroscopic surveys, with their ability to collect an unprecedented number of redshifts, currently dominate large-scale structure studies. However, physical constraints limit these surveys from successfully collecting redshifts from galaxies too close to each other on the focal plane. This ultimately leads to significant systematic effects on galaxy clustering measurements. Using simulated mock catalogs, we demonstrate that fiber collisions have a significant impact on the power spectrum, $P(k)$, monopole and quadrupole that exceeds sample variance at scales smaller than $k\sim0.1~h/Mpc$. We present two methods to account for fiber collisions in the power spectrum. The first, statistically reconstructs the clustering of fiber collided galaxy pairs by modeling the distribution of the line-of-sight displacements between them. It also properly accounts for fiber collisions in the shot-noise correction term of the $P(k)$ estimator. Using this method, we recover the true $P(k)$ monopole of the mock catalogs with residuals of $<0.5\%$ at $k=0.3~h/Mpc$ and $<4\%$ at $k=0.83~h/Mpc$ -- a significant improvement over existing correction methods. The quadrupole, however, does not improve significantly. The second method models the effect of fiber collisions on the power spectrum as a convolution with a configuration space top-hat function that depends on the physical scale of fiber collisions. It directly computes theoretical predictions of the fiber-collided $P(k)$ multipoles and reduces the influence of smaller scales to a set of nuisance parameters. Using this method, we reliably model the effect of fiber collisions on the monopole and quadrupole down to the scale limits of theoretical predictions. The methods we present in this paper will allow us to robustly analyze galaxy power spectrum multipole measurements to much smaller scales than previously possible.

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.

The uncertain origin of the recently-discovered `changing-looking' quasar phenomenon -- in which a luminous quasar dims significantly to a quiescent state in repeat spectroscopy over ~10 year timescales -- may present unexpected challenges to our understanding of quasar accretion. To better understand this phenomenon, we take a first step to building a sample of changing-look quasars with a systematic but simple archival search for these objects in the Sloan Digital Sky Survey Data Release 12. By leveraging the >10 year baselines for objects with repeat spectroscopy, we uncover two new changing-look quasars, and a third discovered previously. Decomposition of the multi-epoch spectra and analysis of the broad emission lines suggest that the quasar accretion disk emission dims due to rapidly decreasing accretion rates (by factors of >2.5), while disfavoring changes in intrinsic dust extinction for the two objects where these analyses are possible. Broad emission line energetics also support intrinsic dimming of quasar emission as the origin for this phenomenon rather than transient tidal disruption events or supernovae. Although our search criteria included quasars at all redshifts and transitions from either quasar-like to galaxy-like states or the reverse, all of the clear cases of changing-look quasars discovered were at relatively low-redshift (z ~ 0.2 - 0.3) and only exhibit quasar-like to galaxy-like transitions.

We measure the clustering of X-ray, radio, and mid-IR-selected active galactic nuclei (AGN) at 0.2 < z < 1.2 using multi-wavelength imaging and spectroscopic redshifts from the PRIMUS and DEEP2 redshift surveys, covering 7 separate fields spanning ~10 square degrees. Using the cross-correlation of AGN with dense galaxy samples, we measure the clustering scale length and slope, as well as the bias, of AGN selected at different wavelengths. Similar to previous studies, we find that X-ray and radio AGN are more clustered than mid-IR-selected AGN. We further compare the clustering of each AGN sample with matched galaxy samples designed to have the same stellar mass, star formation rate, and redshift distributions as the AGN host galaxies and find no significant differences between their clustering properties. The observed differences in the clustering of AGN selected at different wavelengths can therefore be explained by the clustering differences of their host populations, which have different distributions in both stellar mass and star formation rate. Selection biases inherent in AGN selection, therefore, determine the clustering of observed AGN samples. We further find no significant difference between the clustering of obscured and unobscured AGN, using IRAC or WISE colors or X-ray hardness ratio.

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

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), one of three core programs in the Sloan Digital Sky Survey-IV (SDSS-IV), is an integral-field spectroscopic (IFS) survey of roughly 10,000 nearby galaxies. It employs dithered observations using 17 hexagonal bundles of 2 arcsec fibers to obtain resolved spectroscopy over a wide wavelength range of 3,600-10,300A. To map the internal variations within each galaxy, we need to perform accurate {\it spectral surface photometry}, which is to calibrate the specific intensity at every spatial location sampled by each individual aperture element of the integral field unit. The calibration must correct only for the flux loss due to atmospheric throughput and the instrument response, but not for losses due to the finite geometry of the fiber aperture. This requires the use of standard star measurements to strictly separate these two flux loss factors (throughput versus geometry), a difficult challenge with standard single-fiber spectroscopy techniques due to various practical limitations. Therefore, we developed a technique for spectral surface photometry using multiple small fiber-bundles targeting standard stars simultaneously with galaxy observations. We discuss the principles of our approach and how they compare to previous efforts, and we demonstrate the precision and accuracy achieved. MaNGA's relative calibration between the wavelengths of H$\alpha$ and H$\beta$ has a root-mean-square (RMS) of 1.7%, while that between [NII] $\lambda$6583A and [OII] $\lambda$3727A has an RMS of 4.7%. Using extinction-corrected star formation rates and gas-phase metallicities as an illustration, this level of precision guarantees that flux calibration errors will be sub-dominant when estimating these quantities. The absolute calibration is better than 5% for more than 89% of MaNGA's wavelength range.

We report small-scale clustering measurements from the PRIMUS spectroscopic redshift survey as a function of color and luminosity. We measure the real-space cross-correlations between 62,106 primary galaxies with PRIMUS redshifts and a tracer population of 545,000 photometric galaxies over redshifts from z=0.2 to z=1. We separately fit a power-law model in redshift and luminosity to each of three independent color-selected samples of galaxies. We report clustering amplitudes at fiducial values of z=0.5 and L=1.5 L*. The clustering of the red galaxies is ~3 times as strong as that of the blue galaxies and ~1.5 as strong as that of the green galaxies. We also find that the luminosity dependence of the clustering is strongly dependent on physical scale, with greater luminosity dependence being found between r=0.0625 Mpc/h and r=0.25 Mpc/h, compared to the r=0.5 Mpc/h to r=2 Mpc/h range. Moreover, over a range of two orders of magnitude in luminosity, a single power-law fit to the luminosity dependence is not sufficient to explain the increase in clustering at both the bright and faint ends at the smaller scales. We argue that luminosity-dependent clustering at small scales is a necessary component of galaxy-halo occupation models for blue, star-forming galaxies as well as for red, quenched galaxies.

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

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

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

We utilize $\Lambda$CDM halo occupation models of galaxy clustering to investigate the evolving stellar mass dependent clustering of galaxies in the PRIsm MUlti-object Survey (PRIMUS) and DEEP2 Redshift Survey over the past eight billion years of cosmic time, between $0.2<z<1.2$. These clustering measurements provide new constraints on the connections between dark matter halo properties and galaxy properties in the context of the evolving large-scale structure of the universe. Using both an analytic model and a set of mock galaxy catalogs, we find a strong correlation between central galaxy stellar mass and dark matter halo mass over the range $M_\mathrm{halo}\sim10^{11}$-$10^{13}~h^{-1}M_\odot$, approximately consistent with previous observations and theoretical predictions. However, the stellar-to-halo mass relation (SHMR) and the mass scale where star formation efficiency reaches a maximum appear to evolve more strongly than predicted by other models, including models based primarily on abundance-matching constraints. We find that the fraction of satellite galaxies in haloes of a given mass decreases significantly from $z\sim0.5$ to $z\sim0.9$, partly due to the fact that haloes at fixed mass are rarer at higher redshift and have lower abundances. We also find that the $M_1/M_\mathrm{min}$ ratio, a model parameter that quantifies the critical mass above which haloes host at least one satellite, decreases from $\approx20$ at $z\sim0$ to $\approx13$ at $z\sim0.9$. Considering the evolution of the subhalo mass function vis-\`{a}-vis satellite abundances, this trend has implications for relations between satellite galaxies and halo substructures and for intracluster mass, which we argue has grown due to stripped and disrupted satellites between $z\sim0.9$ and $z\sim0.5$.

We infer stellar metallicity and abundance ratio gradients for a sample of red galaxies in the Sloan Digital Sky Survey (SDSS) Main galaxy sample. Because this sample does not have multiple spectra at various radii in a single galaxy, we measure these gradients statistically. We separate galaxies into stellar mass bins, stack their spectra in redshift bins, and calculate the measured absorption line indices in projected annuli by differencing spectra in neighboring redshift bins. After determining the line indices, we use stellar population modeling from the EZ\_Ages software to calculate ages, metallicities, and abundance ratios within each annulus. Our data covers the central regions of these galaxies, out to slightly higher than $1 R_{e}$. We find detectable gradients in metallicity and relatively shallow gradients in abundance ratios, similar to results found for direct measurements of individual galaxies. The gradients are only weakly dependent on stellar mass, and this dependence is well-correlated with the change of $R_e$ with mass. Based on this data, we report mean equivalent widths, metallicities, and abundance ratios as a function of mass and velocity dispersion for SDSS early-type galaxies, for fixed apertures of 2.5 kpc and of 0.5 $R_e$.

We present a model for spectrophotometric calibration errors in observations of quasars from the third generation of the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) and describe the correction procedure we have developed and applied to this sample. Calibration errors are primarily due to atmospheric differential refraction and guiding offsets during each exposure. The corrections potentially reduce the systematics for any studies of BOSS quasars, including the measurement of baryon acoustic oscillations using the Lyman-$\alpha$ forest. Our model suggests that, on average, the observed quasar flux in BOSS is overestimated by $\sim 19\%$ at 3600 \AA\ and underestimated by $\sim 24\%$ at 10,000 \AA. Our corrections for the entire BOSS quasar sample are publicly available.

We study the baryon content of low mass galaxies selected from the Sloan Digital Sky Survey (SDSS DR8), focusing on galaxies in isolated environments where the complicating physics of galaxy-galaxy interactions are minimized. We measure neutral hydrogen (HI) gas masses and line-widths for 148 isolated galaxies with stellar mass between $10^7$ and $10^{9.5} M_{\odot}$. We compare isolated low mass galaxies to more massive galaxies and galaxies in denser environments by remeasuring HI emission lines from the Arecibo Legacy Fast ALFA (ALFALFA) survey 40% data release. All isolated low mass galaxies either have large atomic gas fractions or large atomic gas fractions cannot be ruled out via their upper limits. We measure a median atomic gas fraction of $f_{\rm gas} = 0.82 \pm 0.13$ for our isolated low mass sample with no systems below 0.30. At all stellar masses, the correlations between galaxy radius, baryonic mass and velocity width are not significantly affected by environment. Finally, we estimate a median baryon to total dynamical mass fraction of $f_{\rm baryon, disk} = 0.15 \pm 0.18$. We also estimate two different median baryon to halo mass fractions using the results of semi-analytic models $(f_{\rm baryon, halo} = 0.04 \pm 0.06)$ and abundance matching $(f_{\rm baryon, halo} = 0.04 \pm 0.02)$. Baryon fractions estimated directly using HI observations appear independent of environment and maximum circular velocity, while baryon fractions estimated using abundance matching show a significant depletion of baryons at low maximum circular velocities.

We study the evidence for a connection between active galactic nuclei (AGN) fueling and star formation by investigating the relationship between the X-ray luminosities of AGN and the star formation rates (SFRs) of their host galaxies. We identify a sample of 309 AGN with $10^{41}<L_\mathrm{X}<10^{44} $ erg s$^{-1}$ at $0.2 < z < 1.2$ in the PRIMUS redshift survey. We find AGN in galaxies with a wide range of SFR at a given $L_X$. We do not find a significant correlation between SFR and the observed instantaneous $L_X$ for star forming AGN host galaxies. However, there is a weak but significant correlation between the mean $L_\mathrm{X}$ and SFR of detected AGN in star forming galaxies, which likely reflects that $L_\mathrm{X}$ varies on shorter timescales than SFR. We find no correlation between stellar mass and $L_\mathrm{X}$ within the AGN population. Within both populations of star forming and quiescent galaxies, we find a similar power-law distribution in the probability of hosting an AGN as a function of specific accretion rate. Furthermore, at a given stellar mass, we find a star forming galaxy $\sim2-3$ more likely than a quiescent galaxy to host an AGN of a given specific accretion rate. The probability of a galaxy hosting an AGN is constant across the main sequence of star formation. These results indicate that there is an underlying connection between star formation and the presence of AGN, but AGN are often hosted by quiescent galaxies.

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.

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is an integral-field spectroscopic survey of 10,000 nearby galaxies that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). MaNGA's 17 pluggable optical fiber-bundle integral field units (IFUs) are deployed across a 3 deg field, they yield spectral coverage 3600-10,300 Ang at a typical resolution R ~ 2000, and sample the sky with 2" diameter fiber apertures with a total bundle fill factor of 56%. Observing over such a large field and range of wavelengths is particularly challenging for obtaining uniform and integral spatial coverage and resolution at all wavelengths and across each entire fiber array. Data quality is affected by the IFU construction technique, chromatic and field differential refraction, the adopted dithering strategy, and many other effects. We use numerical simulations to constrain the hardware design and observing strategy for the survey with the aim of ensuring consistent data quality that meets the survey science requirements while permitting maximum observational flexibility. We find that MaNGA science goals are best achieved with IFUs composed of a regular hexagonal grid of optical fibers with rms displacement of 5 microns or less from their nominal packing position, this goal is met by the MaNGA hardware, which achieves 3 microns rms fiber placement. We further show that MaNGA observations are best obtained in sets of three 15-minute exposures dithered along the vertices of a 1.44 arcsec equilateral triangle, these sets form the minimum observational unit, and are repeated as needed to achieve a combined signal-to-noise ratio of 5 per Angstrom per fiber in the r-band continuum at a surface brightness of 23 AB/arcsec^2. (abbrev.)

We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST.

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is a 6-year SDSS-IV survey that will obtain resolved spectroscopy from 3600 $\AA$ to 10300 $\AA$ for a representative sample of over 10,000 nearby galaxies. In this paper, we derive spatially resolved stellar population properties and radial gradients by performing full spectral fitting of observed galaxy spectra from P-MaNGA, a prototype of the MaNGA instrument. These data include spectra for eighteen galaxies, covering a large range of morphological type. We derive age, metallicity, dust and stellar mass maps, and their radial gradients, using high spectral-resolution stellar population models, and assess the impact of varying the stellar library input to the models. We introduce a method to determine dust extinction which is able to give smooth stellar mass maps even in cases of high and spatially non-uniform dust attenuation. With the spectral fitting we produce detailed maps of stellar population properties which allow us to identify galactic features among this diverse sample such as spiral structure, smooth radial profiles with little azimuthal structure in spheroidal galaxies, and spatially distinct galaxy sub-components. In agreement with the literature, we find the gradients for galaxies identified as early-type to be on average flat in age, and negative (- 0.15 dex / R$_e$ ) in metallicity, whereas the gradients for late-type galaxies are on average negative in age (- 0.39 dex / R$_e$ ) and flat in metallicity. We demonstrate how different levels of data quality change the precision with which radial gradients can be measured. We show how this analysis, extended to the large numbers of MaNGA galaxies, will have the potential to shed light on galaxy structure and evolution.

We present an analysis of the data produced by the MaNGA prototype run (P-MaNGA), aiming to test how the radial gradients in recent star formation histories, as indicated by the 4000AA-break (D4000), Hdelta absorption (EW(Hd_A)) and Halpha emission (EW(Ha)) indices, can be useful for understanding disk growth and star formation cessation in local galaxies. We classify 12 galaxies observed on two P-MaNGA plates as either centrally quiescent (CQ) or centrally star-forming (CSF), according to whether D4000 measured in the central spaxel of each datacube exceeds 1.6. For each galaxy we generate both 2D maps and radial profiles of D4000, EW(Hd_A) and EW(Ha). We find that CSF galaxies generally show very weak or no radial variation in these diagnostics. In contrast, CQ galaxies present significant radial gradients, in the sense that D4000 decreases, while both EW(Hd_A) and EW(Ha) increase from the galactic center outward. The outer regions of the galaxies show greater scatter on diagrams relating the three parameters than their central parts. In particular, the clear separation between centrally-measured quiescent and star-forming galaxies in these diagnostic planes is largely filled in by the outer parts of galaxies whose global colors place them in the green valley, supporting the idea that the green valley represents a transition between blue-cloud and red-sequence phases, at least in our small sample. These results are consistent with a picture in which the cessation of star formation propagates from the center of a galaxy outwards as it moves to the red sequence.

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 investigate the effects of galaxy environment on the evolution of the quiescent fraction ($f_\mathrm{Q}$) from z =0.8 to 0.0 using spectroscopic redshifts and multi-wavelength imaging data from the PRIsm MUlti-object Survey (PRIMUS) and the Sloan Digitial Sky Survey (SDSS). Our stellar mass limited galaxy sample consists of ~14,000 PRIMUS galaxies within z = 0.2-0.8 and ~64,000 SDSS galaxies within z = 0.05-0.12. We classify the galaxies as quiescent or star-forming based on an evolving specific star formation cut, and as low or high density environments based on fixed cylindrical aperture environment measurements on a volume-limited environment defining population. For quiescent and star-forming galaxies in low or high density environments, we examine the evolution of their stellar mass function (SMF). Then using the SMFs we compute $f_\mathrm{Q}(M_{*})$ and quantify its evolution within our redshift range. We find that the quiescent fraction is higher at higher masses and in denser environments. The quiescent fraction rises with cosmic time for all masses and environments. At a fiducial mass of $10^{10.5}M_\odot$, from z~0.7 to 0.1, the quiescent fraction rises by 15% at the lowest environments and by 25% at the highest environments we measure. These results suggest that for a minority of galaxies their cessation of star formation is due to external influences on them. However, in the recent Universe a substantial fraction of the galaxies that cease forming stars do so due to internal processes.

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 a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $\Lambda$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,\sigma$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

We present measurements of the luminosity and color-dependence of galaxy clustering at 0.2<z<1.0 in the PRIsm MUlti-object Survey (PRIMUS). We quantify the clustering with the redshift-space and projected two-point correlation functions, xi(rp,pi) and wp(rp), using volume-limited samples constructed from a parent sample of over 130,000 galaxies with robust redshifts in seven independent fields covering 9 sq. deg. of sky. We quantify how the scale-dependent clustering amplitude increases with increasing luminosity and redder color, with relatively small errors over large volumes. We find that red galaxies have stronger small-scale (0.1<rp<1 Mpc/h) clustering and steeper correlation functions compared to blue galaxies, as well as a strong color dependent clustering within the red sequence alone. We interpret our measured clustering trends in terms of galaxy bias and obtain values between b_gal=0.9-2.5, quantifying how galaxies are biased tracers of dark matter depending on their luminosity and color. We also interpret the color dependence with mock catalogs, and find that the clustering of blue galaxies is nearly constant with color, while redder galaxies have stronger clustering in the one-halo term due to a higher satellite galaxy fraction. In addition, we measure the evolution of the clustering strength and bias, and we do not detect statistically significant departures from passive evolution. We argue that the luminosity- and color-environment (or halo mass) relations of galaxies have not significantly evolved since z=1. Finally, using jackknife subsampling methods, we find that sampling fluctuations are important and that the COSMOS field is generally an outlier, due to having more overdense structures than other fields; we find that 'cosmic variance' can be a significant source of uncertainty for high-redshift clustering measurements.

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.

Many classes of active galactic nuclei (AGN) have been observed and recorded since the discovery of Seyfert galaxies. In this paper, we examine the sample of luminous galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS). We find a potentially new observational class of AGN, one with strong and broad MgII $2799$\AA\ line emission, but very weak emission in other normal indicators of AGN activity, such as the broad line H$\alpha$, H$\beta$, and the near-ultraviolet AGN continuum, leading to an extreme ratio of broad H$\alpha$/MgII flux relative to normal quasars. Meanwhile, these objects' narrow-line flux ratios reveal AGN narrow-line regions with levels of activity consistent with the MgII fluxes and in agreement with that of normal quasars. These AGN may represent an extreme case of the Baldwin effect, with very low continuum and high equivalent width relative to typical quasars, but their ratio of broad MgII to broad Balmer emission remains very unusual. They may also be representative of a class of AGN where the central engine is observed indirectly with scattered light. These galaxies represent a small fraction of the total population of luminous galaxies ($\simeq 0.1$\%), but are more likely (about 3.5 times) to have AGN-like nuclear line emission properties than other luminous galaxiess. Because MgII is usually inaccessible for the population of nearby galaxies, there may exist a related population of broad-line MgII emitters in the local universe which are currently classified as narrow-line emitters (Seyfert 2s) or LINERs.

We have identified a low surface brightness stellar stream from visual inspection of SDSS imaging for the edge-on, spiral galaxy NGC5387. A blue overdensity was also identified in SDSS coincident with the stream intersection with the NGC5387 disk. The overdensity was also detected in the GALEX Deep Imaging Survey and found to contribute 38% of the total FUV integrated flux from NGC5387, which suggests that the region is actively forming stars. Deeper imaging was acquired with the Vatican Advanced Technology Telescope (VATT) in the B, V, and R filters that confirmed the presence of both the stellar stream and the blue overdensity. Analysis of the VATT photometry indicates the stellar stream is red in color, B-V = 0.7, and has a stellar mass of 6x10^8 M$_{\odot}$, which implies a merger ratio of 1:50. Assessment of the stream morphology suggests that the merger event has a circular radius, R~16 kpc, the stream formed in ~400 Myr, and the progenitor had a total mass of ~2x10 M$_{\odot}$. Spectroscopy from LBT+MODS1 was used to determine that the blue overdensity is at the same redshift as NGC5387, consists of young stellar populations (~10 Myr), is metal-poor (12 + log(O/H) = 8.03), and forming stars at an enhanced rate (~1-3 M$_{\odot}$/year) given its total stellar mass (2x10^7 M$_{\odot}$). Several interpretations are posited to explain these observational data, of which the most likely are (i) that the blue overdensity is a region of enhanced star formation in the outer disk of NGC5387 induced by the minor accretion event, and (ii) that the blue overdensity is the progenitor of the stellar stream undergoing a period of enhanced star formation as a result of its interaction with NGC5387. Confirmation and theoretical exploration of these scenarios are presented in a companion paper.

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

We present a study of Spitzer/IRAC and X-ray active galactic nuclei (AGNs) selection techniques in order to quantify the overlap, uniqueness, contamination, and completeness of each. We investigate how the overlap and possible contamination of the samples depends on the IR and X-ray depths. We use Spitzer/IRAC imaging, Chandra and XMM X-ray imaging, and PRism MUlti-object Survey (PRIMUS) spectroscopic redshifts to construct galaxy and AGN samples at 0.2<z<1.2 over 8 deg^2. We construct samples over a wide range of IRAC flux limits (SWIRE to GOODS depth) and X-ray flux limits (10 ks to 2 Ms). We compare IR-AGN samples defined using the IRAC color selection of Stern et al. and Donley et al. with X-ray detected AGN samples. For roughly similar depth IR and X-ray surveys, we find that ~75% of IR-AGN are identified as X-ray AGN. This fraction increases to ~90% when comparing against the deepest X-ray data, indicating that only ~10% of IR-selected AGN may be heavily obscured. The IR-AGN selection proposed by Stern et al. suffers from contamination by star-forming galaxies at various redshifts when using deeper IR data, though the selection technique works well for shallow IR data. While similar overall, the IR-AGN samples preferentially contain more luminous AGN, while the X-ray AGN samples preferentially contain lower specific accretion rate AGN, where the host galaxy light dominates at IR wavelengths. The host galaxy populations of the IR and X-ray AGN samples have similar restframe colors and stellar masses; both selections identify AGN in blue, star-forming and red, quiescent galaxies.

We present an observationally motivated model to connect the AGN and galaxy populations at 0.2<z<1.0 and predict the AGN X-ray luminosity function (XLF). We start with measurements of the stellar mass function of galaxies (from the Prism Multi-object Survey) and populate galaxies with AGNs using models for the probability of a galaxy hosting an AGN as a function of specific accretion rate. Our model is based on measurements indicating that the specific accretion rate distribution is a universal function across a wide range of host stellar mass with slope gamma_1 = -0.65 and an overall normalization that evolves with redshift. We test several simple assumptions to extend this model to high specific accretion rates (beyond the measurements) and compare the predictions for the XLF with the observed data. We find good agreement with a model that allows for a break in the specific accretion rate distribution at a point corresponding to the Eddington limit, a steep power-law tail to super-Eddington ratios with slope gamma_2 = -2.1 +0.3 -0.5, and a scatter of 0.38 dex in the scaling between black hole and host stellar mass. Our results show that samples of low luminosity AGNs are dominated by moderately massive galaxies (M* ~ 10^{10-11} M_sun) growing with a wide range of accretion rates due to the shape of the galaxy stellar mass function rather than a preference for AGN activity at a particular stellar mass. Luminous AGNs may be a severely skewed population with elevated black hole masses relative to their host galaxies and in rare phases of rapid accretion.

We present results on the clustering of 282,068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4<z<0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ~0.5 to ~90 Mpc/h. We compare these estimates with the expectations of the flat LCDM cosmological model with parameters compatible with WMAP7 data. We use the MultiDark cosmological simulation together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the LCDM model gives a reasonable description to the observed correlation functions at z~0.5, which is a remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a deviation (>~10%) in the correlation functions for scales less than ~1 Mpc/h and ~10-40 Mpc/h. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12% of the "galaxies" in the abundance-matched sample are satellites inhabiting central haloes with mass M>~1e14 M_sun/h. Using the MultiDark simulation we also study the real space halo bias b(r) of the matched catalogue finding that b=2.00+/-0.07 at large scales, consistent with the one obtained using the measured BOSS projected correlation function. Furthermore, the linear large-scale bias depends on the number density n of the abundance-matched sample as b=-0.048-(0.594+/-0.02)*log(n/(h/Mpc)^3). Extrapolating these results to BAO scales we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ~2-4% dips at ~3 sigma level for wavenumbers k>~0.1 h/Mpc in the linear large-scale bias.

We measure the luminosity and color dependence and the redshift evolution of galaxy clustering in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, which includes about 260,000 galaxies over ~3,300 sq. deg in the redshift range 0.43<z<0.7. To minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole CMASS sample, if approximated by a power-law, has a correlation length of r_0=7.93\pm0.06Mpc/h and an index of \gamma=1.85\pm0.01. Clear dependences on galaxy luminosity and color are found for the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. At a given luminosity (k+e corrected), no significant evolution of the projected 2PCFs with redshift is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints on galaxy formation and evolution.

We measure the evolution of the stellar mass function (SMF) from z=0-1 using multi-wavelength imaging and spectroscopic redshifts from the PRism MUlti-object Survey (PRIMUS) and the Sloan Digital Sky Survey (SDSS). From PRIMUS we construct an i<23 flux-limited sample of ~40,000 galaxies at z=0.2-1.0 over five fields totaling ~5.5 deg^2, and from the SDSS we select ~170,000 galaxies at z=0.01-0.2 that we analyze consistently with respect to PRIMUS to minimize systematic errors in our evolutionary measurements. We find that the SMF of all galaxies evolves relatively little since z=1, although we do find evidence for mass assembly downsizing; we measure a ~30% increase in the number density of ~10^10 Msun galaxies since z~0.6, and a <10% change in the number density of all >10^11 Msun galaxies since z~1. Dividing the sample into star-forming and quiescent using an evolving cut in specific star-formation rate, we find that the number density of ~10^10 Msun star-forming galaxies stays relatively constant since z~0.6, whereas the space-density of >10^11 Msun star-forming galaxies decreases by ~50% between z~1 and z~0. Meanwhile, the number density of ~10^10 Msun quiescent galaxies increases steeply towards low redshift, by a factor of ~2-3 since z~0.6, while the number of massive quiescent galaxies remains approximately constant since z~1. These results suggest that the rate at which star-forming galaxies are quenched increases with decreasing stellar mass, but that the bulk of the stellar mass buildup within the quiescent population occurs around ~10^10.8 Msun. In addition, we conclude that mergers do not appear to be a dominant channel for the stellar mass buildup of galaxies at z<1, even among massive (>10^11 Msun) quiescent galaxies.

The PRIsm MUti-object Survey (PRIMUS) is a spectroscopic galaxy redshift survey to z~1 completed with a low-dispersion prism and slitmasks allowing for simultaneous observations of ~2,500 objects over 0.18 square degrees. The final PRIMUS catalog includes ~130,000 robust redshifts over 9.1 sq. deg. In this paper, we summarize the PRIMUS observational strategy and present the data reduction details used to measure redshifts, redshift precision, and survey completeness. The survey motivation, observational techniques, fields, target selection, slitmask design, and observations are presented in Coil et al 2010. Comparisons to existing higher-resolution spectroscopic measurements show a typical precision of sigma_z/(1+z)=0.005. PRIMUS, both in area and number of redshifts, is the largest faint galaxy redshift survey completed to date and is allowing for precise measurements of the relationship between AGNs and their hosts, the effects of environment on galaxy evolution, and the build up of galactic systems over the latter half of cosmic history.

We report on the small scale (0.5<r<40h^-1 Mpc) clustering of 78895 massive (M*~10^11.3M_sun) galaxies at 0.2<z<0.4 from the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS), to be released as part of SDSS Data Release 9 (DR9). We describe the sample selection, basic properties of the galaxies, and caveats for working with the data. We calculate the real- and redshift-space two-point correlation functions of these galaxies, fit these measurements using Halo Occupation Distribution (HOD) modeling within dark matter cosmological simulations, and estimate the errors using mock catalogs. These galaxies lie in massive halos, with a mean halo mass of 5.2x10^13 h^-1 M_sun, a large scale bias of ~2.0, and a satellite fraction of 12+/-2%. Thus, these galaxies occupy halos with average masses in between those of the higher redshift BOSS CMASS sample and the original SDSS I/II LRG sample.

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

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

We measure the mass growth rate by merging for a wide range of galaxy types. We present the small-scale (0.014 < r < 11 h70^{-1} Mpc) projected cross-correlation functions w(rp) of galaxy subsamples from the spectroscopic sample of the NYU VAGC (5 \times 10^5 galaxies of redshifts 0.03 < z < 0.15) with galaxy subsamples from the SDSS imaging (4 \times 10^7 galaxies). We use smooth fits to de-project the two-dimensional functions w(rp) to obtain smooth three-dimensional real-space cross-correlation functions \xi(r) for each of several spectroscopic subsamples with each of several imaging subsamples. Because close pairs are expected to merge, the three-space functions and dynamical evolution time estimates provide galaxy accretion rates. We find that the accretion onto massive blue galaxies and onto red galaxies is dominated by red companions, and that onto small-mass blue galaxies, red and blue galaxies make comparable contributions. We integrate over all types of companions and find that at fixed stellar mass, the total fractional accretion rates onto red galaxies (\sim 1.5 h70 percent per Gyr) is greater than that onto blue galaxies (\sim 0.5 h70 percent per Gyr). Although these rates are very low, they are almost certainly over-estimates because we have assumed that all close pairs merge as quickly as dynamical friction permits.

We demonstrate that dwarf galaxies (10^7 < M_stellar < 10^9 Msun) with no active star formation are extremely rare (<0.06%) in the field. Our sample is based on the NASA-Sloan Atlas which is a re-analysis of the Sloan Digital Sky Survey Data Release 8. We examine the relative number of quenched versus star forming dwarf galaxies, defining quenched galaxies as having no Halpha emission (EW_Halpha < 2 AA) and a strong 4000AA-break. The fraction of quenched dwarf galaxies decreases rapidly with increasing distance from a massive host, leveling off for distances beyond 1.5 Mpc. We define galaxies beyond 1.5 Mpc of a massive host galaxy to be in the field. We demonstrate that there is a stellar mass threshold of M_stellar < 1.0x10^9 Msun below which quenched galaxies do not exist in the field. Below this threshold, we find that none of the 2951 field dwarf galaxies are quenched; all field dwarf galaxies show evidence for recent star formation. Correcting for volume effects, this corresponds to a 1-sigma upper limit on the quenched fraction of 0.06%. In more dense environments, quenched galaxies account for 23% of the dwarf population over the same stellar mass range. The majority of quenched dwarf galaxies (often classified as dwarf elliptical galaxies) are within 2 virial radii of a massive galaxy, and only a few percent of quenched dwarf galaxies exist beyond 4 virial radii. Thus, for galaxies with stellar mass less than 1.0x10^9 Msun, ending star-formation requires the presence of a more massive neighbor, providing a stringent constraint on models of star formation feedback.

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 analyze the density field of galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) included in the SDSS Data Release Nine (DR9). DR9 includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of 3,275 deg^2. We identify, characterize, and mitigate the impact of sources of systematic uncertainty on large-scale clustering measurements, both for angular moments of the redshift-space correlation function and the spherically averaged power spectrum, P(k), in order to ensure that robust cosmological constraints will be obtained from these data. A correlation between the projected density of stars and the higher redshift (0.43 < z < 0.7) galaxy sample (the `CMASS' sample) due to imaging systematics imparts a systematic error that is larger than the statistical error of the clustering measurements at scales s > 120h^-1Mpc or k < 0.01hMpc^-1. We find that these errors can be ameliorated by weighting galaxies based on their surface brightness and the local stellar density. We use mock galaxy catalogs that simulate the CMASS selection function to determine that randomly selecting galaxy redshifts in order to simulate the radial selection function of a random sample imparts the least systematic error on correlation function measurements and that this systematic error is negligible for the spherically averaged correlation function. The methods we recommend for the calculation of clustering measurements using the CMASS sample are adopted in companion papers that locate the position of the baryon acoustic oscillation feature (Anderson et al. 2012), constrain cosmological models using the full shape of the correlation function (Sanchez et al. 2012), and measure the rate of structure growth (Reid et al. 2012). (abridged)

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5\sigma in both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7\sigma. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

We present evidence that the incidence of active galactic nuclei (AGNs) and the distribution of their accretion rates do not depend on the stellar masses of their host galaxies, contrary to previous studies. We use hard (2-10 keV) X-ray data from three extragalactic fields (XMM-LSS, COSMOS and ELAIS-S1) with redshifts from the Prism Multi-object Survey to identify 242 AGNs with L_{2-10 keV}=10^{42-44} erg /s within a parent sample of ~25,000 galaxies at 0.2<z<1.0 over ~3.4 deg^2 and to i~23. We find that although the fraction of galaxies hosting an AGN at fixed X-ray luminosity rises strongly with stellar mass, the distribution of X-ray luminosities is independent of mass. Furthermore, we show that the probability that a galaxy will host an AGN can be defined by a universal Eddington ratio distribution that is independent of the host galaxy stellar mass and has a power-law shape with slope -0.65. These results demonstrate that AGNs are prevalent at all stellar masses in the range 9.5<log M_*/M_sun<12 and that the same physical processes regulate AGN activity in all galaxies in this stellar mass range. While a higher AGN fraction may be observed in massive galaxies, this is a selection effect related to the underlying Eddington ratio distribution. We also find that the AGN fraction drops rapidly between z~1 and the present day and is moderately enhanced (factor~2) in galaxies with blue or green optical colors. Consequently, while AGN activity and star formation appear to be globally correlated, we do not find evidence that the presence of an AGN is related to the quenching of star formation or the color transformation of galaxies.

The Sloan Digital Sky Survey (SDSS) surveyed 14,555 square degrees, and delivered over a trillion pixels of imaging data. We present a study of galaxy clustering using 900,000 luminous galaxies with photometric redshifts, spanning between $z=0.45$ and $z=0.65$, constructed from the SDSS using methods described in Ross et al. (2011). This data-set spans 11,000 square degrees and probes a volume of $3h^{-3} \rm{Gpc}^3$, making it the largest volume ever used for galaxy clustering measurements. We present a novel treatment of the observational systematics and its applications to the clustering signals from the data set. In this paper, we measure the angular clustering using an optimal quadratic estimator at 4 redshift slices with an accuracy of ~15% with bin size of delta_l = 10 on scales of the Baryon Acoustic Oscillations (BAO) (at l~40-400). We derive cosmological constraints using the full-shape of the power-spectra. For a flat Lambda CDM model, when combined with Cosmic Microwave Background Wilkinson Microwave Anisotropy Probe 7 (WMAP7) and H_0 constraints from 600 Cepheids observed by HST, we find \Omega_\Lambda = 0.73 +/- 0.019 and H_0 to be 70.5 +/- 1.6 km/s/Mpc. For an open Lambda CDM model, when combined with WMAP7 + HST, we find $\Omega_K = 0.0035 +/- 0.0054, improved over WMAP7+HST alone by 40%. For a wCDM model, when combined with WMAP7+HST+SN, we find w = -1.071 +/- 0.078, and H_0 to be 71.3 +/- 1.7 km/s/Mpc, which is competitive with the latest large scale structure constraints from large spectroscopic surveys such as SDSS Data Release 7 (DR7) (Reid et al. 2010, Percival et al. 2010, Montesano et al. 2011) and WiggleZ (Blake et al. 2011). The SDSS-III Data Release 8 (SDSS-III DR8) Angular Clustering Data allows a wide range of investigations into the cosmological model, cosmic expansion (via BAO), Gaussianity of initial conditions and neutrino masses. (abridged)

Passive red galaxies frequently contain warm ionized gas and have spectra similar to low-ionization nuclear emission-line regions (LINERs). Here we investigate the nature of the ionizing sources powering this emission, by comparing nuclear spectroscopy from the Palomar survey with larger aperture data from the Sloan Digital Sky Survey. We find the line emission in the majority of passive red galaxies is spatially extended; the Halpha surface brightness profile depends on radius (r) as r^(-1.28). We detect strong line ratio gradients with radius in [N II]/Ha, [S II]/Ha, and [O III]/[S II], requiring the ionization parameter to increase outwards. Combined with a realistic gas density profile, this outward increasing ionization parameter convincingly rules out AGN as the dominant ionizing source, and strongly favors distributed ionizing sources. Sources that follow the stellar density profile can additionally reproduce the observed luminosity-dependence of the line ratio gradient. Post-AGB stars provide a natural ionization source candidate, though they have an ionization parameter deficit. Velocity width differences among different emission lines disfavor shocks as the dominant ionization mechanism, and suggest that the interstellar medium in these galaxies contains multiple components. We conclude that the line emission in most LINER-like galaxies found in large aperture (>100pc) spectroscopy is not primarily powered by AGN activity and thus does not trace the AGN bolometric luminosity. However, they can be used to trace warm gas in these red galaxies.

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

We study the optical properties of the host galaxies of nuclear 22 GHz (lambda=1.35 cm) water masers. To do so, we cross-match the galaxy sample surveyed for water maser emission (123 detections and 3806 non-detections) with the SDSS low-redshift galaxy sample (z<0.05). Out of 1636 galaxies with SDSS photometry, we identify 48 detections; out of the 1063 galaxies that also have SDSS spectroscopy, we identify 33 detections. We find that maser detection rate is higher at higher optical luminosity (M_B), larger velocity dispersion ($\sigma$), and higher [OIII]5007 luminosity, with [OIII]5007 being the dominant factor. These detection rates are essentially the result of the correlations of isotropic maser luminosity with all three of these variables. These correlations are natural if maser strength increases with central black hole mass and the level of AGN activity. We also find that the detection rate is higher in galaxies with higher extinction. Based on these results, we propose that maser surveys seeking to efficiently find masers should rank AGN targets by extinction-corrected [OIII]5007 flux when available. This prioritization would improve maser detection efficiency, from an overall ~ 3% without pre-selection to ~ 16% for the strongest intrinsic [OIII]5007 emitters, by a factor of ~ 5.

We present the PRIsm MUlti-object Survey (PRIMUS), a spectroscopic faint galaxy redshift survey to z~1. PRIMUS uses a low-dispersion prism and slitmasks to observe ~2,500 objects at once in a 0.18 deg^2 field of view, using the IMACS camera on the Magellan I Baade 6.5m telescope at Las Campanas Observatory. PRIMUS covers a total of 9.1 deg^2 of sky to a depth of i_AB~23.5 in seven different deep, multi-wavelength fields that have coverage from GALEX, Spitzer and either XMM or Chandra, as well as multiple-band optical and near-IR coverage. PRIMUS includes ~130,000 robust redshifts of unique objects with a redshift precision of dz/(1+z)~0.005. The redshift distribution peaks at z=0.6 and extends to z=1.2 for galaxies and z=5 for broad-line AGN. The motivation, observational techniques, fields, target selection, slitmask design, and observations are presented here, with a brief summary of the redshift precision; a companion paper presents the data reduction, redshift fitting, redshift confidence, and survey completeness. PRIMUS is the largest faint galaxy survey undertaken to date. The high targeting fraction (~80%) and large survey size will allow for precise measures of galaxy properties and large-scale structure to z~1.

We investigate in detail the probability distribution function (pdf) of the proper-motion measurement errors in the SDSS+USNO-B proper-motion catalog of \citet{mun04} using clean quasar samples. The pdf of the errors is well-represented by a Gaussian core with extended wings, plus a very small fraction ($<0.1%$) of "outliers". We find while formally the pdf could be well-fit by a five-parameter fitting function, for many purposes it is also adequately to represent the pdf with a one-parameter approximation to this function. We apply this pdf to the calculation of the confidence intervals on the true proper motion for a SDSS+USNO-B proper motion measurement, and discuss several scientific applications of the SDSS proper motion catalogue. Our results have various applications in studies of the galactic structure and stellar kinematics. Specifically, they are crucial for searching hyper-velocity stars in the Galaxy.

We investigate the effects of potential sources of systematic error on the angular and photometric redshift, z_phot, distributions of a sample of redshift 0.4 < z < 0.7 massive galaxies whose selection matches that of the Baryon Oscillation Spectroscopic Survey (BOSS) constant mass sample. Utilizing over 112,778 BOSS spectra as a training sample, we produce a photometric redshift catalog for the galaxies in the SDSS DR8 imaging area that, after masking, covers nearly one quarter of the sky (9,913 square degrees). We investigate fluctuations in the number density of objects in this sample as a function of Galactic extinction, seeing, stellar density, sky background, airmass, photometric offset, and North/South Galactic hemisphere. We find that the presence of stars of comparable magnitudes to our galaxies (which are not traditionally masked) effectively remove area. Failing to correct for such stars can produce systematic errors on the measured angular auto-correlation function, w, that are larger than its statistical uncertainty. We describe how one can effectively mask for the presence of the stars, without removing any galaxies from the sample, and minimize the systematic error. Additionally, we apply two separate methods that can be used to correct the systematic errors imparted by any parameter that can be turned into a map on the sky. We find that failing to properly account for varying sky background introduces a systematic error on w. We measure w, in four z_phot slices of width 0.05 between 0.45 < z_phot < 0.65 and find that the measurements, after correcting for the systematic effects of stars and sky background, are generally consistent with a generic LambdaCDM model, at scales up to 60 degrees. At scales greater than 3 degrees and z_phot > 0.5, the magnitude of the corrections we apply are greater than the statistical uncertainty in w.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

We describe a procedure for background subtracting Sloan Digital Sky Survey (SDSS) imaging that improves the resulting detection and photometry of large galaxies on the sky. Within each SDSS drift scan run, we mask out detected sources and then fit a smooth function to the variation of the sky background. This procedure has been applied to all SDSS-III Data Release 8 images, and the results are available as part of that data set. We have tested the effect of our background subtraction on the photometry of large galaxies by inserting fake galaxies into the raw pixels, reanalyzing the data, and measuring them after background subtraction. Our technique results in no size-dependent bias in galaxy fluxes up to half-light radii of 100 arcsec; in contrast, for galaxies of that size the standard SDSS photometric catalog underestimates fluxes by about 1.5 mag. Our results represent a substantial improvement over the standard SDSS catalog results and should form the basis of any analysis of nearby galaxies using the SDSS imaging data.

The baryonic acoustic feature in galaxy clustering is a promising tool for constraining the nature of the cosmic acceleration, through measurements of expansion rates H and angular diameter distances D_A. Angle-averaged measurements of clustering yield constraints on the quantity D_A^2/H. However, to break the degeneracy between these two parameters one must measure the anisotropic correlation function as a function of both line-of-sight (radial) and transverse separations. Here we investigate how to most effectively do so, using analytic techniques and mock catalogues. In particular, we examine multipole expansions of the correlation function as well as "clustering wedges" xi(Delta mu, s), where mu = s_||/s and s_|| is the radial component of separation s. Both techniques allow strong constraints on H and D_A, as expected. The radial wedges strongly depend on H and the transverse wedges are sensitive to D_A. Analyses around the region of the acoustic peak constrain H ~20% better when using the wedge statistics than when using the monopole-quadrupole combination. However, we show that the hexadecapole allows substantially stronger constraints than the monopole and quadrupole alone. Our findings here demonstrate that wedge statistics provide a practical alternative technique to multipoles, that should be useful to test systematics and will provide comparable or better constraints. Finally, we predict the constraints from galaxy clustering that will be possible with a completed version of the ongoing Baryonic Oscillation Spectroscopic Survey.

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

We measure the luminosity and color dependence of galaxy clustering in the SDSS DR7 main galaxy sample, focusing on the projected correlation function w_p(r_p) of volume-limited samples. We interpret our measurements using halo occupation distribution (HOD) modeling assuming a Lambda-CDM cosmology. The amplitude of w_p(r_p) grows slowly with luminosity for L < L_* and increases sharply at higher luminosities, with bias factor b(>L)=1.06+0.23(L/L_*)^{1.12}. At fixed luminosity, redder galaxies have a stronger and steeper w_p(r_p), a trend that runs steadily from the bluest galaxies to the reddest galaxies. The individual luminosity trends for the red and blue galaxy populations are strikingly different. Blue galaxies show a slow but steady increase of w_p(r_p) with luminosity, at all scales. The large-scale clustering of red galaxies shows little luminosity dependence until a sharp increase at L > 4L_*, but the lowest luminosity red galaxies (0.04-0.25 L_*) show very strong clustering on scales r_p < 2 Mpc/h. Most of the observed trends can be naturally understood within the LCDM+HOD framework. The growth of w_p(r_p) with luminosity reflects an overall shift in the halo mass scale, in particular an increase in the minimum host halo mass Mmin. The mass at which a halo has, on average, one satellite galaxy brighter than L is M_1 ~ 17 Mmin(L) over most of the luminosity range. The growth and steepening of w_p(r_p) for redder galaxies reflects the increasing fraction of galaxies that are satellite systems in high mass halos instead of central systems in low mass halos, a trend that is especially marked at low luminosities. Our extensive measurements, provided in tabular form, will allow detailed tests of theoretical models of galaxy formation, a firm grounding of semi-empirical models of the galaxy population, and new cosmological tests.

We place constraints on the average density (Omega_m) and clustering amplitude (sigma_8) of matter using a combination of two measurements from the Sloan Digital Sky Survey: the galaxy two-point correlation function, w_p, and the mass-to-galaxy-number ratio within galaxy clusters, M/N, analogous to cluster M/L ratios. Our w_p measurements are obtained from DR7 while the sample of clusters is the maxBCG sample, with cluster masses derived from weak gravitational lensing. We construct non-linear galaxy bias models using the Halo Occupation Distribution (HOD) to fit both w_p and M/N for different cosmological parameters. HOD models that match the same two-point clustering predict different numbers of galaxies in massive halos when Omega_m or sigma_8 is varied, thereby breaking the degeneracy between cosmology and bias. We demonstrate that this technique yields constraints that are consistent and competitive with current results from cluster abundance studies, even though this technique does not use abundance information. Using w_p and M/N alone, we find Omega_m^0.5*sigma_8=0.465+/-0.026, with individual constraints of Omega_m=0.29+/-0.03 and sigma_8=0.85+/-0.06. Combined with current CMB data, these constraints are Omega_m=0.290+/-0.016 and sigma_8=0.826+/-0.020. All errors are 1-sigma. The systematic uncertainties that the M/N technique are most sensitive to are the amplitude of the bias function of dark matter halos and the possibility of redshift evolution between the SDSS Main sample and the maxBCG sample. Our derived constraints are insensitive to the current level of uncertainties in the halo mass function and in the mass-richness relation of clusters and its scatter, making the M/N technique complementary to cluster abundances as a method for constraining cosmology with future galaxy surveys.

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

We calculate the real- and redshift-space clustering of massive galaxies at z~0.5 using the first semester of data by the Baryon Oscillation Spectroscopic Survey (BOSS). We study the correlation functions of a sample of 44,000 massive galaxies in the redshift range 0.4<z<0.7. We present a halo-occupation distribution modeling of the clustering results and discuss the implications for the manner in which massive galaxies at z~0.5 occupy dark matter halos. The majority of our galaxies are central galaxies living in halos of mass 10^{13}Msun/h, but 10% are satellites living in halos 10 times more massive. These results are broadly in agreement with earlier investigations of massive galaxies at z~0.5. The inferred large-scale bias (b~2) and relatively high number density (nbar=3e-4 h^3 Mpc^{-3}) imply that BOSS galaxies are excellent tracers of large-scale structure, suggesting BOSS will enable a wide range of investigations on the distance scale, the growth of large-scale structure, massive galaxy evolution and other topics.

Tidal interactions between galaxies can trigger star formation, which contributes to the global star formation rate density of the universe and could be a factor in the transformation of blue, star-forming galaxies to red, quiescent galaxies over cosmic time. We investigate tidally-triggered star formation in isolated close galaxy pairs drawn from the Prism Multi-Object Survey (PRIMUS), a low-dispersion prism redshift survey that has measured ~120,000 robust galaxy redshifts over 9.1 deg^2 out to z ~ 1. We select a sample of galaxies in isolated galaxy pairs at redshifts 0.25 < z < 0.75, with no other objects within a projected separation of 300 h^-1 kpc and dz/(1+z) = 0.01, and compare them to a control sample of isolated galaxies to test for systematic differences in their rest-frame FUV-r and NUV-r colors as a proxy for relative specific SFR. We find that galaxies in r_p < 50 h^-1 kpc pairs have bluer dust-corrected UV-r colors on average than the control galaxies by -0.134 +/- 0.045 magnitudes in FUV-r and -0.075 +/- 0.038 magnitudes in NUV-r, corresponding to a ~15-20% increase in SSFR. This indicates an enhancement in SSFR due to tidal interactions. We also find that this relative enhancement is greater for a subset of r_p < 30 h^-1 kpc pair galaxies, for which the average colors offsets are -0.193 +/- 0.065 magnitudes in FUV-r and -0.159 +/- 0.048 magnitudes in NUV-r, corresponding to a ~25-30% increase in SSFR. We test for evolution in the enhancement of tidally-triggered star formation with redshift across our sample redshift range and find marginal evidence for a decrease in SSFR enhancement from 0.25 < z < 0.5 to 0.5 < z < 0.75. This indicates that a change in enhanced star formation triggered by tidal interactions in low density environments is not a contributor to the decline in the global star formation rate density across this redshift range.

We quantify the fraction of galaxies at moderate redshifts (0.1<z<0.5) that appear red-and-dead in the optical, but in fact contain obscured star formation detectable in the infrared (IR), with the PRIsm MUlti-object Survey (PRIMUS). PRIMUS has measured ~120,000 robust redshifts with a precision of sigma_z/(1+z)~0.5% over 9.1 square degrees of the sky to the depth of i~23 (AB), up to redshift z~1. We specifically targeted 6.7 square degree fields with existing deep IR imaging from the Spitzer Space Telescope from the SWIRE and S-COSMOS surveys. We select in these fields an i band flux-limited sample (i<20 mag in the SWIRE fields and i<21 mag in the S-COSMOS field) of 3310 red-sequence galaxies at 0.1<z<0.5 for which we can reliably classify obscured star-forming and quiescent galaxies using IR color. Our sample constitutes the largest galaxy sample at intermediate redshift to study obscured star formation on the red sequence, and we present the first quantitative analysis of the fraction of obscured star-forming galaxies as a function of luminosity. We find that on average, at L ~ L*, about 15% of red-sequence galaxies have IR colors consistent with star-forming galaxies. The percentage of obscured star-forming galaxies increases by ~8% per mag with decreasing luminosity from the highest luminosities to L~0.2L*. Our results suggest that a significant fraction of red-sequence galaxies have ongoing star formation and that galaxy evolution studies based on optical color therefore need to account for this complication.

We measure the evolution of the [OII]lambda 3727 luminosity function at 0.75<z<1.45 using high-resolution spectroscopy of ~14,000 galaxies observed by the DEEP2 galaxy redshift survey. We find that brighter than L_{OII}=10^{42} erg s^(-1) the luminosity function is well-represented by a power law dN/dL ~ L^{\alpha} with slope \alpha ~ -3. The number density of [OII] emitting galaxies above this luminosity declines by a factor of >~2.5 between z ~ 1.35 and z ~ 0.84. In the limit of no number-density evolution, the characteristic [OII] luminosity, L^*_[OII], defined as the luminosity where the space density equals 10^{-3.5} dex^{-1} Mpc^{-3}, declines by a factor of ~1.8 over the same redshift interval. Assuming that L_[OII] is proportional to the star-formation rate (SFR), and negligible change in the typical dust attenuation in galaxies at fixed [OII] luminosity, the measured decline in L^*_[OII] implies a ~25% per Gyr decrease in the amount of star formation in galaxies during this epoch. Adopting a faint-end power-law slope of -1.3\pm0.2, we derive the comoving SFR density in four redshift bins centered around z~1 by integrating the observed [OII] luminosity function using a local, empirical calibration between L_[OII] and SFR, which statistically accounts for variations in dust attenuation and metallicity among galaxies. We find that our estimate of the SFR density at z~1 is consistent with previous measurements based on a variety of independent SFR indicators.

We study the stellar populations of 1,923 elliptical galaxies at z<0.05 selected from the Sloan Digital Sky Survey as a function of velocity dispersion, sigma, and environment. We construct average, high signal-to-noise spectra and find the following: (1) lower-sigma galaxies have a bluer optical continuum and stronger (but still weak) emission lines; (2) at fixed sigma, field ellipticals have a slightly bluer stellar continuum, especially at wavelengths \lesssim 4000 \AA, and have stronger (but still weak) emission lines compared to their group counterparts, although this environmental dependence is strongest for low-sigma ellipticals. Based on Lick indices measured from both the individual and average spectra, we find that: (1) at a given sigma, elliptical galaxies in groups have systematically weaker Balmer absorption than their field counterparts, although this environmental dependence is most pronounced at low sigma; (2) there is no clear environmental dependence of <Fe>, while the alpha-element absorption indices such as Mgb are only slightly stronger in galaxies belonging to rich groups. An analysis based on simple stellar populations (SSPs) reveals that more massive elliptical galaxies are older, more metal-rich and more strongly alpha-enhanced. We also find that: (1) the SSP-equivalent ages of galaxies in rich groups are, on average, ~1 Gyr older than in the field, although once again this effect is strongest at low sigma; (2) galaxies in rich groups have slightly lower [Fe/H] and are marginally more strongly alpha-enhanced; and (3) there is no significant environmental dependence of total metallicity, [Z/H]. Our results are generally consistent with stronger low-level recent star formation in field ellipticals at low sigma, similar to recent results based on ultraviolet and infrared observations. (Abridged)

We use measurements of the HI content, stellar mass and star formation rates in ~190 massive galaxies with stellar masses greater than 10^10 Msun, obtained from the Galex Arecibo SDSS Survey (GASS) described in Paper I (Catinella et al. 2010) to explore the global scaling relations associated with the bin-averaged ratio of the star formation rate over the HI mass, which we call the HI-based star formation efficiency (SFE). Unlike the mean specific star formation rate, which decreases with stellar mass and stellar mass surface density, the star formation efficiency remains relatively constant across the sample with a value close to SFE = 10^-9.5 yr^-1 (or an equivalent gas consumption timescale of ~3 Gyr). Specifically, we find little variation in SFE with stellar mass, stellar mass surface density, NUV-r color and concentration. We interpret these results as an indication that external processes or feedback mechanisms that control the gas supply are important for regulating star formation in massive galaxies. An investigation into the detailed distribution of SFEs reveals that approximately 5% of the sample shows high efficiencies with SFE > 10^-9 yr^-1, and we suggest that this is very likely due to a deficiency of cold gas rather than an excess star formation rate. Conversely, we also find a similar fraction of galaxies that appear to be gas-rich for their given specific star-formation rate, although these galaxies show both a higher than average gas fraction and lower than average specific star formation rate. Both of these populations are plausible candidates for "transition" galaxies, showing potential for a change (either decrease or increase) in their specific star formation rate in the near future. We also find that 36+/-5% of the total HI mass density and 47+/-5% of the total SFR density is found in galaxies with stellar mass greater than 10^10 Msun. [abridged]

The attenuation of starlight by interstellar dust is investigated in a sample of low redshift, disk-dominated star-forming galaxies using photometry from GALEX and SDSS. By considering broadband colors as a function of galaxy inclination we are able to confidently separate trends arising from increasing dust opacity from possible differences in stellar populations, since stellar populations do not correlate with inclination. All commonly employed dust attenuation curves (such as the Calzetti curve for starbursts, or a power-law curve) provide poor fits to the ultraviolet colors for moderately and highly inclined galaxies. This conclusion rests on the fact that the average FUV-NUV color varies little from face-on to edge-on galaxies, while other colors such as NUV-u and u-r vary strongly with inclination. After considering a number of model variations, we are led to speculate that the presence of the strong dust extinction feature at 2175A seen in the Milky Way (MW) extinction curve is responsible for the observed trends. Independent of our interpretation, these results imply that the modeling of dust attenuation in the ultraviolet is significantly more complicated than traditionally assumed. These results also imply a very weak dependence of the FUV-NUV color on total FUV attenuation, and we conclude from this that it is extremely difficult to use only the observed UV spectral slope to infer the total UV dust attenuation, as is commonly done. We propose several simple tests that might finally identify the grain population responsible for the 2175A feature.

We analyze the line-of-sight baryonic acoustic feature in the two-point correlation function {\xi} of the Sloan Digital Sky Survey (SDSS) luminous red galaxy (LRG) sample (0.16 < z < 0.47). By defining a narrow line-of-sight region, rp < 5.5 Mpc/h, where rp is the transverse separation component, we measure a strong excess of clustering at ~ 110 Mpc/h, as previously reported in the literature. We also test these results in an alternative coordinate system, by defining the line-of-sight as {\theta} < 3{\deg}, where {\theta} is the opening angle. This clustering excess appears much stronger than the feature in the better-measured monopole. A fiducial {\Lambda}CDM non-linear model in redshift-space predicts a much weaker signature. We use realistic mock catalogs to model the expected signal and noise. We find that the line-of-sight measurements can be explained well by our mocks as well as by a featureless {\xi} = 0. We conclude that there is no convincing evidence that the strong clustering measurement is the line-of-sight baryonic acoustic feature. We also evaluate how detectable such a signal would be in the upcoming Baryon Oscillation Spectroscopic Survey LRG volume (BOSS). Mock LRG catalogs (z < 0.6) suggest that: (i) the narrow line- of-sight cylinder and cone defined above probably will not reveal a detectable acoustic feature in BOSS; (ii) a clustering measurement as high as that in the current sample can be ruled out (or confirmed) at a high confidence level using a BOSS-sized data set; and (iii) an analysis with wider angular cuts, which provide better signal-to-noise ratios, can nevertheless be used to compare line-of-sight and transverse distances, and thereby constrain the expansion rate H(z) and diameter distance DA(z).

We present the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog, which is based upon the SDSS Seventh Data Release. The catalog, which contains 105,783 spectroscopically confirmed quasars, represents the conclusion of the SDSS-I and SDSS-II quasar survey. The catalog consists of the SDSS objects that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s or have interesting/complex absorption features, are fainter than i > 15.0 and have highly reliable redshifts. The catalog covers an area of 9380 deg^2. The quasar redshifts range from 0.065 to 5.46, with a median value of 1.49; the catalog includes 1248 quasars at redshifts greater than four, of which 56 are at redshifts greater than five. The catalog contains 9210 quasars with i < 18; slightly over half of the entries have i< 19. For each object the catalog presents positions accurate to better than 0.1" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800-9200 Ang. at a spectral resolution R = 2000 the spectra can be retrieved from the SDSS public database using the information provided in the catalog. Over 96% of the objects in the catalog were discovered by the SDSS. We also include a supplemental list of an additional 207 quasars with SDSS spectra whose archive photometric information is incomplete.

We introduce the GALEX Arecibo SDSS Survey (GASS), an on-going large program that is gathering high quality HI-line spectra using the Arecibo radio telescope for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05, selected from the SDSS spectroscopic and GALEX imaging surveys. The galaxies are observed until detected or until a low gas mass fraction limit (1.5-5%) is reached. This paper presents the first Data Release, consisting of ~20% of the final GASS sample. We use this data set to explore the main scaling relations of HI gas fraction with galaxy structure and NUV-r colour. A large fraction (~60%) of the galaxies in our sample are detected in HI. We find that the atomic gas fraction decreases strongly with stellar mass, stellar surface mass density and NUV-r colour, but is only weakly correlated with galaxy bulge-to-disk ratio (as measured by the concentration index of the r-band light). We also find that the fraction of galaxies with significant (more than a few percent) HI decreases sharply above a characteristic stellar surface mass density of 10^8.5 Msun kpc^-2. The fraction of gas-rich galaxies decreases much more smoothly with stellar mass. One of the key goals of GASS is to identify and quantify the incidence of galaxies that are transitioning between the blue, star-forming cloud and the red sequence of passively-evolving galaxies. Likely transition candidates can be identified as outliers from the mean scaling relations between gas fraction and other galaxy properties. [abridged]

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We have built a reliable and robust system that takes as input an astronomical image, and returns as output the pointing, scale, and orientation of that image (the astrometric calibration or WCS information). The system requires no first guess, and works with the information in the image pixels alone; that is, the problem is a generalization of the "lost in space" problem in which nothing--not even the image scale--is known. After robust source detection is performed in the input image, asterisms (sets of four or five stars) are geometrically hashed and compared to pre-indexed hashes to generate hypotheses about the astrometric calibration. A hypothesis is only accepted as true if it passes a Bayesian decision theory test against a background hypothesis. With indices built from the USNO-B Catalog and designed for uniformity of coverage and redundancy, the success rate is 99.9% for contemporary near-ultraviolet and visual imaging survey data, with no false positives. The failure rate is consistent with the incompleteness of the USNO-B Catalog; augmentation with indices built from the 2MASS Catalog brings the completeness to 100% with no false positives. We are using this system to generate consistent and standards-compliant meta-data for digital and digitized imaging from plate repositories, automated observatories, individual scientific investigators, and hobbyists. This is the first step in a program of making it possible to trust calibration meta-data for astronomical data of arbitrary provenance.

We review the physical properties of nearby, relatively luminous galaxies, using results from newly available massive data sets together with more detailed observations. First, we present the global distribution of properties, including the optical and ultraviolet luminosity, stellar mass, and atomic gas mass functions. Second, we describe the shift of the galaxy population from "late" galaxy types in underdense regions to "early" galaxy types in overdense regions. We emphasize that the scaling relations followed by each galaxy type change very little with environment, with the exception of some minor but detectable effects. The shift in the population is apparent even at the densities of small groups and therefore cannot be exclusively due to physical processes operating in rich clusters. Third, we divide galaxies into four crude types -- spiral, lenticular, elliptical, and merging systems -- and describe some of their more detailed properties. We attempt to put these detailed properties into the global context provided by large surveys.

We examine the correlation function \xi of the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxy sample (LRG) at large scales (60<s<400 Mpc/h) using the final data release (DR7; 105,831 LRGs between 0.16<z<0.47). Using mock catalogs, we demonstrate that the observed baryonic acoustic peak and larger scale signal are consistent with LCDM at the 1.5\sigma level. The signal at 155<s<200 Mpc/h tends to be high relative to theoretical expectations; this slight deviation can be attributed to a bright subsample of the LRGs. Fitting data to a non-linear, redshift-space, template based-model, we constrain the peak position at s_p=103.6+3.6-2.4 Mpc/h when fitting the range 60<s<150 Mpc/h (1\sigma uncertainties measured from the mocks. This redshift-space distance s_p is related to the comoving sound horizon scale r_s after taking into account matter clustering non-linearities, redshift distortions and galaxy clustering bias. Mock catalogs show that the probability that a DR7-sized sample would not have an identifiable peak is at least 10%. As a consistency check of a fiducial cosmology, we use the observed s_p to obtain the distance D_V=[(1+z)^2D_A^2cz/H(z)]^(1/3) relative to the acoustic scale. We find r_s/D_V(z=0.278)=0.1394+-0.0049. This result is in excellent agreement with Percival et. al (2009), who examine roughly the same data set, but using the power spectrum. Comparison with other determinations in the literature are also in very good agreement. We have tested our results against a battery of possible systematic effects, finding all effects are smaller than our estimated sample variance.

Astronomical software is now a fact of daily life for all hands-on members of our community. Purpose-built software for data reduction and modeling tasks becomes ever more critical as we handle larger amounts of data and simulations. However, the writing of astronomical software is unglamorous, the rewards are not always clear, and there are structural disincentives to releasing software publicly and to embedding it in the scientific literature, which can lead to significant duplication of effort and an incomplete scientific record. We identify some of these structural disincentives and suggest a variety of approaches to address them, with the goals of raising the quality of astronomical software, improving the lot of scientist-authors, and providing benefits to the entire community, analogous to the benefits provided by open access to large survey and simulation datasets. Our aim is to open a conversation on how to move forward. We advocate that: (1) the astronomical community consider software as an integral and fundable part of facility construction and science programs; (2) that software release be considered as integral to the open and reproducible scientific process as are publication and data release; (3) that we adopt technologies and repositories for releasing and collaboration on software that have worked for open-source software; (4) that we seek structural incentives to make the release of software and related publications easier for scientist-authors; (5) that we consider new ways of funding the development of grass-roots software; (6) and that we rethink our values to acknowledge that astronomical software development is not just a technical endeavor, but a fundamental part of our scientific practice.

Wide-angle surveys have been an engine for new discoveries throughout the modern history of astronomy, and have been among the most highly cited and scientifically productive observing facilities in recent years. This trend is likely to continue over the next decade, as many of the most important questions in astrophysics are best tackled with massive surveys, often in synergy with each other and in tandem with the more traditional observatories. We argue that these surveys are most productive and have the greatest impact when the data from the surveys are made public in a timely manner. The rise of the "survey astronomer" is a substantial change in the demographics of our field; one of the most important challenges of the next decade is to find ways to recognize the intellectual contributions of those who work on the infrastructure of surveys (hardware, software, survey planning and operations, and databases/data distribution), and to make career paths to allow them to thrive.

We present measurements of the excess mass-to-light ratio measured aroundMaxBCG galaxy clusters observed in the SDSS. This red sequence cluster sample includes objects from small groups with masses ranging from ~5x10^{12} to ~10^{15} M_{sun}/h. Using cross-correlation weak lensing, we measure the excess mass density profile above the universal mean \Delta \rho(r) = \rho(r) - \bar{\rho} for clusters in bins of richness and optical luminosity. We also measure the excess luminosity density \Delta l(r) = l(r) - \bar{l} measured in the z=0.25 i-band. For both mass and light, we de-project the profiles to produce 3D mass and light profiles over scales from 25 kpc/ to 22 Mpc/h. From these profiles we calculate the cumulative excess mass M(r) and excess light L(r) as a function of separation from the BCG. On small scales, where \rho(r) >> \bar{\rho}, the integrated mass-to-light profile may be interpreted as the cluster mass-to-light ratio. We find the M/L_{200}, the mass-to-light ratio within r_{200}, scales with cluster mass as a power law with index 0.33+/-0.02. On large scales, where \rho(r) ~ \bar{\rho}, the M/L approaches an asymptotic value independent of cluster richness. For small groups, the mean M/L_{200} is much smaller than the asymptotic value, while for large clusters it is consistent with the asymptotic value. This asymptotic value should be proportional to the mean mass-to-light ratio of the universe <M/L>. We find <M/L>/b^2_{ml} = 362+/-54 h (statistical). There is additional uncertainty in the overall calibration at the ~10% level. The parameter b_{ml} is primarily a function of the bias of the L <~ L_* galaxies used as light tracers, and should be of order unity. Multiplying by the luminosity density in the same bandpass we find \Omega_m/b^2_{ml} = 0.02+/-0.03, independent of the Hubble parameter.

Differences in clustering properties between galaxy subpopulations complicate the cosmological interpretation of the galaxy power spectrum, but can also provide insights about the physics underlying galaxy formation. To study the nature of this relative clustering, we perform a counts-in-cells analysis of galaxies in the Sloan Digital Sky Survey (SDSS) in which we measure the relative bias between pairs of galaxy subsamples of different luminosities and colours. We use a generalized chi-squared test to determine if the relative bias between each pair of subsamples is consistent with the simplest deterministic linear bias model, and we also use a maximum likelihood technique to further understand the nature of the relative bias between each pair. We find that the simple, deterministic model is a good fit for the luminosity-dependent bias on scales above ~2 Mpc/h, which is good news for using magnitude-limited surveys for cosmology. However, the colour-dependent bias shows evidence for stochasticity and/or non-linearity which increases in strength toward smaller scales, in agreement with previous studies of stochastic bias. Also, confirming hints seen in earlier work, the luminosity-dependent bias for red galaxies is significantly different from that of blue galaxies: both luminous and dim red galaxies have higher bias than moderately bright red galaxies, whereas the biasing of blue galaxies is not strongly luminosity-dependent. These results can be used to constrain galaxy formation models and also to quantify how the colour and luminosity selection of a galaxy survey can impact measurements of the cosmological matter power spectrum.

Galaxies are seen from different viewing angles and their observed properties change as a function of viewing angle. In many circumstances we would rather know the intrinsic properties of galaxies - those properties that do not depend on viewing angle. For a large sample of galaxies it is possible to recover the intrinsic properties of galaxies, statistically, by looking for correlations with galaxy inclination, and then applying a correction to remove those dependencies. Studying the intrinsic properties of galaxies can give a different impression of the galaxy population and help avoid the mistake of connecting observed properties to quantities that don't depend on inclination like halo mass.

The observed properties of galaxies vary with inclination; for most applications we would rather have properties that are independent of inclination, intrinsic properties. One way to determine inclination corrections is to consider a large sample of galaxies, study how the observed properties of these galaxies depend on inclination and then remove this dependence to recover the intrinsic properties. We perform such an analysis for galaxies selected from the Sloan Digital Sky Survey which have been matched to galaxies from the Two-Micron All Sky Survey. We determine inclination corrections for these galaxies as a function of galaxy luminosity and Sersic index. In the g-band these corrections reach as as high as 1.2 mag and have a median value of 0.3 mag for all galaxies in our sample. We find that the corrections show little dependence on galaxy luminosity, except in the $u$ band, but are strongly dependent on galaxy Sersic index. We find that the ratio of red-to-blue galaxies changes from 1:1 to 1:2 when going from observed to intrinsic colors for galaxies in the range -22.75 < M_K < -17.75. We also discuss how survey completeness and photometric redshifts should be determined when taking into account that observed and intrinsic properties differ. Finally, we examine whether previous determinations of stellar mass give an intrinsic quantity or one that depends on galaxy inclination.

We present an algorithm to photometrically calibrate wide field optical imaging surveys, that simultaneously solves for the calibration parameters and relative stellar fluxes using overlapping observations. The algorithm decouples the problem of "relative" calibrations, from that of "absolute" calibrations; the absolute calibration is reduced to determining a few numbers for the entire survey. We pay special attention to the spatial structure of the calibration errors, allowing one to isolate particular error modes in downstream analyses. Applying this to the Sloan Digital Sky Survey imaging data, we achieve ~1% relative calibration errors across 8500 sq.deg. in griz; the errors are ~2% for the u band. These errors are dominated by unmodelled atmospheric variations at Apache Point Observatory. These calibrations, dubbed "ubercalibration", are now public with SDSS Data Release 6, and will be a part of subsequent SDSS data releases.

We study the role of major and minor mergers in the mass growth of luminous red galaxies. We present small-scale ($0.01<r<8\,\hMpc$) projected cross-correlation functions of $23043$ luminous early-type galaxies from the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxy (LRG) sample ($0.16<z<0.30$, $\MMi\approx -22.75\,\mag$) with all their companions in the SDSS imaging sample, split into color and luminosity subsamples with $\MMi<-18\,\mag$. We de-project the two-dimensional functions to obtain three-dimensional real-space LRG--galaxy cross-correlation functions for each companion subsample. We find that the cross-correlation functions are not purely power-law and that there is a clear ``one-halo'' to ``two-halo'' transition near $1\,\hMpc$. We convert these results into close pair statistics and estimate the LRG accretion rate from each companion galaxy subsample using timescales from dynamical friction arguments for each subsample of the companions. We find that the accretion onto LRGs is dominated by dry mergers of galaxies more luminous than $\Lstar$. We integrate the luminosity accretion rate from mergers over all companion galaxy subsamples and find that LRGs are growing by $[1.7\pm 0.1]$ percent per $\Gyr$, on average, from merger activity at redshift $z\sim 0.25$. This rate is almost certainly an over-estimate because we have assumed that all close pairs are merging as quickly as dynamical friction allows; nonetheless it is on the low side of the panoply of measurements in the literature, and lower than any rate predicted from theory.

The galaxy circular velocity function at small masses is related to the matter power spectrum on small scales. Although this function is well-studied for Local Group dwarfs, theoretical predictions and observational measurements are difficult for satellite galaxies, because of ram pressure and tidal stripping. By contrast, isolated dwarf galaxies are less affected by these processes, and almost always have enough 21cm emission to trace their dynamics robustly. Here, we test cold dark matter cosmology using isolated low mass dwarf galaxies from the SDSS with measured 21cm widths. We find consistency between the predicted and observed number density of isolated galaxies down to circular velocities of 50 km/s. Our technique yields a direct test of small-scale cosmology independent of the Lyman-alpha forest power spectrum, but our sample is currently statistically less powerful: warm dark matter particles heavier than 0.5 keV cannot be ruled out. Our major systematic uncertainty is the surface brightness limit of the SDSS. Blind HI surveys, such as the ALFALFA survey on Arecibo, will uncover a larger number of isolated low mass galaxies and increase the power of our constraints. With our sample, we also find that the Tully-Fisher relation for dwarf galaxies is a strong function of environment, and that the baryonic fraction is only a weak function of mass. These results suggest that for dwarf galaxies, gas loss is dominated by external, not internal, processes. [abridged]

We measure the relation between galaxy luminosity and disk circular velocity (the Tully-Fisher [TF] relation), in the g, r, i, and z-bands, for a broadly selected sample of galaxies from the Sloan Digital Sky Survey, with the goal of providing well defined observational constraints for theoretical models of galaxy formation. The input sample of 234 galaxies has a roughly flat distribution of absolute magnitudes in the range -18.5 > Mr > -22, and our only morphological selection is an axis-ratio cut b/a < 0.6 to allow accurate inclination corrections. Long-slit spectroscopy yields usable H-alpha rotation curves for 162 galaxies. Observational errors, including distance errors due to peculiar velocities, are small compared to the intrinsic scatter of the TF relation. The slope of the forward TF relation steepens from -5.5 +/- 0.2 mag/log(km/s) in the g-band to -6.6 +/- 0.2 mag/log(km/s) in the z-band. The intrinsic scatter is approximately 0.4 mag in all bands. The scatter is not dominated by rare outliers or by any particular class of galaxies, though it drops slightly, to 0.36 mag, if we restrict the sample to nearly bulgeless systems. Correlations of TF residuals with other galaxy properties are weak: bluer galaxies are significantly brighter than average in the g-band but only marginally brighter in the i-band; more concentrated galaxies are slightly fainter than average; and the TF residual is virtually independent of half-light radius, contrary to the trend expected for gravitationally dominant disks. The observed residual correlations do not account for most of the intrinsic scatter, implying that this scatter is instead driven largely by variations in the ratio of dark to luminous matter within the disk galaxy population.

We present the fourth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog. The catalog contains 77,429 objects; this is an increase of over 30,000 entries since the previous edition. The catalog consists of the objects in the SDSS Fifth Data Release that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s, or have interesting/complex absorption features, are fainter than i=15.0, and have highly reliable redshifts. The area covered by the catalog is 5740 sq. deg. The quasar redshifts range from 0.08 to 5.41, with a median value of 1.48; the catalog includes 891 quasars at redshifts greater than four, of which 36 are at redshifts greater than five. Approximately half of the catalog quasars have i < 19; nearly all have i < 21. For each object the catalog presents positions accurate to better than 0.2 arcsec. rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains basic radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800--9200A at a spectral resolution of ~2000. The spectra can be retrieved from the public database using the information provided in the catalog. The average SDSS colors of quasars as a function of redshift, derived from the catalog entries, are presented in tabular form. Approximately 96% of the objects in the catalog were discovered by the SDSS.

We investigate the dependence of physical properties of galaxies brighter than M_r=-18.0 in the SDSS on environment, as measured by local density using an adaptive smoothing kernel. We find that variations of galaxy properties with environment are almost entirely due to the dependence of morphology and luminosity on environment. Because galaxy properties depend not only on luminosity but also on morphology, it is clear that galaxy properties cannot be determined solely by dark halo mass. When morphology and luminosity are fixed, other physical properties, such as color, color-gradient, concentration, size, velocity dispersion, and star formation rate, are nearly independent of local density. The only feature is the sharp decrease of late type fraction above the critical luminosity of about M_r=-21.3 in the morphology versus luminosity relation. Weak residual dependences on environment include that of the color of late-types (bluer at lower density) and of the L-sigma relation of early-types (larger dispersion at higher density for bright galaxies). The early-type fraction is a monotonically increasing function of local density and luminosity. The morphology-density- luminosity relation should be a key constraint on galaxy formation models. We demonstrate that the dependence on environment of the morphology of galaxie originates from variations in density on effective Gaussian smoothing scales much smaller than 12 h^{-1}Mpc. We find that galaxy morphology varies both with density measured on an effective Gaussian smoothing scale of 4.7 h^{-1} Mpc and on distance to the nearest bright galaxy (particularly, for about 0.2 h^{-1}Mpc). We propose a mechanism that the morphology of galaxies in galaxy systems is transformed by the tidal force.

In addition to optical photometry of unprecedented quality, the Sloan Digital Sky Survey (SDSS) is also producing a massive spectroscopic database. We discuss determination of stellar parameters, such as effective temperature, gravity and metallicity from SDSS spectra, describe correlations between kinematics and metallicity, and study their variation as a function of the position in the Galaxy. We show that stellar parameter estimates by Beers et al. show a good correlation with the position of a star in the g-r vs. u-g color-color diagram, thereby demonstrating their robustness as well as a potential for photometric parameter estimation methods. Using Beers et al. parameters, we find that the metallicity distribution of the Milky Way stars at a few kpc from the galactic plane is bimodal with a local minimum at [Z/Zo]~ -1.3. The median metallicity for the low-metallicity [Z/Zo]<-1.3 subsample is nearly independent of Galactic cylindrical coordinates R and z, while it decreases with z for the high-metallicity [Z/Zo]> -1.3 sample. We also find that the low-metallicity sample has ~2.5 times larger velocity dispersion and that it does not rotate (at the ~10 km/s level), while the rotational velocity of the high-metallicity sample decreases smoothly with the height above the galactic plane.

We investigate the dependences of galaxy star-formation history and galaxy morphology on environment, using color and H-alpha equivalent width as star-formation history indicators, using concentration and central surface brightness as morphological indicators, and using clustocentric distance as an environment indicator. Clustocentric distance has the virtue that it can be measured with very high precision over a large dynamic range. We find the following asymmetry between morphological and star-formation history parameters: star-formation history parameters relate directly to the clustocentric distance while morphological parameters relate to the clustocentric distance only indirectly through their relationships with star-formation history. This asymmetry has important implications for the role that environment plays in shaping galaxy properties and it places strong constraints on theoretical models of galaxy formation. Current semi-analytic models do not reproduce this effect.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

We investigate the clustering of galaxy groups and clusters in the SDSS using the Berlind et al. (2006) group sample, which is designed to identify galaxy systems that each occupy a single dark matter halo. We estimate group masses from their abundances, and measure their relative large-scale bias as a function of mass. Our measurements are in agreement with the theoretical halo bias function, given a standard LCDM cosmological model, and they tend to favor a low value of the power spectrum amplitude sigma_8. We search for a residual dependence of clustering on other group properties at fixed mass, and find the strongest signal for central galaxy color in high mass groups. Massive groups with less red central galaxies are more biased on large scales than similar mass groups with redder central galaxies. We show that this effect is unlikely to be caused by errors in our mass estimates, and is most likely observational evidence of recent theoretical findings that halo bias depends on a ``second parameter'' other than mass, such as age or concentration. To compare with the data, we study the bias of massive halos in N-body simulations and quantify the strength of the relation between halo bias and concentration at fixed mass. In addition to confirming a non-trivial prediction of the LCDM cosmological model, these results have important implications for the role that environment plays in shaping galaxy properties.

We determine what aspects of the density field surrounding galaxies most affect their properties. For Sloan Digital Sky Survey galaxies, we measure the group environment, meaning the host group luminosity and the distance from the group center (hereafter, ``groupocentric distance''). For comparison, we measure the surrounding density field on scales ranging from 100 kpc/h to 10 Mpc/h. We use the relationship between color and group environment to test the null hypothesis that only the group environment matters, searching for a residual dependence of properties on the surrounding density. Generally, red galaxies are slightly more clustered on small scales (about 100--300 kpc/h) than the null hypothesis predicts, possibly indicating that substructure within groups has some importance. At large scales (> 1 Mpc/h), the actual projected correlation functions of galaxies are biased at less than the 5% level with respect to the null hypothesis predictions. We exclude strongly the converse null hypothesis, that only the surrounding density (on any scale) matters. These results generally encourage the use of the halo model description of galaxy bias, which models the galaxy distribution as a function of host halo mass alone. We compare these results to proposed galaxy formation scenarios within the Cold Dark Matter cosmological model.

We investigate the gas content and baryonic Tully-Fisher relationship for extremely low luminosity dwarf galaxies in the absolute magnitude range -13.5 > Mr > -16. The sample is selected from the Sloan Digital Sky Survey and consists of 101 galaxies for which we have obtained follow-up HI observations using the Arecibo Observatory and Green Bank Telescope. This represents the largest homogeneous sample of dwarfs at low luminosities with well-measured HI and optical properties. The sample spans a range of environments, from dense groups to truly isolated galaxies. The average neutral gas fraction is f_gas=0.6, significantly exceeding that of typical gas-rich galaxies at higher luminosities. Dwarf galaxies are therefore less efficient at turning gas into stars over their lifetimes. The strong environmental dependence of the gas fraction distribution demonstrates that while internal processes can reduce the gas fractions to roughly f_gas=0.4, external processes are required to fully remove gas from a dwarf galaxy. The average rotational velocity of our sample is vrot=50 km/s. Including more massive galaxies from the literature, we fit a baryonic Tully-Fisher slope of M_baryon \propto vrot^(3.70+/- 0.15). This slope compares well with CDM models that assume an equal baryon to dark matter ratio at all masses. While gas stripping or other processes may modify the baryon to dark matter ratio for dwarfs in the densest environments, the majority of dwarf galaxies in our sample have not preferentially lost significant baryonic mass relative to more massive galaxies.

We identify galaxy groups and clusters in volume-limited samples of the SDSS redshift survey, using a redshift-space friends-of-friends algorithm. We optimize the friends-of-friends linking lengths to recover galaxy systems that occupy the same dark matter halos, using a set of mock catalogs created by populating halos of N-body simulations with galaxies. Extensive tests with these mock catalogs show that no combination of perpendicular and line-of-sight linking lengths is able to yield groups and clusters that simultaneously recover the true halo multiplicity function, projected size distribution, and velocity dispersion. We adopt a linking length combination that yields, for galaxy groups with ten or more members: a group multiplicity function that is unbiased with respect to the true halo multiplicity function; an unbiased median relation between the multiplicities of groups and their associated halos; a spurious group fraction of less than ~1%; a halo completeness of more than ~97%; the correct projected size distribution as a function of multiplicity; and a velocity dispersion distribution that is ~20% too low at all multiplicities. These results hold over a range of mock catalogs that use different input recipes of populating halos with galaxies. We apply our group-finding algorithm to the SDSS data and obtain three group and cluster catalogs for three volume-limited samples that cover 3495.1 square degrees on the sky. We correct for incompleteness caused by fiber collisions and survey edges, and obtain measurements of the group multiplicity function, with errors calculated from realistic mock catalogs. These multiplicity function measurements provide a key constraint on the relation between galaxy populations and dark matter halos.

We present results of a measurement of the shape of the density profile of galaxy groups and clusters traced by 43 335 Luminous Red Galaxies (LRGs) with spectroscopic redshifts from the Sloan Digital Sky Survey (SDSS). The galaxies are selected such that they are the brightest within a cylindrical aperture, split into two luminosity samples, and modeled as the sum of stellar and dark matter components. We present a detailed investigation of many possible systematic effects that could contaminate our signal and develop methods to remove them, including a detected intrinsic alignment for galaxies within 100 kpc/h of LRGs which we remove using photometric redshift information. The resulting lensing signal is consistent with NFW profile dark matter halos; the SIS profile is ruled out at the 96 (conservatively) and 99.96 per cent confidence level (CL) for the fainter and brighter lens samples (respectively) when we fit using lensing data between 40 kpc/h and 2 Mpc/h with total signal-to-noise of 19 and 25 for the two lens samples. The lensing signal amplitude suggests that the faint and bright sample galaxies typically reside in haloes of mass (2.9+/-0.4)x10^{13} M_{sun}/h and (6.7+/-0.8)x10^{13} M_{sun}/h respectively, in good agreement with predictions based on halo spatial density with normalization lower than the 'concordance' sigma_8=0.9. When fitting for the concentration parameter in the NFW profile, we find c = 5.0+/-0.6 (stat) +/-1 (sys), and c = 5.6+/-0.6 (stat) +/-1 (sys) for the faint and bright samples, consistent with Lambda-CDM simulations. We also split the bright sample further to determine masses and concentrations for cluster-mass halos, finding mass (1.3+/-0.2)x10^{14} M_{sun}/h for the sample of LRGs brighter than -22.6 in r. (Abridged)

(abridged) There are good observational reasons to believe that the progenitors of red galaxies have undergone starbursts, followed by a post-starburst phase. We investigate the environments of post-starburst galaxies by measuring \textsl{(1)} number densities in $8 h^{-1} \mathrm{Mpc}$ radius comoving spheres, \textsl{(2)} transverse distances to nearest Virgo-like galaxy clusters, and \textsl{(3)} transverse distances to nearest luminous-galaxy neighbors. We compare the post-starburst galaxies to currently star-forming galaxies identified solely by A-star excess or $\Halpha$ emission. We find that post-starburst galaxies are in the same kinds of environments as star-forming galaxies; this is our ``null hypothesis''. More importantly, we find that at each value of the A-star excess, the star-forming and post-starburst galaxies lie in very similar distributions of environment. The only deviations from our null hypothesis are barely significant: a slight deficit of post-starburst galaxies (relative to the star-forming population) in very low-density regions, a small excess inside the virial radii of clusters, and a slight excess with nearby neighbors. None of these effects is strong enough to make the post-starburst galaxies a high-density phenomenon, or to argue that the starburst events are primarily triggered by external tidal impulses (e.g., from close passages of massive galaxies). The small excess inside cluster virial radii suggests that some post-starbursts are triggered by interactions with the intracluster medium, but this represents a very small fraction of all post-starburst galaxies.

We discuss the panchromatic properties of 99,088 galaxies selected from the SDSS Data Release 1 spectroscopic sample (a flux-limited sample for 1360 deg^2). These galaxies are positionally matched to sources detected by ROSAT, GALEX, 2MASS, IRAS, GB6, FIRST, NVSS and WENSS. We find strong correlations between the detection fraction at other wavelengths and optical properties such as flux, colors, and emission-line strengths. Using GALEX, SDSS, and 2MASS data, we construct the UV-IR broad-band spectral energy distributions for various types of galaxies, and find that they form a nearly one-parameter family. For example, based on SDSS u- and r-band data, supplemented with redshift, the K-band 2MASS magnitudes can be "predicted" with an rms scatter of only 0.2 mag. When a dust content estimate determined from SDSS data by Kauffmann et al. (2003) is also utilized, this scatter decreases to 0.1 mag. We demonstrate that this dust content is indeed higher for galaxies detected by IRAS and that it can be used to "predict" measured IRAS 60 micron flux density within a factor of two using only SDSS data. We also show that the position of a galaxy in the emission-line-based Baldwin-Phillips-Terlevich diagram is correlated with the optical light concentration index and u-r color determined from the SDSS broad-band imaging data, and discuss changes in the morphology of this diagram induced by requiring detections at other wavelengths. We study the IR-radio correlation and find evidence that its slope may be different for AGN and star-forming galaxies and related to the H_alpha/H_beta line strength ratio.

Template fits to observed galaxy fluxes allow calculation of K-corrections and conversions among observations of galaxies at various wavelengths. We present a method for creating model-based template sets given a set of heterogeneous photometric and spectroscopic galaxy data. Our technique, non-negative matrix factorization, is akin to principle component analysis (PCA), except that it is constrained to produce nonnegative templates, it can use a basis set of models (rather than the delta function basis of PCA), and it naturally handles uncertainties, missing data, and heterogeneous data (including broad-band fluxes at various redshifts). The particular implementation we present here is suitable for ultraviolet, optical, and near-infrared observations in the redshift range 0 < z < 1.5. Since we base our templates on stellar population synthesis models, the results are intepretable in terms of approximate stellar masses and star-formation histories. We present templates fit with this method to data from GALEX, Sloan Digital Sky Survey spectroscopy and photometry, the Two-Micron All Sky Survey, the Deep Extragalactic Evolutionary Probe and the Great Observatories Origins Deep Survey. In addition, we present software for using such data to estimate K-corrections and stellar masses.

We present the 3D real space clustering power spectrum of a sample of \~600,000 luminous red galaxies (LRGs) measured by the Sloan Digital Sky Survey (SDSS), using photometric redshifts. This sample of galaxies ranges from redshift z=0.2 to 0.6 over 3,528 deg^2 of the sky, probing a volume of 1.5 (Gpc/h)^3, making it the largest volume ever used for galaxy clustering measurements. We measure the angular clustering power spectrum in eight redshift slices and combine these into a high precision 3D real space power spectrum from k=0.005 (h/Mpc) to k=1 (h/Mpc). We detect power on gigaparsec scales, beyond the turnover in the matter power spectrum, on scales significantly larger than those accessible to current spectroscopic redshift surveys. We also find evidence for baryonic oscillations, both in the power spectrum, as well as in fits to the baryon density, at a 2.5 sigma confidence level. The statistical power of these data to constrain cosmology is ~1.7 times better than previous clustering analyses. Varying the matter density and baryon fraction, we find \Omega_M = 0.30 \pm 0.03, and \Omega_b/\Omega_M = 0.18 \pm 0.04, The detection of baryonic oscillations also allows us to measure the comoving distance to z=0.5; we find a best fit distance of 1.73 \pm 0.12 Gpc, corresponding to a 6.5% error on the distance. These results demonstrate the ability to make precise clustering measurements with photometric surveys (abridged).

We use the Sloan Digital Sky Survey (SDSS) spectroscopic sample to constrain the projected radial distribution of satellites around isolated ~ L* galaxies. We employ mock galaxy catalogs derived from high-resolution cosmological simulations to investigate the effects of interloper contamination and show that interlopers significantly bias the estimated slope of the projected radial distribution of satellites. We also show that the distribution of interlopers around galaxies is expected to be non-uniform in velocity space because galaxies are clustered and reside in crowded environments. Successful methods of interloper contamination correction should therefore take into account environments of the host galaxies. Two such new methods are presented and the most reliable of them is used to correct for interloper contamination in analyses of the SDSS galaxy sample. The best fit power-law slope of the interloper-corrected surface density distribution of satellites, Sigma(R) ~ R^alpha, in the volume-limited SDSS sample is alpha ~ -1.7 +/- 0.1, independent of the galaxy and satellite luminosities. Comparison with Lambda-CDM simulations shows that the radial distribution of the SDSS satellites is more concentrated than that of subhalos around galaxy-sized halos, especially at R < 100 kpc/h. The predicted dark matter radial distribution is somewhat more concentrated than the profile of the SDSS satellites, but the difference is not statistically significant for our sample.

SDSSJ1049+5103, commonly known as Willman 1, is an extremely low-luminosity Milky Way companion whose properties are intermediate between those of globular clusters and dwarf spheroidals. In this paper, we use deep photometry to show that this object is old and moderately metal-poor, has a distance of 38 +/- 7 kpc, has an M_V of -2.5 mag, and has a half-light radius of 21 +/- 7 pc, consistent with previous estimates. The spatial distribution of Willman 1's main sequence stars shows 1) its total spatial extent exceeds its tidal radius for a range of assumptions about its total mass and its orbit and 2) the presence of prominent multi-directional stellar tails. The tidal interactions causing these tail features may explain the large physical size of Willman 1 relative to low-luminosity globular clusters. It is the most distant Galactic object yet known to display prominent tails, and is the only distant satellite to display multi-directional tails. Although we cannot at present determine the cause of this unusual morphology, preliminary comparisons between the morphology of Willman 1 and published simulations suggest that it may be near the apocenter of its orbit, or that it may have interacted with another halo object. We find a significant difference between the luminosity functions of stars in the center and in the tails of Willman 1, strongly suggesting mass segregation much like that seen in Palomar 5. Although Willman 1 has more pronounced tidal tails than most confirmed Milky Way dwarf galaxies, because of its very low stellar mass we cannot at present rule out the possibility that it has a dark matter halo. (Abridged)

We investigate the extent to which galaxies' star-formation histories and morphologies are determined by their clustocentric distance and their nearest cluster's richness. We look at three tracers of star-formation history (i-band absolute magnitude M_i, [g-r] color, and Halpha emission line equivalent width) and two indicators of galaxy morphology (surface brightness and radial concentration) for 52,569 galaxies in the redshift range of 0.015<z<0.068. We find that our morphology indicators relate to the clustocentric distance only indirectly through their relationships with stellar population and starformation rate. Galaxies that are near the cluster center tend to be more luminous, redder and have lower Halpha EW (ie, lower star-formation rates) than those that lie near or outside the virial radius of the cluster. The detailed relationships between these galaxy properties and clustocentric distance depend on cluster richness. For richer clusters, we find that (i) the transition in color and Halpha EW from cluster center to field values is more abrupt and occurs closer to the cluster virial radius, and (ii) the color and Halpha EW distributions are overall narrower than in less rich clusters. We also find that the radial gradient seen in the luminosity distribution is strongest around the smaller clusters and decreases as cluster richness increases. We find there is a `characteristic distance' at around one virial radius (the infall region) where the change with radius of galaxy property distributions is most dramatic, but we find no evidence for infall-triggered star bursts. These results suggest that galaxies `know' the distance to, and the size of, their nearest cluster and they express this information in their star-formation histories.

We present measurements of the normalised redshift-space three-point correlation function (Q_z) of galaxies from the Sloan Digital Sky Survey (SDSS) main galaxy sample. We have applied our "npt" algorithm to both a volume-limited (36738 galaxies) and magnitude-limited sample (134741 galaxies) of SDSS galaxies, and find consistent results between the two samples, thus confirming the weak luminosity dependence of Q_z recently seen by other authors. We compare our results to other Q_z measurements in the literature and find it to be consistent within the full jack-knife error estimates. However, we find these errors are significantly increased by the presence of the ``Sloan Great Wall'' (at z ~ 0.08) within these two SDSS datasets, which changes the 3-point correlation function (3PCF) by 70% on large scales (s>=10h^-1 Mpc). If we exclude this supercluster, our observed Q_z is in better agreement with that obtained from the 2dFGRS by other authors, thus demonstrating the sensitivity of these higher-order correlation functions to large-scale structures in the Universe. This analysis highlights that the SDSS datasets used here are not ``fair samples'' of the Universe for the estimation of higher-order clustering statistics and larger volumes are required. We study the shape-dependence of Q_z(s,q,theta) as one expects this measurement to depend on scale if the large scale structure in the Universe has grown via gravitational instability from Gaussian initial conditions. On small scales (s <= 6h^-1 Mpc), we see some evidence for shape-dependence in Q_z, but at present our measurements are consistent with a constant within the errors (Q_z ~ 0.75 +/- 0.05). On scales >10h^-1 Mpc, we see considerable shape-dependence in Q_z.

In this paper, we present Sloan Digital Sky Survey (SDSS) photometry and spectroscopy in the fields of 24 gamma-ray bursts (GRBs) observed by Swift, including bursts localized by Swift, HETE-2, and INTEGRAL, after December 2004. After this bulk release, we plan to provide individual releases of similar data shortly after the localization of future bursts falling in the SDSS survey area. These data provide a solid basis for the astrometric and photometric calibration of follow-up afterglow searches and monitoring. Furthermore, the images provided with this release will allow observers to find transient objects up to a magnitude fainter than possible with Digitized Sky Survey image comparisons.

We present the small-scale (0.01<r<8 h^{-1} Mpc) projected correlation function w_p(r_p) and real space correlation function xi(r) of 24520 luminous early-type galaxies from the Sloan Digital Sky Survey Luminous Red Galaxy (LRG) sample (0.16<z<0.36). ``Fiber collision'' incompleteness of the SDSS spectroscopic sample at scales smaller than 55 arcseconds prevents measurements of the correlation function for LRGs on scales smaller than ~0.3 Mpc by the usual methods. In this work, we cross-correlate the spectroscopic sample with the imaging sample, with a weighting scheme to account for the collisions, extensively tested against mock catalogs. We correct for photometric biases in the SDSS imaging of close galaxy pairs. We find that the correlation function xi(r) is surprisingly close to a r^{-2} power law over more than 4 orders of magnitude in separation r. This result is too steep at small scales to be explained in current versions of the halo model for galaxy clustering. We infer an LRG-LRG merger rate of less than or similar to 0.6 Gyr^{-1} Gpc^{-3} for this sample. This result suggests that the LRG-LRG mergers are not the main mode of mass growth for LRGs at z<0.36.

In the six billion years between redshifts z=1 and z=0.1, galaxies change due to the aging of their stellar populations, the formation of new stars, and mergers with other galaxies. Here I explore the relative importance of these various effects, finding that while mergers are likely to be important for the red galaxy sequence they are unlikely to affect more than 10% of the blue galaxy sequence. I compare the galaxy population at redshift z=0.1 from the Sloan Digital Sky Survey to that at z=1 from the Deep Extragalactic Evolutionary Probe 2. Galaxies are bluer at z=1: the blue sequence by about 0.3 mag and the red sequence by about 0.1 mag, in redshift z=0.1 (u-g) color. I evaluate the change in color and in the luminosity functions of the two sequences using some simplistic stellar population synthesis models. These models indicate that the luminous end of the red sequence fades less than passive evolution allows by about 0.2 mag. Due to a lack of luminous blue progenitors, ``dry'' mergers betweeen red galaxies then must create the luminous red population at z=0.1, if stellar population models are correct. The blue sequence colors and luminosity function are consistent with a reduction in the star-formation rate since redshift z=1 by a factor of about three, with no change in the number density to within 10%. These results restrict the number of blue galaxies that can fall onto the red sequence by any process, and in particular suggest that if mergers are catastrophic events they must be rare for blue galaxies.

We study the luminosity and color dependence of the galaxy 2-point correlation function in the Sloan Digital Sky Survey, starting from a sample of 200,000 galaxies over 2500 deg^2. We concentrate on the projected correlation function w(r_p), which is directly related to the real space \xi(r). The amplitude of w(r_p) grows continuously with luminosity, rising more steeply above the characteristic luminosity L_*. Redder galaxies exhibit a higher amplitude and steeper correlation function at all luminosities. The correlation amplitude of blue galaxies increases continuously with luminosity, but the luminosity dependence for red galaxies is less regular, with bright red galaxies more strongly clustered at large scales and faint red galaxies more strongly clustered at small scales. We interpret these results using halo occupation distribution (HOD) models assuming concordance cosmological parameters. For most samples, an HOD model with two adjustable parameters fits the w(r_p) data better than a power-law, explaining inflections at r_p ~ 1-3 Mpc/h as the transition between the 1-halo and 2-halo regimes of \xi(r). The implied minimum mass for a halo hosting a central galaxy above a luminosity threshold L grows as M_min ~ L at low luminosities and more steeply above L_*. The mass at which an average halo has one satellite galaxy brighter than L is M_1 ~ 23 M_min(L). These results imply a conditional luminosity function (at fixed halo mass) in which central galaxies lie far above a Schechter function extrapolation of the satellite population. HOD models nicely explain the joint luminosity-color dependence of w(r_p) in terms of the color fractions of central and satellite populations as a function of halo mass. The inferred HOD properties are in good qualitative agreement with theoretical predictions.

We measure the projected correlation function w_p(r_p) from the Sloan Digital Sky Survey for a flux-limited sample of 118,000 galaxies and for a volume limited subset of 22,000 galaxies with absolute magnitude M_r<-21. Both correlation functions show subtle but systematic departures from the best-fit power law, in particular a change in slope at r_p~1-2 Mpc/h. These departures are stronger for the volume-limited sample, which is restricted to relatively luminous galaxies. We show that the inflection point in w_p(r_p) can be naturally explained by contemporary models of galaxy clustering, according to which it marks the transition from a large scale regime dominated by galaxy pairs in separate dark matter halos to a small scale regime dominated by galaxy pairs in the same dark matter halo. For example, given the dark halo population predicted by an inflationary cold dark matter scenario, the projected correlation function of the volume-limited sample can be well reproduced by a model in which the mean number of M_r<-21 galaxies in a halo of mass M>M_1=4.74 X 10^{13}\msun/h is <N>_M=(M/M_1)^{0.89}, with 75% of the galaxies residing in less massive, single-galaxy halos. This physically motivated model has the same number of free parameters as a power law, and it fits the w_p(r_p) data better, with a \chi^2/d.o.f.=0.93 compared to 6.12 (for 10 degrees of freedom, incorporating the covariance of the correlation function errors). Departures from a power-law correlation function encode information about the relation between galaxies and dark matter halos. Higher precision measurements of these departures for multiple classes of galaxies will constrain galaxy bias and provide new tests of the theory of galaxy formation.

Weak lensing is emerging as a powerful observational tool to constrain cosmological models, but is at present limited by an incomplete understanding of many sources of systematic error. Many of these errors are multiplicative and depend on the population of background galaxies. We show how the commonly cited geometric test, which is rather insensitive to cosmology, can be used as a ratio test of systematics in the lensing signal at the 1 per cent level. We apply this test to the galaxy-galaxy lensing analysis of the Sloan Digital Sky Survey (SDSS), which at present is the sample with the highest weak lensing signal to noise and has the additional advantage of spectroscopic redshifts for lenses. This allows one to perform meaningful geometric tests of systematics for different subsamples of galaxies at different mean redshifts, such as brighter galaxies, fainter galaxies and high-redshift luminous red galaxies, both with and without photometric redshift estimates. We use overlapping objects between SDSS and the DEEP2 and 2SLAQ spectroscopic surveys to establish accurate calibration of photometric redshifts and to determine the redshift distributions for SDSS. We use these redshift results to compute the projected surface density contrast DeltaSigma around 259 609 spectroscopic galaxies in the SDSS; by measuring DeltaSigma with different source samples we establish consistency of the results at the 10 per cent level (1-sigma). We also use the ratio test to constrain shear calibration biases and other systematics in the SDSS survey data to determine the overall galaxy-galaxy weak lensing signal calibration uncertainty. We find no evidence of any inconsistency among many subsamples of the data.

In this Letter, we report the discovery of a new dwarf satellite to the Milky Way, located at ($\alpha_{2000}, \delta_{2000}$) $=$ (158.72,51.92) in the constellation of Ursa Major. This object was detected as an overdensity of red, resolved stars in Sloan Digital Sky Survey data. The color-magnitude diagram of the Ursa Major dwarf looks remarkably similar to that of Sextans, the lowest surface brightness Milky Way companion known, but with approximately an order of magnitude fewer stars. Deeper follow-up imaging confirms this object has an old and metal-poor stellar population and is $\sim$ 100 kpc away. We roughly estimate M$_V =$ -6.75 and $r_{1/2} =$ 250 pc for this dwarf. Its luminosity is several times fainter than the faintest known Milky Way dwarf. However, its physical size is typical for dSphs. Even though its absolute magnitude and size are presently quite uncertain, Ursa Major is likely the lowest luminosity and lowest surface brightness galaxy yet known.

[abridged] A "missing data" algorithm is developed to model (ie, reconstruct) the three-dimensional velocity distribution function of a sample of stars using data (velocity measurements) every one of which has one dimension unmeasured (the radial direction). It also accounts for covariant measurement uncertainties on the tangential velocity components. The algorithm is applied to tangential velocities measured in a kinematically unbiased sample of 11,865 stars taken from the Hipparcos catalog. The local stellar velocity distribution function of each of a set of 20 color-selected subsamples is modeled as a mixture of two three-dimensional Gaussian ellipsoids of arbitrary relative responsibility. In the fitting, one Gaussian (the "halo") is fixed at the known mean velocity and velocity variance tensor of the Galaxy halo, and the other (the "disk") is allowed to take arbitrary mean and arbitrary variance tensor. The mean and variance tensor (commonly the "velocity ellipsoid") of the disk velocity distribution are both found to be strong functions of stellar color, with long-lived populations showing larger velocity dispersion, slower mean rotation velocity, and smaller vertex deviation than short-lived populations. The local standard of rest (LSR) is inferred in the usual way and the Sun's motion relative to the LSR is found to be (U,V,W)_{\odot}=(10.1,4.0,6.7)+/-(0.5,0.8,0.2) km/s. Artificial data sets are made and analyzed, with the same error properties as the Hipparcos data, to demonstrate that the analysis is unbiased. The results are shown to be insensitive to the assumption that the velocity distributions are Gaussian.

We examine a sample of low redshift (10 < d < 150 Mpc) galaxies including galaxies with r-band absolute magnitudes as faint as -12.5 (for h=1), selected from the Sloan Digital Sky Survey Data Release 2 (SDSS). The sample is unique in containing galaxies of extremely low luminosities in a wide range of environments, selected with uniform and well-understood criteria. We present the luminosity function as well as the broad-band properties of low luminosity galaxies in this sample. A Schechter function is an insufficient parameterization of the r-band luminosity function; there is an upturn in the slope at low luminosities. The resulting slope at low luminosities in this sample is around -1.3. However, we almost certainly miss a large number of galaxies at very low luminosities due to low surface brightness selection effects, and we estimate that the true low luminosity slope may be as steep or steeper than -1.5. The results here are consistent with previous SDSS results and, in the g-band, roughly consistent with the results of the Two degree Field Galaxy Redshift Survey. Extremely low luminosity galaxies are predominantly blue, low surface brightness, exponential disks.

We present the third edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog. The catalog consists of the 46,420 objects in the SDSS Third Data Release that have luminosities larger than M_i = -22 (in a cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega_Lambda = 0.7), have at least one emission line with FWHM larger than 1000 km/s or are unambiguously broad absorption line quasars, are fainter than i = 15.0, and have highly reliable redshifts. The area covered by the catalog is 4188 sq. deg. The quasar redshifts range from 0.08 to 5.41, with a median value of 1.47; the high-redshift sample includes 520 quasars at redshifts greater than four, of which 17 are at redshifts greater than five. For each object the catalog presents positions accurate to better than 0.2 arcsec. rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800--9200A at a spectral resolution about 2000; the spectra can be retrieved from the public database using the information provided in the catalog. A total of 44,221 objects in the catalog were discovered by the SDSS; 28,400 of the SDSS discoveries are reported here for the first time.

We report the discovery of SDSSJ1049+5103, an overdensity of resolved blue stars at (\alpha_{2000}, \delta_{2000}) = (162.343, 51.051). This object appears to be an old, metal-poor stellar system at a distance of 45 +/- 10 kpc, with a half-light radius of 23$\pm 10$ pc and an absolute magnitude of M_V = -3.0^{+2.0}_{-0.7}. One star that is likely associated with this companion has an SDSS spectrum confirming it as a blue horizontal branch star at 48 kpc. The color-magnitude diagram of SDSSJ1049+5103 contains few, if any, horizontal or red giant branch stars, similar to the anomalously faint globular cluster AM 4. The size and luminosity of SDSSJ1049+5103 places it at the intersection of the size-luminosity relationships followed by known globular clusters and by Milky Way dwarf spheroidals. If SDSSJ1049+5103 is a globular cluster, then its properties are consistent with the established trend that the largest radius Galactic globular clusters are all in the outer halo. However, the five known globular clusters with similarly faint absolute magnitudes all have half-mass radii that are smaller than SDSSJ1049+5103 by a factor of $\gtrsim$ 5. If it is a dwarf spheroidal, then it is the faintest yet known by two orders of magnitude, and is the first example of the ultra-faint dwarfs predicted by some theories. The uncertain nature of this new system underscores the sometimes ambiguous distinction between globular clusters and dwarf spheroidals. A simple friends-of-friends search for similar blue, small scalesize star clusters detected all known globulars and dwarfs closer than 50 kpc in the SDSS area, but yielded no other candidates as robust as SDSSJ1049+5103.

Here we present the New York University Value-Added Galaxy Catalog (NYU-VAGC), a catalog of local galaxies (mostly below a redshift of about 0.3) based on a set of publicly-released surveys (including the 2dFGRS, 2MASS, PSCz, FIRST, and RC3) matched to the Sloan Digital Sky Survey (SDSS) Data Release 2. Excluding areas masked by bright stars, the photometric sample covers 3514 square degrees and the spectroscopic sample covers 2627 square degrees (with about 85% completeness). Earlier, proprietary versions of this catalog have formed the basis of many SDSS investigations of the power spectrum, correlation function, and luminosity function of galaxies. We calculate and compile derived quantities (for example, K-corrections and structural parameters for galaxies). The SDSS catalog presented here is photometrically recalibrated, reducing systematic calibration errors across the sky from about 2% to about 1%. We include an explicit description of the geometry of the catalog, including all imaging and targeting information as a function of sky position. Finally, we have performed eyeball quality checks on a large number of objects in the catalog in order to flag deblending and other errors. This catalog is complementary to the SDSS Archive Servers, in that NYU-VAGC's calibration, geometrical description, and conveniently small size are specifically designed for studying galaxy properties and large-scale structure statistics using the SDSS spectroscopic catalog.

We combine the measurements of luminosity dependence of bias with the luminosity dependent weak lensing analysis of dark matter around galaxies to derive the galaxy bias and constrain nonlinear mass and cosmological parameters. We take advantage of theoretical and simulation predictions that predict that while halo bias is rapidly increasing with mass for high masses, it is nearly constant in low mass halos. We use a new weak lensing analysis around the same SDSS galaxies to determine their halo mass probability distribution. These halo mass probability distributions are used to predict the bias for each luminosity subsample and we find an excellent agreement with observed values. The required nonlinear mass suggests slightly lower matter density than usually assumed, Omegam=0.25+/- 0.03 for the simplest models. We combine the bias constraints with those from the WMAP and the SDSS power spectrum analysis to derive new constraints on bias and sigma_8. For the most general parameter space we find sigma_8=0.88+/- 0.06 and b_*=0.99+/- 0.07. In the context of spatially flat models we improve the limit on the neutrino mass for the case of 3 degenerate families from m_nu<0.6eV without bias to m_nu<0.18eV with bias (95% c.l.), which is weakened to m_nu<0.24eV if running is allowed. The corresponding limit for 3 massless + 1 massive neutrino is 1.37eV.

We use the projected correlation function w_p(r_p) of a volume-limited subsample of the Sloan Digital Sky Survey (SDSS) main galaxy redshift catalogue to measure the halo occupation distribution (HOD) of the galaxies of the sample. Simultaneously, we allow the cosmology to vary within cosmological constraints imposed by cosmic microwave background experiments in a Lambda-CDM model. We find that combining w_p(r_p) for this sample alone with the observations by WMAP, ACBAR, CBI and VSA can provide one of the most precise techniques available to measure cosmological parameters. For a minimal flat six-parameter Lambda-CDM model with an HOD with three free parameters, we find Omega_m=0.278^{+0.027}_{-0.027}, sigma_8=0.812^{+0.028}_{-0.027}, and H_0=69.8^{+2.6}_{-2.6}km s^{-1} Mpc^{-1}; these errors are significantly smaller than from CMB alone and similar to those obtained by combining CMB with the large-scale galaxy power spectrum assuming scale-independent bias. The corresponding HOD parameters describing the minimum halo mass and the normalization and cut-off of the satellite mean occupation are M_min=(3.03^{+0.36}_{-0.36})x 10^{12} h^{-1} M_sun, M_1 = (4.58^{+0.60}_{-0.60})x 10^{13} h^{-1} M_sun, and kappa=4.44^{+0.51}_{-0.69}. When more parameters are added to the HOD model, the error bars on the HOD parameters increase because of degeneracies, but the error bars on the cosmological parameters do not increase greatly. Similar modeling for other galaxy samples could reduce the statistical errors on these results, while more thorough investigations of the cosmology dependence of nonlinear halo bias and halo mass functions are needed to eliminate remaining systematic uncertainties, which may be comparable to statistical uncertainties.

We present the large-scale correlation function measured from a spectroscopic sample of 46,748 luminous red galaxies from the Sloan Digital Sky Survey. The survey region covers 0.72 h^{-3} Gpc^3 over 3816 square degrees and 0.16<z<0.47, making it the best sample yet for the study of large-scale structure. We find a well-detected peak in the correlation function at 100h^{-1} Mpc separation that is an excellent match to the predicted shape and location of the imprint of the recombination-epoch acoustic oscillations on the low-redshift clustering of matter. This detection demonstrates the linear growth of structure by gravitational instability between z=1000 and the present and confirms a firm prediction of the standard cosmological theory. The acoustic peak provides a standard ruler by which we can measure the ratio of the distances to z=0.35 and z=1089 to 4% fractional accuracy and the absolute distance to z=0.35 to 5% accuracy. From the overall shape of the correlation function, we measure the matter density Omega_mh^2 to 8% and find agreement with the value from cosmic microwave background (CMB) anisotropies. Independent of the constraints provided by the CMB acoustic scale, we find Omega_m = 0.273 +- 0.025 + 0.123 (1+w_0) + 0.137 Omega_K. Including the CMB acoustic scale, we find that the spatial curvature is Omega_K=-0.010+-0.009 if the dark energy is a cosmological constant. More generally, our results provide a measurement of cosmological distance, and hence an argument for dark energy, based on a geometric method with the same simple physics as the microwave background anisotropies. The standard cosmological model convincingly passes these new and robust tests of its fundamental properties.

We construct and analyze a u-band selected galaxy sample from the SDSS Southern Survey, which covers 275 sq. deg. The sample includes 43223 galaxies with spectroscopic redshifts in the range 0.005<z<0.3 and with 14.5<u<20.5. The S/N in the u-band Petrosian aperture is improved by coadding multiple epochs of imaging data and by including sky-subtraction corrections. Luminosity functions for the near-UV ^{0.1}u band (322+-26 nm) are determined in redshift slices of width 0.02, which show a highly significant evolution in M^{star} of -0.8+-0.1 mag between z=0 and 0.3; with M^{star} = -18.84+-0.05 (AB mag), log phi^{star} = -2.06+-0.03 (Mpc^{-3}) and log rho_L = 19.11+-0.02 (W Hz^{-1} Mpc^{-3}) at z=0.1. The faint-end slope determined for z<0.06 is given by alpha=-1.05+-0.08. This is in agreement with recent determinations from GALEX at shorter wavelengths. Comparing our z<0.3 luminosity density measurements with 0.2<z<1.2 from COMBO-17, we find that the 280-nm density evolves as rho_L proportional to (1+z)^{beta} with beta=2.1+-0.2; and find no evidence for any change in slope over this redshift range. By comparing with other measurements of cosmic star formation history, we estimate that the effective dust attenuation at 280 nm has increased by 0.8+-0.3 mag between z=0 and 1.

We have compiled a large sample of low-redshift active galactic nuclei (AGN) identified via their emission line characteristics from the spectroscopic data of the Sloan Digital Sky Survey. Since emission lines are often contaminated by stellar absorption lines, we developed an objective and efficient method of subtracting the stellar continuum from every galaxy spectrum before making emission line measurements. The distribution of the measured H$\alpha$ Full Width at Half Maxima values of emission line galaxies is strongly bimodal, with two populations separated at about 1,200km s$^{-1}$. This feature provides a natural separation between narrow-line and broad-line AGN. The narrow-line AGN are identified using standard emission line ratio diagnostic diagrams. 1,317 broad-line and 3,074 narrow-line AGN are identified from about 100,000 galaxy spectra selected over 1151 square degrees. This sample is used in a companion paper to determine the emission-line luminosity function of AGN.

The emission line luminosity function of active galactic nuclei (AGN) is measured from about 3000 AGN included in the main galaxy sample of the Sloan Digital Sky Survey within a redshift range of $0<z<0.15$. The $\Ha$ and [OIII]$\lambda 5007$ luminosity functions for Seyferts cover luminosity range of $10^{5-9}$$L_\odot$ in H$\alpha$ and the shapes are well fit by broken power laws, without a turnover at fainter nuclear luminosities. Assuming a universal conversion from emission line strength to continuum luminosity, the inferred B band magnitude luminosity function is comparable both to the AGN luminosity function of previous studies and to the low redshift quasar luminosity function derived from the 2dF redshift survey. The inferred AGN number density is approximately 1/5 of all galaxies and about $6\times 10^{-3}$ of the total light of galaxies in the $r$-band comes from the nuclear activity. The numbers of Seyfert 1s and Seyfert 2s are comparable at low luminosity, while at high luminosity, Seyfert 1s outnumber Seyfert 2s by a factor of 2-4. In making the luminosity function measurements, we assumed that the nuclear luminosity is independent of the host galaxy luminosity, an assumption we test {\it a posteriori}, and show to be consistent with the data. Given the relationship between black hole mass and host galaxy bulge luminosity, the lack of correlation between nuclear and host luminosity suggests that the main variable that determines the AGN luminosity is the Eddington ratio, not the black hole mass. This appears to be different from luminous quasars, which are most likely to be shining near the Eddington limit.

We report the intermediate-scale (0.3 to 40 Mpc/h) clustering of 35,000 luminous early-type galaxies at redshifts 0.16 to 0.44 from the Sloan Digital Sky Survey. We present the redshift-space two-point correlation function \xi(s), the projected correlation function w_p(r_p), and the deprojected real-space correlation function \xi(r), for approximately volume-limited samples. As expected, the galaxies are highly clustered, with the correlation length varying from 9.8 +/- 0.2 Mpc/h to 11.2 +/- 0.2 Mpc/h, dependent on the specific luminosity range. For the -23.2 < Mg < -21.2 sample, the inferred bias relative to that of L* galaxies is 1.84 +/- 0.11 for 1 Mpc/h < r_p < 10 Mpc/h, with yet stronger clustering on smaller scales. We detect luminosity-dependent bias within the sample but see no evidence for redshift evolution between z=0.2 and z=0.4. We find a clear indication for deviations from a power-law in the real-space correlation function, with a dip at ~ 2 Mpc/h scales and an upturn on smaller scales. The precision measurements of these clustering trends offer new avenues for the study of the formation and evolution of these massive galaxies.

We present the small-scale (0.2h^{-1} to 7h^{-1} Mpc) cross-correlations between 32,000 luminous early-type galaxies and a reference sample of 16 million normal galaxies from the Sloan Digital Sky Survey. Our method allows us to construct the spherically averaged, real-space cross-correlation function between the spectroscopic LRG sample and galaxies from the SDSS imaging. We report the cross-correlation as a function of scale, luminosity, and redshift. We find very strong luminosity dependences in the clustering amplitudes, up to a factor of 4 over a factor of 4 in luminosity, and measure this dependence with high signal-to-noise ratio. The luminosity dependence of bias is found to depend on scale, with more variation on smaller scales. The clustering as a function of scale is not a power law, but instead has a dip at 1h^{-1} Mpc and an excess on small scales. The fraction of red galaxies within the L* sample surrounding LRGs is a strong function of scale, as expected. However, the fraction of red galaxies evolves in redshift similarly on small and large scales, suggesting that cluster and field populations are changing in the same manner. The results highlight the advantage on small scales of using cross-correlation methods as a means of avoiding shot noise in samples of rare galaxies.

We estimate the Mass Function of void galaxies in the second public data release of the Sloan Digital Sky Survey from a sample of 1000 galaxies with local density contrasts of delta_v < -0.6. The galaxy sample is split into ellipticals and spirals using a color-Sersic index criteria. We estimate the virial masses of ellipticals using the measured spectral line-widths along with the observed size. Projection effects and uncertainties in halo properties make mass estimates of spirals more difficult. We use an inversion of the Tully-Fisher relation to estimate the isothermal rotational velocity, and introduce a scaling factor to estimate the halo extent. We then fit the measured mass function against a theoretical Press-Schechter model, and find that the distribution of galaxies in voids appears to be nearly unbiased compared to the mass.

We test the homogeneity of the Universe at $z\sim 0.3$ with the Luminous Red Galaxy (LRG) spectroscopic sample of the Sloan Digital Sky Survey. First, the mean number $N(R)$ of LRGs within completely surveyed LRG-centered spheres of comoving radius $R$ is shown to be proportional to $R^3$ at radii greater than $R\sim 70 h^{-1} \mathrm{Mpc}$. The test has the virtue that it does not rely on the assumption that the LRG sample has a finite mean density; its results show, however, that there \emph{is} such a mean density. Secondly, the survey sky area is divided into 10 disjoint solid angular regions and the fractional rms density variations of the LRG sample in the redshift range $0.2<z<0.35$ among these ($\sim 2\times10^7 h^{-3} \mathrm{Mpc^3}$) regions is found to be 7 percent of the mean density. This variance is consistent with typical biased \lcdm models and puts very strong constraints on the quality of SDSS photometric calibration.

We investigate the relationship between the colors, luminosities, and environments of galaxies in the Sloan Digital Sky Survey spectroscopic sample, using environmental measurements on scales ranging from 0.2 to 6 Mpc/h. We find: (1) that the relationship between color and environment persists even to the lowest luminosities we probe (absolute magnitude in the r band of about -14 for h=1); (2) at luminosities and colors for which the galaxy correlation function has a large amplitude, it also has a steep slope; and (3) in regions of a given overdensity on small scales (1 Mpc/h), the overdensity on large scales (6 Mpc/h) does not appear to relate to the recent star formation history of the galaxies. Of these results, the last has the most immediate application to galaxy formation theory. In particular, it lends support to the notion that a galaxy's properties are related only to the mass of its host dark matter halo, and not to the larger scale environment.

We present H-alpha rotation curves of two low luminosity field galaxies with r-band absolute magnitudes M_r=-13.9^{+0.8}_{-0.5} and M_r=-14.7^{+0.3}_{-0.2} (for h=0.7; the large error bars reflect distance uncertainties). Most previously studied galaxies in this luminosity range are members of groups defined by brighter galaxies, but these two systems, selected from Blanton et al.'s (2004) sample of low luminosity galaxies in the Sloan Digital Sky Survey (SDSS), appear to have no bright companions. The measured rotation speeds at the outer extent of the H-alpha rotation curves, 34.8+/-3.8 km/s and 30.9+/-7.2 km/s, are larger than the values of 16.0^{+6.1}_{-5.4} km/s and 20.9^{+6.2}_{-5.2} km/s predicted by extrapolating the inverse Tully-Fisher relation of luminous SDSS galaxies to these faint luminosities. However, a previous HI measurement of the first galaxy shows that it has a gas mass similar to its stellar mass, and the total baryonic mass is consistent with that predicted by McGaugh et al.'s (2000) "baryonic Tully-Fisher relation." We find r-band dynamical mass-to-light ratios within the radii of the last H-alpha data points (about 1.8 disk scale lengths in each case) of 12.6^{+4.7}_{-4.5} Msun/Lsun and 4.8^{+2.5}_{-2.1} Msun/Lsun, much higher than the values ~ 1 Msun/Lsun expected for the stellar populations. The dynamical properties of these galaxies, including the rotation speeds and evidence for high gas fractions and dark matter domination within the luminous extent of the galaxy, are consistent with those of previously studied faint galaxies in nearby groups. Further studies of the SDSS sample will allow characterization of low luminosity galaxies over the full range of environments in which they reside.

The Spitzer Space Telescope First Look Survey Infrared Array Camera (IRAC) near and mid-infrared imaging data partially overlaps the Sloan Digital Sky Survey (SDSS), with 313 visually selected (r<17.6 mag) SDSS Main Sample galaxies in the overlap region. The 3.5 and 7.8 um properties of the galaxies are investigated in the context of their visual properties, where the IRAC [3.5] magnitude primarily measures starlight, and the [7.8] magnitude primarily measures PAH emission from the interstellar medium. As expected, we find a strong inverse correlation between [3.5]-[7.8] and visual color; galaxies red in visual colors (`red galaxies') tend to show very little dust and molecular emission (low `PAH-to-star' ratios), and galaxies blue in visual colors (`blue galaxies,' ie, star-forming galaxies) tend to show large PAH-to-star ratios. Red galaxies with high PAH-to-star ratios tend to be edge-on disks reddened by dust lanes. Simple, visually inferred attenuation corrections bring the visual colors of these galaxies in line with those of face-on disks; ie, PAH emission is closely related to attenuation-corrected, optically inferred star-formation rates. Blue galaxies with anomalously low PAH-to-star ratios are all low-luminosity star-forming galaxies. There is some weak evidence in this sample that the deficiency in PAH emission for these low-luminosity galaxies may be related to emission-line metallicity.

We study the relationship between galaxy luminosity, color, and environment in a cosmological simulation of galaxy formation. We compare the predicted relationship with that found for SDSS galaxies and find that the model successfully predicts most of the qualitative features seen in the data, but also shows some interesting differences. Specifically, the simulation predicts that the local density around bright red galaxies is a strong increasing function of luminosity, but does not depend much on color at fixed luminosity. Moreover, we show that this is due to central galaxies in dark matter halos whose baryonic masses correlate strongly with halo mass. The simulation also predicts that the local density around blue galaxies is a strong increasing function of color, but does not depend much on luminosity at fixed color. We show that this is due to satellite galaxies in halos whose stellar ages correlate with halo mass. Finally, the simulation fails to predict the luminosity dependence of environment observed around low luminosity red galaxies. However, we show that this is most likely due to the simulation's limited resolution. A study of a higher resolution, smaller volume simulation suggests that this dependence is caused by the fact that all low luminosity red galaxies are satellites in massive halos, whereas intermediate luminosity red galaxies are a mixture of satellites in massive halos and central galaxies in less massive halos.

The Hubble Space Telescope has obtained some of its well-deserved impact by producing stunning three-color (RGB) pictures from three-band imaging data. Here we produce a new RGB representation of the $I$, $V$, and $B$ bandpass images of the HST Ultra Deep Field (HUDF). Our representation is based on principles set forth elsewhere (Lupton et al 2004, PASP, 116, 133--137). The principal difference between our RGB representation of the data and the more traditional representation provided by the Space Telescope Science Institute is that the (necessarily) nonlinear transformation between data values and RGB values is done in a color-preserving (i.e., hue- and saturation-preserving) way. For example, if one of the image pixel values saturates the dynamic range of the RGB representation, all three of the R, G, and B values in the representation are truncated such that the hue and saturation of the pixel is the same as it would have been if the pixel had had lower flux but the same astronomical color. This, in effect, makes the bright parts of the representation an informative color map, not "whited out" as they are in traditional representations. For the HUDF, this difference is seen best in the centers of bright galaxies, which show significant detail not visible in the more traditional representation.

We summarize the detection rates at wavelengths other than optical for \~99,000 galaxies from the Sloan Digital Sky Survey (SDSS) Data Release 1 ``main'' spectroscopic sample. The analysis is based on positional cross-correlation with source catalogs from ROSAT, 2MASS, IRAS, GB6, FIRST, NVSS and WENSS surveys. We find that the rest-frame UV-IR broad-band galaxy SEDs form a remarkably uniform, nearly one parameter, family. As an example, the SDSS u and r band data, supplemented with redshift, can be used to predict K band magnitudes measured by 2MASS with an rms scatter of only 0.2 mag; when measurement uncertainties are taken into account, the astrophysical scatter appears not larger than ~0.1 mag.

We present estimates of cosmological parameters from the application of the Karhunen-Loeve transform to the analysis of the 3D power spectrum of density fluctuations using Sloan Digital Sky Survey galaxy redshifts. We use Omega_m*h and f_b = Omega_b/Omega_m to describe the shape of the power spectrum, sigma8 for the (linearly extrapolated) normalization, and beta to parametrize linear theory redshift space distortions. On scales k < 0.16 h/Mpc, our maximum likelihood values are Omega_m*h = 0.264 +/-0.043, f_b = 0.286 +/- 0.065, sigma8 = 0.966 +/- 0.048, and beta = 0.45 +/- 0.12. When we take a prior on Omega_b from WMAP, we find Omega_m*h = 0.207 +/- 0.030, which is in excellent agreement with WMAP and 2dF. This indicates that we have reasonably measured the gross shape of the power spectrum but we have difficulty breaking the degeneracy between Omega_m*h and f_b because the baryon oscillations are not resolved in the current spectroscopic survey window function.

We measure cosmological parameters using the three-dimensional power spectrum P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and other data. Our results are consistent with a ``vanilla'' flat adiabatic Lambda-CDM model without tilt (n=1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than halves the WMAP-only error bars on some parameters, tightening 1 sigma constraints on the Hubble parameter from h~0.74+0.18-0.07 to h~0.70+0.04-0.03, on the matter density from Omega_m~0.25+/-0.10 to Omega_m~0.30+/-0.04 (1 sigma) and on neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical origin of the constraints, i.e., what we do and do not know when using different data sets and prior assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only constraint on the measured age of the Universe tightens from t0~16.3+2.3-1.8 Gyr to t0~14.1+1.0-0.9 Gyr by adding SDSS and SN Ia data. Including tensors, running tilt, neutrino mass and equation of state in the list of free parameters, many constraints are still quite weak, but future cosmological measurements from SDSS and other sources should allow these to be substantially tightened.

Compact groups (CGs) of galaxies -- relatively poor groups of galaxies in which the typical separations between members is of the order of a galaxy diameter -- offer an exceptional laboratory for the study of dense galaxian environments with short (<1Gyr) dynamical time-scales. In this paper, we present an objectively defined catalog of CGs in 153 sq deg of the Sloan Digital Sky Survey Early Data Release (SDSS EDR). To identify CGs, we applied a modified version of Hickson's (1982) criteria aimed at finding the highest density CGs and thus reducing the number of chance alignments. Our catalog contains 175 CGs down to a limiting galaxy magnitude of r* = 21. The resulting catalog has a median depth of approximately z = 0.13, substantially deeper than previous CG catalogs. Since the SDSS will eventually image up to one quarter of the celestial sphere, we expect our final catalog, based upon the completed SDSS, will contain on the order of 5,000 - 10,000 CGs. This catalog will be useful for conducting studies of the general characteristics of CGs, their environments, and their component galaxies.