Issue 
A&A
Volume 542, June 2012



Article Number  A34  
Number of page(s)  11  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201118017  
Published online  01 June 2012 
Wavelet analysis of baryon acoustic structures in the galaxy distribution^{⋆}
^{1} Observatori Astronòmic, Universitat de València, Apartat de Correus 22085, 46071 València, Spain
^{2} Institut de Física Corpuscular (CSICUVEG), València, Spain
^{3} Institute for Computational Cosmology, Department of Physics, Durham University, South Road, Durham DH1 3LE, UK
email: pablo.arnaltemur@durham.ac.uk
^{4} Laboratoire AIM (UMR 7158), CEA/DSMCNRSUniversité Paris Diderot, IRFU, SEDISAP, Service d’Astrophysique, Centre de Saclay, 91191 GifSurYvette Cedex, France
^{5} Departament d’Astronomia i Astrofísica, Universitat de València, 46100Burjassot, València, Spain
^{6} Astroparticule et Cosmologie (APC), CNRSUMR 7164, Université Paris 7, Denis Diderot 10, rue Alice Domon et Léonie Duquet, 75205 Paris Cedex 13, France
^{7} Tartu Observatoorium, 61602 Tõravere, Estonia
^{8} Departamento de Matemática Aplicada y Estadística, Universidad Politécnica de Cartagena, C/Dr. Fleming s/n, 30203 Cartagena, Spain
Received: 5 September 2011
Accepted: 7 March 2012
Context. Baryon acoustic oscillations (BAO) are imprinted in the density field by acoustic waves travelling in the plasma of the early universe. Their fixed scale can be used as a standard ruler to study the geometry of the universe.
Aims. The BAO have been previously detected using correlation functions and power spectra of the galaxy distribution. We present a new method to detect the realspace structures associated with BAO. These baryon acoustic structures are spherical shells of relatively small density contrast, surrounding high density central regions.
Methods. We design a specific wavelet adapted to search for shells, and exploit the physics of the process by making use of two different mass tracers, introducing a specific statistic to detect the BAO features. We show the effect of the BAO signal in this new statistic when applied to the Λ – cold dark matter (ΛCDM) model, using an analytical approximation to the transfer function. We confirm the reliability and stability of our method by using cosmological Nbody simulations from the MareNostrum Institut de Ciències de l’Espai (MICE).
Results. We apply our method to the detection of BAO in a galaxy sample drawn from the Sloan Digital Sky Survey (SDSS). We use the “main” catalogue to trace the shells, and the luminous red galaxies (LRG) as tracers of the high density central regions. Using this new method, we detect, with a high significance, that the LRG in our sample are preferentially located close to the centres of shelllike structures in the density field, with characteristics similar to those expected from BAO. We show that stacking selected shells, we can find their characteristic density profile.
Conclusions. We delineate a new feature of the cosmic web, the BAO shells. As these are real spatial structures, the BAO phenomenon can be studied in detail by examining those shells.
Key words: largescale structure of Universe / distance scale / galaxies: cluster: general / methods: data analysis / methods: statistical
Full Table 1 is only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/542/A34
© ESO, 2012
1. Introduction
Before recombination, the energy of photons is high enough to avoid the formation of neutral hydrogen atoms. This means that baryons and photons are coupled through Compton scattering and electromagnetic interaction between protons and electrons, forming a plasma. In this fluid, two phenomena act in opposite directions: gravitational forces tend to compress the plasma around high density regions, while radiation pressure tends to dilute any such overdensity. The combination of both in the presence of any initial inhomogeneity gives rise to acoustic waves propagating in the baryonphoton plasma. This phenomenon ends abruptly at the epoch of recombination, when the temperature drops sufficiently to allow hydrogen atoms to form, and therefore radiation decouples from the baryons.
Baryon acoustic oscillations (BAO) are therefore due to the propagation of these sound waves in the baryonphoton plasma in the early universe (Peebles & Yu 1970; Hu et al. 1997; Eisenstein & Hu 1998; Bassett & Hlozek 2010). Any primordial overdensity in the early universe produces a spherical acoustic wave in the baryonphoton plasma, travelling outwards: the radiation pressure drags the baryons that are coupled to the photons, and compensates the gravity force that pulls all matter towards the centre. Dark matter, however, is totally decoupled from the photons, and therefore its density at the centre continues growing. About 380 000 years after the Big Bang, the temperature drops so that photons and baryons decouple, and the scale of the baryon shells freezes. After this time, both the central overdensity and the shell grow gravitationally, accreting both dark matter and baryons. The result at late times is a large overdensity at the position of the original perturbation, surrounded by a faint spherical shell on a fixed comoving scale (Eisenstein et al. 2007).
The BAO scale is fixed by the sound horizon at decoupling because it is the distance that the expanding acoustic shells can travel before decoupling. It has been accurately measured by the study of the anisotropies in the cosmic microwave background (CMB) to be (Komatsu et al. 2009) r_{s} = 153.3 ± 2.0 Mpc = 110.4 ± 1.4 h^{1} Mpc (where we take h = 0.72, Freedman et al. 2001)^{1}. Therefore, this scale, once measured, could be used as a standard ruler to measure the Hubble expansion rate H(z) and the angular diameter distance D_{A}(z) as function of redshift (Cooray et al. 2001; Blake & Glazebrook 2003; Seo & Eisenstein 2005).
The BAO should appear as a series of damping wiggles in the matter power spectrum, with the locations of the peaks and troughs in kspace being a function of r_{s} and other cosmological parameters (Eisenstein & Hu 1998). All the harmonics sum up to the same peak in the galaxy correlation function ξ(r) on the scale r_{s}, and therefore it seems more appropriate to use this statistic to detect the BAO feature in the available galaxy redshift surveys encompassing large volumes of the universe (Sánchez et al. 2008).
The first detection of BAO (at about the 3σ level) was reported in an analysis of the correlation function (Eisenstein et al. 2005) of the Sloan Digital Sky Survey (SDSS, York et al. 2000) luminous red galaxies (LRG) sample (Eisenstein et al. 2001), and later in the power spectrum (Cole et al. 2005) of the Twodegree Field Galaxy Redshift Survey (2dFGRS) (Colless et al. 2001). However, this is certainly a controversial topic. Cabré & Gaztañaga (2011) are not finding such level of detection using a data set twice as large in both volume and number of galaxies. They do not claim that this result contradicts the standard Λ – cold dark matter (ΛCDM) model, but instead that it is a consequence of insufficient data. One of the arguments in Cabré & Gaztañaga (2011) is that mixing model selection with parameter determination can lead to some confusion in interpreting the results and their significance. Different authors use different criteria to assess the significance of their BAO detection. For example, when Eisenstein et al. (2005) affirm that the baryon signature was detected at 3.4σ (or at 3.0σ when including only data points between 60 and 180 h^{1} Mpc), they compare their results of the SDSSLRG correlation function with those expected for the bestfit pure CDM model and different BAO models. The most robust BAO detection up to now (Percival et al. 2010) was obtained by studying the combined power spectrum of the LRG and the “main” (Strauss et al. 2002) samples of SDSS, together with the 2dFGRS sample, and is at the ~ 3.6σ level. The authors explicitly state that since this number is obtained by comparing to an arbitrary smooth model, the significance cannot be directly compared with the one reported in Eisenstein et al. (2005). This is a clear example of different authors using different ways to assess the significance of their results that in practice are not comparable. Hütsi (2006) calculated the redshift space power spectrum of the SDSSLRG sample drawn from the data release 4. He concludes that BAO models are more likely by 3.3σ over the corresponding models without any oscillatory behaviour in the power spectrum.
Percival et al. (2007) detected BAO in the clustering of the combined 2dFGRS and SDSS galaxy samples, and used their measurements to constrain cosmological models, in particular a given combination of the angular diameter distance D_{A}(z) and the Hubble parameter H(z). Cabré & Gaztañaga (2009a,b) studied the LRG anisotropic redshiftspace correlation function ξ(σ,π), where π is the lineofsight or radial separation and σ is the transverse separation. Moreover, Gaztañaga et al. (2009) illustrated how H(z) could be constrained using the correlations in the radial direction. Kazin et al. (2010a) found similar results for the correlation measurements and uncertainties, but manifest disagreement in the interpretation of the results regarding the detection of a lineofsight baryon acoustic feature.
More recent studies (Martínez et al. 2009; Cabré & Gaztañaga 2009a; Sánchez et al. 2009; Kazin et al. 2010b) have confirmed this detection in the last data release (DR7, Abazajian et al. 2009) of the SDSSLRG, containing twice as many galaxies as the original sample, although the observed peak is in these cases wider than that observed in the original detection – an issue that needs further explanation. These measurements of the BAO scale at low redshift, combined with other cosmological probes, have been used to place stringent constraints on the values of cosmological parameters (Tegmark et al. 2006; Percival et al. 2007; Sánchez et al. 2009; Percival et al. 2010; Reid et al. 2010; Kazin et al. 2010a).
While Bassett & Afshordi (2010) argue that lowlevel detections may be insufficient to robustly estimate the cosmological parameters, Cabré & Gaztañaga (2011) show instead that it is still possible – assuming a model – to locate the BAO position with data providing very low significant BAO detection.
It is important, therefore, to find evidence of BAO in the galaxy distribution based on complementary methods. A step further is to search for real structures in the galaxy distribution that are responsible for the BAO feature in these secondorder statistics. The detection of these structures would confirm the existence of the baryon acoustic phenomenon. Moreover, if we were able to localise these structures in configuration space, this would allow us to study in more detail the properties of the BAO.
In this paper, we introduce a new method for the detection of BAO that is closely tied to the underlying physics of the process, and apply it to a sample drawn from the SDSS catalogue. This method (described in Sect. 2) is based on analysing directly the threedimensional (3D) galaxy distribution using a very specific wavelet function (which we called “BAOlet”), which is especially wellsuited to searching for BAO features. The method makes use of two different tracers, one to map the overall density field (including the BAO shells), and the other to locate the position of the largest overdensities, which should correspond to the centres of the shells. As we directly study the galaxy distribution in configuration space, this method also allows us to identify regions of space where the BAO signal is either stronger or fainter. We describe the expected signal in the ΛCDM model in Sect. 3, using both an analytical prediction and a Nbody simulation catalogue. We describe the samples used in the case of SDSS in Sect. 4. In Sect. 5, we show the results obtained in this case. We also perform a test to assess the significance of these results, and explore the implications of this analysis regarding the localisation of BAO structures. Finally, we summarise our conclusions and discuss possibilities for future work in Sect. 6.
2. The wavelet detection method
The basis of the new BAO detection method is to focus on the positions of massive dark matter haloes, which correspond to the location of large initial perturbations, and to search for structures resembling the acoustic shells around these. Once we locate the positions of the large overdensities, we need to study the density field to identify the structures corresponding to the acoustic shells around these centres. An appropriate method for identifying structures in continuous fields is a wavelet analysis (Martínez et al. 1993; Starck & Murtagh 2006; Jones 2009). Wavelet transforms are widely used in many areas, especially in image analysis (Mallat 2008; Starck et al. 2010). They are especially suited to the analysis of data on different scales, and the identification of characteristic patterns or structures. Wavelets have been used in cosmology to analyse largescale structure and CMB anisotropies (Martínez et al. 1993; Rauzy et al. 1993; Vielva et al. 2004; Starck et al. 2006; Saar 2009).
Standard wavelet functions such as the Mexican hat are, however, unsuitable for the detection of shells. We need instead a family of wavelets whose shape matches the type of structures we wish to find in our data. We use therefore a especially designed wavelet (the “BAOlet”), which is well adapted to the search of BAO features – shelllike structures around our selected centres. We design this new family of wavelet functions as a transformation of the wideused Bspline wavelets (Saar 2009). These ψ_{R,s}(x) functions are spherically symmetric, and their radial profiles are defined as (1)where R and s are the two parameters that define the scale and width of the BAOlet function, α_{R,s} is the normalisation constant defined so that (2)and B_{3}(x) is the box spline of the third degree, defined by The BAOlet function is shown in Fig. 1. It can be thought of as a spherical shell of radius R and width s, with zero amplitude at its centre and therefore adapted to the detection of spherical shells of a given radius. This specific choice is motivated by the integrated profile being the widely used onedimensional (1D) “Bspline” wavelet function that has a null mean and compact support [ − 2,2] . These properties directly translate into the BAOlet, which also has a null mean – a requirement for any wavelet function – if R > 2s, and takes nonzero values only for R − 2s ≤ x ≤ R + 2s.
Fig. 1 The BAOlet function. Here we show a twodimensional (2D) plot (bottom) of the wavelet ψ_{R,s}(x) used in the analysis, as defined by Eq. (1). The top panel shows a 1D slice along the dasheddotted axis. The wavelet is plotted here for R = 105 h^{1} Mpc, s = 30 h^{1} Mpc. The red dot marks the centre of the wavelet. This function has a null mean (provided that R > 2s), and compact support. It takes nonzero values only for R − 2s ≤ x ≤ R + 2s. 

Open with DEXTER 
We describe the density field using the density contrast δ(x), defined as where ρ(x) is the density field, and ρ_{0} is its mean. Then, given a density contrast map δ(x), properly normalised as in Eq. (2), we can construct, for each point in the parameter space (R,s), a BAOlet coefficient map as the convolution of our density field with the corresponding wavelet (3)The BAOlet acts as a matched filter, which is sensitive to data containing shells of different radii and different widths. Its property of zero mean is also of high importance since it makes the statistics derived from the BAOlet coefficients independent of the background level. It is indeed obvious that any constant added to the input data would not change the BAOlet coefficients. In comparison, the estimation of this baseline level is a very delicate aspect of the BAO detection in the twopoint correlation function.
Owing to the properties of the wavelet, the coefficient maps W_{R,s}(x) should have a null mean when averaged over all points in the volume considered. Equivalently, if we sampled these maps at N random points uniformly distributed in the volume (), the expected value of the average of the coefficients would be zero (4)This condition holds even in the presence of shelllike structures in the density field. Of course, for such structures the value of W_{R,s}(x_{c}) (x_{c} is the centre of the shell) is positive, and remains positive in nearby points. For an ideal δ(r − R) density shell, the radius of the region around the centre where the wavelet amplitude is positive, is s; the positive signal in this region is compensated by negative amplitudes around x = R. However, if we are able to identify the positions of N massive haloes in the same volume (), we can define a new statistic B(R,s) as the mean value of the coefficients W_{R,s}(x) at these positions (5)If there are indeed shelllike structures around the selected density maxima , as expected for baryon acoustic structures, we should find positive values of B(R,s) with the maximum of B at the (R,s) values characterising these shells.
We can obtain further information from the wavelet coefficients W_{R,s}(x), as we have information on the actual dependence of the signal picked up by the BAOlet function on the position. In particular, fixing a set of parameters of interest (R_{i},s_{i}), we could use the coefficients W_{Ri,si}(x_{c}) to identify which of the selected massive haloes produce the strongest signal for these shell characteristics. In the context of BAO, the parameters R_{i}, s_{i} can be chosen a priori using a theoretical model, or a posteriori using the parameters for which the function B(R,s) attains its maximum. In this way, we can localise in configuration space the structures responsible for the largest BAO signal in a given sample.
For our calculation of B(R,s), we sample the (R,s) parameter space on a grid. For each point (R,s), we calculate the coefficient map W_{R,s}(x) as the convolution of the BAOlet with the density field (Eq. (3)). We perform the convolution in Fourier space using a fast Fourier transform (FFT) technique. To avoid problems with the FFT, we zeropad a large region around our density cube. To obtain B(R,s), we sample W_{R,s}(x) at the position of the N selected centres, and calculate the average value (Eq. (5)).
Therefore, to apply this method, we need a way to map the overall density field δ(x), but also to locate the positions of the massive matter haloes . We have to use two different populations of mass tracers, so that they play the appropriate role in the detection algorithm. The idea of using two different tracer sets, one for the small perturbations and another for the high peaks, in a crosscorrelation analysis was anticipated by Eisenstein et al. (2007). We implement here a similar idea, but applying a wavelet tool directly to the density field. As detailed below, we use galaxies from the “main” and LRG samples of SDSS in this case. However, this choice depends on the kind of data available in each case.
3. Prediction from ΛCDM
To illustrate our method more clearly, we show here the results we would expect according to the ΛCDM model, and the effect of BAO in our new statistic B(R,s). For this aim, we use both the analytical approximation to the transfer function of Eisenstein & Hu (1998), and the results from the MareNostrum Institut de Ciències de l’Espai (MICE) simulation (Fosalba et al. 2008).
In the first place, we use the ΛCDM transfer function, which allow us to study directly the effect of the BAO. However, in this case, we must make a series of approximations in order to derive a prediction for B(R,s). We wish to predict the typical result for the wavelet coefficient W_{R,s} at the position of massive matter haloes x_{c}, as a function of R,s. From Eq. (3), we see that this is equivalent to studying the typical density profile around these haloes, δ(y − x_{c}). The ΛCDM transfer function allows us to calculate this profile, provided that we know the initial perturbation corresponding to the selected haloes. We make here the simple approximation of considering that these initial perturbations are pointlike and spherically symmetric, and can thus be simply described by a Dirac delta function in configuration space. This corresponds to a constant value in Fourier space. As the transfer function T(k) describes the relative evolution of the different Fourier modes, the presentday radial density profile corresponding to this initial perturbation is given simply by (Eisenstein et al. 2007) (6)where is the Fourier transform of the transfer function T(k), and C is a normalisation constant that depends on both the details of the initial perturbation and the cosmic growth function D_{1}(z). From Eqs. (3) and (5), we see that the effect of C is just to change the normalisation of our statistic B(R,s).
We use the fitting formulae to the transfer function T(k) of Eisenstein & Hu (1998), and obtain the expected W_{R,s} at a large overdensity using Eqs. (6) and (3). To highlight the particular signature of BAO, we also calculate W_{R,s} using the “no wiggle” transfer function formula, in which the BAO have been removed. We use here the values Ω_{M} = 0.25, Ω_{Λ} = 0.75, Ω_{b} = 0.044, and h = 0.7 for the cosmological parameters, to allow a direct comparison with the MICE simulation. Following Eisenstein & Hu (1998), the sound horizon scale in this case is r_{s} = 109.3 h^{1} Mpc. The results in both cases are shown in Fig. 2. In the plot, we mask the region R < 2s, as for these values of the parameters our BAOlet is not compensated (its mean is different from 0). Comparing both panels of the figure, we see clearly the effect of the presence of BAO in our statistic B(R,s). In the case without BAO, W_{R,s} is always negative, and has a smooth gradient across the (R,s) plane. This gradient is due to the overall shape of the radial profile (Eq. (6)). However, in the presence of BAO, W_{R,s} has a prominent peak with positive values. This clearly illustrates the concept behind the B(R,s) statistic. The BAOlet ψ_{R,s} acts as a matched filter with a shape adapted to detecting BAO shells. Therefore, the positive values in the coefficients W_{R,s} correspond to the cases in which the radial profile is matched by the BAOlet shape. The values at which W_{R,s} attains its absolute maximum, R_{max} = 110 h^{1} Mpc and s_{max} = 22 h^{1} Mpc, thus correspond to the characteristics of the shell that most closely match the observed profile around the selected centres.
Fig. 2 Values of the BAOlet coefficients W_{R,s} expected at the positions of large initial pointlike perturbations, as a function of the BAOlet parameters (R,s). The bottom panel shows the result using a standard ΛCDM transfer function, while the top panel shows the result using a transfer function with the BAO wiggles edited out (Eisenstein & Hu 1998). The normalisation is arbitrary. The contours are drawn in steps of 1000 for W_{R,s} < 0 (dotted), W_{R,s} = 0 (solid), and W_{R,s} > 0 (dashed). The map attains a maximum at R = 110 h^{1} Mpc, s = 22 h^{1} Mpc. 

Open with DEXTER 
To test the reliability of the method, and this ΛCDM prediction, we calculate the B(R,s) for a halo catalogue drawn from the MICE simulations. We use the publicly available halo catalogue from the “MICE3072” run (Crocce et al. 2010). This particular run contains 2048^{3} particles in a box of side 3072 h^{1} Mpc, therefore covers a volume of 29 h^{3} Gpc^{3}. Crocce et al. performed the simulation with the GADGET2 code (Springel 2005), assuming a ΛCDM model with the parameters mentioned above. The haloes in the simulation were selected using a friendsoffriends (FoF) algorithm.
We use the resulting halo catalogue at z = 0, which contains a total of 2 819 031 haloes containing 143 or more particles. This corresponds to haloes with masses ≥ 3.35 × 10^{13} h^{1} M_{⊙}. The halo number density is thus 9.72 × 10^{5} h^{3} Mpc^{3}. We use the full halo catalogue as a tracer of the overall density field. We then select as centres for the calculation of B(R,s) in Eq. (5) only the haloes with a mass ≥ 1.76 × 10^{14} h^{1} M_{⊙}. We choose this mass threshold in order to select approximately the 10% most massive haloes in the simulation box. This choice is somewhat arbitrary, but serves for the purpose of testing the BAOlet method and illustrating the expected result.
Figure 3 shows the BAOlet result B(R,s) for these MICE samples, compared to the theoretical results obtained above from the Eisenstein & Hu (1998) transfer functions. We obtain a result very similar to that of Fig. 2, as B(R,s) has a clear peak, and attains its absolute maximum for R_{max} = 108 h^{1} Mpc and s_{max} = 28 h^{1} Mpc. This indicates that our BAOlet method can be applied to two sets of mass tracers, although the details of the tracers used here are very different from the ones we use later on with the SDSS samples. This also confirms the expected effect of the presence of BAO in the B(R,s) function, namely a large peak with positive values of B, located approximately at the values of R and s corresponding to the radius and width of the acoustic shells. That we obtain here slightly different values for R_{max} and s_{max} than those predicted above may be due to nonlinear evolution effects, which slightly reduce the radius and increase the width of the shells. A similar effect is present in the correlation function (see e.g. Crocce & Scoccimarro 2008).
Fig. 3 The BAOlet statistic B calculated for the MICE simulation sample described in the text as a function of the parameters (R,s) (bottom panel). The contours are drawn in steps of 5 for B < 0 (dotted), B = 0 (solid), and B > 0 (dashed). This function attains its maximum for R = 108 h^{1} Mpc, s = 28 h^{1} Mpc. The top two panels show cuts at the values s = 28 h^{1} Mpc (top) and s = 22 h^{1} Mpc (middle), marked with grey horizontal lines in the 2D panel. In each case, the solid blue line corresponds to the value obtained from MICE, the dashed red line corresponds to the theoretical expectation from the Eisenstein & Hu (1998) transfer function (bottom panel of Fig. 2), and the dotted green line to the theoretical expectation using the “no wiggle” transfer function (top panel of Fig. 2). These theoretical predictions have been renormalised to get the same value at the maximum in B(R,s). 

Open with DEXTER 
We also use this halo catalogue from MICE to make a qualitative estimation of how different observational effects would affect the BAOlet result. In the first place, we study the effect of redshiftspace distortions. To this end, we calculate the redshiftspace positions of all haloes taking into account their peculiar velocities, as given by the simulation, and considering an observer located in one of the vertices of the simulation cube. The result for B(R,s) in this case is shown in the top panel of Fig. 4, where it is compared to the realspace result discussed above. As can be seen from the figure, although small differences appear between the real and redshiftspace results, the main features of the B(R,s) prediction remain the same, with the position of the maximum changing by only ~1 h^{1} Mpc.
In the second case, we add the effect of a decreasing radial selection function across the sample. We model this selection as an exponential decay function, such that the final number of haloes used to trace the overall density field is ~20% of the total. In our calculations, we then weight each halo by the inverse of the mean density at its redshift, as we do later for the SDSS data. We do not apply any selection function to the centres. The results for B(R,s) obtained in this case (including also the redshiftspace effects) are shown in the bottom panel of Fig. 4. As above, these observational effects do not significantly change the overall behaviour of B(R,s), or the location of the maximum of the peak. Overall, although the MICE catalogue used does not mimic the characteristics of our SDSS samples, we can be confident that neither redshiftspace distortions nor a radial selection function (when it is taken into account in the calculation) should significantly bias our results.
Fig. 4 The BAOlet statistic B for the MICE simulation when some observational effects are taken into account. In the top panel, we show the B(R,s) obtained when redshiftspace distortions are included in the simulation. In the bottom panel, we show the result when a radial selection function is applied to the halo catalogue. In both cases, the contours are drawn in steps of 5 in B. Solid contours correspond to the results with the observational effects included. The dashed contours correspond to the original realspace result without selection, i.e., they are identical to those in the bottom panel of Fig. 3. 

Open with DEXTER 
4. SDSS samples used
We use data from two different samples of the DR7 of the spectroscopic SDSS. On the one hand, we use the “main” galaxy sample (Strauss et al. 2002) as mass tracers to reconstruct the overall density field δ(x). On the other hand, we use the LRG as tracers of the central overdensities, and therefore use them as the selected centres to compute B(R,s).
LRG were selected by the SDSS team using several colour and magnitude cuts to obtain a highly biased sample extending to high redshift (Eisenstein et al. 2001). The galaxies selected in this way are known to reside near the centres of massive dark matter haloes (Zheng et al. 2009) and are thus adequate tracers for the centres of baryon acoustic structures. We apply an extra cut in the Kcorrected, evolved, gband absolute magnitude of −23.2 < M_{g} < −21.2, as in the previous BAO analysis by Eisenstein et al. (2005). This results in an approximately volumelimited sample in the redshift range 0.15 < z < 0.30.
“Main” galaxies in the SDSS constitute a much denser sample, and are therefore more suitable to map small density changes such as BAO shells. We use the “main” sample from the ValueAdded Galaxy Catalogue (Blanton et al. 2005), which constitutes a magnitudelimited sample in the r band, with r < 17.6. We apply an extra simple cut, M_{r} < −20.
For the conversion of angles and redshifts into comoving distances, we use a fiducial cosmology with the parameters Ω_{M} = 0.25, Ω_{Λ} = 0.75. In our analysis, we use distances in units of h^{1} Mpc, so that they do not depend on the specific value of h. We convert the distribution of the “main” galaxies into a density field δ(x) binning it into a grid with cubic pixels of 3 h^{1} Mpc side. We correct for the selection effects by weighting each galaxy by the inverse of the average density at its redshift. As explained below, we performed some tests by slightly changing this weighting scheme. Although this weighting may not be optimal, it should not significantly affect our results, given that the wavelet method does not depend on the local background level. We used the density field constructed in this way for the calculation of the BAOlet coefficients following Eq. (3).
In our calculations, we could only use the region in which these two samples overlap, which corresponds to the redshift limits 0.15 < z < 0.26. To minimise border effects in the B(R,s) calculation, we defined a buffer region of r_{buff} = 175 h^{1} Mpc from any of the borders of the “main” sample volume. As centres, we use only the LRG in the inner volume. This allows us to use the density field, as traced by the “main” sample galaxies, from z > 0.09. To minimise angular selection effects and border effects, we use a compact area of the sky where the angular completeness is nearly uniform. This area covers 5511deg^{2} and is defined, in the SDSS survey coordinates (Stoughton et al. 2002), by the limits −31.25° < η < 28.75°, −54.8° < λ < 51.8°. This results in finally using the density field in a volume of 2.2 × 10^{8} h^{3} Mpc^{3}, as traced by N_{Main} = 198 342 galaxies. The number of LRG used as centres (avoiding the buffer region) is N_{LRG} = 1599.
In Fig. 5, we show a slice of this survey showing both the “main” galaxies and the LRG. We show how, given the buffer used, the LRG used as centres are located in only an inner volume of the larger “main” sample. To illustrate the idea of the method, we show a zoom around a given LRG galaxy. Even for this single centre, a slight overdensity of “main” galaxies is seen at the radii of 105 − 110 h^{1} Mpc.
Fig. 5 The SDSS catalogues used. We show a sixdegree thick slice drawn from the catalogues used in our analysis. The red points are the LRG galaxies, which form a nearly volumelimited sample at z > 0.15. The “main” galaxies of the SDSS are depicted in blue. As shown, we use the “main” galaxies from z > 0.09, and we use only the LRG in an inner volume, allowing for a buffer region of r_{buff} = 175 h^{1} Mpc from any of the borders of the “main” sample volume. For our analysis, we use the samples covering a total area of 5511deg^{2}. The radius of a typical BAO shell is shown as a segment. At the top insets, we show two orthogonal slices of width 20 h^{1} Mpc, centred on a particular LRG. The BAOlet coefficient map W_{R,s}(x) has a large value at the position of this centre (for R_{max},s_{max}), thus we expect to find a strong BAO signal. The two circles have the radii of 100 and 110 h^{1} Mpc, respectively. This is a single BAO shell where the overdensity can be appreciated by eye on the right scale. 

Open with DEXTER 
As the structures we search for are huge, with radii of about 100 h^{1} Mpc, we have to consider the effect of the assumed cosmology (different comoving distances) on our result. To estimate the distance differences, we compare the distances in our adopted MICE cosmology (Ω_{M} = 0.25, Ω_{Λ} = 0.75) with those in the WMAP sevenyear cosmological model (Komatsu et al. 2011), Ω_{M} = 0.271, Ω_{Λ} = 0.729. We fixed the redshift difference δz = 0.07 that corresponds approximately to our shell diameter of 200 h^{1} Mpc, and found that this gives distance differences of only a 0.3 and 0.8 per cent at the near and far borders of our sample (the MICE distances are larger than the WMAP7 ones in each case). Hence, for our nearby volume, the effect is small, and does not influence our results given that the statistical uncertainties are much larger (see next section). However, this effect will be significant for deep samples.
5. Results for the SDSS samples
We perform the calculation of B(R,s) for the SDSS in an analogous way to the case of the MICE simulation, using the samples defined in Sect. 4. Our results are shown in Fig. 6. As above, we mask the region R < 2s. As we do not introduce any border correction when calculating the B(R,s) statistic, we also mask the region corresponding to the values R > r_{buf} − s. Values obtained at such large values of R could contain some spurious signals, as the calculation of W_{R,s} would rely on the density field in regions outside the survey boundaries.
Fig. 6 The BAOlet statistic B calculated for SDSS data as a function of the parameters (R,s). The bottom panel shows the results in the full parameter space considered, where we sampled both R and s at intervals of 1 h^{1} Mpc. We mask two areas, at the upper right and left corners, where our results are not reliable (see details in the text). The contours are drawn in steps of 5 for B < 0 (dotted), B = 0 (solid), and B > 0 (dashed). The top two panels show cuts at the arbitrarily chosen values s = 36 h^{1} Mpc (top) and s = 20 h^{1} Mpc (middle), marked with grey horizontal lines in the 2D panel. In these panels, the blue line is B(R,s), while the green line and the red band show the mean () and 1σ^{MC} interval for the Monte Carlo realisations of random centres. We obtain a clear significant peak at different values of s, with a maximum for R = 116 h^{1} Mpc, s = 36 h^{1} Mpc. 

Open with DEXTER 
The resulting B(R,s) map is qualitatively very similar to that expected, either using an analytical ΛCDM model (Fig. 2) or the MICE simulation (Figs. 3 and 4). This indicates that the observed pattern does not originate from spurious features in the SDSS but is closely related to the largescale structure and more specifically the BAO. The B(R,s) attains a maximum at R_{max} = 116 h^{1} Mpc, s_{max} = 36 h^{1} Mpc. This maximum is clearly related to the characteristics of the BAO structures present in our samples. We studied the robustness of this result by changing the weighting scheme applied to construct the density map (see Sect. 4). We do so by capping at different maximum values the possible weights associated with each galaxy, and repeating the calculation of B(R,s) in each case. The results are qualitatively similar, a peak in B(R,s) being obtained in all cases. However, the position of the peak changes in each case, with maximum changes of the order of ± 5 h^{1} Mpc in R_{max}, and ± 10 h^{1} Mpc in s_{max}. Therefore, the difference between the position of the peak obtained from the SDSS data and that given by the MICE simulation is not significant. In any case, we cannot use the scale and the width of the observed maximum of B(R,s) as direct estimates of the radius or width of the shells, especially given that our analysis of the possible observational biases (Fig. 4) was only qualitative.
To assess the significance of the BAO detection with this method, we focused on the value of B(R,s) obtained at the maximum, B_{max} = B(R_{max},s_{max}) = 22.9 ± 3.7^{2}. A more thorough analysis would model the B(R,s) statistic in the full parameter space. However, given the large covariances between measurements at different values of (R,s) we do not expect there to be a large difference from the simple case we consider. We assess the probability of finding such a maximum in the case in which there are no baryon acoustic structures present in our sample. We model this null hypothesis by using randomly distributed centres for the calculation of B(R,s) in Eq. (5), instead of LRG. Even using the W_{R,s}(x) coefficients from the observed density field (traced by SDSS “main” galaxies), the expected value of B(R,s) in this case is 0 (see Eq. (4)), and we expect to obtain a significantly stronger signal in the data. In this way, we test the null hypothesis that, either there are no shelllike structures in the density field traced by the “main” sample, or these shelllike structures are not found preferentially around LRG centres. In either case, that would mean that there are no BAOlike structures present in our sample.
To perform the significance test, we generate 10^{5} random realisations of a Poisson process, with the mean number of points N_{LRG}, in the same volume as the LRG considered in the calculation (i.e. taking into account the buffer zone). For a realisation j, we use the generated points as our centres to compute the B(R,s) statistic following Eq. (5), using the W_{R,s}(x) coefficients obtained from the data. We can then obtain the mean value , and the standard deviation σ^{MC}(R,s) of the MonteCarlo realisations of the centres. We show and a band of 1σ^{MC}(R,s) around it in the top panels of Fig. 6.
We now calculate our signaltonoise ratio at the maximum as , and assess the probability of finding such a large value of SNR_{max} anywhere in the parameter space for the Monte Carlo realisations. We use SNR_{max} instead of directly using B_{max} because for some regions of parameter space, especially at low s, σ^{MC}(R,s) is extremely large. Therefore, if we used B_{max}, we would need to arbitrarily restrict the parameter space studied, thus possibly introducing a posteriori bias. When using SNR_{max}, we sample the full parameter space considered in the calculations (as shown in Fig. 6). We compute the maximum value of SNR for each realisation j in the full (R,s) range, . The distribution of the values of is shown in Fig. 7, where it is compared to the value of SNR_{max} obtained in the real data. We find that only one of the realisations give a value of larger than SNR_{max}. Thus, the probability of obtaining such a large maximum value of SNR in the absence of baryon acoustic structures (our null hypothesis) is p ≃ 10^{5}, equivalent to a ~4.4σ detection in the Gaussian case.
However, we emphasise that the significance found in this work cannot be compared directly to other detection levels found in the literature, as it has been stated in the introduction. In particular, we do not compare our results with an analytical noBAO model of B(R,s) (such as that shown in the top panel of Fig. 2), since to do so would require the detailed modelling of all the selection effects affecting the two samples used.
Fig. 7 Histogram showing the distribution of the maximum SNR values obtained, in the full (R,s) space, for the 10^{5} Monte Carlo realisations of Poissondistributed centres (). This histogram has a mean of 3.72 and a standard deviation of 0.46. We show as a dashed vertical line the value obtained from the data (using the LRG as centres), SNR_{max} = 6.60. Only one of the Monte Carlo realisations give a maximum value larger than SNR_{max}. 

Open with DEXTER 
Fig. 8 A 2D projection of the distribution in redshift space of the LRG used in the analysis (i.e. inside our buffer region)^{3}. 

Open with DEXTER 
Basic characteristics (J2000 sky coordinates, redshift and redshift uncertainty) of the 10 LRG with the larger values of the BAOlet coefficient at the maximum W_{max}.
As explained in Sect. 2, we can extract more information about the BAO phenomenon in our samples by making further use of the BAOlet coefficient maps W_{R,s}(x). Here, we use the coefficient values at the positions of the LRG for the parameters R_{max}, s_{max}, which correspond to the characteristics of the BAO shells present in our samples. In this way, the values W_{max} ≡ W_{Rmax,smax} are a measure of the strength of the signal coming from a BAO shell around a given point, and in particular, a given LRG. Therefore, using W_{max} we can localise in configuration space the regions of the volume covered by our samples where the BAO signal mostly originates.
We illustrate this idea in Fig. 8, where we plot a 2D projection of the distribution of the LRG used as centres in our analysis. We also show the value of W_{max} for each of them according to a colour gradient, the highest values of W_{max} corresponding to the red points in the plot.
In Table 1, we provide the 10 LRG used as centres with the larger values of W_{max}. The whole catalogue of the N_{LRG} = 1599 LRG used as our centres, and the value of W_{max} obtained for each of them is available online. This catalogue could be used to study the relation of the BAO signal at a given LRG to its properties or the environment. It could also be used to select LRG centres with strong signal, and use them to refine the measurements of the BAO characteristics.
As an illustration of this later use, we show a simple way to study the overall properties of the BAO structures, its shape and scale. We select the centres that we know display a prominent acoustic feature, i.e., those for which W_{max} > 0. This leaves us with N_{r} = 809 centres. To improve the signaltonoise ratio in this illustration when studying the BAO structures, we stack together the 3D density maps around the N_{r} selected LRG. In doing so, we keep the lineofsight direction aligned for all the centres, as this direction defines the possible anisotropies in the distribution. We show a 3D view and a 2D cut of this stacked density map in Fig. 9. Thanks to this selection, the characteristic elements of the BAO are amplified, namely on the one hand, a central bump of high density corresponding to the massive halo traced by the LRG, and on the other hand, the shell surrounding it of scale ~109 h^{1} Mpc, which corresponds to a fainter overdensity. We also observe the anisotropic nature of these structures. This is a combination of our working in redshift space, and the redshiftdependent selection function for “main” galaxies.
Fig. 9 Stacked 3D and 2D density fields. In the top three panels, we show the density field after stacking the N_{r} centres with W_{Rmax,smax} > 0. We show surfaces encompassing the regions above different thresholds in density after an isotropic Gaussian smoothing with σ = 10 h^{1} Mpc. The density threshold decreases from top to bottom, with values of δ = 1.24,1.18,1.13 (δ is the density relative to the average density of the sample). We show only the bottom half of the density field for clarity. It can be seen that the acoustic shell clearly appears around the central overdensity at the detected horizon scale. A 2D slice of this density field is shown in the bottom panel. Here, the dotted line is a circle whose radius corresponds to the one we measure for the BAO shells, r_{max} = 109.5 h^{1} Mpc, and the arrow marks the direction of the line of sight. 

Open with DEXTER 
A simpler view can be obtained by calculating the average radial density profile ρ(r) around the N_{r} centres. The resulting profile, shown in Fig. 10, has the same features as the 3D view, namely a high bump on short scales and a clear peak on about the acoustic scale, with a maximum at r_{max} = 109.5 ± 3.9 h^{1} Mpc. The error in r_{max} was estimated using bootstrap realisations (Lupton 1993). This scale gives the radius of the baryon acoustic shells, and is therefore a good estimator of the acoustic scale in the sample. We also show in Fig. 10 the radial profiles restricted to different regions of the sphere, to better characterise the anisotropy of the distribution. We define two cones with a width of 45° with respect to the lineofsight in each direction (we call these “near” and “far” regions), and a “transverse” region covering the belt between the cones. We obtain for each of these regions qualitatively similar results. As expected, we see how the “near” and “far” subsamples are more strongly affected by observational effects, such as redshiftspace distortions, which are more significant along the line of sight. In contrast, the result for the “transverse” subsample matches, within the errors, that for the full sphere. We note that the value of r_{max} is slightly larger for the “far” sample than for the “near” one.
We emphasise that this approach would be impossible with any statistical BAO detection method used this far, since the spatial localisation of the shells is completely lost in either the correlation function or the power spectrum, while the local nature of the wavelet approach has allowed us to identify the positions of the most representative structures in our sample. Moreover, we measure the acoustic scale at positions selected for their low contamination from other structures, which is not the case when averaging over the full sample. In this way, we maximise the BAO signal, while minimising the effect of signals coming from different largescale structures.
Fig. 10 Average radial profiles. We show the radial density profile averaged over the N_{r} centres with W_{Rmax,smax} > 0. We plot , where ρ_{0} is the average density of the sample. The continuous line with a red error band shows the radial profile for the full sphere. We also show the error bands for the radial profile restricted to regions of the sphere, as defined in the text: the “near” region (green), the “far” region (orange), and the “transverse” region (blue). The arrow signals the location of the maximum, r_{max} = 109.5 h^{1} Mpc. The error band corresponds to the 1σ dispersion of 1000 bootstrap realisations. The profiles were estimated using a B_{3} kernel of width h = 4 h^{1} Mpc in the radial coordinate, but similar results are obtained when using slightly different widths or shapes of the kernel. 

Open with DEXTER 
6. Discussion and conclusions
We have designed a new method to detect BAO in the galaxy distribution and to localise, in configuration space, the structures responsible for them. This method is based on the use of a specially designed wavelet applied directly to the density field. Our approach also relies on the use of two different tracers: one for the overall density distribution, and the other for the central overdensities of the baryon acoustic structures.
After testing the method with simulations, we have applied this method to the detection of baryon acoustic structures in a sample drawn from the SDSS. In this case, we have used galaxies from the “main” catalogue to trace the overall density field, and galaxies from the LRG catalogue to trace the location of massive dark matter haloes. We clearly detected BAO in the sample providing a confirmation of the detection obtained previously using general twopoint statistics (the power spectrum and correlation function). Our approach provides an independent method for the detection. Finally, we showed how this method allows us, through the use of W_{max}(x), to localise in configuration space the actual structures responsible for the BAO signal obtained. This is a consequence of using a wavelet acting directly on the density field. We illustrated the utility of this approach by showing the density distribution stacked around a set of centres known to show the BAO feature given their W_{max} value.
Recent works have proposed alternative methods to study the BAO based on wavelets (Xu et al. 2010; Tian et al. 2011). In particular, Tian et al. use a Mexican hat wavelet function with two parameters, conceptually similar to ours. They use it to search for a peak in the two point correlation function of the “main” SDSS sample, obtaining a detection with a pvalue p = 0.002 (equivalent to 3.1σ in the Gaussian case). As in our case, this shows the utility of using the “main” sample to reduce the shot noise in the calculation and obtain significant detections. However, these works apply the wavelet to the measured twopoint correlation function, instead of directly to the density field. In this way, these studies take advantage of the capabilities of the wavelets to characterise accurately the BAO signal (in terms of radius and width), but are unable to get any information about the localisation of these structures in space.
The use of wavelets directly in the density field isolates valuable information about the baryon acoustic structures that is hidden in the standard twopoint statistics. In particular, it gives us information, in terms of the coefficients W_{R,s}(x), to localise regions in the sampled volume containing either the strongest or weakest signal. We expect that this new method for studying BAO will be of much use to ongoing and planned surveys, such as the WiggleZ Survey (Drinkwater et al. 2010), the Baryon Oscillation Spectroscopic Survey (BOSS, Eisenstein et al. 2011), or the Physics of the Accelerating Universe (PAU) Survey (Benítez et al. 2009), which will cover a much larger volume than studied here and explore higher redshifts.
Acknowledgments
This work has been supported by the European Research Council grant SparseAstro (ERC 228261), by the Spanish CONSOLIDER projects AYA200614056, CSD200700060 and AYA201022111C0302, including FEDER contributions, by the Generalitat Valenciana project of excellence PROMETEO/2009/064, and by the Estonian grants SF0060067s08 and ETF8005. P.A.M. was supported by the Spanish Ministerio de Educación through a FPU grant, and by an ERC StG Grant (DEGAS259586). We acknowledge the use of data from the MICE simulations, publicly available at http://www.ice.cat/mice. We also acknowledge the use of public data from SDSS. Funding for the SDSS and SDSSII has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the US Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the MaxPlanckInstitute for Astronomy (MPIA), the MaxPlanckInstitute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington.
References
 Abazajian, K. N., AdelmanMcCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543 [NASA ADS] [CrossRef] [Google Scholar]
 Bassett, B. A., & Afshordi, N. 2010 [arXiv:1005.1664] [Google Scholar]
 Bassett, B. A., & Hlozek, R. 2010, in Dark Energy: Observational and Theoretical Approaches, ed. P. RuizLapuente (Cambridge: Cambridge University Press), 246 [Google Scholar]
 Benítez, N., Gaztañaga, E., Miquel, R., et al. 2009, ApJ, 691, 241 [NASA ADS] [CrossRef] [Google Scholar]
 Blake, C., & Glazebrook, K. 2003, ApJ, 594, 665 [NASA ADS] [CrossRef] [Google Scholar]
 Blanton, M. R., Schlegel, D. J., Strauss, M. A., et al. 2005, AJ, 129, 2562 [NASA ADS] [CrossRef] [Google Scholar]
 Cabré, A., & Gaztañaga, E. 2009a, MNRAS, 393, 1183 [NASA ADS] [CrossRef] [Google Scholar]
 Cabré, A., & Gaztañaga, E. 2009b, MNRAS, 396, 1119 [NASA ADS] [CrossRef] [Google Scholar]
 Cabré, A., & Gaztañaga, E. 2011, MNRAS, 412, L98 [NASA ADS] [Google Scholar]
 Cole, S., Percival, W. J., Peacock, J. A., et al. 2005, MNRAS, 362, 505 [NASA ADS] [CrossRef] [Google Scholar]
 Colless, M., Dalton, G., Maddox, S., et al. 2001, MNRAS, 328, 1039 [NASA ADS] [CrossRef] [Google Scholar]
 Cooray, A., Hu, W., Huterer, D., & Joffre, M. 2001, ApJ, 557, L7 [NASA ADS] [CrossRef] [Google Scholar]
 Crocce, M., & Scoccimarro, R. 2008, Phys. Rev. D, 77, 023533 [NASA ADS] [CrossRef] [Google Scholar]
 Crocce, M., Fosalba, P., Castander, F. J., & Gaztañaga, E. 2010, MNRAS, 403, 1353 [NASA ADS] [CrossRef] [Google Scholar]
 Drinkwater, M. J., Jurek, R. J., Blake, C., et al. 2010, MNRAS, 401, 1429 [NASA ADS] [CrossRef] [Google Scholar]
 Eisenstein, D. J., & Hu, W. 1998, ApJ, 496, 605 [NASA ADS] [CrossRef] [Google Scholar]
 Eisenstein, D. J., Seo, H., & White, M. 2007, ApJ, 664, 660 [NASA ADS] [CrossRef] [Google Scholar]
 Eisenstein, D. J., Weinberg, D. H., Agol, E., et al. 2011, AJ, 142, 72 [NASA ADS] [CrossRef] [Google Scholar]
 Eisenstein, D. J., Annis, J., Gunn, J. E., et al. 2001, AJ, 122, 2267 [NASA ADS] [CrossRef] [Google Scholar]
 Eisenstein, D. J., Zehavi, I., Hogg, D. W., et al. 2005, ApJ, 633, 560 [NASA ADS] [CrossRef] [Google Scholar]
 Fosalba, P., Gaztañaga, E., Castander, F. J., & Manera, M. 2008, MNRAS, 391, 435 [NASA ADS] [CrossRef] [Google Scholar]
 Freedman, W. L., Madore, B. F., Gibson, B. K., et al. 2001, ApJ, 553, 47 [NASA ADS] [CrossRef] [Google Scholar]
 Gaztañaga, E., Cabré, A., & Hui, L. 2009, MNRAS, 399, 1663 [NASA ADS] [CrossRef] [Google Scholar]
 Hu, W., Sugiyama, N., & Silk, J. 1997, Nature, 386, 37 [NASA ADS] [CrossRef] [Google Scholar]
 Hütsi, G. 2006, A&A, 449, 891 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jones, B. J. T. 2009, in Data Analysis in Cosmology, ed. V. J. Martínez, E. Saar, E. MartínezGonzález, & M. J. PonsBordería (Berlin: Springer Verlag), Lect. Notes Phys., 665, 3 [Google Scholar]
 Kazin, E. A., Blanton, M. R., Scoccimarro, R., McBride, C. K., & Berlind, A. A. 2010a, ApJ, 719, 1032 [NASA ADS] [CrossRef] [Google Scholar]
 Kazin, E. A., Blanton, M. R., Scoccimarro, R., et al. 2010b, ApJ, 710, 1444 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Lupton, R. 1993, Statistics in Theory and Practice (Princeton, New Jersey: Princeton University Press) [Google Scholar]
 Mallat, S. 2008, A Wavelet Tour of Signal Processing, The Sparse Way, 3rd edn. (Academic Press) [Google Scholar]
 Martínez, V. J., Paredes, S., & Saar, E. 1993, MNRAS, 260, 365 [NASA ADS] [Google Scholar]
 Martínez, V. J., ArnalteMur, P., Saar, E., et al. 2009, ApJ, 696, L93; Erratum: 2009, 703, L184 [NASA ADS] [CrossRef] [Google Scholar]
 Peebles, P. J. E., & Yu, J. T. 1970, ApJ, 162, 815 [NASA ADS] [CrossRef] [Google Scholar]
 Percival, W. J., Cole, S., Eisenstein, D. J., et al. 2007, MNRAS, 381, 1053 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Percival, W. J., Reid, B. A., Eisenstein, D. J., et al. 2010, MNRAS, 401, 2148 [NASA ADS] [CrossRef] [Google Scholar]
 Rauzy, S., LachiezeRey, M., & Henriksen, R. N. 1993, A&A, 273, 357 [NASA ADS] [Google Scholar]
 Reid, B. A., Percival, W. J., Eisenstein, D. J., et al. 2010, MNRAS, 404, 60 [NASA ADS] [CrossRef] [Google Scholar]
 Saar, E. 2009, in Data Analysis in Cosmology, ed. V. J. Martínez, E. Saar, E. MartínezGonzález, & M. J. PonsBordería (Berlin: Springer Verlag), Lect. Notes Phys., 665, 523 [Google Scholar]
 Sánchez, A. G., Baugh, C. M., & Angulo, R. 2008, MNRAS, 390, 1470 [NASA ADS] [Google Scholar]
 Sánchez, A. G., Crocce, M., Cabré, A., Baugh, C. M., & Gaztañaga, E. 2009, MNRAS, 400, 1643 [NASA ADS] [CrossRef] [Google Scholar]
 Seo, H.J., & Eisenstein, D. J. 2005, ApJ, 633, 575 [NASA ADS] [CrossRef] [Google Scholar]
 Springel, V. 2005, MNRAS, 364, 1105 [NASA ADS] [CrossRef] [Google Scholar]
 Starck, J.L., & Murtagh, F. 2006, Astronomical Image and Data Analysis (Springer), 2nd edn. [Google Scholar]
 Starck, J., Moudden, Y., Abrial, P., & Nguyen, M. 2006, A&A, 446, 1191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Starck, J.L., Murtagh, F., & Fadili, M. 2010, Sparse Image and Signal Processing (Cambridge University Press) [Google Scholar]
 Stoughton, C., Lupton, R. H., Bernardi, M., et al. 2002, AJ, 123, 485 [NASA ADS] [CrossRef] [Google Scholar]
 Strauss, M. A., Weinberg, D. H., Lupton, R. H., et al. 2002, AJ, 124, 1810 [NASA ADS] [CrossRef] [Google Scholar]
 Tegmark, M., Eisenstein, D. J., Strauss, M. A., et al. 2006, Phys. Rev. D, 74, 123507 [NASA ADS] [CrossRef] [Google Scholar]
 Tian, H. J., Neyrinck, M. C., Budavári, T., & Szalay, A. S. 2011, ApJ, 728, 34 [NASA ADS] [CrossRef] [Google Scholar]
 Vielva, P., MartínezGonzález, E., Barreiro, R. B., Sanz, J. L., & Cayón, L. 2004, ApJ, 609, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Xu, X., White, M., Padmanabhan, N., et al. 2010, ApJ, 718, 1224 [NASA ADS] [CrossRef] [Google Scholar]
 York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 120, 1579 [NASA ADS] [CrossRef] [Google Scholar]
 Zheng, Z., Zehavi, I., Eisenstein, D. J., Weinberg, D. H., & Jing, Y. P. 2009, ApJ, 707, 554 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Basic characteristics (J2000 sky coordinates, redshift and redshift uncertainty) of the 10 LRG with the larger values of the BAOlet coefficient at the maximum W_{max}.
All Figures
Fig. 1 The BAOlet function. Here we show a twodimensional (2D) plot (bottom) of the wavelet ψ_{R,s}(x) used in the analysis, as defined by Eq. (1). The top panel shows a 1D slice along the dasheddotted axis. The wavelet is plotted here for R = 105 h^{1} Mpc, s = 30 h^{1} Mpc. The red dot marks the centre of the wavelet. This function has a null mean (provided that R > 2s), and compact support. It takes nonzero values only for R − 2s ≤ x ≤ R + 2s. 

Open with DEXTER  
In the text 
Fig. 2 Values of the BAOlet coefficients W_{R,s} expected at the positions of large initial pointlike perturbations, as a function of the BAOlet parameters (R,s). The bottom panel shows the result using a standard ΛCDM transfer function, while the top panel shows the result using a transfer function with the BAO wiggles edited out (Eisenstein & Hu 1998). The normalisation is arbitrary. The contours are drawn in steps of 1000 for W_{R,s} < 0 (dotted), W_{R,s} = 0 (solid), and W_{R,s} > 0 (dashed). The map attains a maximum at R = 110 h^{1} Mpc, s = 22 h^{1} Mpc. 

Open with DEXTER  
In the text 
Fig. 3 The BAOlet statistic B calculated for the MICE simulation sample described in the text as a function of the parameters (R,s) (bottom panel). The contours are drawn in steps of 5 for B < 0 (dotted), B = 0 (solid), and B > 0 (dashed). This function attains its maximum for R = 108 h^{1} Mpc, s = 28 h^{1} Mpc. The top two panels show cuts at the values s = 28 h^{1} Mpc (top) and s = 22 h^{1} Mpc (middle), marked with grey horizontal lines in the 2D panel. In each case, the solid blue line corresponds to the value obtained from MICE, the dashed red line corresponds to the theoretical expectation from the Eisenstein & Hu (1998) transfer function (bottom panel of Fig. 2), and the dotted green line to the theoretical expectation using the “no wiggle” transfer function (top panel of Fig. 2). These theoretical predictions have been renormalised to get the same value at the maximum in B(R,s). 

Open with DEXTER  
In the text 
Fig. 4 The BAOlet statistic B for the MICE simulation when some observational effects are taken into account. In the top panel, we show the B(R,s) obtained when redshiftspace distortions are included in the simulation. In the bottom panel, we show the result when a radial selection function is applied to the halo catalogue. In both cases, the contours are drawn in steps of 5 in B. Solid contours correspond to the results with the observational effects included. The dashed contours correspond to the original realspace result without selection, i.e., they are identical to those in the bottom panel of Fig. 3. 

Open with DEXTER  
In the text 
Fig. 5 The SDSS catalogues used. We show a sixdegree thick slice drawn from the catalogues used in our analysis. The red points are the LRG galaxies, which form a nearly volumelimited sample at z > 0.15. The “main” galaxies of the SDSS are depicted in blue. As shown, we use the “main” galaxies from z > 0.09, and we use only the LRG in an inner volume, allowing for a buffer region of r_{buff} = 175 h^{1} Mpc from any of the borders of the “main” sample volume. For our analysis, we use the samples covering a total area of 5511deg^{2}. The radius of a typical BAO shell is shown as a segment. At the top insets, we show two orthogonal slices of width 20 h^{1} Mpc, centred on a particular LRG. The BAOlet coefficient map W_{R,s}(x) has a large value at the position of this centre (for R_{max},s_{max}), thus we expect to find a strong BAO signal. The two circles have the radii of 100 and 110 h^{1} Mpc, respectively. This is a single BAO shell where the overdensity can be appreciated by eye on the right scale. 

Open with DEXTER  
In the text 
Fig. 6 The BAOlet statistic B calculated for SDSS data as a function of the parameters (R,s). The bottom panel shows the results in the full parameter space considered, where we sampled both R and s at intervals of 1 h^{1} Mpc. We mask two areas, at the upper right and left corners, where our results are not reliable (see details in the text). The contours are drawn in steps of 5 for B < 0 (dotted), B = 0 (solid), and B > 0 (dashed). The top two panels show cuts at the arbitrarily chosen values s = 36 h^{1} Mpc (top) and s = 20 h^{1} Mpc (middle), marked with grey horizontal lines in the 2D panel. In these panels, the blue line is B(R,s), while the green line and the red band show the mean () and 1σ^{MC} interval for the Monte Carlo realisations of random centres. We obtain a clear significant peak at different values of s, with a maximum for R = 116 h^{1} Mpc, s = 36 h^{1} Mpc. 

Open with DEXTER  
In the text 
Fig. 7 Histogram showing the distribution of the maximum SNR values obtained, in the full (R,s) space, for the 10^{5} Monte Carlo realisations of Poissondistributed centres (). This histogram has a mean of 3.72 and a standard deviation of 0.46. We show as a dashed vertical line the value obtained from the data (using the LRG as centres), SNR_{max} = 6.60. Only one of the Monte Carlo realisations give a maximum value larger than SNR_{max}. 

Open with DEXTER  
In the text 
Fig. 8 A 2D projection of the distribution in redshift space of the LRG used in the analysis (i.e. inside our buffer region)^{3}. 

Open with DEXTER  
In the text 
Fig. 9 Stacked 3D and 2D density fields. In the top three panels, we show the density field after stacking the N_{r} centres with W_{Rmax,smax} > 0. We show surfaces encompassing the regions above different thresholds in density after an isotropic Gaussian smoothing with σ = 10 h^{1} Mpc. The density threshold decreases from top to bottom, with values of δ = 1.24,1.18,1.13 (δ is the density relative to the average density of the sample). We show only the bottom half of the density field for clarity. It can be seen that the acoustic shell clearly appears around the central overdensity at the detected horizon scale. A 2D slice of this density field is shown in the bottom panel. Here, the dotted line is a circle whose radius corresponds to the one we measure for the BAO shells, r_{max} = 109.5 h^{1} Mpc, and the arrow marks the direction of the line of sight. 

Open with DEXTER  
In the text 
Fig. 10 Average radial profiles. We show the radial density profile averaged over the N_{r} centres with W_{Rmax,smax} > 0. We plot , where ρ_{0} is the average density of the sample. The continuous line with a red error band shows the radial profile for the full sphere. We also show the error bands for the radial profile restricted to regions of the sphere, as defined in the text: the “near” region (green), the “far” region (orange), and the “transverse” region (blue). The arrow signals the location of the maximum, r_{max} = 109.5 h^{1} Mpc. The error band corresponds to the 1σ dispersion of 1000 bootstrap realisations. The profiles were estimated using a B_{3} kernel of width h = 4 h^{1} Mpc in the radial coordinate, but similar results are obtained when using slightly different widths or shapes of the kernel. 

Open with DEXTER  
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.