Source author record

Daniel J. Eisenstein

Daniel J. Eisenstein appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

135works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

135 published item(s)

preprint2026arXiv

On the origins of oxygen: ALMA and JWST characterise the multi-phase, metal-enriched, star-bursting medium within a 'normal' $z > 11$ galaxy

The unexpectedly high abundance of galaxies at $z > 11$ revealed by JWST has sparked a debate on the nature of early galaxies and the physical mechanisms regulating their formation. The Atacama Large Millimeter/submillimeter Array (ALMA) has begun to provide vital insights on their gas and dust content, but so far only for extreme 'blue monsters'. Here we present new, deep ALMA observations of JADES-GS-z11-0, a more typical (sub-$L^*$) $z > 11$ galaxy that bridges the discovery space of JWST and the Hubble Space Telescope. These data confirm the presence of the [O III] 88 $μ$m line at $4.5σ$ significance, precisely at the redshift of several faint emission lines previously seen with JWST/NIRSpec, while the underlying dust continuum remains undetected ($F_ν< 9.0 \, \mathrm{μJy}$), implying an obscured star formation rate (SFR) of $\text{SFR}_\text{IR} \lesssim 6 \, \mathrm{M_\odot \, yr^{-1}}$ and dust mass of $M_\text{dust} \lesssim 1.0 \times 10^{6} \, \mathrm{M_\odot}$ (all $3σ$). The accurate ALMA redshift of $z_\text{[O III]} = 11.1221 \pm 0.0006$ ($\gtrsim \! 5\times$ refined over NIRSpec) helps confirm that redshifts measured purely from the Lyman-$α$ break, even spectroscopically, should properly take into account the effects of potential damped Lyman-$α$ absorption (DLA) systems to avoid systematic overestimates of up to $Δz \approx 0.5$. The [O III] 88 $μ$m luminosity of $L_\text{[O III]} = (1.1 \pm 0.3) \times 10^{8} \, \mathrm{L_\odot}$, meanwhile, agrees well with the scaling relation for local metal-poor dwarfs given the SFR measured by NIRCam, NIRSpec, and MIRI. The spatially resolved MIRI and ALMA emission also underscores that JADES-GS-z11-0 is likely to consist of two low-mass components that are undergoing strong bursts of star formation yet are already pre-enriched in oxygen (~20-30% solar), only 400 Myr after the Big Bang.

preprint2025arXiv

Galaxy-Multiplet Clustering from DESI DR2

We present an efficient estimator for higher-order galaxy clustering using small groups of nearby galaxies, or multiplets. Using the Luminous Red Galaxy (LRG) sample from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2, we identify galaxy multiplets as discrete objects and measure their cross-correlations with the general galaxy field. Our results show that the multiplets exhibit stronger clustering bias as they trace more massive dark matter halos than individual galaxies. When comparing the observed clustering statistics with the mock catalogs generated from the N-body simulation AbacusSummit, we find that the mocks underpredict multiplet clustering despite reproducing the galaxy two-point auto-correlation reasonably well. This discrepancy indicates that the standard Halo Occupation Distribution (HOD) model is insufficient to describe the properties of galaxy multiplets, revealing the greater constraining power of this higher-order statistic on galaxy-halo connection and the possibility that multiplets are specific to additional assembly bias. We demonstrate that incorporating secondary biases into the HOD model improves agreement with the observed multiplet statistics, specifically by allowing galaxies to preferentially occupy halos in denser environments. Our results highlight the potential of utilizing multiplet clustering, beyond traditional two-point correlation measurements, to break degeneracies in models describing the galaxy-dark matter connection.

preprint2024arXiv

Inside the bubble: exploring the environments of reionisation-era Lyman-$α$ emitting galaxies with JADES and FRESCO

We present a study of the environments of 17 Lyman-$α$ (Ly$α$) emitting galaxies (LAEs) in the reionisation era ($5.8 < z < 8$) identified by JWST/NIRSpec as part of the JWST Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Ly$α$ emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Ly$α$ velocity offsets ($Δv_\text{Ly$α$}<300\,\mathrm{km\,s^{-1}}$) with moderately high Ly$α$ escape fractions ($f_\mathrm{esc,\,Lyα}>5\%$) observed in our sample of LAEs, suggesting the presence of ionised hydrogen along the line of sight towards at least eight out of 17 LAEs. We find minimum physical `bubble' sizes of the order of $R_\text{ion}\sim0.1$-$1\,\mathrm{pMpc}$ are required in a patchy reionisation scenario where ionised bubbles containing the LAEs are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at $z \sim 5.8$-$5.9$ and $z\sim7.3$, suggesting Ly$α$ transmission is strongly enhanced in such overdense regions, and underlining the importance of LAEs as tracers of the first large-scale ionised bubbles. Considering only spectroscopically confirmed galaxies, we find our sample of UV-faint LAEs ($M_\text{UV}\gtrsim-20\,\mathrm{mag}$) and their direct neighbours are generally not able to produce the required ionised regions based on the Ly$α$ transmission properties, suggesting lower-luminosity sources likely play an important role in carving out these bubbles. These observations demonstrate the combined power of JWST multi-object and slitless spectroscopy in acquiring a unique view of the early Universe during cosmic reionisation via the most distant LAEs.

preprint2022arXiv

Constructing high-fidelity halo merger trees in AbacusSummit

Tracking the formation and evolution of dark matter haloes is a critical aspect of any analysis of cosmological $N$-body simulations. In particular, the mass assembly of a halo and its progenitors, encapsulated in the form of its merger tree, serves as a fundamental input for constructing semi-analytic models of galaxy formation and, more generally, for building mock catalogues that emulate galaxy surveys. We present an algorithm for constructing halo merger trees from AbacusSummit, the largest suite of cosmological $N$-body simulations performed to date consisting of nearly 60 trillion particles, and which has been designed to meet the Cosmological Simulation Requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. Our method tracks the cores of haloes to determine associations between objects across multiple timeslices, yielding lists of halo progenitors and descendants for the several tens of billions of haloes identified across the entire suite. We present an application of these merger trees as a means to enhance the fidelity of AbacusSummit halo catalogues by flagging and "merging" haloes deemed to exhibit non-monotonic past merger histories. We show that this cleaning technique identifies portions of the halo population that have been deblended due to choices made by the halo finder, but which could have feasibly been part of larger aggregate systems. We demonstrate that by cleaning halo catalogues in this post-processing step, we remove potentially unphysical features in the default halo catalogues, leaving behind a more robust halo population that can be used to create highly-accurate mock galaxy realisations from AbacusSummit.

preprint2022arXiv

Illustrating galaxy-halo connection in the DESI era with IllustrisTNG

We employ the hydrodynamical simulation IllustrisTNG to inform the galaxy-halo connection of the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples of the Dark Energy Spectroscopic Instrument (DESI) survey at redshift z ~ 0.8. Specifically, we model the galaxy colors of IllustrisTNG and apply sliding DESI color-magnitude cuts, matching the DESI target densities. We study the halo occupation distribution model (HOD) of the selected samples by matching them to their corresponding dark matter halos in the IllustrisTNG dark matter run. We find the HOD of both the LRG and ELG samples to be consistent with their respective baseline models, but also we find important deviations from common assumptions about the satellite distribution, velocity bias, and galaxy secondary biases. We identify strong evidence for concentration-based and environment-based occupational variance in both samples, an effect known as "galaxy assembly bias". The central and satellite galaxies have distinct dependencies on secondary halo properties, showing that centrals and satellites have distinct evolutionary trajectories and should be modelled separately. These results serve to inform the necessary complexities in modeling galaxy-halo connection for DESI analyses and also prepare for building high-fidelity mock galaxies. Finally, we present a shuffling-based clustering analysis that reveals a 10-15% excess in the LRG clustering of modest statistical significance due to secondary galaxy biases. We also find a similar excess signature for the ELGs, but with much lower statistical significance. When a larger hydrodynamical simulation volume becomes available, we expect our analysis pipeline to pinpoint the exact sources of such excess clustering signatures.

preprint2022arXiv

Stringent $σ_8$ constraints from small-scale galaxy clustering using a hybrid MCMC+emulator framework

We present a novel simulation-based hybrid emulator approach that maximally derives cosmological and Halo Occupation Distribution (HOD) information from non-linear galaxy clustering, with sufficient precision for DESI Year 1 (Y1) analysis. Our hybrid approach first samples the HOD space on a fixed cosmological simulation grid to constrain the high-likelihood region of cosmology+HOD parameter space, and then constructs the emulator within this constrained region. This approach significantly reduces the parameter volume emulated over, thus achieving much smaller emulator errors with fixed number of training points. We demonstrate that this combined with state-of-the-art simulations result in tight emulator errors comparable to expected DESI Y1 LRG sample variance. We leverage the new AbacusSummit simulations and apply our hybrid approach to CMASS non-linear galaxy clustering data. We infer constraints on $σ_8 = 0.762\pm0.024$ and $fσ_8 (z_{eff} = 0.52) = 0.444\pm0.016$, the tightest among contemporary galaxy clustering studies. We also demonstrate that our $fσ_8$ constraint is robust against secondary biases and other HOD model choices, a critical first step towards showcasing the robust cosmology information accessible in non-linear scales. We speculate that the additional statistical power of DESI Y1 should tighten the growth rate constraints by at least another 50-60%, significantly elucidating any potential tension with Planck. We also address the "lensing is low" tension, where we find that the combined effect of a lower $fσ_8$ and environment-based bias lowers the predicted lensing signal by 15%, accounting for approximately 50% of the discrepancy between the lensing measurement and clustering-based predictions.

preprint2022arXiv

Target Selection and Validation of DESI Luminous Red Galaxies

The Dark Energy Spectroscopic Instrument (DESI) is carrying out a 5-year survey that aims to measure the redshifts of tens of millions of galaxies and quasars, including 8 million luminous red galaxies (LRGs) in the redshift range of $0.4<z<{\sim}\,1.0$. Here we present the selection of the DESI LRG sample and assess its spectroscopic performance using data from Survey Validation (SV) and the first 2 months of the Main Survey. The DESI LRG sample, selected using $g$, $r$, $z$, and $W1$ photometry from the DESI Legacy Imaging Surveys, is highly robust against imaging systematics. The sample has a target density of 605 deg$^{-2}$ and a comoving number density of $5\times10^{-4}\ h^3\mathrm{Mpc}^{-3}$ in $0.4<z<0.8$; this is a significantly higher density than previous LRG surveys (such as SDSS, BOSS and eBOSS) while also extending to $z \sim 1$. After applying a bright star veto mask developed for the sample, $98.9\%$ of the observed LRG targets yield confident redshifts (with a catastrophic failure rate of $0.2\%$ in the confident redshifts), and only $0.5\%$ of the LRG targets are stellar contamination. The LRG redshift efficiency varies with source brightness and effective exposure time, and we present a simple model that accurately characterizes this dependence. In the appendices, we describe the extended LRG samples observed during SV.

preprint2022arXiv

The DESI $N$-body Simulation Project -- II. Suppressing sample variance with fast simulations

Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise three-dimensional map of our Universe. The survey effective volume reaches $\sim20\Gpchcube$. It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. \textsc{AbacusSummit} is a suite of high-resolution dark-matter-only simulations designed for this purpose, with $200\Gpchcube$ (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide a more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi-$N$-body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variance than the expected DESI statistical variance at the scales we are interested in, e.g. $k < 0.3\hMpc$ for the halo power spectrum. In addition, it can significantly suppress the sample variance of the halo bispectrum. We further generalize the method for other cosmologies with only one realization in \textsc{AbacusSummit} suite to extend the effective volume $\sim 20$ times. In summary, our proposed strategy of combining high-fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data.

preprint2021arXiv

Checkpointing with cp: the POSIX Shared Memory System

We present the checkpointing scheme of Abacus, an $N$-body simulation code that allocates all persistent state in POSIX shared memory, or ramdisk. Checkpointing becomes as simple as copying files from ramdisk to external storage. The main simulation executable is invoked once per time step, memory mapping the input state, computing the output state directly into ramdisk, and unmapping the input state. The main executable remains unaware of the concept of checkpointing, with the top-level driver code launching a file-system copy between executable invocations when a checkpoint is needed. Since the only information flow is through files on ramdisk, the checkpoint must be correct so long as the simulation is correct. However, we find that with multi-GB of state, there is a significant overhead to unmapping the shared memory. This can be partially mitigated with multithreading, but ultimately, we do not recommend shared memory for use with a large state.

preprint2021arXiv

Evidence for galaxy assembly bias in BOSS CMASS redshift-space galaxy correlation function

Building accurate and flexible galaxy-halo connection models is crucial in modeling galaxy clustering on non-linear scales. Recent studies have found that halo concentration by itself cannot capture the full galaxy assembly bias effect and that the local environment of the halo can be an excellent indicator of galaxy assembly bias. In this paper, we propose an extended halo occupation distribution model (HOD) that includes both a concentration-based assembly bias term and an environment-based assembly bias term. We use this model to achieve a good fit (chi^2/DoF = 1.35) on the 2D redshift-space 2-point correlation function (2PCF) of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy sample. We find that the inclusion of both assembly bias terms is strongly favored by the data and the standard 5-parameter HOD is strongly rejected. More interestingly, the redshift-space 2PCF drives the assembly bias parameters in a way that preferentially assigns galaxies to lower mass halos. This results in galaxy-galaxy lensing predictions that are within 1sigma agreement with the observation, alleviating the perceived tension between galaxy clustering and lensing. We also showcase a consistent 3-5sigma preference for a positive environment-based assembly bias that persists over variations in the fit. We speculate that the environmental dependence might be driven by underlying processes such as mergers and feedback, but might also be indicative of a larger halo boundaries such as the splashback radius. Regardless, this work highlights the importance of building flexible galaxy-halo connection models and demonstrates the extra constraining power of the redshift-space 2PCF.

preprint2021arXiv

The galaxy-halo connection of emission-line galaxies in IllustrisTNG

We employ the hydrodynamical simulation IllustrisTNG-300-1 to explore the halo occupation distribution (HOD) and environmental dependence of luminous star-forming emission-line galaxies (ELGs) at $z \sim 1$. Such galaxies are key targets for current and upcoming cosmological surveys. We select model galaxies through cuts in colour-colour space allowing for a direct comparison with the Extended Baryon Oscillation Spectroscopic Survey and the Dark Energy Spectroscopic Instrument (DESI) surveys and then compare them with galaxies selected based on specific star-formation rate (sSFR) and stellar mass. We demonstrate that the ELG populations are twice more likely to reside in lower-density regions (sheets) compared with the mass-selected populations and twice less likely to occupy the densest regions of the cosmic web (knots). We also show that the colour-selected and sSFR-selected ELGs exhibit very similar occupation and clustering statistics, finding that the agreement is best for lower redshifts. In contrast with the mass-selected sample, the occupation of haloes by a central ELG peaks at $\sim$20\%. We furthermore explore the dependence of the HOD and the auto-correlation on environment, noticing that at fixed halo mass, galaxies in high-density regions cluster about 10 times more strongly than low-density ones. This result suggests that we should model carefully the galaxy-halo relation and implement assembly bias effects into our models (estimated at $\sim$4\% of the clustering of the DESI colour-selected sample at $z = 0.8$). Finally, we apply a simple mock recipe to recover the clustering on large scales ($r \gtrsim 1 \ {\rm Mpc}/h$) to within 1\% by augmenting the HOD model with an environment dependence, demonstrating the power of adopting flexible population models.

preprint2020arXiv

Baryon Acoustic Oscillations in the projected cross-correlation function between the eBOSS DR16 quasars and photometric galaxies from the DESI Legacy Imaging Surveys

We search for the Baryon Acoustic Oscillations in the projected cross-correlation function binned into transverse comoving radius between the SDSS-IV DR16 eBOSS quasars and a dense photometric sample of galaxies selected from the DESI Legacy Imaging Surveys. We estimate the density of the photometric sample of galaxies in this redshift range to be about 2900 deg$^{-2}$, which is deeper than the official DESI ELG selection, and the density of the spectroscopic sample is about 20 deg$^{-2}$. In order to mitigate the systematics related to the use of different imaging surveys close to the detection limit, we use a neural network approach that accounts for complex dependencies between the imaging attributes and the observed galaxy density. We find that we are limited by the depth of the imaging surveys which affects the density and purity of the photometric sample and its overlap in redshift with the quasar sample, which thus affects the performance of the method. When cross-correlating the photometric galaxies with quasars in $0.6 \leq z \leq 1.2$, the cross-correlation function can provide better constraints on the comoving angular distance, $D_{\rm M}$ (6\% precision) compared to the constraint on the spherically-averaged distance $D_{\rm V}$ (9\% precision) obtained from the auto-correlation. Although not yet competitive, this technique will benefit from the arrival of deeper photometric data from upcoming surveys which will enable it to go beyond the current limitations we have identified in this work.

preprint2020arXiv

Can Assembly Bias Explain the Lensing Amplitude of the BOSS CMASS Sample in a Planck Cosmology?

In this paper, we investigate whether galaxy assembly bias can reconcile the 20-40% disagreement between the observed galaxy projected clustering signal and the galaxy-galaxy lensing signal in the BOSS CMASS galaxy sample reported in Leauthaud et al. (2017). We use the suite of AbacusCosmos Lamda-CDM simulations at Planck best-fit cosmology and two flexible implementations of extended halo occupation distribution (HOD) models that incorporate galaxy assembly bias to build forward models and produce joint fits of the observed galaxy clustering signal and the galaxy-galaxy lensing signal. We find that our models using the standard HODs without any assembly bias generalizations continue to show a 20-40% over-prediction of the observed galaxy-galaxy lensing signal. We find that our implementations of galaxy assembly bias do not reconcile the two measurements at Planck best-fit cosmology. In fact, despite incorporating galaxy assembly bias, the satellite distribution parameter, and the satellite velocity bias parameter into our extended HOD model, our fits still strongly suggest a 31-34% discrepancy between the observed projected clustering and galaxy-galaxy lensing measurements. It remains to be seen whether a combination of other galaxy assembly bias models, alternative cosmological parameters, or baryonic effects can explain the amplitude difference between the two signals.

preprint2020arXiv

HST Grism-derived Forecasts for Future Galaxy Redshift Surveys

The mutually complementary Euclid and Roman galaxy redshift surveys will use Halpha- and [OIII]-selected emission line galaxies as tracers of the large scale structure at $0.9 \lesssim z \lesssim 1.9$ (Halpha) and $1.5 \lesssim z \lesssim 2.7$ ([OIII]). It is essential to have a reliable and sufficiently precise knowledge of the expected numbers of Halpha-emitting galaxies in the survey volume in order to optimize these redshift surveys for the study of dark energy. Additionally, these future samples of emission-line galaxies will, like all slitless spectroscopy surveys, be affected by a complex selection function that depends on galaxy size and luminosity, line equivalent width, and redshift errors arising from the misidentification of single emission-line galaxies. Focusing on the specifics of the Euclid survey, we combine two slitless spectroscopic WFC3-IR datasets -- 3D-HST+AGHAST and the WISP survey -- to construct a Euclid-like sample that covers an area of 0.56 deg$^2$ and includes 1277 emission line galaxies. We detect 1091 ($\sim$3270 deg$^{-2}$) Halpha+[NII]-emitting galaxies in the range $0.9\leq z \leq 1.6$ and 162 ($\sim$440 deg$^{-2}$) [OIII]$λ$5007-emitters over $1.5\leq z \leq 2.3$ with line fluxes $\geq 2 \times 10^{-16}$ erg s$^{-1}$ cm$^{-2}$. The median of the Halpha+[NII] equivalent width distribution is $\sim$250Å, and the effective radii of the continuum and Halpha+[NII] emission are correlated with a median of $\sim$0.38" and significant scatter ($σ\sim $0.2"$-$0.35"). Finally, we explore the prevalence of redshift misidentification in future Euclid samples, finding potential contamination rates of $\sim$14-20% and $\sim$6% down to $2\times 10^{-16}$ and $6 \times 10^{-17}$ erg s$^{-1}$ cm$^{-2}$, respectively, though with increased wavelength coverage these percentages drop to nearly zero.

preprint2020arXiv

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: BAO and RSD measurements from anisotropic clustering analysis of the Quasar Sample in configuration space between redshift 0.8 and 2.2

We measure the anisotropic clustering of the quasar sample from Data Release 16 (DR16) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (eBOSS). A sample of $343,708$ spectroscopically confirmed quasars between redshift $0.8<z<2.2$ are used as tracers of the underlying dark matter field. In comparison with DR14 sample, the final sample doubles the number of objects as well as the survey area. In this paper, we present the analysis in configuration space by measuring the two-point correlation function and decompose using the Legendre polynomials. For the full-shape analysis of the Legendre multipole moments, we measure the BAO distance and the growth rate of the cosmic structure. At an effective redshift of $z_{\rm eff}=1.48$, we measure the comoving angular diameter distance $D_{\rm M}(z_{\rm eff})/r_{\rm drag} = 30.66\pm0.88$, the Hubble distance $D_{\rm H}(z_{\rm eff})/r_{\rm drag} = 13.11\pm0.52$, and the growth rate $fσ_8(z_{\rm eff}) = 0.439\pm0.048$. The accuracy of these measurements is confirmed using an extensive set of mock simulations developed for the quasar sample. The uncertainties on the distance and growth rate measurements have been reduced substantially ($\sim 45\%$ and $\sim30\%$) with respect to the DR14 results. We also perform a BAO-only analysis to cross check the robustness of the methodology of the full-shape analysis. Combining our analysis with the Fourier space analysis, we arrive at $D^{\bf{c}}_{\rm M}(z_{\rm eff})/r_{\rm drag} = 30.22 \pm 0.79$, $D^{\bf{c}}_{\rm H}(z_{\rm eff})/r_{\rm drag} = 13.26 \pm 0.47$, and $fσ_8^{\bf{c}}(z_{\rm eff}) = 0.464 \pm 0.045$.

preprint2020arXiv

Validation of Emission-Line Galaxies Target Selection Algorithms for the Dark Energy Spectroscopic Instrument Using the MMT Binospec

The forthcoming Dark Energy Spectroscopic Instrument (DESI) experiment plans to measure the effects of dark energy on the expansion of the Universe and create a $3$D map of the Universe using galaxies up to $z \sim 1.6$ and QSOs up to $z \sim 3.5$. In order to create this map, DESI will obtain spectroscopic redshifts of over $30$ million objects; among them, a majority are \oii emitting star-forming galaxies known as emission-line galaxies (ELGs). These ELG targets will be pre-selected by drawing a selection region on the $g - r$ vs. $r - z$ colour-colour plot, where high redshift ELGs form a separate locus from the lower redshift ELGs and interlopers. In this paper, we study the efficiency of three ELG target selection algorithms -- the final design report (FDR) cut based on the DEEP2 photometry, Number Density Modelling and Random Forest -- to determine how the combination of these three algorithms can be best used to yield a simple selection boundary that will be best suited to meet DESI's science goals. To do this, we selected $17$ small patches in the DESI footprint where we run the three target selection algorithms to pre-select ELGs based on their photometry. We observed the pre-selected ELGs using the MMT Binospec, which is similar in functionality to the DESI instrument, to obtain their spectroscopic redshifts and fluxes of $1054$ ELGs. By analysing the redshift and fluxing distribution of these galaxies, we find that although NDM performed the best, simple changes in the FDR definition would also yield sufficient performance.

preprint2019arXiv

A Hybrid Deep Learning Approach to Cosmological Constraints From Galaxy Redshift Surveys

We present a deep machine learning (ML)-based technique for accurately determining $σ_8$ and $Ω_m$ from mock 3D galaxy surveys. The mock surveys are built from the AbacusCosmos suite of $N$-body simulations, which comprises 40 cosmological volume simulations spanning a range of cosmological models, and we account for uncertainties in galaxy formation scenarios through the use of generalized halo occupation distributions (HODs). We explore a trio of ML models: a 3D convolutional neural network (CNN), a power-spectrum-based fully connected network, and a hybrid approach that merges the two to combine physically motivated summary statistics with flexible CNNs. We describe best practices for training a deep model on a suite of matched-phase simulations and we test our model on a completely independent sample that uses previously unseen initial conditions, cosmological parameters, and HOD parameters. Despite the fact that the mock observations are quite small ($\sim0.07h^{-3}\,\mathrm{Gpc}^3$) and the training data span a large parameter space (6 cosmological and 6 HOD parameters), the CNN and hybrid CNN can constrain $σ_8$ and $Ω_m$ to $\sim3\%$ and $\sim4\%$, respectively.

preprint2019arXiv

Cosmology with galaxy-galaxy lensing on non-perturbative scales: Emulation method and application to BOSS LOWZ

We describe our nonlinear emulation (i.e., interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with $N$-body simulations of nonlinear structure formation, designed to accurately predict the projected clustering and galaxy-galaxy lensing signals from luminous red galaxies (LRGs) in the redshift range $0.16 < z < 0.36$ on comoving scales $0.6 < r_p < 30$ \hMpc. The interpolation accuracy is $\lesssim 1-2$ per cent across the entire physically plausible range of parameters for all scales considered. We correctly recover the true value of the cosmological parameter $S_8 = ({σ_8}/{0.8228}) ({Ω_{\text{m}}}/{0.3107})^{0.6}$ from mock measurements produced via subhalo abundance matching (SHAM)-based lightcones designed to approximately match the properties of the SDSS LOWZ galaxy sample. Applying our model to Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 14 (DR14) LOWZ galaxy clustering and galaxy-shear cross-correlation measurements made with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) imaging, we perform a prototype cosmological analysis marginalizing over $w$CDM cosmological parameters and galaxy HOD parameters. We obtain a 4.4 per cent measurement of $S_8 = 0.847 \pm 0.037$, in $3.5σ$ tension with the Planck cosmological results of $1.00 \pm 0.02$. We discuss the possibility of underestimated systematic uncertainties or astrophysical effects that could explain this discrepancy.

preprint2019arXiv

Galaxy-Galaxy Lensing in HSC: Validation Tests and the Impact of Heterogeneous Spectroscopic Training Sets

Although photometric redshifts (photo-z's) are crucial ingredients for current and upcoming large-scale surveys, the high-quality spectroscopic redshifts currently available to train, validate, and test them are substantially non-representative in both magnitude and color. We investigate the nature and structure of this bias by tracking how objects from a heterogeneous training sample contribute to photo-z predictions as a function of magnitude and color, and illustrate that the underlying redshift distribution at fixed color can evolve strongly as a function of magnitude. We then test the robustness of the galaxy-galaxy lensing signal in 120 deg$^2$ of HSC-SSP DR1 data to spectroscopic completeness and photo-z biases, and find that their impacts are sub-dominant to current statistical uncertainties. Our methodology provides a framework to investigate how spectroscopic incompleteness can impact photo-z-based weak lensing predictions in future surveys such as LSST and WFIRST.

preprint2019arXiv

RascalC: A Jackknife Approach to Estimating Single and Multi-Tracer Galaxy Covariance Matrices

To make use of clustering statistics from large cosmological surveys, accurate and precise covariance matrices are needed. We present a new code to estimate large scale galaxy two-point correlation function (2PCF) covariances in arbitrary survey geometries that, due to new sampling techniques, runs $\sim 10^4$ times faster than previous codes, computing finely-binned covariance matrices with negligible noise in less than 100 CPU-hours. As in previous works, non-Gaussianity is approximated via a small rescaling of shot-noise in the theoretical model, calibrated by comparing jackknife survey covariances to an associated jackknife model. The flexible code, RascalC, has been publicly released, and automatically takes care of all necessary pre- and post-processing, requiring only a single input dataset (without a prior 2PCF model). Deviations between large scale model covariances from a mock survey and those from a large suite of mocks are found to be be indistinguishable from noise. In addition, the choice of input mock are shown to be irrelevant for desired noise levels below $\sim 10^5$ mocks. Coupled with its generalization to multi-tracer data-sets, this shows the algorithm to be an excellent tool for analysis, reducing the need for large numbers of mock simulations to be computed.

preprint2016arXiv

Clustering properties of $g$-selected galaxies at $z\sim0.8$

Current and future large redshift surveys, as the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), will use emission-line galaxies (ELG) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range $0.6 < z < 1.7$. With current data, we explore the halo-galaxy connection by measuring three clustering properties of $g$-selected ELGs as matter tracers in the redshift range $0.6 < z < 1$: (i) the redshift-space two-point correlation function using spectroscopic redshifts from the BOSS ELG sample and VIPERS; (ii) the angular two-point correlation function on the footprint of the CFHT-LS; (iii) the galaxy-galaxy lensing signal around the ELGs using the CFHTLenS. We interpret these observations by mapping them onto the latest high-resolution MultiDark Planck N-body simulation, using a novel (Sub)Halo-Abundance Matching technique that accounts for the ELG incompleteness. ELGs at $z\sim0.8$ live in halos of $(1\pm 0.5)\times10^{12}\,h^{-1}$M$_{\odot}$ and 22.5$\pm2.5$% of them are satellites belonging to a larger halo. The halo occupation distribution of ELGs indicates that we are sampling the galaxies in which stars form in the most efficient way, according to their stellar-to-halo mass ratio.

preprint2016arXiv

Improving Initial Conditions for Cosmological $N$-Body Simulations

In cosmological $N$-body simulations, the representation of dark matter as discrete "macroparticles" suppresses the growth of structure, such that simulations no longer reproduce linear theory on small scales near $k_{\rm Nyquist}$. Marcos et al. demonstrate that this is due to sparse sampling of modes near $k_{\rm Nyquist}$ and that the often-assumed continuum growing modes are not proper growing modes of the particle system. We develop initial conditions that respect the particle linear theory growing modes and then rescale the mode amplitudes to account for growth suppression. These ICs also allow us to take advantage of our very accurate $N$-body code Abacus to implement 2LPT in configuration space. The combination of 2LPT and rescaling improves the accuracy of the late-time power spectra, halo mass functions, and halo clustering. In particular, we achieve 1% accuracy in the power spectrum down to $k_{\rm Nyquist}$, versus $k_{\rm Nyquist}/4$ without rescaling or $k_{\rm Nyquist}/13$ without 2LPT, relative to an oversampled reference simulation. We anticipate that our 2LPT will be useful for large simulations where FFTs are expensive and that rescaling will be useful for suites of medium-resolution simulations used in cosmic emulators and galaxy survey mock catalogs. Code to generate initial conditions is available at https://github.com/lgarrison/zeldovich-PLT

preprint2016arXiv

PRIMUS + DEEP2: Clustering of X-ray, Radio and IR-AGN at z~0.7

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.

preprint2016arXiv

Signatures of the Primordial Universe from Its Emptiness: Measurement of Baryon Acoustic Oscillations from Minima of the Density Field

Sound waves from the primordial fluctuations of the Universe imprinted in the large-scale structure, called baryon acoustic oscillations (BAOs), can be used as standard rulers to measure the scale of the Universe. These oscillations have already been detected in the distribution of galaxies. Here we propose to measure BAOs from the troughs (minima) of the density field. Based on two sets of accurate mock halo catalogues with and without BAOs in the seed initial conditions, we demonstrate that the BAO signal cannot be obtained from the clustering of classical disjoint voids, but is clearly detected from overlapping voids. The latter represent an estimate of all troughs of the density field. We compute them from the empty circumsphere centers constrained by tetrahedra of galaxies using Delaunay triangulation. Our theoretical models based on an unprecedented large set of detailed simulated void catalogues are remarkably well confirmed by observational data. We use the largest recently publicly available sample of luminous red galaxies from SDSS-III BOSS DR11 to unveil for the first time a >3$σ$ BAO detection from voids in observations. Since voids are nearly isotropically expanding regions, their centers represent the most quiet places in the Universe, keeping in mind the cosmos origin and providing a new promising window in the analysis of the cosmological large-scale structure from galaxy surveys.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in Fourier-space

We analyse the Baryon Acoustic Oscillation (BAO) signal of the final Baryon Oscillation Spectroscopic Survey (BOSS) data release (DR12). Our analysis is performed in Fourier-space, using the power spectrum monopole and quadrupole. The dataset includes $1\,198\,006$ galaxies over the redshift range $0.2 < z < 0.75$. We divide this dataset into three (overlapping) redshift bins with the effective redshifts $\zeff = 0.38$, $0.51$ and $0.61$. We demonstrate the reliability of our analysis pipeline using N-body simulations as well as $\sim 1000$ MultiDark-Patchy mock catalogues, which mimic the BOSS-DR12 target selection. We apply density field reconstruction to enhance the BAO signal-to-noise ratio. By including the power spectrum quadrupole we can separate the line-of-sight and angular modes, which allows us to constrain the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$ separately. We obtain two independent $1.6\%$ and $1.5\%$ constraints on $D_A(z)$ and $2.9\%$ and $2.3\%$ constraints on $H(z)$ for the low ($\zeff=0.38$) and high ($\zeff=0.61$) redshift bin, respectively. We obtain two independent $1\%$ and $0.9\%$ constraints on the angular averaged distance $D_V(z)$, when ignoring the Alcock-Paczynski effect. The detection significance of the BAO signal is of the order of $8σ$ (post-reconstruction) for each of the three redshift bins. Our results are in good agreement with the Planck prediction within $Λ$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in~\citet{Alam2016} to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: combining correlated Gaussian posterior distributions

The cosmological information contained in anisotropic galaxy clustering measurements can often be compressed into a small number of parameters whose posterior distribution is well described by a Gaussian. We present a general methodology to combine these estimates into a single set of consensus constraints that encode the total information of the individual measurements, taking into account the full covariance between the different methods. We illustrate this technique by applying it to combine the results obtained from different clustering analyses, including measurements of the signature of baryon acoustic oscillations (BAO) and redshift-space distortions (RSD), based on a set of mock catalogues of the final SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). Our results show that the region of the parameter space allowed by the consensus constraints is smaller than that of the individual methods, highlighting the importance of performing multiple analyses on galaxy surveys even when the measurements are highly correlated. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The methodology presented here is used in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological implications of the configuration-space clustering wedges

We explore the cosmological implications of anisotropic clustering measurements in configuration space of the final galaxy samples from Data Release 12 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We implement a new detailed modelling of the effects of non-linearities, galaxy bias and redshift-space distortions that can be used to extract unbiased cosmological information from our measurements for scales $s \gtrsim 20\,h^{-1}{\rm Mpc}$. We combined the galaxy clustering information from BOSS with the latest cosmic microwave background (CMB) observations and Type Ia supernovae samples and found no significant evidence for a deviation from the $Λ$CDM cosmological model. In particular, these data sets can constrain the dark energy equation of state parameter to $w_{\rm DE}=-0.996\pm0.042$ when assumed time-independent, the curvature of the Universe to $Ω_{k}=-0.0007\pm 0.0030$ and the sum of the neutrino masses to $\sum m_ν < 0.25\,{\rm eV}$ at 95 per cent CL. We explore the constraints on the growth rate of cosmic structures assuming $f(z)=Ω_{\rm m}(z)^γ$ and obtain $γ= 0.609\pm 0.079$, in good agreement with the predictions of general relativity of $γ=0.55$. We compress the information of our clustering measurements into constraints on the parameter combinations $D_{\rm V}(z)/r_{\rm d}$, $F_{\rm AP}(z)$ and $fσ_8(z)$ at the effective redshifts of $z=0.38$, $0.51$ and $0.61$ with their respective covariance matrices and find good agreement with the predictions for these parameters obtained from the best-fitting $Λ$CDM model to the CMB data from the Planck satellite. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: double-probe measurements from BOSS galaxy clustering \& Planck data -- towards an analysis without informative priors

We develop a new methodology called double-probe analysis with the aim of minimizing informative priors in the estimation of cosmological parameters. We extract the dark-energy-model-independent cosmological constraints from the joint data sets of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy sample and Planck cosmic microwave background (CMB) measurement. We measure the mean values and covariance matrix of $\{R$, $l_a$, $Ω_b h^2$, $n_s$, $log(A_s)$, $Ω_k$, $H(z)$, $D_A(z)$, $f(z)σ_8(z)\}$, which give an efficient summary of Planck data and 2-point statistics from BOSS galaxy sample, where $R=\sqrt{Ω_m H_0^2}\,r(z_*)$, and $l_a=πr(z_*)/r_s(z_*)$, $z_*$ is the redshift at the last scattering surface, and $r(z_*)$ and $r_s(z_*)$ denote our comoving distance to $z_*$ and sound horizon at $z_*$ respectively. The advantage of this method is that we do not need to put informative priors on the cosmological parameters that galaxy clustering is not able to constrain well, i.e. $Ω_b h^2$ and $n_s$. Using our double-probe results, we obtain $Ω_m=0.304\pm0.009$, $H_0=68.2\pm0.7$, and $σ_8=0.806\pm0.014$ assuming $Λ$CDM; and $Ω_k=0.002\pm0.003$ and $w=-1.00\pm0.07$ assuming o$w$CDM. The results show no tension with the flat $Λ$CDM cosmological paradigm. By comparing with the full-likelihood analyses with fixed dark energy models, we demonstrate that the double-probe method provides robust cosmological parameter constraints which can be conveniently used to study dark energy models. We extend our study to measure the sum of neutrino mass and obtain $Σm_ν<0.10/0.22$ (68\%/95\%) assuming $Λ$CDM and $Σm_ν<0.26/0.52$ (68\%/95\%) assuming $w$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Observational systematics and baryon acoustic oscillations in the correlation function

We present baryon acoustic oscillation (BAO) scale measurements determined from the clustering of 1.2 million massive galaxies with redshifts 0.2 < z < 0.75 distributed over 9300 square degrees, as quantified by their redshift-space correlation function. In order to facilitate these measurements, we define, describe, and motivate the selection function for galaxies in the final data release (DR12) of the SDSS III Baryon Oscillation Spectroscopic Survey (BOSS). This includes the observational footprint, masks for image quality and Galactic extinction, and weights to account for density relationships intrinsic to the imaging and spectroscopic portions of the survey. We simulate the observed systematic trends in mock galaxy samples and demonstrate that they impart no bias on baryon acoustic oscillation (BAO) scale measurements and have a minor impact on the recovered statistical uncertainty. We measure transverse and radial BAO distance measurements in 0.2 < z < 0.5, 0.5 < z < 0.75, and (overlapping) 0.4 < z < 0.6 redshift bins. In each redshift bin, we obtain a precision that is 2.7 per cent or better on the radial distance and 1.6 per cent or better on the transverse distance. The combination of the redshift bins represents 1.8 per cent precision on the radial distance and 1.1 per cent precision on the transverse distance. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the correlation function of LOWZ and CMASS galaxies in Data Release 12

We present distance scale measurements from the baryon acoustic oscillation signal in the CMASS and LOWZ samples from the Data Release 12 of the Baryon Oscillation Spectroscopic Survey (BOSS). The total volume probed is 14.5 Gpc$^3$, a 10 per cent increment from Data Release 11. From an analysis of the spherically averaged correlation function, we infer a distance to $z=0.57$ of $D_V(z)r^{\rm fid}_{\rm d}/r_ {\rm d}=2028\pm21$ Mpc and a distance to $z=0.32$ of $D_V(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1264\pm22$ Mpc assuming a cosmology in which $r^{\rm fid}_{\rm d}=147.10$ Mpc. From the anisotropic analysis, we find an angular diameter distance to $z=0.57$ of $D_{\rm A}(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1401\pm21$ Mpc and a distance to $z=0.32$ of $981\pm20$ Mpc, a 1.5 per cent and 2.0 per cent measurement respectively. The Hubble parameter at $z=0.57$ is $H(z)r_{\rm d}/r^{\rm fid}_{\rm d}=100.3\pm3.7$ km s$^{-1}$ Mpc$^{-1}$ and its value at $z=0.32$ is $79.2\pm5.6$ km s$^{-1}$ Mpc$^{-1}$, a 3.7 per cent and 7.1 per cent measurement respectively. These cosmic distance scale constraints are in excellent agreement with a $Λ$CDM model with cosmological parameters released by the recent Planck 2015 results.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: single-probe measurements from CMASS anisotropic galaxy clustering

With the largest spectroscopic galaxy survey volume drawn from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), we can extract cosmological constraints from the measurements of redshift and geometric distortions at quasi-linear scales (e.g. above 50 $h^{-1}$Mpc). We analyze the broad-range shape of the monopole and quadrupole correlation functions of the BOSS Data Release 12 (DR12) CMASS galaxy sample, at the effective redshift $z=0.59$, to obtain constraints on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalized growth rate $f(z)σ_8(z)$, and the physical matter density $Ω_mh^2$. We obtain robust measurements by including a polynomial as the model for the systematic errors, and find it works very well against the systematic effects, e.g., ones induced by stars and seeing. We provide accurate measurements $\{D_A(0.59)r_{s,fid}/r_s$ $\rm Mpc$, $H(0.59)r_s/r_{s,fid}$ $km s^{-1} Mpc^{-1}$, $f(0.59)σ_8(0.59)$, $Ω_m h^2\}$ = $\{1427\pm26$, $97.3\pm3.3$, $0.488 \pm 0.060$, $0.135\pm0.016\}$, where $r_s$ is the comoving sound horizon at the drag epoch and $r_{s,fid}=147.66$ Mpc is the sound scale of the fiducial cosmology used in this study. The parameters which are not well constrained by our galaxy clustering analysis are marginalized over with wide flat priors. Since no priors from other data sets, e.g., cosmic microwave background (CMB), are adopted and no dark energy models are assumed, our results from BOSS CMASS galaxy clustering alone may be combined with other data sets, i.e., CMB, SNe, lensing or other galaxy clustering data to constrain the parameters of a given cosmological model. The uncertainty on the dark energy equation of state parameter, $w$, from CMB+CMASS is about 8 per cent. The uncertainty on the curvature fraction, $Ω_k$, is 0.3 per cent. We do not find deviation from flat $Λ$CDM.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2016arXiv

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data

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.

preprint2015arXiv

A new look at lines of sight: using Fourier methods for the wide-angle anisotropic 2-point correlation function

The anisotropic 2-point correlation function (2PCF) of galaxies measures pairwise clustering as a function of the pair separation's angle to the line of sight. The latter is often defined as either the angle bisector of the observer-galaxy-pair triangle or the vector from the observer to the separation midpoint. Here we show how to accelerate either of these measurements with Fourier Transforms, using a slight generalization of the Yamamoto et al. (2006) estimator in which each member of the pair is used successively as the line of sight. We also present perturbation theory predictions for our generalized estimator including wide-angle corrections.

preprint2015arXiv

A Simple Analytic Treatment of Linear Growth of Structure with Baryon Acoustic Oscillations

In linear perturbation theory, all information about the growth of structure is contained in the Green's function, or equivalently, transfer function. These functions are generally computed using numerical codes or by phenomenological fitting formula anchored in accurate analytic results in the limits of large and small scale. Here we present a framework for analytically solving all scales, in particular the intermediate scales relevant for the baryon acoustic oscillations (BAO). We solve for the Green's function and transfer function using spherically-averaged overdensities and the approximation that the density of the coupled baryon-photon fluid is constant interior to the sound horizon.

preprint2015arXiv

Abundances, Stellar Parameters, and Spectra From the SDSS-III/APOGEE Survey

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.

preprint2015arXiv

Accelerating the 2-point and 3-point galaxy correlation functions using Fourier transforms

Though Fourier Transforms (FTs) are a common technique for finding correlation functions, they are not typically used in computations of the anisotropy of the two-point correlation function (2PCF) about the line of sight in wide-angle surveys because the line-of-sight direction is not constant on the Cartesian grid. Here we show how FTs can be used to compute the multipole moments of the anisotropic 2PCF. We also show how FTs can be used to accelerate the 3PCF algorithm of Slepian & Eisenstein (2015). In both cases, these FT methods allow one to avoid the computational cost of pair counting, which scales as the square of the number density of objects in the survey. With the upcoming large datasets of DESI, Euclid, and LSST, FT techniques will therefore offer an important complement to simple pair or triplet counts.

preprint2015arXiv

ASPCAP: The Apogee Stellar Parameter and Chemical Abundances Pipeline

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

preprint2015arXiv

Building a Better Understanding of the High Redshift BOSS Galaxies as Tools for Cosmology

We explore the bluer star-forming population of the Sloan Digital Sky Survey (SDSS) III/BOSS CMASS DR11 galaxies at $z>0.55$ to quantify their differences, in terms of redshift-space distortions and large-scale bias, with respect to the luminous red galaxy sample. We perform a qualitative analysis to understand the significance of these differences and whether we can model and reproduce them in mock catalogs. Specifically, we measure galaxy clustering in CMASS on small and intermediate scales ($r\lesssim 50\,h^{-1}$Mpc) by computing the two-point correlation function $-$ both projected and redshift-space $-$ of these galaxies, and a new statistic, $Σ(π)$, able to provide robust information about redshift-space distortions and large-scale bias. We interpret our clustering measurements by adopting a Halo Occupation Distribution (HOD) scheme that maps them onto high-resolution N-body cosmological simulations to produce suitable mock galaxy catalogs. The traditional HOD prescription can be applied to the red and the blue samples, independently, but this approach is unphysical since it allows the same mock galaxies to be either red or blue. To overcome this failure, we modify the standard formulation and infer the red and the blue mock catalogs directly from the full one, so that they are complementary and non-overlapping. This separation is performed by matching the observed CMASS red and blue galaxy fractions and produces reliable and accurate models.

preprint2015arXiv

Computing the Three-Point Correlation Function of Galaxies in $\mathcal{O}(N^2)$ Time

We present an algorithm that computes the multipole coefficients of the galaxy three-point correlation function (3PCF) without explicitly considering triplets of galaxies. Rather, centering on each galaxy in the survey, it expands the radially-binned density field in spherical harmonics and combines these to form the multipoles without ever requiring the relative angle between a pair about the central. This approach scales with number and number density in the same way as the two-point correlation function, allowing runtimes that are comparable, and 500 times faster than a naive triplet count. It is exact in angle and easily handles edge correction. We demonstrate the algorithm on the LasDamas SDSS-DR7 mock catalogs, computing an edge corrected 3PCF out to $90\;{\rm Mpc}/h$ in under an hour on modest computing resources. We expect this algorithm will render it possible to obtain the large-scale 3PCF for upcoming surveys such as Euclid, LSST, and DESI.

preprint2015arXiv

Cosmological implications of baryon acoustic oscillation (BAO) measurements

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 $Ω_m=0.301 \pm 0.008$ and curvature $Ω_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 $χ^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)

preprint2015arXiv

Dark Matter Halo Models of Stellar Mass-Dependent Galaxy Clustering in PRIMUS+DEEP2 at 0.2<z<1.2

We utilize $Λ$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-à-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$.

preprint2015arXiv

Deriving Photometric Redshifts using Fuzzy Archetypes and Self-Organizing Maps. I. Methodology

We propose a method to substantially increase the flexibility and power of template fitting-based photometric redshifts by transforming a large numbers of galaxy spectral templates into a corrresponding collection of "fuzzy archetypes" using a suitable set of perturbative priors designed to account for empirical variation in dust attenuation and emission line strengths. To bypass widely seperated degeneracies in parameter space (e.g., the redshift-reddening degeneracy), we train Self-Organizing Maps (SOMs) on a large "model catalogs" generated from appropriate Monte Carlo sampling of our fuzzy archetypes to cluster the predicted observables in a topologically smooth fashion. Subsequent sampling over the SOM then allows full reconstruction of the relevant probability distribution functions (PDFs) using the associated set of inverse mappings from the SOM to the underlying model parameters. This combined approach enables the multi-modal exploration of known variation among galaxy spectral energy distributions (SEDs) using large numbers of archetypes with minimal modeling assumptions. We demonstrate the power of this approach to recover full redshift PDFs using discrete Markov Chain Monte Carlo (MCMC) sampling methods combined with SOMs constructed from model catalogs based on LSST $ugrizY$ and Euclid $YJH$ mock photometry.

preprint2015arXiv

Deriving Photometric Redshifts using Fuzzy Archetypes and Self-Organizing Maps. II. Comparing Sampling Techniques Using Mock Data

In a companion paper, we proposed combining large numbers of "fuzzy archetypes" with Self-Organizing Maps (SOMs) to derive photometric redshifts in a data-driven way. In this paper, we investigate the performance of several sampling approaches that build on this general idea using a mock catalog designed to approximately simulate LSST ($ugrizY$) and Euclid ($YJH$) data from $z=0-6$ at fixed LSST $Y=24$ mag. We test eight different approaches: two brute-force methods, two Markov Chain Monte Carlo (MCMC)-based methods, two hierarchical sampling methods, and two "quick-search" methods based on quantities derived during the initial SOM training process. We find most methods perform reasonably well with small catastrophic outlier fractions and are able to robustly identify redshift probability distribution functions that are multi-modal and/or poorly constrained. Once these insecure objects are removed, the results are generally in good agreement with the strict accuracy requirements necessary to meet Euclid weak lensing goals for most redshifts above $z \sim 0.8$. These results demonstrate the utility of our data clustering-based approach and highlight its effectiveness to derive quick and accurate photo-z's using large numbers of templates.

preprint2015arXiv

Exploring Photometric Redshifts as an Optimization Problem: An Ensemble MCMC and Simulated Annealing-Driven Template-Fitting Approach

Using a grid of $\sim 2$ million elements ($Δz = 0.005$) adapted from COSMOS photometric redshift (photo-z) searches, we investigate the general properties of template-based photo-z likelihood surfaces. We find these surfaces are filled with numerous local minima and large degeneracies that generally confound rapid but "greedy" optimization schemes, even with additional stochastic sampling methods. In order to robustly and efficiently explore these surfaces, we develop BAD-Z [Brisk Annealing-Driven Redshifts (Z)], which combines ensemble Markov Chain Monte Carlo (MCMC) sampling with simulated annealing to sample arbitrarily large, pre-generated grids in approximately constant time. Using a mock catalog of 384,662 objects, we show BAD-Z samples $\sim 40$ times more efficiently compared to a brute-force counterpart while maintaining similar levels of accuracy. Our results represent first steps toward designing template-fitting photo-z approaches limited mainly by memory constraints rather than computation time.

preprint2015arXiv

Measuring the luminosity and virial black hole mass dependence of quasar-galaxy clustering at z ~ 0.8

We study the dependence of quasar clustering on quasar luminosity and black hole mass by measuring the angular overdensity of photometrically selected galaxies imaged by WISE about z $\sim$ 0.8 quasars from SDSS. By measuring the quasar-galaxy cross-correlation function and using photometrically selected galaxies, we achieve a higher density of tracer objects and a more sensitive detection of clustering than measurements of the quasar autocorrelation function. We test models of quasar formation and evolution by measuring the luminosity dependence of clustering amplitude. We find a significant overdensity of WISE galaxies about z $\sim$ 0.8 quasars at 0.2--6.4 h$^{-1}$ Mpc in projected comoving separation. We find no appreciable increase in clustering amplitude with quasar luminosity across a decade in luminosity, and a power-law fit between luminosity and clustering amplitude gives an exponent of $-$0.01 $\pm$ 0.06 (1 $σ$ errorbar). We also fail to find a significant relationship between clustering amplitude and black hole mass, although our dynamic range in true mass is suppressed due to the large uncertainties in virial black hole mass estimates. Our results indicate that a small range in host dark matter halo mass maps to a large range in quasar luminosity.

preprint2015arXiv

Modeling Galactic Conformity with the Color-Halo Age Relation in the Illustris Simulation

Comparisons between observational surveys and galaxy formation models find that the mass of dark matter haloes can largely explain galaxies' stellar mass. However, it remains uncertain whether additional environmental variables, generally referred to as assembly bias, are necessary to explain other galaxy properties. We use the Illustris Simulation to investigate the role of assembly bias in producing galactic conformity by considering 18,000 galaxies with $M_{stellar}$ > $2 \times 10^9$ $M_{\odot}$. We find a significant signal of galactic conformity: out to distances of about 10 Mpc, the mean red fraction of galaxies around redder galaxies is higher than around bluer galaxies at fixed stellar mass. Dark matter haloes exhibit an analogous conformity signal, in which the fraction of haloes formed at earlier times (old haloes) is higher around old haloes than around younger ones at fixed halo mass. A plausible interpretation of galactic conformity can be given as a combination of the halo conformity signal with the galaxy color-halo age relation: at fixed stellar mass, particularly toward the low-mass end, Illustris' galaxy colors correlate with halo age, with the reddest galaxies (often satellites) being preferentially found in the oldest haloes. In fact, we can explain the galactic conformity effect with a simple semi-empirical model, by assigning stellar mass based on halo mass (abundance matching) and by assigning galaxy color based on halo age (age matching). We investigate other interpretations for the galactic conformity, particularly its dependence on the isolation criterion and on the central-satellite information. Regarding comparison to observations, we conclude that the adopted selection/isolation criteria, projection effects, and stacking techniques can have a significant impact on the measured amplitude of the conformity signal.

preprint2015arXiv

New white dwarf and subdwarf stars in the Sloan Digital Sky Survey Data Release 12

We report the discovery of 6576 new spectroscopically confirmed white dwarf and subdwarf stars in the Sloan Digital Sky Survey Data Release 12. We obtain Teff, log g and mass for hydrogen atmosphere white dwarf stars (DAs) and helium atmosphere white dwarf stars (DBs), estimate the calcium/helium abundances for the white dwarf stars with metallic lines (DZs) and carbon/helium for carbon dominated spectra DQs. We found one central star of a planetary nebula, one ultra-compact helium binary (AM CVn), one oxygen line dominated white dwarf, 15 hot DO/PG1159s, 12 new cataclysmic variables, 36 magnetic white dwarf stars, 54 DQs, 115 helium dominated white dwarfs, 148 white dwarf+main sequence star binaries, 236 metal polluted white dwarfs, 300 continuum spectra DCs, 230 hot subdwarfs, 2936 new hydrogen dominated white dwarf stars, and 2675 cool hydrogen dominated subdwarf stars. We calculate the mass distribution of all 5883 DAs with S/N>15 in DR12, including the ones in DR7 and DR10, with an average S/N=26, corrected to the 3D convection scale, and also the distribution after correcting for the observed volume, using 1/Vmax.

preprint2015arXiv

PRIMUS: The Effect of Physical Scale on the Luminosity-Dependence of Galaxy Clustering via Cross-Correlations

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.

preprint2015arXiv

Robust New Statistic for fitting the Baryon Acoustic Feature

We investigate the utility and robustness of a new statistic, $ω_{\ell}\left(r_{c}\right)$, for analyzing Baryon Acoustic Oscillations (BAO). We apply $ω_{\ell}\left(r_{c}\right)$, introduced in Xu et al. (2010), to mocks and data from the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) included in the SDSS Data Release Eleven (DR11). We fit the anisotropic clustering using the monopole and quadrupole of the $ω_{\ell}\left(r_{c}\right)$ statistic in a manner similar to conventional multipole fitting methods using the correlation function as detailed in (Xu et al. 2012). To test the performance of the $ω_{\ell}\left(r_{c}\right)$ statistic we compare our results to those obtained using the multipoles. The results are in agreement. We also conduct a brief investigation into some of the possible advantages of using the $ω_{\ell}\left(r_{c}\right)$ statistic for BAO analysis. The $ω_{\ell}\left(r_{c}\right)$ analysis matches the stability of the multipoles analysis in response to artificially introduced distortions in the data, without using extra nuisance parameters to improve the fit. When applied to data with systematics, the $ω_{\ell}\left(r_{c}\right)$ statistic again matches the performance of fitting the multipoles without using nuisance parameters. In all the analyzed circumstances, we find that fitting the $ω_{\ell}\left(r_{c}\right)$ statistic removes the requirement for extra nuisance parameters.

preprint2015arXiv

sCOLA: The N-body COLA Method Extended to the Spatial Domain

We present sCOLA -- an extension of the N-body COmoving Lagrangian Acceleration (COLA) method to the spatial domain. Similar to the original temporal-domain COLA, sCOLA is an N-body method for solving for large-scale structure in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory. Incorporating the sCOLA method in an N-body code allows one to gain computational speed by capturing the gravitational potential from the far field using perturbative techniques, while letting the N-body code solve only for the near field. The far and near fields are completely decoupled, effectively localizing gravity for the N-body side of the code. Thus, running an N-body code for a small simulation volume using sCOLA can reproduce the results of a standard N-body run for the same small volume embedded inside a much larger simulation. We demonstrate that sCOLA can be safely combined with the original temporal-domain COLA. sCOLA can be used as a method for performing zoom-in simulations. It also allows N-body codes to be made embarrassingly parallel, thus allowing for efficiently tiling a volume of interest using grid computing. Moreover, sCOLA can be useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering. Surveys that will benefit the most are ones with large aspect ratios, such as pencil-beam surveys, where sCOLA can easily capture the effects of large-scale transverse modes without the need to substantially increase the simulated volume. As an illustration of the method, we present proof-of-concept zoom-in simulations using a freely available sCOLA-based N-body code.

preprint2015arXiv

SDSS-III Baryon Oscillation Spectroscopic Survey Data Release 12: galaxy target selection and large scale structure catalogues

The Baryon Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey (SDSS) III project, has provided the largest survey of galaxy redshifts available to date, in terms of both the number of galaxy redshifts measured by a single survey, and the effective cosmological volume covered. Key to analysing the clustering of these data to provide cosmological measurements is understanding the detailed properties of this sample. Potential issues include variations in the target catalogue caused by changes either in the targeting algorithm or properties of the data used, the pattern of spectroscopic observations, the spatial distribution of targets for which redshifts were not obtained, and variations in the target sky density due to observational systematics. We document here the target selection algorithms used to create the galaxy samples that comprise BOSS. We also present the algorithms used to create large scale structure catalogues for the final Data Release (DR12) samples and the associated random catalogues that quantify the survey mask. The algorithms are an evolution of those used by the BOSS team to construct catalogues from earlier data, and have been designed to accurately quantify the galaxy sample. The code used, designated MKSAMPLE, is released with this paper.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

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.

preprint2014arXiv

Baryon Acoustic Oscillations in the Lyα forest of BOSS DR11 quasars

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

preprint2014arXiv

Characterizing unknown systematics in large scale structure surveys

Photometric large scale structure (LSS) surveys probe the largest volumes in the Universe, but are inevitably limited by systematic uncertainties. Imperfect photometric calibration leads to biases in our measurements of the density fields of LSS tracers such as galaxies and quasars, and as a result in cosmological parameter estimation. Earlier studies have proposed using cross-correlations between different redshift slices or cross-correlations between different surveys to reduce the effects of such systematics. In this paper we develop a method to characterize unknown systematics. We demonstrate that while we do not have sufficient information to correct for unknown systematics in the data, we can obtain an estimate of their magnitude. We define a parameter to estimate contamination from unknown systematics using cross-correlations between different redshift slices and propose discarding bins in the angular power spectrum that lie outside a certain contamination tolerance level. We show that this method improves estimates of the bias using simulated data and further apply it to photometric luminous red galaxies in the Sloan Digital Sky Survey as a case study.

preprint2014arXiv

New white dwarf stars in the Sloan Digital Sky Survey Data Release 10

We report the discovery of 9 089 new spectroscopically confirmed white dwarfs and subdwarfs in the Sloan Digital Sky Survey Data Release 10. We obtain Teff, log g and mass for hydrogen atmosphere white dwarf stars (DAs) and helium atmosphere white dwarf stars (DBs), and estimate the calcium/helium abundances for the white dwarf stars with metallic lines (DZs) and carbon/helium for carbon dominated spectra DQs. We found 1 central star of a planetary nebula, 2 new oxygen spectra on helium atmosphere white dwarfs, 71 DQs, 42 hot DO/PG1159s, 171 white dwarf+main sequence star binaries, 206 magnetic DAHs, 327 continuum dominated DCs, 397 metal polluted white dwarfs, 450 helium dominated white dwarfs, 647 subdwarfs and 6888 new hydrogen dominated white dwarf stars.

preprint2014arXiv

PRIMUS: Effect of Galaxy Environment on the Quiescent Fraction Evolution at z < 0.8

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.

preprint2014arXiv

PRIMUS: Galaxy Clustering as a Function of Luminosity and Color at 0.2<z<1

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.

preprint2014arXiv

Quasar-Lyman $α$ Forest Cross-Correlation from BOSS DR11 : Baryon Acoustic Oscillations

We measure the large-scale cross-correlation of quasars with the Lyman alpha forest absorption, using over 164,000 quasars from Data Release 11 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We extend the previous study of roughly 60,000 quasars from Data Release 9 to larger separations, allowing a measurement of the Baryonic Acoustic Oscillation (BAO) scale along the line of sight $c/(H(z=2.36) ~ r_s) = 9.0 \pm 0.3$ and across the line of sight $D_A(z=2.36) / ~ r_s = 10.8 \pm 0.4$, consistent with CMB and other BAO data. Using the best fit value of the sound horizon from Planck data ($r_s=147.49 Mpc$), we can translate these results to a measurement of the Hubble parameter of $H(z=2.36) = 226 \pm 8 km/s / Mpc$ and of the angular diameter distance of $D_A(z=2.36) = 1590 \pm 60 Mpc$. The measured cross-correlation function and an update of the code to fit the BAO scale (baofit) are made publicly available.

preprint2014arXiv

The Clustering of Galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS): measuring growth rate and geometry with anisotropic clustering

We use the observed anisotropic clustering of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 CMASS sample to measure the linear growth rate of structure, the Hubble expansion rate and the comoving distance scale. Our sample covers 8498 ${\rm deg}^2$ and encloses an effective volume of 6.0 ${\rm Gpc}^3$ at an effective redshift of $\bar{z} = 0.57$. We find $fσ_8 = 0.441 \pm 0.044$, $H = 93.1 \pm 3.0\ {\mathrm{km}\ \mathrm{s}^{-1} \mathrm{Mpc}^{-1}}$ and $D_{\rm A} = 1380 \pm 23\ {\rm Mpc}$ when fitting the growth and expansion rate simultaneously. When we fix the background expansion to the one predicted by spatially-flat $Λ$CDM model in agreement with recent Planck results, we find $fσ_8 = 0.447 \pm 0.028$ (6 per cent accuracy). While our measurements are generally consistent with the predictions of $Λ$CDM and General Relativity, they mildly favor models in which the strength of gravitational interactions is weaker than what is predicted by General Relativity. Combining our measurements with recent cosmic microwave background data results in tight constraints on basic cosmological parameters and deviations from the standard cosmological model. Separately varying these parameters, we find $w = -0.983 \pm 0.075$ (8 per cent accuracy) and $γ= 0.69 \pm 0.11$ (16 per cent accuracy) for the effective equation of state of dark energy and the growth rate index, respectively. Both constraints are in good agreement with the standard model values of $w=-1$ and $γ= 0.554$.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 10 and 11 galaxy samples

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 $Λ$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\,σ$ 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.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: galaxy clustering measurements in the low redshift sample of Data Release 11

We present the distance measurement to z = 0.32 using the 11th data release of the Sloan Digital Sky Survey-III Baryon Acoustic Oscillation Survey (BOSS). We use 313,780 galaxies of the low-redshift (LOWZ) sample over 7,341 square-degrees to compute $D_V = (1264 \pm 25)(r_d/r_{d,fid})$ - a sub 2% measurement - using the baryon acoustic feature measured in the galaxy two-point correlation function and power-spectrum. We compare our results to those obtained in DR10. We study observational systematics in the LOWZ sample and quantify potential effects due to photometric offsets between the northern and southern Galactic caps. We find the sample to be robust to all systematic effects found to impact on the targeting of higher-redshift BOSS galaxies, and that the observed north-south tensions can be explained by either limitations in photometric calibration or by sample variance, and have no impact on our final result. Our measurement, combined with the baryonic acoustic scale at z = 0.57, is used in Anderson et al. (2013a) to constrain cosmological parameters.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: modeling of the luminosity and colour dependence in the Data Release 10

We investigate the luminosity and colour dependence of clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The halo occupation distribution framework is adopted to model the projected two-point correlation function measurements on small and intermediate scales (from $0.02$ to $60\,h^{-1}{\rm {Mpc}}$) and to interpret the observed trends and infer the connection of galaxies to dark matter halos. We find that luminous red galaxies reside in massive halos of mass $M{\sim}10^{13}$--$10^{14}\,h^{-1}{\rm M_\odot}$ and more luminous galaxies are more clustered and hosted by more massive halos. The strong small-scale clustering requires a fraction of these galaxies to be satellites in massive halos, with the fraction at the level of 5--8 per cent and decreasing with luminosity. The characteristic mass of a halo hosting on average one satellite galaxy above a luminosity threshold is about a factor $8.7$ larger than that of a halo hosting a central galaxy above the same threshold. At a fixed luminosity, progressively redder galaxies are more strongly clustered on small scales, which can be explained by having a larger fraction of these galaxies in the form of satellites in massive halos. Our clustering measurements on scales below $0.4\,h^{-1}{\rm {Mpc}}$ allow us to study the small-scale spatial distribution of satellites inside halos. While the clustering of luminosity-threshold samples can be well described by a Navarro-Frenk-White (NFW) profile, that of the reddest galaxies prefers a steeper or more concentrated profile. Finally, we also use galaxy samples of constant number density at different redshifts to study the evolution of luminous galaxies, and find the clustering to be consistent with passive evolution in the redshift range of $0.5 \lesssim z \lesssim 0.6$.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Testing gravity with redshift-space distortions using the power spectrum multipoles

We analyse the anisotropic clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS Data Release 11 (DR11) sample, which consists of $690\,827$ galaxies in the redshift range $0.43 < z < 0.7$ and has a sky coverage of $8\,498\,\text{deg}^2$. We perform our analysis in Fourier space using a power spectrum estimator suggested by Yamamoto et al. (2006). We measure the multipole power spectra in a self-consistent manner for the first time in the sense that we provide a proper way to treat the survey window function and the integral constraint, without the commonly used assumption of an isotropic power spectrum and without the need to split the survey into sub-regions. The main cosmological signals exploited in our analysis are the Baryon Acoustic Oscillations and the signal of redshift space distortions, both of which are distorted by the Alcock-Paczynski effect. Together, these signals allow us to constrain the distance ratio $D_V(z_{\rm eff})/r_s(z_d) = 13.89\pm 0.18$, the Alcock-Paczynski parameter $F_{\rm AP}(z_{\rm eff}) = 0.679\pm0.031$ and the growth rate of structure $f(z_{\rm eff})σ_8(z_{\rm eff}) = 0.419\pm0.044$ at the effective redshift $z_{\rm eff}=0.57$. We did not find significant systematic uncertainties for $D_V/r_s$ or $F_{\rm AP}$ but include a systematic error for $fσ_8$ of $3.1\%$. Combining our dataset with Planck to test General Relativity (GR) through the simple $γ$-parameterisation, reveals a $\sim 2σ$ tension between the data and the prediction by GR. The tension between our result and GR can be traced back to a tension in the clustering amplitude $σ_8$ between CMASS and Planck.

preprint2014arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Technical Overview

The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul dark/grey time, with an average cadence of ~4 days, totaling more than 30 epochs. Supporting photometric monitoring in the g and i bands was conducted at multiple facilities including the CFHT and the Steward Observatory Bok telescopes in 2014, with a cadence of ~2 days and covering all lunar phases. The RM field (RA, DEC=14:14:49.00, +53:05:00.0) lies within the CFHT-LS W3 field, and coincides with the Pan-STARRS 1 (PS1) Medium Deep Field MD07, with three prior years of multi-band PS1 light curves. The SDSS-RM 6-month baseline program aims to detect time lags between the quasar continuum and broad line region (BLR) variability on timescales of up to several months (in the observed frame) for ~10% of the sample, and to anchor the time baseline for continued monitoring in the future to detect lags on longer timescales and at higher redshift. SDSS-RM is the first major program to systematically explore the potential of RM for broad-line quasars at z>0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science.

preprint2014arXiv

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

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.

preprint2013arXiv

Baryon Acoustic Oscillations in the Ly-α forest of BOSS quasars

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

preprint2013arXiv

Detection of Lyβauto-correlations and Lyα-Lyβcross-correlations in BOSS Data Release 9

The Lyman-$β$ forest refers to a region in the spectra of distant quasars that lies between the rest-frame Lyman-$β$ and Lyman-$γ$ emissions. The forest in this region is dominated by a combination of absorption due to resonant Ly$α$ and Ly$β$ scattering. When considering the 1D Ly$β$ forest in addition to the 1D Ly$α$ forest, the full statistical description of the data requires four 1D power spectra: Ly$α$ and Ly$β$ auto-power spectra and the Ly$α$-Ly$β$ real and imaginary cross-power spectra. We describe how these can be measured using an optimal quadratic estimator that naturally disentangles Ly$α$ and Ly$β$ contributions. Using a sample of approximately 60,000 quasar sight-lines from the BOSS Data Release 9, we make the measurement of the one-dimensional power spectrum of fluctuations due to the Ly$β$ resonant scattering. While we have not corrected our measurements for resolution damping of the power and other systematic effects carefully enough to use them for cosmological constraints, we can robustly conclude the following: i) Ly$β$ power spectrum and Ly$α$-Ly$β$ cross spectra are detected with high statistical significance; ii) the cross-correlation coefficient is $\approx 1$ on large scales; iii) the Ly$β$ measurements are contaminated by the associated OVI absorption, which is analogous to the SiIII contamination of the Ly$α$ forest. Measurements of the Ly$β$ forest will allow extension of the usable path-length for the Ly$α$ measurements while allowing a better understanding of the physics of intergalactic medium and thus more robust cosmological constraints.

preprint2013arXiv

Discovery of a Dynamical Cold Point in the Heart of the Sagittarius dSph Galaxy with Observations from the APOGEE Project

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

preprint2013arXiv

Measuring D_A and H at z=0.35 from the SDSS DR7 LRGs using baryon acoustic oscillations

We present measurements of the angular diameter distance D_A(z) and the Hubble parameter H(z) at z=0.35 using the anisotropy of the baryon acoustic oscillation (BAO) signal measured in the galaxy clustering distribution of the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) Luminous Red Galaxies (LRG) sample. Our work is the first to apply density-field reconstruction to an anisotropic analysis of the acoustic peak. Reconstruction partially removes the effects of non-linear evolution and redshift-space distortions in order to sharpen the acoustic signal. We present the theoretical framework behind the anisotropic BAO signal and give a detailed account of the fitting model we use to extract this signal from the data. Our method focuses only on the acoustic peak anisotropy, rather than the more model-dependent anisotropic information from the broadband power. We test the robustness of our analysis methods on 160 LasDamas DR7 mock catalogues and find that our models are unbiased at the ~0.2% level in measuring the BAO anisotropy. After reconstruction we measure D_A(z=0.35)=1050+/-38 Mpc and H(z=0.35)=84.4+/-7.0 km/s/Mpc assuming a sound horizon of r_s=152.76 Mpc. Note that these measurements are correlated with a correlation coefficient of 0.58. This represents a factor of 1.4 improvement in the error on D_A relative to the pre-reconstruction case; a factor of 1.2 improvement is seen for H.

preprint2013arXiv

Observational Probes of Cosmic Acceleration

The accelerating expansion of the universe is the most surprising cosmological discovery in many decades, implying that the universe is dominated by some form of "dark energy" with exotic physical properties, or that Einstein's theory of gravity breaks down on cosmological scales. The profound implications of cosmic acceleration have inspired ambitious experimental efforts to measure the history of expansion and growth of structure with percent-level precision or higher. We review in detail the four most well established methods for making such measurements: Type Ia supernovae, baryon acoustic oscillations (BAO), weak gravitational lensing, and galaxy clusters. We pay particular attention to the systematic uncertainties in these techniques and to strategies for controlling them at the level needed to exploit "Stage IV" dark energy facilities such as BigBOSS, LSST, Euclid, and WFIRST. We briefly review a number of other approaches including redshift-space distortions, the Alcock-Paczynski test, and direct measurements of H_0. We present extensive forecasts for constraints on the dark energy equation of state and parameterized deviations from GR, achievable with Stage III and Stage IV experimental programs that incorporate supernovae, BAO, weak lensing, and CMB data. We also show the level of precision required for other methods to provide constraints competitive with those of these fiducial programs. We emphasize the value of a balanced program that employs several of the most powerful methods in combination, both to cross-check systematic uncertainties and to take advantage of complementary information. Surveys to probe cosmic acceleration produce data sets with broad applications, and they continue the longstanding astronomical tradition of mapping the universe in ever greater detail over ever larger scales.

preprint2013arXiv

PRIMUS: An observationally motivated model to connect the evolution of the AGN and galaxy populations out to z~1

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.

preprint2013arXiv

PRIMUS: Constraints on Star Formation Quenching and Galaxy Merging, and the Evolution of the Stellar Mass Function From z=0-1

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.

preprint2013arXiv

PRIMUS: Infrared and X-ray AGN Selection Techniques at 0.2<z<1.2

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.

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological constraints from the full shape of the clustering wedges

We explore the cosmological implications of the clustering wedges, xi_perp(s) and xi_para(s), of the CMASS Data Release 9 (DR9) sample of the Baryon Oscillation Spectroscopic Survey (BOSS). These clustering wedges are defined by averaging the full two-dimensional correlation function, xi(mu,s), over the ranges 0<mu<0.5 and 0.5<mu<1, respectively. These measurements allow us to constrain the parameter combinations D_A(z)/r_s(z_d)=9.03 +- 0.21 and cz/(r_s(z_d)H(z)) = 12.14 +- 0.43 at the mean redsfhit of the sample, z=0.57. We combine the information from the clustering wedges with recent measurements of CMB, BAO and type Ia supernovae to obtain constraints on the cosmological parameters of the standard LCDM model and a number of potential extensions. The information encoded in the clustering wedges is most useful when the dark energy equation of state is allowed to deviate from its standard LCDM value. The combination of all datasets shows no evidence of a deviation from a constant dark energy equation of state, in which case we find w_DE = -1.013 +- 0.064, in complete agreement with a cosmological constant. We explore potential deviations from general relativity by constraining the growth rate f(z)=d ln D(a)/ d ln a, in which case the combination of the CMASS clustering wedges with CMB data implies f(z=0.57)=0.719 +- 0.094, in accordance with the predictions of GR. Our results clearly illustrate the additional constraining power of anisotropic clustering measurements with respect to that of angle-averaged quantities.

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological implications of the full shape of the clustering wedges in the data release 10 and 11 galaxy samples

We explore the cosmological implications of the angle-averaged correlation function, xi(s), and the clustering wedges, xi_perp(s) and xi_para(s), of the LOWZ and CMASS galaxy samples from Data Release 10 and 11 of the SDSS-III Baryon Oscillation Spectroscopic Survey. Our results show no significant evidence for a deviation from the standard LCDM model. The combination of the information from our clustering measurements with recent data from the cosmic microwave background is sufficient to constrain the curvature of the Universe to Omega_k = 0.0010 +- 0.0029, the total neutrino mass to Sum m_nu < 0.23 eV (95% confidence level), the effective number of relativistic species to N_eff=3.31 +- 0.27, and the dark energy equation of state to w_DE = -1.051 +- 0.076. These limits are further improved by adding information from type Ia supernovae and baryon acoustic oscillations from other samples. In particular, this data set combination is completely consistent with a time-independent dark energy equation of state, in which case we find w_DE=-1.024 +- 0.052. We explore the constraints on the growth-rate of cosmic structures assuming f(z)=Omega_m(z)^gamma and obtain gamma=0.69 +- 0.15, in agreement with the predictions from general relativity of gamma=0.55.

preprint2013arXiv

The Clustering of Galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Including covariance matrix errors

We present improved methodology for including covariance matrices in the error budget of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy clustering measurements, revisiting Data Release 9 (DR9) analyses, and describing a method that is used in DR10/11 analyses presented in companion papers. The precise analysis method adopted is becoming increasingly important, due to the precision that BOSS can now reach: even using as many as 600 mock catalogues to estimate covariance of 2-point clustering measurements can still lead to an increase in the errors of ~20%, depending on how the cosmological parameters of interest are measured. In this paper we extend previous work on this contribution to the error budget, deriving formulae for errors measured by integrating over the likelihood, and to the distribution of recovered best-fit parameters fitting the simulations also used to estimate the covariance matrix. Both are situations that previous analyses of BOSS have considered. We apply the formulae derived to Baryon Acoustic Oscillation (BAO) and Redshift-Space Distortion (RSD) measurements from BOSS in our companion papers. To further aid these analyses, we consider the optimum number of bins to use for 2-point measurements using the monopole power spectrum or correlation function for BAO, and the monopole and quadrupole moments of the correlation function for anisotropic-BAO and RSD measurements.

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Luminosity and Color Dependence and Redshift Evolution

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

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Measuring D_A and H at z=0.57 from the Baryon Acoustic Peak in the Data Release 9 Spectroscopic Galaxy Sample

We present measurements of the angular diameter distance to and Hubble parameter at z=0.57 from the measurement of the baryon acoustic peak in the correlation of galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey. Our analysis is based on a sample from Data Release 9 of 264,283 galaxies over 3275 square degrees in the redshift range 0.43<z<0.70. We use two different methods to provide robust measurement of the acoustic peak position across and along the line of sight in order to measure the cosmological distance scale. We find D_A(0.57) = 1408 +/- 45 Mpc and H(0.57) = 92.9 +/- 7.8 km/s/Mpc for our fiducial value of the sound horizon. These results from the anisotropic fitting are fully consistent with the analysis of the spherically averaged acoustic peak position presented in Anderson et al, 2012. Our distance measurements are a close match to the predictions of the standard cosmological model featuring a cosmological constant and zero spatial curvature.

preprint2013arXiv

The Clustering of Galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Measuring H(z) and D_A(z) at z = 0.57 with Clustering Wedges

We analyze the 2D correlation function of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of massive galaxies of the ninth data release to measure cosmic expansion H and the angular diameter distance D_A at a mean redshift of <z> = 0.57. We apply, for the first time, a new correlation function technique called clustering wedges. Using a physically motivated model, the anisotropic baryonic acoustic feature in the galaxy sample is detected at a significance level of 4.7 sigma compared to a featureless model. The baryonic acoustic feature is used to obtain model independent constraints cz/H/r_s = 12.28 +- 0.82 (6.7 per-cent accuracy) and D_A/r_s = 9.05 +- 0.27 (3.0 per-cent) with a correlation coefficient of -0.5, where r_s is the sound horizon scale at the end of the baryonic drag era. We conduct thorough tests on the data and 600 simulated realizations, finding robustness of the results regardless of the details of the analysis method. Combining with r_s constraints from the Cosmic Microwave Background we obtain H(0.57) = 90.8 +- 6.2 kms-1Mpc-1 and D_A(0.57) = 1386 +- 45 Mpc. We use simulations to forecast results of the final BOSS CMASS data set. We apply the reconstruction technique on the simulations demonstrating that the sharpening of the anisotropic baryonic acoustic feature should improve the detection as well as tighten constraints of H and D_A by 30 per-cent on average.

preprint2013arXiv

The Clustering of Galaxies in the SDSS-III DR10 Baryon Oscillation Spectroscopic Survey: No Detectable Colour Dependence of Distance Scale or Growth Rate Measurements

We study the clustering of galaxies, as a function of their colour, from Data Release Ten (DR10) of the SDSS-III Baryon Oscillation Spectroscopic Survey. We select 122,967 galaxies with 0.43 < z < 0.7 into a "Blue" sample and 131,969 into a "Red" sample based on k+e corrected (to z=0.55) r-i colours and i band magnitudes. The samples are chosen to each contain more than 100,000 galaxies, have similar redshift distributions, and maximize the difference in clustering amplitude. The Red sample has a 40% larger bias than the Blue (b_Red/b_Blue = 1.39+-0.04), implying the Red galaxies occupy dark matter halos with an average mass that is 0.5 log Mo greater. Spherically averaged measurements of the correlation function, ξ0, and the power spectrum are used to locate the position of the baryon acoustic oscillation (BAO) feature of both samples. Using ξ0, we obtain distance scales, relative to our reference LCDM cosmology, of 1.010+-0.027 for the Red sample and 1.005+-0.031 for the Blue. After applying reconstruction, these measurements improve to 1.013+/-0.020 for the Red sample and 1.008+-0.026 for the Blue. For each sample, measurements of ξ0 and the second multipole moment, ξ2, of the anisotropic correlation function are used to determine the rate of structure growth, parameterized by fσ8. We find fσ8,Red = 0.511+-0.083, fσ8,Blue = 0.509+/-0.085, and fσ8,Cross = 0.423+-0.061 (from the cross-correlation between the Red and Blue samples). We use the covariance between the bias and growth measurements obtained from each sample and their cross-correlation to produce an optimally-combined measurement of fσ8,comb = 0.443+-0.055. In no instance do we detect significant differences in distance scale or structure growth measurements obtained from the Blue and Red samples.

preprint2013arXiv

The PRIsm MUlti-object Survey (PRIMUS). II. Data Reduction and Redshift Fitting

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.

preprint2013arXiv

Very Metal-poor Stars in the Outer Galactic Bulge Found by the Apogee Survey

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

preprint2012arXiv

A 2% Distance to z = 0.35 by Reconstructing Baryon Acoustic Oscillations - III : Cosmological Measurements and Interpretation

We use the 2% distance measurement from our reconstructed baryon acoustic oscillations (BAOs) signature using the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) Luminous Red Galaxies (LRGs) from Padmanabhan et al. (2012) and Xu et al. (2012) combined with cosmic microwave background (CMB) data from Wilkinson Microwave Anisotropy Probe (WMAP7) to measure parameters for various cosmological models. We find a 1.7% measurement of H_0 = 69.8 +/- 1.2 km/s/Mpc and a 5.0% measurement of Omega_m = 0.280 +/- 0.014 for a flat Universe with a cosmological constant. These measurements of H_0 and Omega_m are robust against a range of underlying models for the expansion history. We measure the dark energy equation of state parameter w = -0.97 +/- 0.17, which is consistent with a cosmological constant. If curvature is allowed to vary, we find that the Universe is consistent with a flat geometry (Omega_K = -0.004 +/- 0.005). We also use a combination of the 6 Degree Field Galaxy Survey BAO data, WiggleZ Dark Energy Survey data, Type Ia supernovae (SN) data, and a local measurement of the Hubble constant to explore cosmological models with more parameters. Finally, we explore the effect of varying the energy density of relativistic particles on the measurement of H_0.

preprint2012arXiv

A 2% Distance to z=0.35 by Reconstructing Baryon Acoustic Oscillations - I : Methods and Application to the Sloan Digital Sky Survey

We apply the reconstruction technique to the clustering of galaxies from the SDSS DR7 LRG sample, sharpening the baryon acoustic oscillation (BAO) feature and achieving a 1.9% measurement of the distance to z=0.35. This is the first application of reconstruction of the BAO feature in a galaxy redshift survey. We update the reconstruction algorithm of Eisenstein et al, 2007 to account for the effects of survey geometry as well as redshift-space distortions and validate it on 160 LasDamas simulations. We demonstrate that reconstruction sharpens the BAO feature in the angle averaged galaxy correlation function, reducing the nonlinear smoothing scale Σ_nl from 8.1 Mpc/h to 4.4 Mpc/h. Reconstruction also significantly reduces the effects of redshift-space distortions at the BAO scale, isotropizing the correlation function. This sharpened BAO feature yields an unbiased distance estimate (< 0.2%) and reduces the scatter from 3.3% to 2.1%. We demonstrate the robustness of these results to the various reconstruction parameters, including the smoothing scale, the galaxy bias and the linear growth rate. Applying this reconstruction algorithm to the SDSS LRG DR7 sample improves the significance of the BAO feature in these data from 3.3 sigma for the unreconstructed correlation function, to 4.2 sigma after reconstruction. We estimate a relative distance scale D_V/r_s to z=0.35 of 8.88+/-0.17, where r_s is the sound horizon and D_V = (D_A^2/H)^{1/3} is a combination of the angular diameter distance D_A and Hubble parameter H. Assuming a sound horizon of 154.25 Mpc, this translates into a distance measurement D_V (z=0.35) = 1.356+/-0.025 Gpc. We find that reconstruction reduces the distance error in the DR7 sample from 3.5% to 1.9%, equivalent to a survey with three times the volume of SDSS.

preprint2012arXiv

A 2% Distance to z=0.35 by Reconstructing Baryon Acoustic Oscillations - II: Fitting Techniques

We present results from fitting the baryon acoustic oscillation (BAO) signal in the correlation function obtained from the first application of reconstruction to a galaxy redshift survey, namely, the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) luminous red galaxy (LRG) catalogue. We also introduce more careful approaches for deriving a suitable covariance matrix and fitting model for galaxy correlation functions. These all aid in obtaining a more accurate measurement of the acoustic scale and its error. We validate our reconstruction, covariance matrix and fitting techniques on 160 mock catalogues derived from the LasDamas simulations in real and redshift space. We then apply these techniques to the DR7 LRG sample and find that the error on the acoustic scale decreases from ~3.5% before reconstruction to ~1.9% after reconstruction. This factor of 1.8 reduction in the error is equivalent to the effect of tripling the survey volume. We also see an increase in our BAO detection confidence from ~3-sigma to ~4-sigma after reconstruction with our confidence level in measuring the correct acoustic scale increasing from ~3-sigma to ~5-sigma. Using the mean of the acoustic scale probability distributions produced from our fits, we find D_v/r_s = 8.89 +/- 0.31 before reconstruction and 8.88 +/- 0.17 after reconstruction.

preprint2012arXiv

Clustering of Sloan Digital Sky Survey III Photometric Luminous Galaxies: The Measurement, Systematics and Cosmological Implications

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 Ω_Λ= 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 $Ω_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)

preprint2012arXiv

Evidence of Galaxy Cluster Motions with the Kinematic Sunyaev-Zel'dovich Effect

Using high-resolution microwave sky maps made by the Atacama Cosmology Telescope, we for the first time present strong evidence for motions of galaxy clusters and groups via microwave background temperature distortions due to the kinematic Sunyaev-Zel'dovich effect. Galaxy clusters are identified by their constituent luminous galaxies observed by the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. We measure the mean pairwise momentum of clusters, with a probability of the signal being due to random errors of 0.002, and the signal is consistent with the growth of cosmic structure in the standard model of cosmology.

preprint2012arXiv

PRIMUS: The dependence of AGN accretion on host stellar mass and color

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.

preprint2012arXiv

Spectral Classification and Redshift Measurement for the SDSS-III Baryon Oscillation Spectroscopic Survey

(abridged) We describe the automated spectral classification, redshift determination, and parameter measurement pipeline in use for the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III (SDSS-III) as of Data Release 9, encompassing 831,000 moderate-resolution optical spectra. We give a review of the algorithms employed, and describe the changes to the pipeline that have been implemented for BOSS relative to previous SDSS-I/II versions, including new sets of stellar, galaxy, and quasar redshift templates. For the color-selected CMASS sample of massive galaxies at redshift 0.4 <~ z <~ 0.8 targeted by BOSS for the purposes of large-scale cosmological measurements, the pipeline achieves an automated classification success rate of 98.7% and confirms 95.4% of unique CMASS targets as galaxies (with the balance being mostly M stars). Based on visual inspections of a subset of BOSS galaxies, we find that ~0.2% of confidently reported CMASS sample classifications and redshifts are incorrect, and ~0.4% of all CMASS spectra are objects unclassified by the current algorithm which are potentially recoverable. The BOSS pipeline confirms that ~51.5% of the quasar targets have quasar spectra, with the balance mainly consisting of stars. Statistical (as opposed to systematic) redshift errors propagated from photon noise are typically a few tens of km/s for both galaxies and quasars, with a significant tail to a few hundreds of km/s for quasars. We test the accuracy of these statistical redshift error estimates using repeat observations, finding them underestimated by a factor of 1.19 to 1.34 for galaxies, and by a factor of 2 for quasars. We assess the impact of sky-subtraction quality, S/N, and other factors on galaxy redshift success. Finally, we document known issues, and describe directions of ongoing development.

preprint2012arXiv

The Apache Point Observatory Galactic Evolution Experiment: First Detection of High Velocity Milky Way Bar Stars

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

preprint2012arXiv

The Baryon Oscillation Spectroscopic Survey of SDSS-III

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.

preprint2012arXiv

The BOSS Emission-Line Lens Survey. II. Investigating Mass-Density Profile Evolution in the SLACS+BELLS Strong Gravitational Lens Sample

We present an analysis of the evolution of the central mass-density profile of massive elliptical galaxies from the SLACS and BELLS strong gravitational lens samples over the redshift interval z ~ 0.1-0.6, based on the combination of strong-lensing aperture mass and stellar velocity-dispersion constraints. We find a significant trend towards steeper mass profiles (parameterized by the power-law density model with rho ~ r^[-gamma]) at later cosmic times, with magnitude d<gamma>/dz = -0.60 +/- 0.15. We show that the combined lens-galaxy sample is consistent with a non-evolving distribution of stellar velocity dispersions. Considering possible additional dependence of <gamma> on lens-galaxy stellar mass, effective radius, and Sersic index, we find marginal evidence for shallower mass profiles at higher masses and larger sizes, but with a significance that is sub-dominant to the redshift dependence. Using the results of published Monte Carlo simulations of spectroscopic lens surveys, we verify that our mass-profile evolution result cannot be explained by lensing selection biases as a function of redshift. Interpreted as a true evolutionary signal, our result suggests that major dry mergers involving off-axis trajectories play a significant role in the evolution of the average mass-density structure of massive early-type galaxies over the past 6 Gyr. We also consider an alternative non-evolutionary hypothesis based on variations in the strong-lensing measurement aperture with redshift, which would imply the detection of an "inflection zone" marking the transition between the baryon-dominated and dark-matter halo-dominated regions of the lens galaxies. Further observations of the combined SLACS+BELLS sample can constrain this picture more precisely, and enable a more detailed investigation of the multivariate dependences of galaxy mass structure across cosmic time.

preprint2012arXiv

The BOSS Lyman-alpha Forest Sample from SDSS Data Release 9

We present the BOSS Lyman-alpha (Lya) Forest Sample from SDSS Data Release 9, comprising 54,468 quasar spectra with zqso > 2.15 suitable for Lya forest analysis. This data set probes the intergalactic medium with absorption redshifts 2.0 < z_alpha < 5.7 over an area of 3275 square degrees, and encompasses an approximate comoving volume of 20 h^-3 Gpc^3. With each spectrum, we have included several products designed to aid in Lya forest analysis: improved sky masks that flag pixels where data may be unreliable, corrections for known biases in the pipeline estimated noise, masks for the cores of damped Lya systems and corrections for their wings, and estimates of the unabsorbed continua so that the observed flux can be converted to a fractional transmission. The continua are derived using a principal component fit to the quasar spectrum redwards of restframe Lya (lambda > 1216 Ang), extrapolated into the forest region and normalized by a linear function to fit the expected evolution of the Lya forest mean-flux. The estimated continuum errors are ~5% rms. We also discuss possible systematics arising from uncertain spectrophotometry and artifacts in the flux calibration; global corrections for the latter are provided. Our sample provides a convenient starting point for users to analyze clustering in BOSS Lya forest data, and it provides a fiducial data set that can be used to compare results from different analyses of baryon acoustic oscillations in the Lya forest. The full data set is available from the SDSS-III DR9 web site.

preprint2012arXiv

The Clustering of Galaxies in SDSS-III DR9 Baryon Oscillation Spectroscopic Survey: Constraints on Primordial Non-Gaussianity

We analyze the density field of 264,283 galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) and included in the SDSS data release nine (DR9). In total, the SDSS DR9 BOSS data includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of more than 3,000 deg^2. We measure the power spectrum of these galaxies with redshifts 0.43 < z < 0.7 in order to constrain the amount of local non-Gaussianity, f_NL,local, in the primordial density field, paying particular attention to the impact of systematic uncertainties. The BOSS galaxy density field is systematically affected by the local stellar density and this influences the ability to accurately measure f_NL,local. In the absence of any correction, we find (erroneously) that the probability that f_NL,local is greater than zero, P(f_NL,local >0), is 99.5%. After quantifying and correcting for the systematic bias and including the added uncertainty, we find -45 < f_NL,local < 195 at 95% confidence, and P(f_NL,local >0) = 91.0%. A more conservative approach assumes that we have only learned the k-dependence of the systematic bias and allows any amplitude for the systematic correction; we find that the systematic effect is not fully degenerate with that of f_NL,local, and we determine that -82 < f_NL,local < 178 (at 95% confidence) and P(f_NL,local >0) = 68%. This analysis demonstrates the importance of accounting for the impact of Galactic foregrounds on f_NL,local measurements. We outline the methods that account for these systematic biases and uncertainties. We expect our methods to yield robust constraints on f_NL,local for both our own and future large-scale-structure investigations.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Analysis of potential systematics

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)

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 9 Spectroscopic Galaxy Sample

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σin both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7σ. 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.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: measurements of the growth of structure and expansion rate at z=0.57 from anisotropic clustering

We analyze the anisotropic clustering of massive galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) sample, which consists of 264,283 galaxies in the redshift range 0.43 < z < 0.7 spanning 3,275 square degrees. Both peculiar velocities and errors in the assumed redshift-distance relation ("Alcock-Paczynski effect") generate correlations between clustering amplitude and orientation with respect to the line-of-sight. Together with the sharp baryon acoustic oscillation (BAO) standard ruler, our measurements of the broadband shape of the monopole and quadrupole correlation functions simultaneously constrain the comoving angular diameter distance (2190 +/- 61 Mpc) to z=0.57, the Hubble expansion rate at z=0.57 (92.4 +/- 4.5 km/s/Mpc), and the growth rate of structure at that same redshift (d sigma8/d ln a = 0.43 +/- 0.069). Our analysis provides the best current direct determination of both DA and H in galaxy clustering data using this technique. If we further assume a LCDM expansion history, our growth constraint tightens to d sigma8/d ln a = 0.415 +/- 0.034. In combination with the cosmic microwave background, our measurements of DA, H, and growth all separately require dark energy at z > 0.57, and when combined imply Ω_Λ = 0.74 +/- 0.016, independent of the Universe's evolution at z<0.57. In our companion paper (Samushia et al. prep), we explore further cosmological implications of these observations.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: the low redshift sample

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.

preprint2012arXiv

The Galaxy Optical Luminosity Function from the AGN and Galaxy Evolution Survey (AGES)

We present the galaxy optical luminosity function for the redshift range 0.05<z<0.75 from the AGN and Galaxy Evolution Survey (AGES), a spectroscopic survey of 7.6 sq. deg. in the Bootes field of the NOAO Deep Wide-Field Survey. Our statistical sample is comprised of 12,473 galaxies with known redshifts down to I=20.4 (AB). Our results at low redshift are consistent with those from SDSS; at higher redshift, we find strong evidence for evolution in the luminosity function, including differential evolution between blue and red galaxies. We find that the luminosity density evolves as (1+z)^(0.54+/-0.64) for red galaxies and (1+z)^(1.64+/-0.39) for blue galaxies.

preprint2012arXiv

The Milky Way's circular velocity curve between 4 and 14 kpc from APOGEE data

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

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

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.

preprint2012arXiv

The Sloan Digital Sky Survey quasar catalog: ninth data release

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

preprint2011arXiv

Evolution of the Stellar Mass-Metallicity Relation Since z=0.75

We measure the gas-phase oxygen abundances of ~3000 star-forming galaxies at z=0.05-0.75 using optical spectrophotometry from the AGN and Galaxy Evolution Survey (AGES), a spectroscopic survey of I_AB<20.45 galaxies over 7.9 deg^2 in the NOAO Deep Wide Field Survey (NDWFS) Bootes field. We use state-of-the-art techniques to measure the nebular emission lines and stellar masses, and explore and quantify several potential sources of systematic error, including the choice of metallicity diagnostic, aperture bias, and contamination from unidentified active galactic nuclei (AGN). Combining volume-limited AGES samples in six independent redshift bins and ~75,000 star-forming galaxies with r_AB<17.6 at z=0.05-0.2 selected from the Sloan Digital Sky Survey (SDSS) that we analyze in the identical manner, we measure the evolution of the stellar mass-metallicity (M-Z) between z=0.05 and z=0.75. We find that at fixed stellar mass galaxies at z~0.7 have just 30%-60% the metal content of galaxies at the present epoch, where the uncertainty is dominated by the strong-line method used to measure the metallicity. Moreover, we find no statistically significant evidence that the M-Z relation evolves in a mass-dependent way for M=10^9.8-10^11 Msun star-forming galaxies. Thus, for this range of redshifts and stellar masses the M-Z relation simply shifts toward lower metallicity with increasing redshift without changing its shape.

preprint2011arXiv

Galaxy Bias and its Effects on the Baryon Acoustic Oscillations Measurements

The baryon acoustic oscillation (BAO) feature in the clustering of matter in the universe serves as a robust standard ruler and hence can be used to map the expansion history of the universe. We use high force resolution simulations to analyze the effects of galaxy bias on the measurements of the BAO signal. We apply a variety of Halo Occupation Distributions (HODs) and produce biased mass tracers to mimic different galaxy populations. We investigate whether galaxy bias changes the non-linear shifts on the acoustic scale relative to the underlying dark matter distribution presented by Seo et al (2009). For the less biased HOD models (b < 3), we do not detect any shift in the acoustic scale relative to the no-bias case, typically 0.10% \pm 0.10%. However, the most biased HOD models (b > 3) show a shift at moderate significance (0.79% \pm 0.31% for the most extreme case). We test the one-step reconstruction technique introduced by Eisenstein et al. (2007) in the case of realistic galaxy bias and shot noise. The reconstruction scheme increases the correlation between the initial and final (z = 1) density fields achieving an equivalent level of correlation at nearly twice the wavenumber after reconstruction. Reconstruction reduces the shifts and errors on the shifts. We find that after reconstruction the shifts from the galaxy cases and the dark matter case are consistent with each other and with no shift. The 1-sigma systematic errors on the distance measurements inferred from our BAO measurements with various HODs after reconstruction are about 0.07% - 0.15%.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

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)

preprint2011arXiv

The BOSS Emission-Line Lens Survey (BELLS). I. A large spectroscopically selected sample of Lens Galaxies at redshift ~ 0.5

We present a catalog of 25 definite and 11 probable strong galaxy-galaxy gravitational lens systems with redshifts 0.4 \lesssim z \lesssim 0.7, discovered spectroscopically by the presence of higher redshift emission-lines within the Baryon Oscillation Spectroscopic Survey (BOSS) of luminous galaxies, and confirmed with high-resolution Hubble Space Telescope (HST) images of 44 candidates. Our survey extends the methodology of the Sloan Lens ACS Survey (SLACS: Bolton et al. 2006; 2008) to higher redshift. We describe the details of the BOSS spectroscopic candidate detections, our HST Adanced Camera for Surveys (ACS) image processing and analysis methods, and our strong gravitational lens modeling procedure. We report BOSS spectroscopic parameters and ACS photometric parameters for all candidates, and mass-distribution parameters for the best-fit singular isothermal ellipsoid models of definite lenses. Our sample to date was selected using only the first six months of BOSS survey-quality spectroscopic data. The full 5-year BOSS database should produce a sample of several hundred strong galaxy-galaxy lenses and in combination with SLACS lenses at lower redshift, strongly constrain the redshift evolution of the structure of elliptical, bulge-dominated galaxies as a function of luminosity, stellar mass, and rest-frame color, thereby providing a powerful test for competing theories of galaxy formation and evolution.

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

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.

preprint2011arXiv

The PRIsm MUlti-Object Survey (PRIMUS) I: Survey Overview and Characteristics

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.

preprint2011arXiv

The SDSS-III Baryon Oscillation Spectroscopic Survey: Quasar Target Selection for Data Release Nine

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]

preprint2010arXiv

High-precision predictions for the acoustic scale in the non-linear regime

We measure shifts of the acoustic scale due to nonlinear growth and redshift distortions to a high precision using a very large volume of high-force-resolution simulations. We compare results from various sets of simulations that differ in their force, volume, and mass resolution. We find a consistency within 1.5-sigma for shift values from different simulations and derive shift alpha(z) -1 = (0.300\pm 0.015)% [D(z)/D(0)]^{2} using our fiducial set. We find a strong correlation with a non-unity slope between shifts in real space and in redshift space and a weak correlation between the initial redshift and low redshift. Density-field reconstruction not only removes the mean shifts and reduces errors on the mean, but also tightens the correlations: after reconstruction, we recover a slope of near unity for the correlation between the real and redshift space and restore a strong correlation between the low and the initial redshifts. We derive propagators and mode-coupling terms from our N-body simulations and compared with Zeldovich approximation and the shifts measured from the chi^2 fitting, respectively. We interpret the propagator and the mode-coupling term of a nonlinear density field in the context of an average and a dispersion of its complex Fourier coefficients relative to those of the linear density field; from these two terms, we derive a signal-to-noise ratio of the acoustic peak measurement. We attempt to improve our reconstruction method by implementing 2LPT and iterative operations: we obtain little improvement. The Fisher matrix estimates of uncertainty in the acoustic scale is tested using 5000 (Gpc/h)^3 of cosmological PM simulations from Takahashi et al. (2009). (abridged)

preprint2010arXiv

MARVELS-1b: A Short-Period, Brown Dwarf Desert Candidate from the SDSS-III MARVELS Planet Search

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.

preprint2010arXiv

Morphology and Size Differences between Local & High Redshift Luminous Infrared Galaxies

We show that the star-forming regions in high-redshift luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) and submillimeter galaxies (SMGs) have similar physical scales to those in local normal star-forming galaxies. To first order, their higher infrared (IR) luminosities result from higher luminosity surface density. We also find a good correlation between the IR luminosity and IR luminosity surface density in starburst galaxies across over five orders of magnitude of IR luminosity from local normal galaxies to z ~ 2 SMGs. The intensely star-forming regions of local ULIRGs are significantly smaller than those in their high-redshift counterparts and hence diverge significantly from this correlation, indicating that the ULIRGs found locally are a different population from the high-redshift ULIRGs and SMGs. Based on this relationship, we suggest that luminosity surface density should serve as a more accurate indicator for the IR emitting environment, and hence the observable properties, of star-forming galaxies than their IR luminosity. We demonstrate this approach by showing that ULIRGs at z ~ 1 and a lensed galaxy at z ~ 2.5 exhibit aromatic features agreeing with local LIRGs that are an order of magnitude less luminous, but have similar IR luminosity surface density. A consequence of this relationship is that the aromatic emission strength in star-forming galaxies will appear to increase at z > 1 for a given IR luminosity compared to their local counterparts.

preprint2010arXiv

PRIMUS: Enhanced Specific Star Formation Rates In Close Galaxy Pairs

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.

preprint2010arXiv

PRIMUS: Obscured Star Formation on the Red Sequence

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.

preprint2010arXiv

The Evolution of the Star Formation Rate of Galaxies at 0.0 < z < 1.2

We present the 24 micron rest-frame luminosity function (LF) of star-forming galaxies in the redshift range 0.0 < z < 0.6 constructed from 4047 spectroscopic redshifts from the AGN and Galaxy Evolution Survey of 24 micron selected sources in the Bootes field of the NOAO Deep Wide-Field Survey. This sample provides the best available combination of large area (9 deg^2), depth, and statistically complete spectroscopic observations, allowing us to probe the evolution of the 24 micron LF of galaxies at low and intermediate redshifts while minimizing the effects of cosmic variance. In order to use the observed 24 micron luminosity as a tracer for star formation, active galactic nuclei (AGNs) that could contribute significantly at 24 micron are identified and excluded from our star-forming galaxy sample based on their mid-IR spectral energy distributions or the detection of X-ray emission. The evolution of the 24 micron LF of star-forming galaxies for redshifts of z < 0.65 is consistent with a pure luminosity evolution where the characteristic 24 micron luminosity evolves as (1+z)^(3.8+/-0.3). We extend our evolutionary study to encompass 0.0 < z < 1.2 by combining our data with that of the Far-Infrared Deep Extragalactic Legacy Survey. Over this entire redshift range the evolution of the characteristic 24 micron luminosity is described by a slightly shallower power law of (1+z)^(3.4+/-0.2). We find a local star formation rate density of (1.09+/-0.21) x 10^-2 Msun/yr/Mpc^-3, and that it evolves as (1+z)^(3.5+/-0.2) over 0.0 < z < 1.2. These estimates are in good agreement with the rates using optical and UV fluxes corrected for the effects of intrinsic extinction in the observed sources. This agreement confirms that star formation at z <~ 1.2 is robustly traced by 24 micron observations and that it largely occurs in obscured regions of galaxies. (Abridged)

preprint2009arXiv

Baryon Acoustic Oscillations in the Sloan Digital Sky Survey Data Release 7 Galaxy Sample

The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample represents the final set of galaxies observed using the original SDSS target selection criteria. We analyse the clustering of galaxies within this sample, including both the Luminous Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift Survey (2dFGRS) data. Baryon Acoustic Oscillations are observed in power spectra measured for different slices in redshift; this allows us to constrain the distance--redshift relation at multiple epochs. We achieve a distance measure at redshift z=0.275, of r_s(z_d)/D_V(0.275)=0.1390+/-0.0037 (2.7% accuracy), where r_s(z_d) is the comoving sound horizon at the baryon drag epoch, D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3), D_A(z) is the angular diameter distance and H(z) is the Hubble parameter. We find an almost independent constraint on the ratio of distances D_V(0.35)/D_V(0.2)=1.736+/-0.065, which is consistent at the 1.1sigma level with the best fit Lambda-CDM model obtained when combining our z=0.275 distance constraint with the WMAP 5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, but the discrepancy is now of lower significance, a change caused by a revised error analysis and a change in the methodology adopted, as well as the addition of more data. Using WMAP5 constraints on Omega_bh^2 and Omega_ch^2, and combining our BAO distance measurements with those from the Union Supernova sample, places a tight constraint on Omega_m=0.286+/-0.018 and H_0 = 68.2+/-2.2km/s/Mpc that is robust to allowing curvature and non-Lambda dark energy. This result is independent of the behaviour of dark energy at redshifts greater than those probed by the BAO and supernova measurements. (abridged)

preprint2009arXiv

Cosmological Constraints from the Clustering of the Sloan Digital Sky Survey DR7 Luminous Red Galaxies

We present the power spectrum of the reconstructed halo density field derived from a sample of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Release (DR7). The halo power spectrum has a direct connection to the underlying dark matter power for k <= 0.2 h/Mpc, well into the quasi-linear regime. This enables us to use a factor of ~8 more modes in the cosmological analysis than an analysis with kmax = 0.1 h/Mpc, as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 2006). The observed halo power spectrum for 0.02 < k < 0.2 h/Mpc is well-fit by our model: chi^2 = 39.6 for 40 degrees of freedom for the best fit LCDM model. We find Ω_m h^2 * (n_s/0.96)^0.13 = 0.141^{+0.009}_{-0.012} for a power law primordial power spectrum with spectral index n_s and Ω_b h^2 = 0.02265 fixed, consistent with CMB measurements. The halo power spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to an effective distance to z=0.35: r_s/D_V(0.35) = 0.1097^{+0.0039}_{-0.0042}. Combining the halo power spectrum measurement with the WMAP 5 year results, for the flat LCDM model we find Ω_m = 0.289 +/- 0.019 and H_0 = 69.4 +/- 1.6 km/s/Mpc. Allowing for massive neutrinos in LCDM, we find \sum m_ν < 0.62 eV at the 95% confidence level. If we instead consider the effective number of relativistic species Neff as a free parameter, we find Neff = 4.8^{+1.8}_{-1.7}. Combining also with the Kowalski et al. (2008) supernova sample, we find Ω_{tot} = 1.011 +/- 0.009 and w = -0.99 +/- 0.11 for an open cosmology with constant dark energy equation of state w.

preprint2009arXiv

Halo Occupation Distribution Modeling of Clustering of Luminous Red Galaxies

We perform Halo Occupation Distribution (HOD) modeling to interpret small-scale and intermediate-scale clustering of 35,000 luminous early-type galaxies and their cross-correlation with a reference imaging sample of normal L* galaxies in the Sloan Digital Sky Survey. The modeling results show that most of these luminous red galaxies (LRGs) are central galaxies residing in massive halos of typical mass M ~ a few times 10^13 to 10^14 Msun/h, while a few percent of them have to be satellites within halos in order to produce the strong auto-correlations exhibited on smaller scales. The mean luminosity Lc of central LRGs increases with the host halo mass, with a rough scaling relation of Lc \propto M^0.5. The halo mass required to host on average one satellite LRG above a luminosity threshold is found to be about 10 times higher than that required to host a central LRG above the same threshold. We find that in massive halos the distribution of L* galaxies roughly follows that of the dark matter and their mean occupation number scales with halo mass as M^1.5. The HOD modeling results also allows for an intuitive understanding of the scale-dependent luminosity dependence of the cross-correlation between LRGs and L_* galaxies. Constraints on the LRG HOD provide tests to models of formation and evolution of massive galaxies, and they are also useful for cosmological parameter investigations. In one of the appendices, we provide LRG HOD parameters with dependence on cosmology inferred from modeling the two-point auto-correlation functions of LRGs.

preprint2009arXiv

The Star Formation and Nuclear Accretion Histories of Normal Galaxies in the AGES Survey

We combine IR, optical and X-ray data from the overlapping, 9.3 square degree NOAO Deep Wide-Field Survey (NDWFS), AGN and Galaxy Evolution Survey (AGES), and XBootes Survey to measure the X-ray evolution of 6146 normal galaxies as a function of absolute optical luminosity, redshift, and spectral type over the largely unexplored redshift range 0.1 < z < 0.5. Because only the closest or brightest of the galaxies are individually detected in X-rays, we use a stacking analysis to determine the mean properties of the sample. Our results suggest that X-ray emission from spectroscopically late-type galaxies is dominated by star formation, while that from early-type galaxies is dominated by a combination of hot gas and AGN emission. We find that the mean star formation and supermassive black hole accretion rate densities evolve like (1+z)^3, in agreement with the trends found for samples of bright, individually detectable starburst galaxies and AGN. Our work also corroborates the results of many previous stacking analyses of faint source populations, with improved statistics.

preprint2007arXiv

Improved forecasts for the baryon acoustic oscillations and cosmological distance scale

We present the cosmological distance errors achievable using the baryon acoustic oscillations as a standard ruler. We begin from a Fisher matrix formalism that is upgraded from Seo & Eisenstein (2003). We isolate the information from the baryonic peaks by excluding distance information from other less robust sources. Meanwhile we accommodate the Lagrangian displacement distribution into the Fisher matrix calculation to reflect the gradual loss of information in scale and in time due to nonlinear growth, nonlinear bias, and nonlinear redshift distortions. We then show that we can contract the multi-dimensional Fisher matrix calculations into a 2-dimensional or even 1-dimensional formalism with physically motivated approximations. We present the resulting fitting formula for the cosmological distance errors from galaxy redshift surveys as a function of survey parameters and nonlinearity, which saves us going through the 12-dimensional Fisher matrix calculations. Finally, we show excellent agreement between the distance error estimates from the revised Fisher matrix and the precision on the distance scale recovered from N-body simulations.

preprint2006arXiv

A Catalog of Spectroscopically Confirmed White Dwarfs from the Sloan Digital Sky Survey Data Release 4

We present a catalog of 9316 spectroscopically confirmed white dwarfs from the Sloan Digital Sky Survey Data Release 4. We have selected the stars through photometric cuts and spectroscopic modeling, backed up by a set of visual inspections. Roughly 6000 of the stars are new discoveries, roughly doubling the number of spectroscopically confirmed white dwarfs. We analyze the stars by performing temperature and surface gravity fits to grids of pure hydrogen and helium atmospheres. Among the rare outliers are a set of presumed helium-core DA white dwarfs with estimated masses below 0.3 Msun, including two candidates that may be the lowest masses yet found. We also present a list of 928 hot subdwarfs.

preprint2006arXiv

Black Hole Masses and Eddington Ratios at 0.3<z<4

We study the distribution of Eddington luminosity ratios, L_bol/L_edd, of active galactic nuclei (AGNs) discovered in the AGN and Galaxy Evolution Survey (AGES). We combine H-beta, MgII, and CIV line widths with continuum luminosities to estimate black hole (BH) masses in 407 AGNs, covering the redshift range z~0.3-4 and the bolometric luminosity range L_bol~10^45-10^47 erg/s. The sample consists of X-ray or mid-infrared (24 micron) point sources with optical magnitude R<=21.5 mag and optical emission line spectra characteristic of AGNs. For the range of luminosity and redshift probed by AGES, the distribution of estimated Eddington ratios is well described as log-normal with a peak at L_bol/L_edd ~ 1/4 and a dispersion of 0.3 dex. Since additional sources of scatter are minimal, this dispersion must account for contributions from the scatter between estimated and true BH mass and the scatter between estimated and true bolometric luminosity. Therefore, we conclude that: (1) neither of these sources of error can contribute more than \~0.3 dex rms; and (2) the true Eddington ratios of optically luminous AGNs are even more sharply peaked. Because the mass estimation errors must be smaller than ~0.3 dex, we can also investigate the distribution of Eddington ratios at fixed BH mass. We show for the first time that the distribution of Eddington ratios at fixed BH mass is peaked, and that the dearth of AGNs at a factor ~10 below Eddington is real and not an artifact of sample selection. These results provide strong evidence that supermassive BHs gain most of their mass while radiating close to the Eddington limit, and they suggest that the fueling rates in luminous AGNs are ultimately determined by BH self-regulation of the accretion flow rather than galactic scale dynamical disturbances.

preprint2005arXiv

Discovery of a Fifth Image of the Large Separation Gravitationally Lensed Quasar SDSS J1004+4112

We report the discovery of a fifth image in the large separation lensed quasar system SDSS J1004+4112. A faint point source located 0.2'' from the center of the brightest galaxy in the lensing cluster is detected in images taken with the Advanced Camera for Surveys (ACS) and the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on the Hubble Space Telescope. The flux ratio between the point source and the brightest lensed component in the ACS image is similar to that in the NICMOS image. The location and brightness of the point source are consistent with lens model predictions for a lensed image. We therefore conclude that the point source is likely to be a fifth image of the source quasar. In addition, the NICMOS image reveals the lensed host galaxy of the source quasar, which can strongly constrain the structure of the lensing critical curves and thereby the mass distribution of the lensing cluster.

preprint2005arXiv

SDSSJ102111.02+491330.4: A Newly Discovered Gravitationally Lensed Quasar

We report follow-up observations of two gravitational lens candidates identified in the Sloan Digital Sky Survey (SDSS) dataset. We have confirmed that SDSS J102111.02+491330.4 is a previously unknown gravitationally lensed quasar. This lens system exhibits two images of a $z = 1.72$ quasar, with an image separation of $1{\farcs}14 \pm 0.04$. Optical and near-IR imaging of the system reveals the presence of the lensing galaxy between the two quasar images. Observations of SDSS J112012.12+671116.0 indicate that it is more likely a binary quasar than a gravitational lens. This system has two quasars at a redshift of $z = 1.49$, with an angular separation of $1{\farcs}49 \pm 0.02$. However, the two quasars have markedly different SEDs and no lens galaxy is apparent in optical and near-IR images of this system. We also present a list of 31 SDSS lens candidates which follow-up observations have confirmed are \textit{not} gravitational lenses.

preprint2005arXiv

WFMOS - Sounding the Dark Cosmos

Vast sound waves traveling through the relativistic plasma during the first million years of the universe imprint a preferred scale in the density of matter. We now have the ability to detect this characteristic fingerprint in the clustering of galaxies at various redshifts and use it to measure the acceleration of the expansion of the Universe. The Wide-Field Multi-Object Spectrograph (WFMOS) would use this test to shed significant light on the true nature of dark energy, the mysterious source of this cosmic acceleration. WFMOS would also revolutionise studies of the kinematics of the Milky Way and provide deep insights into the clustering of galaxies at redshifts up to z~4. In this article we discuss the recent progress in large galaxy redshift surveys and detail how WFMOS will help unravel the mystery of dark energy.

preprint2004arXiv

A Catalog of Spectroscopically Identified White Dwarf Stars in the First Data Release of the Sloan Digital Sky Survey

We present the full spectroscopic white dwarf and hot subdwarf sample from the SDSS first data release, DR1. We find 2551 white dwarf stars of various types, 240 hot subdwarf stars, and an additional 144 objects we have identified as uncertain white dwarf stars. Of the white dwarf stars, 1888 are non-magnetic DA types and 171, non-magnetic DBs. The remaining (492) objects consist of all different types of white dwarf stars: DO, DQ, DC, DH, DZ, hybrid stars like DAB, etc., and those with non-degenerate companions. We fit the DA and DB spectra with a grid of models to determine the Teff and log(g) for each object. For all objects, we provide coordinates, proper motions, SDSS photometric magnitudes, and enough information to retrieve the spectrum/image from the SDSS public database. This catalog nearly doubles the known sample of spectroscopically-identified white dwarf stars. In the DR1 imaged area of the sky, we increase the known sample of white dwarf stars by a factor of 8.5. We also comment on several particularly interesting objects in this sample.

preprint2004arXiv

SDSS J1335+0118: A New Two-Image Gravitational Lens

We report the discovery of the two-image gravitationally lensed quasar SDSS J1335+0118. The object was selected as a lens candidate from the Sloan Digital Sky Survey. The imaging and spectroscopic follow-up observations confirm that the system exhibits two gravitationally lensed images of a quasar at z=1.57. The image separation is 1.56''. We also detect an extended component between the two quasar images, likely the lensing galaxy. Preliminary mass modeling predicts the differential time delay Δt ~ 30 h^{-1} day assuming the redshift of the lens galaxy is 0.5.

preprint2003arXiv

SDSS J115517.35+634622.0: A Newly Discovered Gravitationally Lensed Quasar

We report the discovery of SDSSJ115517.35+634622.0, a previously unknown gravitationally lensed quasar. The lens system exhibits two images of a $z = 2.89$ quasar, with an image separation of $1{\farcs}832 \pm 0.007$ . Near-IR imaging of the system reveals the presence of the lensing galaxy between the two quasar images. Based on absorption features seen in the Sloan Digital Sky Survey (SDSS) spectrum, we determine a lens galaxy redshift of $z = 0.1756$. The lens is rather unusual in that one of the quasar images is only $0{\farcs}22\pm0{\farcs}07$ ($\sim 0.1 R_{\rm eff}$) from the center of the lens galaxy and photometric modeling indicates that this image is significantly brighter than predicted by a SIS model. This system was discovered in the course of an ongoing search for strongly lensed quasars in the dataset from the SDSS.

preprint1998arXiv

Cosmic Complementarity: H_0 and Omega_m from Combining CMB Experiments and Redshift Surveys

We show that upcoming CMB satellite experiments and large redshift surveys can be used together to yield 5% determinations of H_0 and Omega_m, an order of magnitude improvement over CMB data alone. CMB anisotropies provide the sound horizon at recombination as a standard ruler. For reasonable baryon fractions, this scale is imprinted on the galaxy power spectrum as a series of spectral features. Measuring these features in redshift space determines the Hubble constant, which in turn yields Omega_m once combined with CMB data. Since the oscillations in both power spectra are frozen in at recombination, this test is insensitive to low-redshift cosmology.

preprint1998arXiv

The Structure of Structure Formation Theories

We study the general structure of models for structure formation, with applications to the reverse engineering of the model from observations. Through a careful accounting of the degrees of freedom in covariant gravitational instability theory, we show that the evolution of structure is completely specified by the stress history of the dark sector. The study of smooth, entropic, sonic, scalar anisotropic, vector anisotropic, and tensor anisotropic stresses reveals the origin, robustness, and uniqueness of specific model phenomenology. We construct useful and illustrative analytic solutions that cover cases with multiple species of differing equations of state relevant to the current generation of models, especially those with effectively smooth components. We present a simple case study of models with phenomenologies similar to that of a LambdaCDM model to highlight reverse-engineering issues. A critical-density universe dominated by a single type of dark matter with the appropriate stress history can mimic a LambdaCDM model exactly.

preprint1997arXiv

Baryonic Features in the Matter Transfer Function

We provide scaling relations and fitting formulae for adiabatic cold dark matter cosmologies that account for all baryon effects in the matter transfer function to better than 10% in the large-scale structure regime. They are based upon a physically well-motivated separation of the effects of acoustic oscillations, Compton drag, velocity overshoot, baryon infall, adiabatic damping, Silk damping, and cold-dark-matter growth suppression. We also find a simpler, more accurate, and better motivated form for the zero baryon transfer function than previous works. These descriptions are employed to quantify the amplitude and location of baryonic features in linear theory. While baryonic oscillations are prominent if the baryon fraction exceeds $Ω_0 h^2 + 0.2$, the main effect in more conventional cosmologies is a sharp suppression in the transfer function below the sound horizon. We provide a simple but accurate description of this effect and stress that it is not well approximated by a change in the shape parameter $Γ$.

preprint1997arXiv

HOP: A New Group-Finding Algorithm for N-body Simulations

We describe a new method (HOP) for identifying groups of particles in N-body simulations. Having assigned to every particle an estimate of its local density, we associate each particle with the densest of the N_hop particles nearest to it. Repeating this process allows us to trace a path, within the particle set itself, from each particle in the direction of increasing density. The path ends when it reaches a particle that is its own densest neighbor; all particles reaching the same such particle are identified as a group. Combined with an adaptive smoothing kernel for finding the densities, this method is spatially adaptive, coordinate-free, and numerically straight-forward. One can proceed to process the output by truncating groups at a particular density contour and combining groups that share a (possibly different) density contour. While the resulting algorithm has several user-chosen parameters, we show that the results are insensitive to most of these, the exception being the outer density cutoff of the groups.