Issue 
A&A
Volume 552, April 2013



Article Number  A96  
Number of page(s)  18  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201220724  
Published online  04 April 2013 
Baryon acoustic oscillations in the Lyα forest of BOSS quasars^{⋆}
^{1}
APC, Université Paris DiderotParis 7, CNRS/IN2P3, CEA, Observatoire de
Paris, 10 rue A. Domon & L.
Duquet, Paris,
France
email:
ngbusca@apc.univparis7.fr
^{2}
CEA, Centre de Saclay, IRFU, 91191
GifsurYvette,
France
^{3}
Lawrence Berkeley National Laboratory,
1 Cyclotron Road, Berkeley, CA
94720,
USA
^{4}
Department of Physics and Astronomy, University of
California, Irvine,
CA
92697,
USA
^{5}
Institute of Cosmology and Gravitation, Dennis Sciama Building,
University of Portsmouth, Portsmouth, PO1
3FX, UK
^{6}
Bldg 510 Brookhaven National Laboratory,
Upton, NY
11973,
USA
^{7}
Apache Point Observatory, PO Box 59, Sunspot, NM
88349,
USA
^{8}
Department of Physics and Astronomy, University of
Utah, 115 S 1400 E,
Salt Lake City, UT
84112,
USA
^{9}
Institute for Advanced Study, Einstein Drive, Princeton, NJ
08540,
USA
^{10}
Bruce and Astrid McWilliams Center for Cosmology, Carnegie Mellon
University, Pittsburgh, PA
15213,
USA
^{11}
HarvardSmithsonian Center for Astrophysics, Harvard
University, 60 Garden
St., Cambridge
MA
02138,
USA
^{12}
Center for Cosmology and Particle Physics, New York
University, New York,
NY
10003,
USA
^{13}
Department of Physics and Center for Cosmology and AstroParticle
Physics, Ohio State University, Columbus, OH
43210,
USA
^{14}
MaxPlanckInstitut für Astronomie, Königstuhl 17, 69117
Heidelberg,
Germany
^{15}
Department of Astronomy, University of Wisconsin,
475 North Charter Street,
Madison, WI
53706,
USA
^{16}
Steward Observatory, University of Arizona,
933 N. Cherry Ave.,
Tucson, AZ
85121,
USA
^{17} Institució Catalana de Recerca i Estudis Avançats,
Barcelona, Catalonia
^{18}
Institut de Ciències del Cosmos, Universitat de
Barcelona/IEEC, 08028
Barcelona, Catalonia
^{19}
Department of Physics and Astronomy, University of
Wyoming, Laramie,
WY
82071,
USA
^{20}
Université Paris 6 et CNRS, Institut d’Astrophysique de
Paris, 98bis blvd.
Arago, 75014
Paris,
France
^{21}
Department of Astronomy and Astrophysics, The Pennsylvania State
University, University
Park, PA
16802,
USA
^{22}
Institute for Gravitation and the Cosmos, ThePennsylvania State
University, University
Park, PA
16802,
USA
^{23}
INAF, Osservatorio Astronomico di Trieste, via G. B. Tiepolo
11, 34131
Trieste,
Italy
^{24} INFN/National Institute for Nuclear Physics, via Valerio 2,
34127 Trieste, Italy
^{25}
Department of Astronomy, Ohio State University,
140 West 18th Avenue,
Columbus, OH
43210,
USA
^{26}
Department of Astronomy and Astrophysics and the Enrico Fermi
Institute, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois, 60615, USA
^{27}
Institute of Theoretical Physics, University of
Zurich, 8057
Zurich,
Switzerland
^{28}
Departamento de Astronomía, Universidad de Chile,
Casilla 36D, Santiago, Chile
Received:
12
November
2012
Accepted:
4
February
2013
We report a detection of the baryon acoustic oscillation (BAO) feature in the threedimensional correlation function of the transmitted flux fraction in the Lyα forest of highredshift quasars. The study uses 48 640 quasars in the redshift range 2.1 ≤ z ≤ 3.5 from the Baryon Oscillation Spectroscopic Survey (BOSS) of the third generation of the Sloan Digital Sky Survey (SDSSIII). At a mean redshift z = 2.3, we measure the monopole and quadrupole components of the correlation function for separations in the range 20 h^{1} Mpc < r < 200 h^{1} Mpc. A peak in the correlation function is seen at a separation equal to (1.01 ± 0.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 D_{A} 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 observations, this result implies H(z = 2.3) = (224 ± 8) km s^{1} Mpc^{1}, indicating that the time derivative of the cosmological scale parameter ȧ = H(z = 2.3)/(1 + z) is significantly greater than that measured with BAO at z ~ 0.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 darkenergy dominated universe.
Key words: cosmology: observations / dark energy / largescale structure of Universe / cosmological parameters
Appendices are available in electronic form at http://www.aanda.org
© ESO, 2013
1. Introduction
Baryon acoustic oscillations (BAO) in the prerecombination universe have striking effects on the anisotropies of the cosmic microwave background (CMB) and on the large scale structure (LSS) of matter (Weinberg et al. 2012, and references therein). The BAO effects were first seen in the series of peaks in the CMB angular power spectrum (de Bernardis et al. 2000). Subsequently, the BAO relic at redshift z ~ 0.3 was seen (Eisenstein et al. 2005; Cole et al. 2005) as a peak in the galaxygalaxy correlation function at a comoving distance corresponding to the sound horizon at recombination. For the WMAP7 cosmological parameters (Komatsu et al. 2011), the expected comoving scale of the BAO peak is r_{s} = 153 Mpc, with an uncertainty of ≈1%.
The BAO peak in the correlation function at a redshift z appears at an angular separation Δθ = r_{s}/(1 + z)D_{A}(z) and at a redshift separation Δz = r_{s}H(z)/c, where D_{A} and H are the angular distances and expansion rates. Measurement of the peak position at any redshift thus constrains the combinations of cosmological parameters that determine r_{s}H and r_{s}/D_{A}. While the possibility of measuring both combinations is beginning to be exploited (Chuang & Wang 2012; Xu et al. 2012), most present measurements have concentrated on the combination , which determines the peak position for an isotropic distribution of galaxy pairs and an isotropic clustering strength. The “BAO Hubble diagram”, D_{V}/r_{s} vs. z, now includes the Sloan Digital Sky Survey (SDSS) measurement (Eisenstein et al. 2005) updated to the DR7 (Abazajian et al. 2009) sample and combined with 2dF data (Percival et al. 2010), the 6dF point at z = 0.1 (Beutler et al. 2011), the WiggleZ points at (0.4 < z < 0.8) (Blake et al. 2011a), and a reanalysis of the SDSS DR7 sample that uses reconstruction (Eisenstein 2007; Padmanabhan et al. 2009) to sharpen the precision of the BAO measurement (Padmanabhan et al. 2012; Mehta et al. 2012). Recently, the Baryon Oscillation Spectroscopic Survey (BOSS; Dawson et al. 2013) of SDSSIII (Eisenstein et al. 2011) has added a precise measurement at z ~ 0.57 (Anderson et al. 2012). BOSS has also reported a measurement of D_{A}(z = 0.55)/r_{s} based on the BAO structure in the angular power spectrum of galaxies (Seo et al. 2012).
In this paper, we present an observation of the BAO peak at z ~ 2.3 found in the flux correlation function of the Lyα forest of BOSS quasars. This is the first such observation at a redshift where the expansion dynamics is matterdominated, z > 0.8. The possibility of such a measurement was suggested by McDonald (2003) and White (2003) and first studied in detail by McDonald & Eisenstein (2007). While the galaxy BAO measurements are most sensitive to , the Lyα flux transmission is more sensitive to peculiar velocity gradient effects, which enhance redshift distortions and shift our sensitivity to the expansion rate. As we shall show below, the most accurately measured combination from the Lyα forest BAO peak is ∝, and the present BOSS data set allows us to determine its value to a precision of 3.5%. Combining this result with constraints from CMB observations allows us to deduce the value of H(z = 2.3) accurate to 4%. Comparing our results with measurements of H_{0} and of H(0.2 < z < 0.8) reveals the expected sequence of deceleration and acceleration in models with dark energy.
The last decade has seen increasing use of Lyα absorption to investigate LSS. The number of quasars in early studies (Croft et al. 1999, 2002; McDonald et al. 2000, 2006; Viel et al. 2004) was enough only to determine the Lyα absorption correlation along individual lines of sight. With the BOSS project the surface density of quasars is sufficient to probe the full threedimensional distribution of neutral hydrogen. A study using the first 10 000 BOSS quasars was presented by Slosar et al. (2011). This sample provided clear evidence for the expected longrange correlations, including the redshiftspace distortions due to the gravitational growth of structure. With the SDSS data release DR9 (Ahn et al. 2012), we now have ~60 000 quasars at z ~ 2.3 (Pâris et al. 2012), with a high enough surface density to observe the BAO peak.
The use of Lyα absorption to trace matter has certain interesting differences from the use of galaxies. Galaxy surveys provide a catalog of positions in redshift space that correspond to points of high overdensities. On the other hand, the forest region of a quasar spectrum provides a complete mapping of the absorption over a ~400 h^{1} Mpc (comoving) range starting about 100 h^{1} Mpc in front of the quasar (so as to avoid the necessity of modeling the quasar’s Lyα emission line). To the extent that quasar lines of sight are random, a large collection of quasars can provide a nearly unbiased sample of points where the absorption is measured. Cosmological simulations (Cen et al. 1994; Petitjean et al. 1995; Zhang et al. 1995; Hernquist et al. 1996; MiraldaEscudé et al. 1996; Theuns et al. 1998) indicate that most of the Lyα absorption is due to cosmic filamentary structures with overdensities of order one to ten, much lower than the overdensities of virialized halos sampled by galaxies. These simulations have also indicated that, on large scales, the mean Lyα absorption is a linear tracer of the mass overdensity (Croft et al. 1997, 1998; Weinberg et al. 1998; McDonald et al. 2000; McDonald 2003), implying a relation of the power spectrum of the measured absorption to that of the underlying mass fluctuations. Finally, the forest is observable in a redshift range inaccessible to current large galaxy surveys and where theoretical modeling is less dependent on nonlinear effects in cosmological structure formation. These factors combine to make Lyα absorption a promising tracer of mass that is complementary to galaxy tracers.
With Lyα forest measurements along multiple sightlines, one can attempt to reconstruct the underlying 3dimensional mass density field (Nusser & Haehnelt 1999; Pichon et al. 2001; Gallerani et al. 2011), from which one can investigate topological characteristics (Caucci et al. 2008) or the power spectrum (Kitaura et al. 2012). However, the BOSS sample is fairly sparse, with a typical transverse sightline separation ~15 h^{1} Mpc (comoving), and the signaltonoise ratio in individual spectra is low (see Fig. 2 below), which makes it poorly suited to such reconstruction techniques. In this paper we take the more direct approach of measuring the BAO feature in the correlation function of transmitted flux fraction (the ratio of the observed flux to the flux expected in the absence of absorption).
The use of the Lyα forest is handicapped by the fact that not all fluctuations of the transmitted flux fraction are due to fluctuations of the density of hydrogen. Because the neutral hydrogen density is believed to be determined by photoionization equilibrium with the flux of ultraviolet (UV) photons from stars and quasars, variations of the UV flux may contribute to fluctuations of the transmitted flux (Worseck & Wisotzki 2006). In addition, hydrogen can be selfshielded to ionizing photons, leading to much higher neutral fractions and damping wings in strong absorption systems (FontRibera & MiraldaEscud 2012). Metals present in the intergalactic medium provide additional absorption superimposed on the Lyα absorption (Pieri et al. 2010). A further complication results from the fact that the flux correlation function uses the fluctuations of the transmitted flux about its mean value; this requires an estimate of the product of the mean absorption (as a function of absorber redshift) and the unabsorbed flux for individual quasars (Le Goff et al. 2011).
All of these complications are important for a complete understanding of the statistics of the Lyα forest. Fortunately, they are not expected to produce a sharp peaklike feature in the correlation function, so the interpretation of a peak as due to BAO should give robust constraints on the cosmological parameters.
The BOSS collaboration has performed three independent analyses to search for BAO in Lyα forest. This paper presents two of them, both of which aim to analyze the forest using simple procedures at the expense of some loss of sensitivity. The third analysis, with the goal of an optimal measurement of the flux correlation function with a more complex method, is described in a separate publication (Slosar et al. 2013).
Our methods are tested extensively on a set of detailed mock catalogs of the BOSS Lyα forest data set. These mock catalogs, which use the method presented by FontRibera et al. (2012), will be described in detail in a forthcoming public release paper (Bailey et al., in prep.). In addition, the BOSS collaboration have also released a fiducial version of the DR9 Lyα forest spectra (Lee et al. 2012b), with various perobject products including masks, continua, and noise correction vectors designed to aid in Lyα forest analysis. While our analyses implement their own sample selection criteria and continuum determination procedures, we have also applied our measurement to this Lee et al. (2012b) sample.
This paper is organized as follows. Section 2 presents the BOSS quasar sample used in this analysis and the procedure used to produce the quasar spectra. Section 3 describes the analysis to measure the correlation function. Section 4 derives the monopole and quadrupole components of the correlation function and determines the significance of the peak observed in these functions at the BAO scale. The cosmological implications of our detection of a BAO peak are discussed in Sect. 5. Finally, Appendix A provides a brief description of the mock spectra used to test our methodology and Appendix B shows the result of our BAO measurement applied to the BOSS Lyα sample of Lee et al. (2012b).
2. The BOSS quasar sample and data reduction
The BOSS project (Dawson et al. 2013) of SDSSIII (Eisenstein et al. 2011) is obtaining the spectra of ~1.6 × 10^{6} luminous galaxies and ~150 000 quasars. The project uses upgraded versions of the SDSS spectrographs (Smee et al. 2012) mounted on the Sloan 2.5m telescope (Gunn et al. 2006) at Apache Point, New Mexico. BOSS galaxy and quasar spectroscopic targets are selected using algorithms based primarily on photometry from the SDSS camera (Gunn et al. 1998; York et al. 2000) in the ugriz bands (Fukugita et al. 1996; Smith et al. 2002) reduced and calibrated as described by Stoughton et al. (2002), Pier et al. (2003), and Padmanabhan et al. (2008). Targets are assigned to fibers appropriately positioned in the 3° diameter focal plane according to a specially designed tiling algorithm (Blanton et al. 2003). Fibers are fixed in place by a pierced metal plate drilled for each observed field and fed to one of two spectrographs. Each exposed plate generates 1000 spectra covering wavelengths of 360 to 1000 nm with a resolving power ranging from 1500 to 3000 (Smee et al. 2012). A median of 631 of these fibers are assigned to galactic targets and 204 to quasar targets. The BOSS spectroscopic targets are observed in dark and gray time, while the brighttime is used by other SDSSIII surveys (see Eisenstein et al. 2011).
The quasar spectroscopy targets are selected from photometric data with a combination of algorithms (Richards et al. 2009; Yeche et al. 2009; Kirkpatrick et al. 2011; Bovy et al. 2011; PalanqueDelabrouille et al. 2011; for a summary, see Ross et al. (2012)). The algorithms use SDSS fluxes and, for SDSS Stripe 82, photometric variability. When available, we also use data from nonoptical surveys (Bovy et al. 2012): the GALEX survey (Martin et al. 2005) in the UV; the UKIDSS survey (Lawrence et al. 2007) in the nearinfrared (NIR), and the FIRST survey (Becker et al. 1995) in the radio.
The quasar spectroscopy targets are divided into two samples “CORE” and “BONUS”. The CORE sample consists of 20 quasar targets per square degree selected from SDSS photometry with a uniform algorithm, for which the selection efficiency for z > 2.1 quasars is ~50%. The selection algorithm for the CORE sample (Bovy et al. 2011) was fixed at the end of the first year of the survey, thus making it useful for studies that require a uniform target selection across the sky. The BONUS sample was chosen from a combination of algorithms with the purpose of increasing the density on the sky of observed quasars beyond that of the CORE sample. The combined samples yield a mean density of identified quasars of 15 deg^{2} with a maximum of 20 deg^{2}, mostly in zones where photometric variability, UV, and/or NIR data are available. The combined BONUS plus CORE sample can be used for Lyα BAO studies, which require the highest possible quasar density in a broad sky area but are insensitive to the uniformity of the quasar selection criteria because the structure being mapped is in the foreground of these quasar backlights.
The data presented here consist of the DR9 data release (Ahn et al. 2012) covering ~3000 deg^{2} of the sky shown in Fig. 1. These data cover about onethird of the ultimate BOSS footprint.
Fig. 1
HammerAitoff projection in equatorial coordinates of the BOSS DR9 footprint. The observations cover ~3000 deg^{2}. 
The data were reduced with the SDSSIII pipeline as described in Bolton et al. (2012). Typically four exposures of 15 min were coadded in pixels of wavelength width ~0.09 nm. Besides providing flux calibrated spectra, the pipeline provided preliminary object classifications (galaxy, quasar, star) and redshift estimates.
The spectra of all quasar targets were visually inspected, as described in Pâris et al. (2012), to correct for misidentifications or inaccurate redshift determinations and to flag broad absorption lines (BAL). Damped Lyα troughs are visually flagged, but also identified and characterized automatically (Noterdaeme et al. 2012). The visual inspection of DR9 confirms 60 369 quasars with 2.1 ≤ z_{q} ≤ 3.5. In order to simplify the analysis of the Lyα forest, we discarded quasars with visually identified BALs and DLAs, leaving 48 640 quasars.
For the measurement of the flux transmission, we use the restframe wavelength interval (1)The range is bracketed by the Lyα and Lyβ emission lines at 121.6 and 102.5 nm. The limits are chosen conservatively to avoid problems of modeling the shapes of the two emission lines and to avoid quasar proximate absorbers. The absorber redshift, z = λ/λ_{Lyα} − 1, is in the range 1.96 < z < 3.38. The lower limit is set by the requirement that the observed wavelength be greater than 360 nm below which the system throughput is less than 10% its peak value. The upper limit comes from the maximum quasar redshift of 3.5, beyond which the BOSS surface density of quasars is not sufficient to be useful. The distribution of absorber redshift is shown in Fig. 2 (top panel). When given the weights used for the calculation of the correlation function (Sect. 3.3), the absorbers have a mean redshift of ⟨ z ⟩ = 2.31.
Fig. 2
Top: weighted distribution of absorber redshifts used in the calculation of the correlation function in the distance range 80 h^{1} Mpc < r < 120 h^{1} Mpc. Bottom: distribution of signaltonoise ratio for analysis pixels (triplets of pipeline pixels) averaged over the forest region. 
For the determination of the correlation function, we use “analysis pixels” that are the flux average over three adjacent pipeline pixels. Throughout the rest of this paper, “pixel” refers to analysis pixels unless otherwise stated. The effective width of these pixels is 210 km s^{1}, i.e. an observedwavelength width ~0.27 nm ~ 2 h^{1} Mpc. The total sample of 48 640 quasars thus provides ~8 × 10^{6} measurements of Lyα absorption over a total volume of ~20 h^{3} Gpc^{3}.
Figure 2 (bottom panel) shows the distribution of the signaltonoise ratio for pixels averaged over the forest region. The relatively modest mean value of 5.17 reflects the exposure times necessary to acquire such a large number of spectra.
In addition to the BOSS spectra, we analyzed 15 sets of mock spectra that were produced by the methods described in Appendix A. These spectra do not yet reproduce all of the characteristics of the BOSS sample, but they are nevertheless useful for a qualitative understanding of the shape of the measured correlation function. More importantly, they are useful for understanding the detectability of a BAOlike peak and the precision of the measurement of its position.
3. Measurement of the correlation function
The flux correlation function can be determined through a simple twostep process. In the first step, for each pixel in the forest region (Eq. (1)) of quasar q, the measured flux f_{q}(λ) at observed wavelength λ is compared with the mean expected flux, , thus defining the “delta field”: (2)Here, C_{q}(λ) is the unabsorbed flux (the socalled “continuum”) and is the mean transmitted fraction at the HI absorber redshift. The quantities λ and z in Eq. (2) are not independent but related via z = λ/λ_{Lyα} − 1.
Figure 3 shows an example of an estimation for C_{q}(λ) (blue line) and (red line). Our two methods for estimating C_{q} and are described in Sects. 3.1 and 3.2.
Fig. 3
An example of a BOSS quasar spectrum of redshift 3.239. The red and blue lines cover the forest region used here, 104.5 < λ_{rf} < 118.0. This region is sandwiched between the quasar’s Lyβ and Lyα emission lines respectively at 435 and 515 nm. The blue line is an estimate of the continuum (unabsorbed flux) by method 2 and the red line is the estimate of the product of the continuum and the mean absorption by method 1. 
In the second step, the correlation function is calculated as a weighted sum of products of the deltas: (3)where the w_{ij} are weights and each i or j indexes a measurement on a quasar q at wavelength λ. The sum over (i,j) is understood to run over all pairs of pixels of all pairs of quasars within A defining a region in space of pixel separations, r_{i} − r_{j}. The region A is generally defined by a range r_{min} < r < r_{max} and μ_{min} < μ < μ_{max} with: (4)where (r_{i} − r_{j})_{∥} is the component along the line of sight. Separations in observational pixel coordinates (RA, Dec, z) are transformed to (r,μ) in units of h^{1} Mpc by using a ΛCDM fiducial cosmology with matter and vacuum densities of (5)In the sum (3), we exclude pairs of pixels from only one quasar to avoid the correlated errors in δ_{i} and δ_{j} coming from the estimate of C_{q}. Note that the weights in Eq. (3) are set to zero for pixels flagged by the pipeline as having problems due, e.g., to sky emission lines or cosmic rays.
A procedure for determining ξ is defined by its method for estimating the expected flux and by its choice of weights, w_{ij}. The two methods described here use the same technique to calculate weights but have different approaches to estimate . We will see that the two methods produce correlation functions that have no significant differences. However, the two independent codes were invaluable for consistency checks throughout the analysis.
The two methods were “blind” to the extent that many of the procedures were defined during tests either with mock data or with the real data in which we masked the region of the peak in the correlation function. Among those aspects fixed in this way were the quasar sample, the continuum determination, the weighting, the extraction of the monopole and quadrupole correlation function and the determination of the peak significance (Sect. 4). This early freezing of procedures resulted in some that are suboptimal but which will be improved in future analyses. We note, however, that the procedures used to extract cosmological information (Sect. 5) were decided on only after demasking the data.
3.1. Continuum fits, method 1
Both methods for estimating the product assume that C_{q} is, to first approximation, proportional to a universal quasar spectrum that is a function of restframe wavelength, λ_{rf} = λ/(1 + z_{q}) (for quasar redshift z_{q}), multiplied by a mean transmission fraction that slowly varies with absorber redshift. Following this assumption, the universal spectrum is found by stacking the appropriately normalized spectra of quasars in our sample, thus averaging out the fluctuating Lyα absorption. The product for individual quasars is then derived from the universal spectrum by normalizing it to account for the quasar’s mean forest flux and then modifying its slope to account for spectralindex diversity and/or photospectroscopic miscalibration.
Method 1 estimates directly the product in Eq. (2). An example is given by the red line in Fig. 3. The estimate is made by modeling each spectrum as (6)where a_{q} is a normalization, b_{q} a “deformation parameter”, and ⟨ λ ⟩ is the mean wavelength in the forest for the quasar q and is the mean normalized flux obtained by stacking spectra in bins of width Δz = 0.1: (7)Here z is the redshift of the absorption line at observed wavelength λ (z = λ/λ_{Lyα} − 1), f_{q} is the observed flux of quasar q at wavelength λ and is the average of the flux of quasar q for 127.5 < λ_{rf} < 128.5 nm. The weight w_{q}(λ) is given by . The quantity ivar is the pipeline estimate of the inverse flux variance in the pixel corresponding to wavelength λ. The quantity is the contribution to the variance in the flux due to the LSS. We approximate it by its value at the typical redshift of the survey, z ~ 2.3: (Sect. 3.3).
Figure 4 shows the resulting mean δ_{i} as a function of observed wavelength. The mean fluctuates about zero with up to 2% deviations with correlated features that include the H and K lines of singly ionized calcium (presumably originating from some combination of solar neighborhood, interstellar medium and the Milky Way halo absorption) and features related to Balmer lines. These Balmer features are a byproduct of imperfect masking of Balmer absorption lines in Fstar spectroscopic standards, which are used to produce calibration vectors (in the conversion of CCD counts to flux) for DR9 quasars. Therefore such Balmer artifacts are constant for all fibers in a plate fed to one of the two spectrographs and so they are approximately constant for every “halfplate”.
If unsubtracted, the artifacts in Fig. 4 would lead to spurious correlations, especially between pairs of pixels with separations that are purely transverse to the line of sight. We have made a global correction by subtracting the quantity ⟨ δ ⟩ (λ) in Fig. 4 (unsmoothed) from individual measurements of δ. This is justified if the variance of the artifacts from halfplatetohalfplate is sufficiently small, as halfplatewide deviations from our global correction could, in principle add spurious correlations.
We have investigated this variance both by measuring the Balmer artifacts in the calibration vectors themselves and by studying continuum regions of all available quasars in the DR9 sample. Both studies yield no detection of excess variance arising from these artifacts, but do provide upper limits. The study of the calibration vectors indicate that the squareroot of the variance is less than 20% of the mean Balmer artifact deviations and the study of quasar spectra indicate that the squareroot of the variance is less than 100% of the mean Balmer artifacts (and less than 50% of the mean calcium line deviations).
We then performed Monte Carlo simulations by adding a random sampling of our measured artifacts to our data to confirm that our global correction is adequate. We found that there is no significant effect on the determination of the BAO peak position, even if the variations are as large as that allowed in our tests.
Fig. 4
Mean of δ(λ) as a function of observed wavelength (method 1). Systematic offsets from zero are seen at the 2% level. The calcium lines (393.4, 396.8 nm) is present. The features around the hydrogen lines Hγ, δ and ϵ (434.1, 410.2, 397.0 nm) are artifacts from the use of Fstars for the photocalibration of the spectrometer. 
3.2. Continuum fits, method 2
Method 1 would be especially appropriate if the fluxes had a Gaussian distribution about the mean absorbed flux, . Since this is not the case, we have developed method 2 which explicitly uses the probability distribution function for the transmitted flux fraction F, P(F,z), where 0 < F < 1. We use the P(F,z) that results from the lognormal model used to generate mock data (see Appendix A).
Using P(F,z), we can construct for each BOSS quasar the PDF of the flux in pixel i, f_{i}, by assuming a continuum C_{q}(λ_{i}) and convolving with the pixel noise, σ_{i}: (8)The continuum is assumed to be of the form (9)where is the mean flux as determined by stacking spectra as follows: (10)as in Eq. (7) except that here there is no redshift binning. The parameters a_{q} and b_{q} are then determined for each quasar by maximizing a likelihood given by (11)Figure 3 shows the C_{q}(λ) estimated for a typical quasar (blue line).
The last element necessary to use Eq. (2) is the mean transmitted flux fraction . If P(F,z) derived from the mocks were the true distribution of the transmitted flux fraction, then could simply be computed from the average of this distribution. Since this is not precisely true, we determine from the data by requiring that the mean of the delta field vanish for all redshifts. The we obtain is shown in Fig. 5. The unphysical wiggles in the derived are associated with the aforementioned residuals in for method 1 (Fig. 4).
There is one inevitable effect of our two continuum estimating procedures. The use of the forest data in fitting the continuum effectively forces each quasar to have a mean absorption near that of the mean for the entire quasar sample. This approach introduces a spurious negative correlation between pixels on a given quasar even when well separated in wavelength. This negative correlation has no direct effect on our measurement of the flux correlation function because we do not use pixel pairs from the same quasar. However, the physical correlation between absorption on neighboring quasars causes the unphysical negative correlation for individual quasars to generate a negative contribution to the correlation measured with quasar pairs. Fortunately, this distortion is a smooth function of scale so it can be expected to have little effect on the observability or position of the BAO peak. This expectation is confirmed by analysis of the mock spectra (Sect. 5).
Fig. 5
Mean transmitted flux fraction as a function of redshift obtained from the continuum fits with method 2. Data are shown in black, mock000 in red and the input mean transmitted flux fraction in blue. 
3.3. Weights
A discussion on the optimal use of weights for the Lyα correlation function is found in McQuinn & White (2011). Here we simply choose the weights w_{ij} so as to approximately minimize the relative error on estimated with Eq. (3). In the approximation of uncorrelated pixels, the variance of is (12)where the pixel variance, ξ_{ii}, includes contributions from both observational noise and LSS. The signaltonoise ratio is: (13)Because of LSS growth and redshift evolution of the mean absorption, the ξ_{ij} depend on redshift and we use the measured dependence of the 1d correlation function (McDonald et al. 2006) (14)Maximizing the signaltonoise ratio with respect to w_{ij} this gives: (15)For this expression to be used, we require a way of estimating the ξ_{ii}. We assume that it can be decomposed into a noise term and a LSS term (σ_{LSS}): (16)where is the pipeline estimate of the noisevariance of pixel i and η is a factor that corrects for a possible misestimate of the variance by the pipeline.
We then organize the data in bins of and redshift. In each such bin, we measure the variance of δ_{i}, which serves as an estimator of ξ_{ii} for the bin in question. The two functions η(z) and can then be determined by fitting Eq. (16).
These fits are shown in Fig. 6. The top panel shows that the measured inverse variance follows the inverse pipeline variance until saturating at the redshiftdependent LSS variance (shown on the bottom left panel). For z > 3, there are not enough pixel pairs to determine η(z) and . In this high redshift range, we assumed η = 1 and extrapolated σ_{LSS}(z) with a seconddegree polynomial fit to the z < 3 data.
Fig. 6
Top panel: inverse total variance in bins of redshift as a function of the pipeline inverse variance. Bottom panel: parameters of the fit: the LSS contribution σ_{LSS} (left) and the pipeline correction factor η (right) as a function of redshift. The lines show fits to the data as explained in the text. 
3.4. ξ(r, μ)
The procedure described above was used to determine ξ(r,μ) through Eq. (3) in rbins of width 4 h^{1} Mpc (centered at 2,6,..., 198 h^{1} Mpc) and in μbins of width 0.02, (centered at 0.01, 0.03, ... 0.99). The 50 × 50 r − μ bins have an average of 6 × 10^{6} terms in the sum (3) with an average nominal variance of ξ for individual bins of (10^{4})^{2} as given by (Eq. (12)).
Figure 7 shows an example of ξ(r,μ) for the r bin centered on 34 h^{1} Mpc. The blue dots are the data and the red dots are the mean of the 15 mocks. The function falls from positive to negative values with increasing μ, as expected from redshift distortions. The effect is enhanced by the deformation due to the continuum subtraction.
Figure 8 presents ξ(r,μ) averaged over three bins in μ. A clear peak at the expected BAO position, r_{s} = 105 h^{1} Mpc, is present in the bin 0.8 < μ < 1.0 corresponding to separation vectors within 37° of the lineofsight. The curves show the best fits for a ΛCDM correlation function, as described in Sect. 5.
The data were divided into various subsamples to search for systematic errors in ξ(r,μ). For example, searches were made for differences between the northern and southern Galactic cap regions and between higher and lower signaltonoise ratio quasars. No significant differences were found in the overall shape and amplitude of the correlation function. We also verified that the BAO peak position does not change significantly when wavelength slices of Lyα forest data are eliminated, in particular slices centered on the Balmer features in Fig. 4. The peak position also does not change significantly if the subtraction of the mean δ (Fig. 4) is suppressed.
Fig. 7
ξ(r,μ) vs. μ^{2} for the bin centered on r = 34 h^{1} Mpc. The red dots are the mean of the 15 mocks and the blue dots are the data. 
Fig. 8
ξ(r,μ) averaged over 0.1 < μ < 0.5, 0.5 < μ < 0.8 and 0.8 < μ < 1. The curves give fits (Sect. 5) to the data imposing concordance ΛCDM cosmology. The BAO peak is most clearly present in the data for μ > 0.8. 
4. The monopole and quadrupole
The analysis of the correlation function was performed in the framework of the standard multipole decomposition (Hamilton 1992). For each bin in r we fit a monopole (ℓ = 0) and quadrupole (ℓ = 2) to the angular dependence: (17)where P_{ℓ} is the ℓLegendre polynomial. We ignore the small and poorly determined ℓ = 4 term. This fit is performed using a simple χ^{2} minimization with the nominal variance (Eq. (12)) and ignoring the correlations between bins. This approach makes the fit slightly suboptimal. (Later, we will correctly take into account correlations between rbins of the monopole and quadrupole.) We also exclude from this fit the portion μ < 0.1 to avoid residual biases due to correlated sky subtraction across quasars; this has a negligible impact on the fits and, at any rate, there is little BAO signal at low μ.
Figure 9 displays the monopole and quadrupole signals found by the two methods. The two methods are slightly offset from one another, but the peak structure is very similar. Figure 9 also shows the combination ξ_{0} + 0.1ξ_{2} which, because of the small monopolequadrupole anticorrelation (Sect. 4.1), is a betterdetermined quantity. The peak structure seen in Fig. 8 is also present in these figures.
Because of the continuum estimation procedure (Sects. 3.1 and 3.2), we can expect that the monopole and quadrupole shown in Fig. 9 are deformed with respect to the true monopole and quadrupole. The most important difference is that the measured monopole is negative for 60 h^{1} Mpc < r < 100 h^{1} Mpc while the true ΛCDM monopole remains positive for all r < 130 h^{1} Mpc. The origin of the deformation in the continuum estimate is demonstrated in Appendix A where both the true and estimated continuum can be used to derive the correlation function (Fig. A.1). As expected, the deformation is a slowly varying function of r so neither the position of the BAO peak nor its amplitude above the slowly varying part of the correlation function are significantly affected.
Fig. 9
Monopole (upper panel) and quadrupole (middle panel) correlation functions found by method 1 (red) and method 2 (black). The bottom panel shows the combination ξ_{0} + 0.1ξ_{2} found by method 1 (red) and method 2 (black). 
4.1. Covariance of the monopole and quadrupole
In order to determine the significance of the peak we must estimate the covariance matrix of the monopole and quadrupole. If the fluctuations δ_{i} in Eq. (3) in different pixels were uncorrelated, the variance of ξ_{A} would simply be the weighted products of the fluctuation variances. This yields a result that is ~30% smaller than the true correlation variance that we compute below. The reason is, of course, that the δpairs are correlated, either from LSS or from correlations induced by instrumental effects or continuum subtraction; this effect reduces the effective number of pairs and introduces correlations between (r,μ) bins.
Rather than determine the full covariance matrix for ξ(r,μ), we determined directly the covariance matrix for ξ_{0}(r) and ξ_{2}(r) by standard techniques of dividing the full quasar sample into subsamples according to position on the sky. In particular we used the subsampling technique described below. We also tried a bootstrap technique (e.g. Efron & Gong 1983) consisting of substituting the entire set of N subdivisions of the data by N of these subdivisions chosen at random (with replacement) to obtain a “bootstrap” sample. The covariances are then measured from the ensemble of bootstrap samples. Both techniques give consistent results.
The adopted covariance matrix for the monopole and quadrupole uses the subsampling technique. We divide the data into angular sectors and calculate a correlation function in each sector. Pairs of pixels belonging to different sectors contribute only to the sector of the pixel with lower right ascension. We investigated two different divisions of the sky data: defining 800 (contiguous but disjoint) sectors of similar solid angle, and taking the plates as defining the sectors (this latter version does not lead to disjoint sectors). The two ways of dividing the data lead to similar covariance matrices.
Fig. 10
rdependence of the product of the monopole (top) and quadrupole (middle) variances and the number of pairs in the rbin. The bottom panel shows the product for the monopolequadrupole covariance (r = r′ elements). The dotted lines show the means for r > 20 h^{1} Mpc. 
Each sector s in each division of the data provides a measurement of ξ_{s}(r,μ) that can be used to derive a monopole and quadrupole, ξ_{ℓs}(r), (ℓ = 0,2). The covariance of the whole BOSS sample can then be estimated from the weighted and rescaled covariances for each sector: (18)The average denoted by ⟨ ⟩ is the simple average over sectors, while denotes the correlation function measured for the whole BOSS sample. The W_{s}(r) are the summed pixelpair weights for the radial bin r for the sector s and W(r) is the same sum for the whole BOSS sample.
The most important terms in the covariance matrix are the r = r′ terms, i.e. the monopole and quadrupole variances. They are shown in Fig. 10 as a function of r. In the figure, they are multiplied by the number N of pixel pairs in the rbin. The product is nearly independent of r, as expected for a variance nearly equal to the pixel variance divided by N. For the monopole, the variances are only about 30% higher than what one would calculate naively assuming uncorrelated pixels and Eq. (12). Figure 10 also displays the monopolequadrupole covariance times number of pairs, which also is nearly independent of r.
Figure 11 displays the monopolemonopole and quadrupolequadrupole covariances. Nearestneighbor covariances are of order 20%. Figure 11 also shows monopolequadrupole covariance.
We used the 15 sets of mock spectra to test our method for calculating the covariance matrix. From the 15 measurements of ξ_{ℓ}(r) one can calculate the average values of ξ_{ℓ}(r)ξ_{ℓ′}(r′) and compare them with those expected from the covariance matrix. Figures 12 shows this comparison for the monopole and quadrupole variance, the monopole and quadrupole covariances between neighboring rbins and the monopolequadrupole covariance. The agreement is satisfactory.
Fig. 11
Monopole and quadrupole covariance matrix. The monopolemonopole and quadrupolequadrupole elements are normalized to the variance: . The monopolequadrupole elements are normalized to the mean of the quadrupole and monopole variances. The first offdiagonal elements of the monopolemonopole and quadrupolequadrupole elements are ~20% of the diagonal elements. The diagonal elements of the quadrupolemonopole covariance are ~− 0.2 times the geometric mean of the monopole and quadrupole variances. 
Fig. 12
Verification of the offdiagonal elements of the covariance matrix with the 15 sets of mock spectra. The black lines show correlations derived from the dispersion of the 15 measurements and the red lines show the correlations expected from the covariance matrix calculated by subsampling. The top and middle panels show the correlation between neighboring bins for monopole and quadrupole respectively. The bottom panel, the correlation between monopole and quadrupole measured at the same distance bin. 
4.2. Detection significance of the BAO peak
In this section, we estimate the significance of our detection of a BAO peak at 105 h^{1} Mpc. At the statistical power of the present data, it is clear that the peak significance will depend to some extent on how we treat the socalled “broadband” correlation function on which the peak is superimposed. In particular, the significance will depend strongly on the rrange over which the correlation functions are fitted. To the extent that the BAO peak is known to be present in the matter correlation function and that the Lyα absorption is known to trace matter, the actual significance is of limited interest for cosmology. Of greater interest is the uncertainty in the derived cosmological parameter constraints (Sect. 5) which will be nonlinear reflections of the peak significance derived here.
A detection of the BAO peak requires comparing the quality of a fit with no peak (the null hypothesis) to that of a fit with a peak. Typically, this exercise would be performed by choosing a test statistic, such as the χ^{2}, computing the distribution of this quality indicator from a large number of peakless simulations and looking at the consistency of the data with this distribution. Since our mock data sets are quite computationally expensive and only a handful are available, we chose a different approach.
Our detection approach uses the following expression to fit the observed monopole and quadrupole. (19)where is a broadband term to describe the LSS correlation function in the absence of a peak, is a peak term, and is a “distortion” term used to model the effects of continuum subtraction. The broadband term is derived from the fiducial ΛCDM cosmology defined by the parameters in Eq. (A.1). It is obtained by fitting the shape of the fiducial correlation function with an 8node spline function masking the region of the peak (80 h^{1} Mpc < r < 120 h^{1} Mpc). The peak term is the difference between the theoretical correlation function and the broadband term. Finally, the distortion term is calculated from simulations, as the difference in the monopole or quadrupole measured using the true continuum and that measured from fitting the continuum as described in Appendix A. The three components are shown in Fig. 13.
Fig. 13
Fitting functions used for the determination of the peak detection significance: r^{2}ξ^{bb}(r) (blue), r^{2} [ξ^{peak}(r) + ξ^{bb}(r)] (black) and r^{2}ξ^{dist}(r) (red) for the monopole (solid lines) and quadrupole (dashed lines). 
Expression (19) contains three parameters each for the monopole and quadrupole (so six in total). We have performed fits leaving all six parameters free and fits where we fix the ratio C_{2}/C_{0} to be equal to its nominal value used to generate our mock spectra (the value given by assuming a “redshift distortion parameter” β = 1.4, see Appendix A). We define the test statistic as the χ^{2} difference between fitting Eq. (19) simultaneously to monopole and quadrupole by fixing C_{0} to zero (a “peakless” four or fiveparameter fit) and fitting for C_{0} (a five or sixparameter fit). In our detection fits we do not fit for the BAO position but fix it to the theoretical prediction. The distribution for this test statistic (“”) under the null hypothesis is a χ^{2} distribution with one degree of freedom. The significance is then given by .
Figure 14 shows the fits to monopole (top panel) and quadrupole (bottom panel) and the corresponding fits with and without peaks and fixing β = 1.4. For method 2, we obtain χ^{2}/d.o.f. = 93.7/85 (111.8/86) with (without) a peak, giving Δχ^{2} = 18.1 for a detection significance of 4.2σ. For method 1, we obtain χ^{2}/d.o.f. = 93.2/85 (102.2/86) for a significance of 3.0σ. Allowing β to be a free parameter gives essentially the same detection significances.
The detection significance of ~4σ is typical of that which we found in the 15 sets of mock spectra. For the mocks, the significances ranged from 0 to 6σ with a mean of 3.5σ.
Our significance depends strongly on the fitting range. For a lower boundary of the range of r_{min} = 20,40,60 h^{1} Mpc we obtain a significance of σ = 4.2, 3.2 and 2.3, respectively (method 2). The reason for this result is illustrated in Fig. 15, where the results of the fits with and without peaks are compared to data for different values of r_{min} = 20,40,60 h^{1} Mpc. Reducing the fitting range poses less stringent constraints on the distortion and broadband terms, thus allowing some of the peak to be attributed to the broadband. In particular, the statistically insignificant bump in the quadrupole at ~65 h^{1} Mpc causes the fitted broadband to increase as r_{min} is increased to 60 h^{1} Mpc, decreasing the amplitude of the BAO peak. For the monopole, the r_{min} = 60 h^{1} Mpc fit predicts a positive slope for ξ_{0}(r) that decreases the amplitude of the peak but predicts a ξ_{0}(r < 50 h^{1} Mpc) to be much less than what is measured.
Fig. 14
Monopole and quadrupole fits with a BAO peak (red line) and without a BAO peak (blue line, Method 2). The fitting range is r_{min} < r < 200 h^{1} Mpc with r_{min} = 20 h^{1} Mpc. 
Fig. 15
Same as Fig. 14 except for different fitting ranges: blue, green and red curves are for fits with r_{min} = 20 h^{1} Mpc, 40 h^{1} Mpc and 60 h^{1} Mpc respectively (Method 2). The solid lines for fits without a BAO peak and the dashed lines with a peak. 
5. Cosmology with the BAO peak
The observed position of the BAO peak in ξ(r,μ) is determined by two sets of cosmological parameters: the “true” parameters and the “fiducial” parameters. Nature uses the true cosmology to create correlations at the true sound horizon, r_{s}. The true cosmology transforms physical separations between Lyα absorbers into angles on the sky and redshift differences: θ_{BAO} = r_{s}/D_{A}(z)(1 + z) and Δz_{BAO} = r_{s}H(z)/c. We, on the other hand, use a “fiducial” cosmology (defined by Eq. (A.1)) to transform angular and redshift differences to local distances at the redshift in question to reconstruct ξ(r,μ). If the fiducial cosmology is the true cosmology, the reconstructed peaks will be at the calculated fiducial sound horizon, r_{s,f}. Limits on the difference between the fiducial and reconstructed peak position can be used to constrain the differences between the fiducial and true cosmological models.
5.1. The peak position
The use of incorrect fiducial D_{A}, H and r_{s} leads to shifts in the BAO peak position in the transverse and radial directions by the multiplicative factors α_{t} and α_{r}: where the subscript f refers to the fiducial model. Following Xu et al. (2012), we will use a fitting function, , for the monopole and quadrupole that follows the expected peak position as a function of (α_{t},α_{r}): (22)Here, the two functions , derived from the power spectrum given in Appendix A, describe the underlying mass correlation function, the linear bias b and redshift distortion parameter β, and the movement of the BAO peak for α_{t},α_{r} ≠ 1. The functions A_{ℓ} take into account distortions, as described below.
For α_{t} = α_{r} ≡ α_{iso}, there is a simple isotropic scaling of the coordinates by α_{iso} and is given by , where ξ_{ℓ,f} are the fiducial monopole and quadrupole and the normalizations f_{ℓ} are the functions of the bias and redshiftdistortion parameter given by Hamilton (1992). For α_{t} ≠ α_{r}, Xu et al. (2012) found an approximate formula for that was good in the limit  α_{t} − α_{r} ≪ 1. We take the more direct route of numerically expanding ξ_{f}(α_{t}r_{t},α_{r}r_{r}) in Legendre polynomials, P_{ℓ}(μ), to directly calculate the .
The functions A_{ℓ}(r) describe broadband distortions due to continuum subtraction and the fact that the broadband correlation function is not expected to change in the same way as the BAO peak position when one deviates from the fiducial model. They correspond to the term in Eq. (19). We have used two forms to represent A_{ℓ}(r): (23)and (24)The observed monopole and quadrupole can then be fit to Eq. (22) with free parameters α_{t}, α_{r}, bias, β, and the nuisance parameters (a_{ℓ},b_{ℓ},c_{ℓ} and d_{ℓ}).
We first fixed (α_{t},α_{r}) = (1,1) to determine if we find reasonable values of (b,β). These two parameters are highly degenerate since both the quadrupole and monopole have amplitudes that are proportional to b^{2} times polynomials in β. A welldetermined combination is b(1 + β), for which we find a value 0.38 ± 0.07; this is in agreement with b(1 + β) = 0.336 ± 0.012 found at r ~ 40 h^{1} Mpc by Slosar et al. (2011). The larger error of our fit reflects the substantial freedom we have introduced with our distortion function.
We next freed all parameters to constrain (α_{t},α_{r}). The contours for the two methods and two broadbands are shown in Fig. 16 and the χ^{2} for the fiducial and bestfit models are given in Table 1. The broadband term in Eq. (24) fits the data better than that in Eq. (23) both for the fiducial parameters and for the best fit. For broadband in Eq. (24), the χ^{2} for the fiducial model is acceptable for both methods: χ^{2}/d.o.f. = 85.0/80 for method 1 and χ^{2}/d.o.f. = 71.5/80 for method 2.
The contours in the figure are elongated along the direction for which the BAO peak position stays approximately fixed at large μ (near the radial direction, where the observations are most sensitive). The best constrained combination of D_{A} and H of the form turns out to have ζ ~ 0.2. This low value of ζ reflects the fact that we are mostly sensitive to the BAO peak in the radial direction. At the one standarddeviation level, the precision on this combination is about 4%. However, even this combination is sensitive to the tails in the contours. A more robust indicator of the statistical accuracy of the peakposition determination comes from fits imposing α_{t} = α_{r} ≡ α_{iso}, as has generally been done in previous BAO studies with the exception of Chuang & Wang (2012) and Xu et al. (2012). This constraint does not correspond to any particular class of cosmological models. It does however eliminate the tails in the contours in a way that is similar to the imposition of outside data sets. The two methods and broadbands give consistent results, as seen in Table 1.
We used the sets of mock spectra to search for biases in our measurement of α_{iso}. The mean value reconstructed for this quantity on individual mocks is 1.002 ± 0.007, suggesting that there are no significant biases in the determination of the BAO scale. Figure 17 shows the values and errors for the individual mocks along with that for the data. Both the measured value and its uncertainty for the data is typical of that found for individual sets of mock spectra.
Fig. 16
Contours for (D_{A}/r_{s},r_{s}H) obtained by fitting the monopole and quadrupole to (22). The broadband distortions are Eq. (23) – dashed lines or Eq. (24) – solid lines. The blue lines are for method 1 and the red lines for method 2. All contours are for Δχ^{2} = 4 except for the interior solid red contour which is for Δχ^{2} = 1. 
Fig. 17
Measurements of α_{iso} (= α_{t} = α_{r}) for the 15 sets of mock spectra and for the data (realization = −1). The large errors for realization 5 and 8 are due to the very low significance of the BAO peak found on these two sets. 
Fig. 18
Constraints on the matter and darkenergy density parameters (Ω_{M},Ω_{Λ}) assuming a darkenergy pressuredensity ratio w = −1. The blue regions are the one and two standard deviation constraints derived from our contours in Fig. 16 (method 2, broadband 24) combined with a measurement of H_{0} (Riess et al. 2011). Also shown are one and two standard deviation contours from lower redshift measurements of D_{V}/r_{s} (also combined with H_{0}) at z = 0.11 (6dF: Beutler et al. 2011), z = 0.35 (LRG: Percival et al. 2010) and z = 0.57 (CMASS: Anderson et al. 2012). All constraints use a WMAP7 (Komatsu et al. 2011) prior on the baryontophoton ratio η but do not otherwise incorporate CMB results. 
Fig. 19
As in Fig. 18 with constraints on the matter density parameter, Ω_{M}, and darkenergy pressuredensity ratio w assuming Ω_{M} + Ω_{Λ} = 1. 
5.2. Constraints on cosmological models
Our constraints on (D_{A}/r_{s},Hr_{s}) can be used to constrain the cosmological parameters. In a ΛCDM cosmology, apart from the prefactors of H_{0} that cancel, D_{A}/r_{s} and Hr_{s} evaluated at z = 2.3 depend primarily on Ω_{M} through r_{s} and on Ω_{Λ} which, with Ω_{M}, determines D_{A} and H. The sound horizon also depends on H_{0} (required to derive Ω_{γ} from T_{CMB}), on the effective number of neutrino species N_{ν} (required to derive the radiation density from the photon density), and on the baryontophoton number ratio, η (required for the speed of sound).
Figure 18 shows the ΛCDM constraints on (Ω_{M},Ω_{Λ}) derived from the contours in Fig. 16 combined with the most recent measurement of H_{0} (Riess et al. 2011). We use the contours for method 2 and the broadband of Eq. (24) which gives better fits to the data than the other method and broadband. The contours also assume N_{ν} = 3 and the WMAP7 value of η (Komatsu et al. 2011). Also shown are constraints from BAO measurements of D_{V}/r_{s} (Percival et al. 2010; Anderson et al. 2012; Beutler et al. 2011).
The Lyα contours are nicely orthogonal to the lower redshift D_{V}/r_{s} measurements, reinforcing the requirement of dark energy from BAO data. In fact, our measurement is the only BAO measurement that by itself requires dark energy: Ω_{Λ} > 0.5. This is because at z = 2.3 the universe is strongly matter dominated and the factor in H partially cancels the in r_{s}, enhancing the importance of the Ω_{Λ} dependence of H.
Figure 19 shows the constraints on (Ω_{M},w; where w is the darkenergy pressuredensity ratio) assuming a flat universe: Ω_{k} = 0. Our result is the only BAO measurement that by itself requires negative w. Our limit w < −0.6 requires matter domination at z = 2.3. (25)
5.3. Constraints on H(z)
The contours in Fig. 16 give the measurements of Hr_{s} given in Table 1. A measurement of the expansion rate deep in the matterdominated epoch can be used to demonstrate the deceleration of the expansion at that time. Unfortunately, our data are not yet precise enough to do this. To make a more precise measurement of H(z = 2.3), we must add further constraints to eliminate the long tails in Fig. 16. These tails correspond to models where 1/H(z = 2.3) is increased (resp. decreased) with respect to the fiducial value while D_{A}(z = 2.3) is decreased (resp. increased). For flat models, this would imply a change in the mean of 1/H (averaged up to z = 2.3) that is opposite to that of the change in 1/H(z = 2.3), which requires a functional form H(z) that strongly differs from the fiducial case. It is possible to construct models with this property by introducing significant nonzero curvature.
Because of the importance of curvature, the tails are eliminated once WMAP7 constraints (Komatsu et al. 2011) are included. This is done in Fig. 20 within the framework of nonflat models where the darkenergy pressuredensity ratio, w(z), is determined by two parameters, w_{0} and w_{a}: w(z) = w_{0} + w_{a}z/(1 + z). As expected, the WMAP7 results in this framework constrain D_{A} and 1/H to migrate in roughly the same direction as one moves away from the fiducial model. Combining WMAP7 constraints with ours gives the values of H(z = 2.3)r_{s} given in the last column of Table 1. For what follows, we adopt the mean of methods 1 and 2 that use the more flexible broadband of Eq. (24): (26)The precision on H is now sufficient to study the redshift evolution of H(z).
The fiducial model has r_{s} = 152.76 Mpc and H(z = 2.3) = 3.23H_{0}, H_{0} = 70 km s^{1} Mpc^{1}. These results produce (27)or equivalently (28)This number can be compared with the measurements of H(z) at lower redshift shown in Table 2 and Fig. 21. Other than those of H_{0}, the measurements that we use can be divided into two classes: those (like ours) that use r_{s} as the standard of length and those that use c/H_{0} as the standard of length.
Fig. 20
Constraints on (D_{A}/r_{s},r_{s}H)_{z = 2.3} within the framework of OwOwaCDM models. The green contours are our 1σ and 2σ constraints using method 2 and broadband (24). The gray contours are the 1σ and 2σ constraints from WMAP7 (Komatsu et al. 2011). The red contours show the combined constraints. 
Recent measurements of H(z)/(1 + z).
The comparison with our measurement is simplest with BAObased measurements that use r_{s} as the standard of length and therefore measure H(z)r_{s} (as is done here). The first attempt at such a measurement was made by Gaztañaga et al. (2009), a result debated in subsequent papers by MiraldaEscudé (2009), Yoo & MiraldaEscudé (2010), Kazin et al. (2010), and Cabré & Gaztañaga (2011). Here, we use four more recent measurements. Chuang & Wang (2012) and Xu et al. (2012) studied the SDSS DR7 LRG sample and decomposed the BAO peak into radial and angular components, thus extracting directly Hr_{s} and D_{A}/r_{s}. Blake et al. (2012) and Reid et al. (2012) took a more indirect route. They first used the angleaveraged peak position to derive . They then studied the AlcockPaczynski effect on the broadband galaxy correlation function to determine D_{A}(z)H(z). Combining the two measurements yielded H(z)r_{s}.
Fig. 21
Measurements of H(z)/(1 + z) vs. z demonstrating the acceleration of the expansion for z < 0.8 and deceleration for z > 0.8. The BAObased measurements are the filled circles: this work: (red), Xu et al. (2012, black)Chuang & Wang (2012, blue), Reid et al. (2012, cyan), and Blake et al. (2012, green). The open green circles are from WiggleZ (Blake et al. 2011b) AlcockPaczynski data combined with supernova data yielding H(z)/H_{0} (without the flatness assumption) plotted here assuming H_{0} = 70 km s^{1} Mpc^{1}. The open blue circle is the H_{0} measurement of Riess et al. (2011). The open black squares with dashed error bars show the results of Riess et al. (2007) which were derived by differentiating the SNIa Hubble diagram and assuming spatial flatness. (For visual clarity, the Riess et al. 2007 point at z = 0.43 has been shifted to z = 0.48.) The line is the ΛCDM prediction for (h,Ω_{M},Ω_{Λ}) = (0.7,0.27,0,73). 
It is evident from comparing our H(z) measurement (filled red circle in Fig. 21) to the other BAObased measurements (other filled circles) that H(z)/(1 + z) decreases between z = 2.3 and z = 0.35−0.8. To demonstrate deceleration quantitatively, we fit the eight BAObased values of H(z) in Table 2 to the oΛCDM form H(z) = H_{0}(Ω_{Λ} + Ω_{M}(1 + z)^{3} + (1−Ω_{Λ} − Ω_{M})(1 + z)^{2})^{1/2}. Marginalizing over Ω_{Λ} and H_{0} we find (29)clearly indicating deceleration between z = 2.3 and z = 0.5. This measurement is in good agreement with the fiducial value of 1.146. We emphasize that this result is independent of r_{s}, assuming only that the BAOpeak position is redshiftindependent in comoving coordinates. The result also does not assume spatial flatness.
To map the expansion rate over the full range 0 < z < 2.3, we must adopt the fiducial value of r_{s} and compare the resulting H(z) with H_{0} and with other BAOfree measurements. Besides the H_{0} measurement of Riess et al. (2011), we use the WiggleZ analysis combining their AlcockPaczynski data with distant supernova data from the Union2 compilation (Amanullah et al. 2010). The supernova analysis does not use the poorly known mean SNIa luminosity, so the SNIa Hubble diagram gives the luminosity distance in units of , D_{L}(z)H_{0}. Combining this result with the AlcockPaczynski measurement of D_{A}(z)H(z) yields H(z)/H_{0}. The values are given in Table 2.
We fit all the data in Table 2 (filled and open circles in Fig. 21) to the ΛCDM form of H(z). This yields an estimate of the redshift of minimum H(z)/(1 + z) (30)which compares well with the fiducial value: z_{d − a} = (2Ω_{Λ}/Ω_{M})^{1/3} − 1 = 0.755.
In this analysis, we have not used two other sources of information on H(z) at high redshift. The first use highredshift type Ia supernovae to probe the era where the universe transitions from deceleration to acceleration (e.g., Riess et al. 2004, 2007). The data of Riess et al. (2007) (plotted as the open squares in Fig. 21) yielded useful measurements up to z ~ 1.1. However, this data yields constraints on H(z) that are weaker than those of BAObased methods because of the need to differentiate the distanceredshift relation. Moreover, these inferences of H(z) assume spatial flatness. Fitting the SNe data to a model with an evolving deceleration parameter q(z) = q_{0} + (dq/dz)_{0}z and assuming flatness, Riess et al. (2007, 2004) were able to demonstrate that (dq/dz)_{0} > 0, i.e. a negative 3rdderivative of a(t). However, we point out that in a more general q(z) model, the demonstration that dq/dz > 0 at low redshift is not equivalent to a demonstration that becomes negative in the past.
Another approach to determining H(z) uses the evolution of stellar populations as a clock to infer dt/dz (Stern et al. 2010; Moresco et al. 2012). This method yields results that are consistent with ΛCDM expectations, but the uncertainties (statistical and systematic) are larger than those of the determinations in Table 2, so we have not plotted them in Fig. 21.
6. Conclusions
In this paper, we have presented the first observation of the BAO peak using the Lyα forest. It represents both the first BAO detection deep in the matter dominated epoch and the first to use a tracer of mass that is not galactic. The results are consistent with concordance ΛCDM, and require, by themselves, the existence of dark energy. Combined with CMB constraints, we deduce the expansion rate at z = 2.3 and demonstrate directly the sequence of deceleration and acceleration expected in darkenergy dominated cosmologies. These results have been confirmed with higher precision by Slosar et al. (2013) using the same underlying DR9 data set but more aggressive data cuts and a more nearly optimal statistical method.
BOSS continues to acquire data and will eventually produce a quasar sample three times larger than DR9. We can thus expect improved precision in our measurements of distances and expansion rates, leading to improved constraints on cosmological parameters. The Lyα forest may well be the most practical method for obtaining precise D_{A}(z) and H(z) measurements at z > 2, thanks to the large number of independent density measurements per quasar. It is reassuring that the first sample large enough to yield a detection of BAO produces a signal in good agreement with expectations. In the context of BAO dark energy constraints, high redshift measurements are especially valuable for breaking the degeneracy between curvature and the equation of state history More generally, however, by probing an epoch largely inaccessible to other methods, BAO in the Lyα forest have the potential to reveal surprises, which could provide critical insights into the origin of cosmic acceleration.
Acknowledgments
We thank Carlos Allende Prieto, Ashley Ross and Uros Seljak for stimulating discussions and Adam Riess for providing the data points of Riess et al. (2007) in Fig. 21. Funding for SDSSIII has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the US Department of Energy Office of Science. The SDSSIII web site is http://www.sdss3.org/. The French Participation Group of SDSSIII was supported by the Agence Nationale de la Recherche under contract ANR08BLAN0222. The research leading to these results has received funding from the European Union Seventh Framework Programme (FP7/20072013) under grant agreement No. PIIFGA2011301665. SDSSIII is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSSIII Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, University of Cambridge, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University.
References
 Abazajian, K. N., AdelmanMcCarthy, J. K., Agüeros, M. S., et al. 2009, ApJS, 182, 543 [NASA ADS] [CrossRef] [Google Scholar]
 Ahn, C. P., Alexandroff, R., Allen de Prieto, C., et al. 2012, ApJS, 203, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Alcock, C., & Paczynski, B. 1979, Nature, 281, 358 [NASA ADS] [CrossRef] [Google Scholar]
 Amanullah, R., Lidman, C., Rubin, D., et al. 2010, ApJ, 716, 712 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, L., Aubourg, E., Bailey, S., et al. 2012, MNRAS, 427, 3435 [Google Scholar]
 Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559 [NASA ADS] [CrossRef] [Google Scholar]
 Beutler, F., Blake, C., Colless, M., et al. 2011, MNRAS, 416, 3017 [NASA ADS] [CrossRef] [Google Scholar]
 Bi, H. G., Börner, G., & Chu, Y. 1992, A&A, 266, 1 [NASA ADS] [Google Scholar]
 Blake, C., Davis, T., Poole, G. B., et al. 2011a, MNRAS, 415, 2892 [NASA ADS] [CrossRef] [Google Scholar]
 Blake, C., Glazebrook, K., Davis, T. M., et al. 2011b, MNRAS, 418, 1725 [NASA ADS] [CrossRef] [Google Scholar]
 Blake, C., Brough, S., Colless, M., et al. 2012, MNRAS, 425, 405 [NASA ADS] [CrossRef] [Google Scholar]
 Blanton, M. R., Lin, H., Luptonn, R. H., et al. 2003, AJ, 125, 2276 [NASA ADS] [CrossRef] [Google Scholar]
 Bolton, A., Schlegel, D. J., Aubourg, É., et al. 2012, AJ, 144, 144 [NASA ADS] [CrossRef] [Google Scholar]
 Bovy, J., Hennawi, J. F., Hogg, D. W., et al. 2011, ApJ, 729, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Bovy, J., Myers, A. D., Hennawi, J. F., et al. 2012, ApJ, 749, 41 [NASA ADS] [CrossRef] [Google Scholar]
 Cabré, A., & Gaztaǹaga, E. 2011, MNRAS, 412, L98 [NASA ADS] [Google Scholar]
 Caucci, S., Colombi, S., Pichon, C., et al. 2008, MNRAS, 386, 211 [NASA ADS] [CrossRef] [Google Scholar]
 Cen, R., MiraldaEscudé, J., Ostriker, J. P., & Rauch, M. 1994, ApJ, 437, L9 [Google Scholar]
 Chuang, ChiaHsun, & Yun Wang 2012, MNRAS, 426, 226 [NASA ADS] [CrossRef] [Google Scholar]
 Cole, S., Percival, W. J., Peacock, J. A., et al. 2005, MNRAS, 362, 505 [NASA ADS] [CrossRef] [Google Scholar]
 Coles, P., & Jones, B. 1991, MNRAS, 248, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Croft, R., Weinberg, D. H., Katz, N., & Hernquist, L. 1997, ApJ, 488, 532 [NASA ADS] [CrossRef] [Google Scholar]
 Croft, R., Weinberg, D. H., Katz, N., & Hernquist, L. 1998, ApJ, 495, 44 [NASA ADS] [CrossRef] [Google Scholar]
 Croft, R. A. C., Weinberg, D. H., Pettini, M., et al. 1999, ApJ, 520, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Croft, R. A. C., Weinberg, D., Bolte, M., et al. 2002, ApJ, 581, 20 [NASA ADS] [CrossRef] [Google Scholar]
 Dawson, K., Schlegel, D., Ahn, C., et al. 2013, AJ, 145, 10 [Google Scholar]
 de Bernardis, P., Ade, P. A. R., Bock, J. J., et al. 2000, Nature, 404, 955 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Eisenstein, D. J. 2007, ApJ, 664, 675 [Google Scholar]
 Eisenstein, D. J., Zehavi, I., Hogg, D. W., et al. 2005, ApJ, 633, 560 [NASA ADS] [CrossRef] [Google Scholar]
 Eisenstein, D. J., Weinberg, D. H., Agol, E., et al. 2011, AJ, 142, 72 [Google Scholar]
 Efron, B., & Gong, G. 1983, The American Statistician, 37, 36 [Google Scholar]
 FontRibera, A., & MiraldaEscud, J. 2012, J. Cosmol. Astropart. Phys., 07, 028 [NASA ADS] [CrossRef] [Google Scholar]
 FontRibera, A., McDonald, P., & Escudé, J. M. 2012, J. Cosmol. Astropart. Phys., 01, 001 [NASA ADS] [CrossRef] [Google Scholar]
 Fukugita, M., Ichikawa, T., Gunn, J. E., et al. 1996, AJ, 111, 1748 [NASA ADS] [CrossRef] [Google Scholar]
 Gallerani, S., Kitaura, F. S., & Ferrara, A. 2011, MNRAS, 413, L6 [NASA ADS] [Google Scholar]
 Gaztañaga, E., Cabré, A., & Hui, L. 2009, MNRAS, 399, 1663 [NASA ADS] [CrossRef] [Google Scholar]
 Gnedin, N., & Hui, L. 1996, ApJ, 472, L73 [NASA ADS] [CrossRef] [Google Scholar]
 Greig, B., Bolton, J., Wyithe, J., & Stuart, B. 2011, MNRAS, 418, 1980 [NASA ADS] [CrossRef] [Google Scholar]
 Gunn, J. E., Carr, M., Rockosi, C., et al. 1998, AJ, 116, 3040 [NASA ADS] [CrossRef] [Google Scholar]
 Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, AJ, 131, 2332 [NASA ADS] [CrossRef] [Google Scholar]
 Hamilton, A. J. S. 1992, ApJ, 385, L5 [NASA ADS] [CrossRef] [Google Scholar]
 Hernquist, L., Katz, N., Weinberg, D. H., & MiraldaEscudé, J. 1996, ApJ, 457, L51 [NASA ADS] [CrossRef] [Google Scholar]
 Kazin, E. A., Blanton, M. R., Scoccimarro, R., et al. 2010, ApJ, 719, 1032 [Google Scholar]
 Kirkpatrick, J. A., Schlegel, D. J., Ross, N. P., et al. 2011, AJ, 743, 125 [Google Scholar]
 Kitaura, F.S., Gallerani, S., & Ferrara, A. 2012, MNRAS, 420, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Komaatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Lawrence, A., Warren, S. J., Almaini, O., et al. 2007, MNRAS, 379, 1599 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Lee, K.G., Suzuki, N., & Spergel, D. N. 2012a, AJ, 143, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Lee, K.G., Bailey, S., Bartsch, L. E., et al. 2012b, AJ, submitted [arXiv:1211.5146] [Google Scholar]
 Le Goff, J. M., Magneville, C., Rollinde, E., et al. 2011, A&A, 534, 135 [Google Scholar]
 Lidz, A., FaucherGiguére, C.A., Dall’Aglio, A., et al. 2010, ApJ, 718, 199 [NASA ADS] [CrossRef] [Google Scholar]
 Martin, D. C., Fansonn, J., Schiminovich, D., et al. 2005, ApJ, 619, L1 [Google Scholar]
 McDonald, P. 2003, ApJ, 585, 34 [NASA ADS] [CrossRef] [Google Scholar]
 McDonald, P., & Eisenstein, D. 2007, Phys. Rev. D, 76, 063009 [NASA ADS] [CrossRef] [Google Scholar]
 McDonald, P., White, M., Cohn, J. D., et al. 2000, ApJ, 543, 1 [NASA ADS] [CrossRef] [Google Scholar]
 McDonald, P., Seljak, U., Schlegel, D. J., et al. 2006, ApJS, 163, 80 [NASA ADS] [CrossRef] [Google Scholar]
 McQuinn, M., & White, M. 2011, MNRAS, 415, 2257 [NASA ADS] [CrossRef] [Google Scholar]
 Mehta, K. T., Cuesta, A., Xu, X., et al. 2012, MNRAS, 427, 2168 [NASA ADS] [CrossRef] [Google Scholar]
 MiraldaEscudé, J. 2009 [arXiv:0901.1219] [Google Scholar]
 MiraldaEscudé, J., Cen, R., Ostriker, J. P., & Rauch, M. 1996, ApJ, 471, 582 [NASA ADS] [CrossRef] [Google Scholar]
 Moresco, M., Cimatti, A., Jimenez, R., et al. 2012, J. Cosmol. Astropart. Phys., 08, 006 [Google Scholar]
 Noterdaeme, P., Petitjean, P., Carithers, W. C., et al. 2012, A&A, 547, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nusser, A., & Haehnelt, M. 1999, MNRAS, 303 179 [Google Scholar]
 Padmanabhan, N., Schlegel, D. J., Finkbeiner, D. P., et al. 2008, ApJ, 674, 1217 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Padmanabhan, N., White, M., Cohn, J. D. 2009, Phys. Rev. D, 79, 3523 [CrossRef] [Google Scholar]
 Padmanabhan, N., Xu, X., Eisenstein, D. J., et al. 2012, MNRAS, 427, 2132 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 PalanqueDelabrouille, N., Yèche, C., Myers, A. D., et al. 2011, A&A, 530, 122 [Google Scholar]
 Pâris, I., Petitjean, P., Rollinde, E., et al. 2011, A&A, 530, 50 [Google Scholar]
 Pâris, I., Petitjean, P., Aubourg, E., et al. 2012, A&A, 548, 66 [Google Scholar]
 Percival, W. J., Reid, B. A., Eisenstein, D. J., et al. 2010, MNRAS, 401, 2148 [NASA ADS] [CrossRef] [Google Scholar]
 Petitjean, P., Mücket, J. P., & Kates, R. E. 1995, A&A, 295, L9 [NASA ADS] [Google Scholar]
 Pichon, C., Vergely, J. L., Rollinde, E., et al. 2001, MNRAS, 326, 597 [NASA ADS] [CrossRef] [Google Scholar]
 Pier, J. R., Munn, J. A., Hindsley, R. B., et al. 2003, AJ, 125, 1559 [NASA ADS] [CrossRef] [Google Scholar]
 Pieri, M. M., Frank, S., Weinberg, D. H., et al. 2010, ApJ, 724, L69 [NASA ADS] [CrossRef] [Google Scholar]
 Reid, B. A., Samushia, L., White, M., et al. 2012, MNRAS, 426, 2719 [NASA ADS] [CrossRef] [Google Scholar]
 Richards, G. T., Myers, A. D., Gray, A. G., et al. 2009, ApJS, 180, 67 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Strolger, L.G., Tonry, J., et al. 2004, ApJ, 607, 665 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Strolger, L.G., Casertano, S., et al. 2007, ApJ, 659, 98 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Macri, L., Casertano, S., et al. 2011, ApJ, 730, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Ross, N. P., Myers, A. D., Sheldon, E. S., et al. 2012, ApJS, 199, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Seo, H.J., Ho, S., White, M., et al. 2012, ApJ, 761, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Slosar, A., FontRibera, A., Pieri, M. M., et al. 2011, J. Cosmol. Astropart. Phys., 09, 001 [Google Scholar]
 Slosar, A., Iršič, V., Kirkby, D., et al. 2013, JCAP, submitted [arXiv:1301.3459] [Google Scholar]
 Smee, S., Gunn, J. E., Uomoto, A., et al. 2012, AJ, submitted [arXiv:1208.2233] [Google Scholar]
 Smith, J. A., Tucker, D. L., Kent, S., et al. 2002, AJ, 123, 2121 [Google Scholar]
 Stern, D., Jimenez, R., Verde, L., et al. 2010, J. Cosmol. Astropart. Phys., 02, 008 [Google Scholar]
 Stoughton, C., Lupton, R. H., Bernardi, M., et al. 2002, AJ, 123, 485 [NASA ADS] [CrossRef] [Google Scholar]
 Suzuki, N. 2006, ApJS, 163, 110 [NASA ADS] [CrossRef] [Google Scholar]
 Suzuki, N., Tytler, D., Kirkman, D., et al. 2005, ApJ, 618, 592 [NASA ADS] [CrossRef] [Google Scholar]
 Theuns, T., Leonard, A., Efstathiou, G., Pearce, F. R., & Thomas, P. A. 1998, MNRAS, 301, 478 [NASA ADS] [CrossRef] [Google Scholar]
 Viel, M., Haehnelt, M. G., & Springel, V. 2004, MNRAS, 354, 684 [NASA ADS] [CrossRef] [Google Scholar]
 Weinberg, D. H., Katz, N., Hernquist, L., et al. 1998, ASPC, 148, 21 [NASA ADS] [Google Scholar]
 Weinberg, D. H., Mortonson, J. J., Eisenstein, D. J., et al. 2012 [arXiv:1201.2434] [Google Scholar]
 White, M. 2003, in The Davis Meeting On Cosmic Inflation 2003, March 22−25, Davis CA., 18 [Google Scholar]
 Worseck, G., & Wisotzki, L. 2006, A&A, 450, 495 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Xu, X., Cuesta, A. J., Padmanabhan, N., et al. 2012, MNRAS, submitted [arXiv:1206.6732] [Google Scholar]
 Yèche, C., Petitjean, P., Rich, J., et al. 2009, A&A, 523, A14 [Google Scholar]
 Yoo, J., & MiraldaEscudé, J. 2010, Phys. Rev. D, 82, 3527 [NASA ADS] [Google Scholar]
 York, D. G., Adelman, J., Anderson, J. E., Jr., et al. (SDSS Collaboration) 2000, AJ, 120, 1579 [Google Scholar]
 Zhang, Y., Anninos, P., & Norman, M. L. 1995, ApJ, 453, L57 [NASA ADS] [CrossRef] [Google Scholar]
Online material
Appendix A: Mock quasar spectra
We have produced mock spectra in order to tune the analysis procedure and to study statistical uncertainties and systematic effects in the measured correlation function.
In some galaxy clustering studies (e.g. Anderson et al. 2012) the covariance matrix of the measured correlation function is obtained from mock data sets. In this case, it is crucial to have very realistic mocks with the right statistics.
In order to do so, we would need to generate several realizations of hydrodynamical simulations, with a large enough box to cover the whole survey (several Gpc^{3}) and at the same time have a good enough resolution to resolve the Jeans mass of the gas (tenths of kpc). This type of simulations are not possible to generate with current technology, but luckily in this study the covariance matrix is obtained from the data itself, and the mock data sets are only used to test our analysis and to study possible systematic effects.
In the last few years there have been several methods proposed to generate simplified mock Lymanα surveys by combining Gaussian fields and nonlinear transformations of the field (Le Goff et al. 2011; Greig et al. 2011; FontRibera et al. 2012). In this study we used a set of mocks generated using the process described in FontRibera et al. (2012), the same method used in the first publication of the Lymanα correlation function from BOSS (Slosar et al. 2011).
The mock quasars were generated at the angular positions and redshifts of the BOSS quasars. The unabsorbed spectra (continua) of the quasars were generated using the Principal Component Analysis (PCA) eigenspectra of Suzuki et al. (2005), with amplitudes for each eigenspectrum randomly drawn from Gaussian distributions with sigma equal to the corresponding eigenvalues as published in Suzuki (2006) Table 1. A detailed description will be provided by Bailey et al. (in prep.), accompanying a public release of the mock catalogs.
We generated the field of transmitted flux fraction, F, that have a ΛCDM power spectrum with the fiducial parameters (A.1)where h = H_{0}/100 km s^{1}Mpc^{1}. These values produce a fiducial sound horizon of (A.2)Here, we use the parametrized fitting formula introduced by McDonald (2003) to fit the results of the power spectrum from several numerical simulations, (A.3)where μ_{k} = k_{∥}/k is the cosine of the angle between k and the line of sight, b_{δ} is the density bias parameter, β is the redshift distortion parameter, P_{L}(k) is the linear matter power spectrum, and D_{F}(k,μ_{k}) is a nonlinear term that approaches unity at small k. This form of P_{F} is the expected one at small k in linear theory, and provides a good fit to the 3D Lyα observations reported in Slosar et al. (2011). We do not generate a density and a velocity field, but directly create the Lyα forest absorption field instead, with the redshift distortions being directly introduced in the input power spectrum model of Eq. (A.3), with the parameter β that measures the strength of the redshift distortion.
To model the evolution of the forest with redshift, b_{δ} varies with redshift according to b_{δ} = 0.14 [(1 + z)/3.25] ^{1.9} (McDonald et al. 2006). The redshift distortion parameter is given a fixed value of β_{F} = 1.4. The nonlinear correction factor D(k,μ_{k}) is taken from McDonald (2003). The flux field was constructed by generating Gaussian random fields g with an appropriately chosen power spectrum (FontRibera et al. 2012) to which the lognormal transformation F = exp(−ae^{υg}) is applied (Coles & Jones 1991; Bi et al. 1992; Gnedin & Hui 1996). Here a and υ are free parameters chosen to reproduce the flux variance and mean transmitted flux fraction (McDonald et al. 2006).
DLA’s were added to the spectra according to the procedure described in FontRibera & MiraldaEscud (2012).
Finally, the spectra were modified to include the effects of the BOSS spectrograph point spread function (PSF), readout noise, photon noise, and flux miscalibration.
Fifteen independent realizations of the BOSS data were produced and analyzed with the same procedures as those for the real data.
Fig. A.1
Effect of the continuum estimation procedure on the correlation function found with the mock spectra. The black dots are the average of monopole and quadrupole obtained with the 15 sets using the exact continua. The blue (red) dots show those obtained with the continuum estimation of method 1 (method 2) as described in Sect. 3.1 (3.2). 
Fig. A.2
Comparison of the correlation function for the mock spectra and that for the data. The red dots show the mean of the 15 sets of mock spectra and the black dots show the data. 
We used the mock spectra to understand how our analysis procedure modifies the correlation function. Figure A.1 shows the average over 15 mocks of the reconstructed quadrupole and monopole using methods 1 and 2 (Sects. 3.1 and 3.2) and that reconstructed with the true continuum. The monopole and quadrupole for the two methods have a general shape that follows that found with the true continuum including the position of the BAO peak. However, both methods produce a monopole that becomes negative for 60 h^{1} Mpc < r < 100 h^{1} Mpc while the true monopole remains positive for all r < 130 h^{1} Mpc. As discussed in Sect. 3.2, this result is due to the continuum estimation of the two methods which introduced negative correlations. For both methods, however, the BAO peak remains visible with a deviation above the “broadband” correlation function that is hardly affected by the distortion.
Figure A.2 presents ξ_{0}(r) and ξ_{2}(r) found with the data, along with the mean of 15 mocks. The figure demonstrates that our mocks do not perfectly reproduce the data. In particular, for r < 80 h^{1} Mpc, the monopole is underestimated and the quadrupole overestimated. Since we use only peak positions to extract cosmological constraints, we only use the mocks qualitatively to search for possible systematic problems in extracting the peak position.
Appendix B: Results for a fiducial BOSS Lyα forest sample
The spectra analyzed here are all available through SDSS DR9 (Ahn et al. 2012), and the DR9 quasar catalog is described
by Pâris et al. (2012). A Lyα forest analysis requires many detailed choices about data selection and continuum determination. To aid community analyses and comparison of results from different groups, Lee et al. (2012b) has presented a fiducial BOSS Lyα forest sample that uses constrained PCA continuum determination (Lee et al. 2012a) and reasonable choices of masks for DLAs, BALs, and data reduction artifacts. Both the data selection and the continuum determination differ from those used here. Figure B.1 compares the Method 2 correlation function from this paper’s analysis to that obtained by applying the Method 2 weights and correlation measurement code directly to the continuumnormalized spectra of the fiducial Lee et al. sample. The good agreement in this figure, together with the good agreement between our Method 1 and Method 2 results, demonstrates the robustness of the BAO measurement, and the more general correlation function measurement, to the BOSS DR9 data.
Fig. B.1
Comparison of the monopole and quadrupole correlation functions for the sample used here (black dots) and for the sample and continua of Lee et al. (2012b) (red dots). 
All Tables
All Figures
Fig. 1
HammerAitoff projection in equatorial coordinates of the BOSS DR9 footprint. The observations cover ~3000 deg^{2}. 

In the text 
Fig. 2
Top: weighted distribution of absorber redshifts used in the calculation of the correlation function in the distance range 80 h^{1} Mpc < r < 120 h^{1} Mpc. Bottom: distribution of signaltonoise ratio for analysis pixels (triplets of pipeline pixels) averaged over the forest region. 

In the text 
Fig. 3
An example of a BOSS quasar spectrum of redshift 3.239. The red and blue lines cover the forest region used here, 104.5 < λ_{rf} < 118.0. This region is sandwiched between the quasar’s Lyβ and Lyα emission lines respectively at 435 and 515 nm. The blue line is an estimate of the continuum (unabsorbed flux) by method 2 and the red line is the estimate of the product of the continuum and the mean absorption by method 1. 

In the text 
Fig. 4
Mean of δ(λ) as a function of observed wavelength (method 1). Systematic offsets from zero are seen at the 2% level. The calcium lines (393.4, 396.8 nm) is present. The features around the hydrogen lines Hγ, δ and ϵ (434.1, 410.2, 397.0 nm) are artifacts from the use of Fstars for the photocalibration of the spectrometer. 

In the text 
Fig. 5
Mean transmitted flux fraction as a function of redshift obtained from the continuum fits with method 2. Data are shown in black, mock000 in red and the input mean transmitted flux fraction in blue. 

In the text 
Fig. 6
Top panel: inverse total variance in bins of redshift as a function of the pipeline inverse variance. Bottom panel: parameters of the fit: the LSS contribution σ_{LSS} (left) and the pipeline correction factor η (right) as a function of redshift. The lines show fits to the data as explained in the text. 

In the text 
Fig. 7
ξ(r,μ) vs. μ^{2} for the bin centered on r = 34 h^{1} Mpc. The red dots are the mean of the 15 mocks and the blue dots are the data. 

In the text 
Fig. 8
ξ(r,μ) averaged over 0.1 < μ < 0.5, 0.5 < μ < 0.8 and 0.8 < μ < 1. The curves give fits (Sect. 5) to the data imposing concordance ΛCDM cosmology. The BAO peak is most clearly present in the data for μ > 0.8. 

In the text 
Fig. 9
Monopole (upper panel) and quadrupole (middle panel) correlation functions found by method 1 (red) and method 2 (black). The bottom panel shows the combination ξ_{0} + 0.1ξ_{2} found by method 1 (red) and method 2 (black). 

In the text 
Fig. 10
rdependence of the product of the monopole (top) and quadrupole (middle) variances and the number of pairs in the rbin. The bottom panel shows the product for the monopolequadrupole covariance (r = r′ elements). The dotted lines show the means for r > 20 h^{1} Mpc. 

In the text 
Fig. 11
Monopole and quadrupole covariance matrix. The monopolemonopole and quadrupolequadrupole elements are normalized to the variance: . The monopolequadrupole elements are normalized to the mean of the quadrupole and monopole variances. The first offdiagonal elements of the monopolemonopole and quadrupolequadrupole elements are ~20% of the diagonal elements. The diagonal elements of the quadrupolemonopole covariance are ~− 0.2 times the geometric mean of the monopole and quadrupole variances. 

In the text 
Fig. 12
Verification of the offdiagonal elements of the covariance matrix with the 15 sets of mock spectra. The black lines show correlations derived from the dispersion of the 15 measurements and the red lines show the correlations expected from the covariance matrix calculated by subsampling. The top and middle panels show the correlation between neighboring bins for monopole and quadrupole respectively. The bottom panel, the correlation between monopole and quadrupole measured at the same distance bin. 

In the text 
Fig. 13
Fitting functions used for the determination of the peak detection significance: r^{2}ξ^{bb}(r) (blue), r^{2} [ξ^{peak}(r) + ξ^{bb}(r)] (black) and r^{2}ξ^{dist}(r) (red) for the monopole (solid lines) and quadrupole (dashed lines). 

In the text 
Fig. 14
Monopole and quadrupole fits with a BAO peak (red line) and without a BAO peak (blue line, Method 2). The fitting range is r_{min} < r < 200 h^{1} Mpc with r_{min} = 20 h^{1} Mpc. 

In the text 
Fig. 15
Same as Fig. 14 except for different fitting ranges: blue, green and red curves are for fits with r_{min} = 20 h^{1} Mpc, 40 h^{1} Mpc and 60 h^{1} Mpc respectively (Method 2). The solid lines for fits without a BAO peak and the dashed lines with a peak. 

In the text 
Fig. 16
Contours for (D_{A}/r_{s},r_{s}H) obtained by fitting the monopole and quadrupole to (22). The broadband distortions are Eq. (23) – dashed lines or Eq. (24) – solid lines. The blue lines are for method 1 and the red lines for method 2. All contours are for Δχ^{2} = 4 except for the interior solid red contour which is for Δχ^{2} = 1. 

In the text 
Fig. 17
Measurements of α_{iso} (= α_{t} = α_{r}) for the 15 sets of mock spectra and for the data (realization = −1). The large errors for realization 5 and 8 are due to the very low significance of the BAO peak found on these two sets. 

In the text 
Fig. 18
Constraints on the matter and darkenergy density parameters (Ω_{M},Ω_{Λ}) assuming a darkenergy pressuredensity ratio w = −1. The blue regions are the one and two standard deviation constraints derived from our contours in Fig. 16 (method 2, broadband 24) combined with a measurement of H_{0} (Riess et al. 2011). Also shown are one and two standard deviation contours from lower redshift measurements of D_{V}/r_{s} (also combined with H_{0}) at z = 0.11 (6dF: Beutler et al. 2011), z = 0.35 (LRG: Percival et al. 2010) and z = 0.57 (CMASS: Anderson et al. 2012). All constraints use a WMAP7 (Komatsu et al. 2011) prior on the baryontophoton ratio η but do not otherwise incorporate CMB results. 

In the text 
Fig. 19
As in Fig. 18 with constraints on the matter density parameter, Ω_{M}, and darkenergy pressuredensity ratio w assuming Ω_{M} + Ω_{Λ} = 1. 

In the text 
Fig. 20
Constraints on (D_{A}/r_{s},r_{s}H)_{z = 2.3} within the framework of OwOwaCDM models. The green contours are our 1σ and 2σ constraints using method 2 and broadband (24). The gray contours are the 1σ and 2σ constraints from WMAP7 (Komatsu et al. 2011). The red contours show the combined constraints. 

In the text 
Fig. 21
Measurements of H(z)/(1 + z) vs. z demonstrating the acceleration of the expansion for z < 0.8 and deceleration for z > 0.8. The BAObased measurements are the filled circles: this work: (red), Xu et al. (2012, black)Chuang & Wang (2012, blue), Reid et al. (2012, cyan), and Blake et al. (2012, green). The open green circles are from WiggleZ (Blake et al. 2011b) AlcockPaczynski data combined with supernova data yielding H(z)/H_{0} (without the flatness assumption) plotted here assuming H_{0} = 70 km s^{1} Mpc^{1}. The open blue circle is the H_{0} measurement of Riess et al. (2011). The open black squares with dashed error bars show the results of Riess et al. (2007) which were derived by differentiating the SNIa Hubble diagram and assuming spatial flatness. (For visual clarity, the Riess et al. 2007 point at z = 0.43 has been shifted to z = 0.48.) The line is the ΛCDM prediction for (h,Ω_{M},Ω_{Λ}) = (0.7,0.27,0,73). 

In the text 
Fig. A.1
Effect of the continuum estimation procedure on the correlation function found with the mock spectra. The black dots are the average of monopole and quadrupole obtained with the 15 sets using the exact continua. The blue (red) dots show those obtained with the continuum estimation of method 1 (method 2) as described in Sect. 3.1 (3.2). 

In the text 
Fig. A.2
Comparison of the correlation function for the mock spectra and that for the data. The red dots show the mean of the 15 sets of mock spectra and the black dots show the data. 

In the text 
Fig. B.1
Comparison of the monopole and quadrupole correlation functions for the sample used here (black dots) and for the sample and continua of Lee et al. (2012b) (red dots). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.