# Coherent Combination of Baryon Acoustic Oscillation Statistics and Peculiar Velocity Measurements from Redshift Survey

## Abstract

New statistical method is proposed to coherently combine Baryon Acoustic Oscillation statistics (BAO) and peculiar velocity measurements exploiting decomposed density–density and velocity–velocity spectra in real space from the observed redshift distortions in redshift space, 1) to achieve stronger dark energy constraints, and , which are enhanced from BAO or velocity measurements alone, and 2) to cross–check consistency of dark energy constraints from two different approaches; BAO as geometrical measurements and peculiar velocity as large scale structure formation observables. In addition to those advantages, as power spectra decomposition procedure is free from uncertainty of galaxy bias, this simultaneous fitting is an optimal method to extract cosmological parameters without any pre–assumption about galaxy bias.

###### pacs:

draft^{1}

## I Introduction

The evolution of large scale structure, as revealed in the clustering of galaxies observed in wide–deep redshift surveys, has been one of the key cosmological probes. Structure formation is driven to grow by a competition between gravitational attraction and the expansion of space-time. This enables us to test our model of gravity at cosmological scales as well as the expansion history of the Universe. Although galaxies are not an unbiased tracer of the large scale matter distribution predictable from linear gravitational theory, various statistical methods using galaxy redshift survey have been developed to probe the nature of cosmic acceleration such as dark energy or modified gravity. Acoustic peak structure imprinted on density field inhomogeniety provides mapping to the large scale power spectra in transverse and longitudinal directions. Before decoupling, the photons and baryons are tightly coupled to form plasma, and overdensity with baryon acoustic oscillation feature (BAO) develops through the competition between outward pressure support of radiation and gravitational attraction. After decoupling, as baryonic matter is influenced only by gravitational force at large scale, the baryons and dark matter instabilities form a configuration including predictable BAO pattern extending from the sound horizon measured by CMB at decoupling epoch. This appears to be a standard ruler tracing the relation between redshift and expansion rate Seo and Eisenstein (2003); Blake and Glazebrook (2003); Linder (2003); Hu and Haiman (2003); Cooray (2004); Matsubara (2004); Amendola et al. (2005); Blake and Bridle (2005); Glazebrook and Blake (2005); Dolney et al. (2006); Seo and Eisenstein (2007); Jeong and Komatsu (2009); Simpson and Peacock (2010); Samushia et al. (2010). Additionally, as galaxies are expected to act nearly as test particles within the cosmological matter flow, the motions of galaxies carry an imprint on the rate of growth of large-scale structure and allows us to constrain cosmological models Wang (2008); Song and Percival (2009); McDonald and Seljak (2009); Bean and Tangmatitham (2010); Jennings et al. (2010); Taruya et al. (2010); Song et al. (2010). Despite the intensive studies of both methods constraining cosmological parameters, the coherent way to combine both has not yet been fully optimized.

Both different statistics can be combined coherently exploiting decomposed density–density and velocity–velocity spectra in real space from the observed redshift distortions in redshift space. The Fisher matrix exposition was developed White et al. (2008), assuming strict functional forms for the power spectra or allowing them to float freely. As expected the constraints are tightest when theoretical investigations can provide good priors for the form of distorted power spectra in redshift space, but even relatively conservative assumptions suggest that percent level of decomposition should be possible with future surveys White et al. (2008); Song and Kayo (2010). The cosmological distance errors are achievable using BAO preserved on decomposed density–density spectra, which is degraded a few factors of order from BAO of full observed power spectra. Decrement is compensated by adding cosmological constraints from the decomposed velocity–velocity spectra utilizing the growth of large-scale structure as imprinted on dynamics of galaxies observed in large redshift surveys. Tighter constraints are achieved by this combination as well as dark energy constraints can be cross–checked by two different alternative approaches; distance measures and large scale structure formation observables.

The detailed formalism is presented in the next section. The Fisher matrix analysis to decompose spectra is briefly reviewed, then BAO statistics is redefined with the decomposed density–density spectra. Coherent formalism to combine BAO and decomposed velocity–velocity spectra is introduced, and advantages are discussed to gain stronger constraints on dark energy models with future surveys. For illustration, a fiducial CDM cosmology is assumed with , , and (in good agreement with a variety of observations) when computing specific predictions for future surveys.

## Ii Coherent approach to combine BAO and peculiar velocity

The formalism is briefly reviewed to forecast errors on galaxy–galalxy density spectra, , and velocity–velocity spectra, ( is divergence of velocity map in unit of , ) from the observed spectra, , of redshift surveys. We explain in detail how to estimate error forecast of BAO from the decomposed density–density spectra as well as how to combine decomposed peculiar velocity measurement coherently.

### ii.1 Decomposition of the observed power spectra in redshift space

The observed two-point correlation function in redshift space is decomposed into spectra of density fluctuations and peculiar velocity fields in real space White et al. (2008). The power spectra in redshift space, , are given by,

(1) | |||||

where denotes the cosine of angle between orientation of two point correlation and the line of sight. The cross-correlation coefficient is defined as . The density and velocity fields are highly correlated for thus we assume that the density and velocities are perfectly correlated, . The density-velocity cross-spectrum becomes the geometric mean of the two auto-spectra to leave only two free functions, and . Uncertainty in the observed redshifts is modeled by a line-of-smearing of the structure using fitting function where denotes one–dimensional velocity dispersion (FoG: Finger-of-God effect).

The accuracy of decomposition of and out of is estimated using Fisher matrix analysis determining the sensitivity of a particular measurement. Fisher matrix for this decomposition, , is written as,

(2) |

where . The effective volume is given by,

(3) |

where denotes galaxy number density, here . Comoving volume, , given by each redshift shell from to 2 with () is written as,

(4) |

where and denote comoving distances of outer and inner shell of the given redshift bin respectively. The weight function is given by

(5) |

where is the finger-of-god suppression factor and is a threshold value indicating our confidence in the accuracy of the FoG model. We fix indicating trust of clearance of non-linear effect up to 10 contamination of the factor beyond Kaiser limit. Derivative terms in Eq. (2) are given by,

(6) |

### ii.2 Variances of power spectra

The detectability of galaxy–galaxy density spectra in redshift space, , is given by,

(7) |

where denotes redshift bin from to 2, and denotes space bin running from to , i.e. integration of in each bin from lower to upper bounds denoted by and respectively.

The detectability of galaxy–galaxy density spectra in real space, , can be estimated from decomposition Fisher matrix . The –component of the inverse matrix of gives the variance of as,

(8) |

In the upper panel of Fig. 1, solid curve represents fractional variance of at the same redshift bin . The decomposition procedure to transfrom measured galaxy–galaxy density spectra in redshift space into real space degrades detectability, but is not very significant about a factor of 1.8.

Variance of decomposed velocity–velocity spectra is derived from –component of the inverse matrix of ,

(9) |

In the lower panel of Fig. 1, solid curve represents fractional variance of at redshift bin of . As is set to be 0.1, the detectability is significantly weakened at . It is observed V–shape of variance pointed around threshold scale of . Due to the short range of detectability, probing BAO structure using proves to be challenge, but, as is not a biased tracer, the measured amplitude can be used to constrain cosmological parameters.

Threshold limit representing cutting–off scale of due to the unpredictable non–linear physics, , reduces detectability of spectra shown in Fig. 1, but differently affects on and . The introduced calibrates Fisher matrix integrand at the combination of not . The decomposition of is less affected, as it is most determined at limit, while is better decomposed at limit. Poorer detectability of at is due to increasing FoG calibration at limit.

### ii.3 Formalism of BAO Fisher matrix

Cosmological constraints are forecasted from a measured position of BAO in the three dimensional power spectra following the formalism of Seo and Eisenstein (2007). The baryonic part of the power spectrum can be modeled as

(10) |

where is the sound horizon at drag epoch and is the Silk damping scale. The Silk damping scale can be accurately fitted by .

As in Seo and Eisenstein (2007), Eq. (10) is multiplied by additional Gaussian functions to account for the erasure of information due to nonlinear evolution. The final is given by

(11) |

where is galaxy power spectrum at , is a normalization factor, models the loss of information due to nonlinear growth, and is defined as,

(12) |

Eq. (10) is compared with the fitting formula of Seo and Eisenstein to estimate . The fiducial cosmology, as chosen in this research, has an analytical estimate of . As in Seo and Eisenstein (2007) the numerical values , where is the growth function normalized as at early time, and for the cosmology with . When the physical value of the sound horizon is known to high precision from CMB measurements the error on is equivalent to the error on .

Centroid approximation is available for , as it is decomposed into spherically symmetric Fourier space. BAO Fisher matrix from is given by,

(13) |

where is the total number of bin upto . The variance of is simply given by of . In Fig. 2, fractional errors of calculated from only are presented.

### ii.4 Tight constraints on dark energy from

The observed power spectra depend not only on dark energy parameters but the entire matter content and primordial power spectrum. We do not assume these quantities to be known. Instead, it is assumed that CMB data are available to constrain them. The CMB power spectra included in our analyses are the spectra of temperature–temperature, temperature–polarization and polarization–polarization. Cosmological parameter space is given by (, , , , , , , ).

With CMB temeprature fixed at , the sound horizon at last scattering surface, , is given in unit by and measured by CMB acoustic peak structure. With the measured angular size of sound horizon at last scattering surface, the angular diameter distance to the last scattering surface is also given in unit. The geometrical factor determination in unit is converted into the Fourier space dimension of in unit. As primordial spectra are tightly constrained by CMB physics, the pivoting scale of primordial spectra is given in unit. In addition to this, as the shape of spectra in terms of (transfer function ) is tightly determined by parameters measured by CMB physics, we express in terms of in unit.

In redshift survey, measured redshift of objects is converted using distance in unit in which mapping is less affected by most uncertain distance parameter of the Hubble constant. The velocity–velocity power spectra are given by in terms of in unit Song et al. (2010),

Variation of dark energy parameter impact on becomes complicated when we consider the combined constraints from CMB and redshift survey in which dependent quantities are measured in different units. First, multiplication factor becomes which changes due to floating of and . Second, primordial tilt factor changes as and are determined by different units. Third, transfer function factor is unmodified in terms of with varying dark energy parameters, as decay due to silk damping is well determined by cosmological parameters measured at last scattering surface with CMB physics. But it is altered due to using unit of instead of for distance measure of redshift survey. Finally, growth factor changes.

In Fig. 3, total change of is presented with varying (upper panel) and (lower panel). If degeneracy is not considered between and , tight constraints on dark energy paramters are observed. Fourier space is split into 10 bins from to . This figure is a representation by one of 10 bins chosen at .

### ii.5 Coherent combination of and

To calculate the expected parameter errors we make a first order Taylor expansion of the parameter dependence including CMB, BAO and decomposed . The expected parameter errors are easily calculated from the inversed covariance matrix of the Fisher matrix using this ‘linear response’ approximation. The linear response approximation is improved and susceptibility to numerical error is reduced with a careful choice of the parameters including . denotes a Fisher matrix of CMB physics equivalent to combination of distance measure at last scattering surface and primordial parameter constraints, thus exploiting unlensed CMB power spectra only. We use Planck experiment specs at three frequencies of 100, 143 and 217 GeV. Each channel has the resolution of 9.1, 7.2, 5.0 arcmin , the temperature sensitivity of 5.5, 6.0 and 13 K and the polarization sensitivity of , 11, 27 K respectively. We calculate using the method developed in our previous paper Kaplinghat et al. (2002).

Constraint on distance measure by BAO is converted into cosmological parameter constraints using the following Fisher matrix formalism,

(14) |

where indices and denote elements in cosmological parameter space. Dark energy constraints from are shown to be the outer dash contour in Fig. 4 and the dash contour in Fig. 5, and , which is weaker than dark energy constraint from BAO of full observed spectra, and (the inner dash contour in Fig. 4). Here error estimation denotes 68 statistical significance.

While constrains dark energy by distance measure, dark energy is probed by signature on the large scale structure growth of peculiar velocity measurements . Shown in the bottom panel of Fig. 1 for fractional errors of , the effective range of is too narrow to observe the baryonic acoustic peak structure. However unlike , it is an unbiased quantity that is predictable with a given set of cosmological parameters. Fisher matrix is given by,

At scales greater than , the uncertainty increases due to the unknown non–linear physics causing the FoG effect and the violation of assumption of perfect cross–correlation between density–density and velocity–velocity spectra Song and Kayo (2010). It was shown that the decomposed starts to be biased as scale approaches to . Cut–off scale needs to be introduced at in order to skim away contributions contaminated from non–linear physics at .

In Fig. 5, dotted contour represents dark energy constraint from . Tight contraints achieved using is worth being investigated. At last redshift bin of , dark energy constraints can be read from Fig. 3 about or which confirms the results from Eq. II.5. If fitting is extended to all redshift bins then or (derived from Eq. II.5). This still agrees with that is observed in Fig. 3. But simultaneous determination of and becomes less tight as and , because of degenerate variation in terms of and . The results shown in Fig. 5 are given after full marginalization with all other cosmological parameters.

Then coherent combination of BAO and peculiar velocity measurements can be simply given by

(15) |

Solid contour of Fig. 5 represents dark energy constraints from , and , which is tighter than any other case considered in this paper.

## Iii Conclusion

The known galaxy bias has always been considered an obstacle utilizing redshift surveys for cosmological purpose. If just peak structure is considered as a standard ruler at the given targeted redshift, cosmological information is extractable with less dependence on bias uncertainty. However enriched information of large scale structure is not fully exploited. Previously, decomposed contributions of density–density and velocity–velocity spectra were studied in the observed redshift distortion using anisotropic feature aligned along the line of sight. In this paper, it is proposed to coherently combine both decomposed spectra for constraining cosmological parameters (still in bias uncertainty independent way) to enhance dark energy constraints.

What is implemented in this new method is to cross–check constraints on dark energy in two different approaches; distance measure by BAO and large scale stucture growth factors by peculiar velocity measurements. There are many different theoretical models proposed to generate cosmic acceleration. It is essential to test those theories in many different observational methods. When we decompose spectra into density–density and velocity–velocity pairs, we can test theoretical models using geometrical tools as well as using large scale structure formation. In addition, dark energy constraints are enhanced from performing BAO alone using density–density spectra before decomposition. In this sense, our approach is optimizing way to exploit redshift surveys for the purpose of constraining dark energy.

It is important to test consistency between geometrical measures and structural observable. As general relativity can be modified at large scales. One of vital mission of dark energy probe is to distinguish two different class of models; modifying energy momentum components or modifying gravity. Coherent motion measurements provide us with a way of tracing Newtonian force sourcing dynamics of gravitational instability. Our new statistical methods assist us testing modified gravity type models too.

It has not yet been fully tested whether the assumption of preservation of baryon acoustic peak structure is true or not, although baryon acoustic structure is not expected to be lost in the decomposition procedure. It is necessary to test this assumption using N-body simulation. For decomposition of velocity measurements, it was tested whether there was velocity bias or not. However at the same time, it was known that velocity–velocity spectra are most sensitive to be extracted without systematic uncertainty. In the future, higher resolution experiments will be launched, and it would be interesting whether velocity bias is not induced yet.

## Acknowledgments

We would like to thank Lado Samuthia for discussion about BAO Fisher matrix, Eric Linder for helpful comments in details, and the referee for improving the manuscript. We thank Korea Institute for Advanced Study for providing computing resources (KIAS linux cluster system) for this work.

### Footnotes

- preprint:

### References

- H.-J. Seo and D. J. Eisenstein, Astrophys. J. 598, 720 (2003), eprint astro-ph/0307460.
- C. Blake and K. Glazebrook, Astrophys. J. 594, 665 (2003), eprint astro-ph/0301632.
- E. V. Linder, Phys. Rev. D68, 083504 (2003), eprint astro-ph/0304001.
- W. Hu and Z. Haiman, Phys. Rev. D68, 063004 (2003), eprint astro-ph/0306053.
- A. Cooray, MNRAS 348, 250 (2004), eprint arXiv:astro-ph/0311515.
- T. Matsubara, Astrophys. J. 615, 573 (2004), eprint astro-ph/0408349.
- L. Amendola, C. Quercellini, and E. Giallongo, Mon. Not. Roy. Astron. Soc. 357, 429 (2005), eprint astro-ph/0404599.
- C. Blake and S. Bridle, Mon. Not. Roy. Astron. Soc. 363, 1329 (2005), eprint astro-ph/0411713.
- K. Glazebrook and C. Blake, Astrophys. J. 631, 1 (2005), eprint astro-ph/0505608.
- D. Dolney, B. Jain, and M. Takada, Mon. Not. Roy. Astron. Soc. 366, 884 (2006), eprint astro-ph/0409445.
- H.-J. Seo and D. J. Eisenstein, Astrophys. J. 665, 14 (2007), eprint astro-ph/0701079.
- D. Jeong and E. Komatsu, Astrophys. J. 691, 569 (2009), eprint 0805.2632.
- F. Simpson and J. A. Peacock, Phys. Rev. D81, 043512 (2010), eprint 0910.3834.
- L. Samushia et al. (2010), eprint 1006.0609.
- Y. Wang, JCAP 0805, 021 (2008), eprint 0710.3885.
- Y.-S. Song and W. J. Percival, JCAP 0910, 004 (2009), eprint 0807.0810.
- P. McDonald and U. Seljak, JCAP 0910, 007 (2009), eprint 0810.0323.
- R. Bean and M. Tangmatitham, Phys. Rev. D81, 083534 (2010), eprint 1002.4197.
- E. Jennings, C. M. Baugh, and S. Pascoli (2010), eprint 1003.4282.
- A. Taruya, T. Nishimichi, and S. Saito (2010), eprint 1006.0699.
- Y.-S. Song, C. G. Sabiu, I. Kayo, and R. C. Nichol (2010), eprint 1006.4630.
- M. White, Y.-S. Song, and W. J. Percival, Mon. Not. Roy. Astron. Soc. 397, 1348 (2008), eprint 0810.1518.
- Y.-S. Song and I. Kayo Mon. Not. Roy. Astron. Soc. 407, 1123 (2010), eprint 1003.2420.
- M. Kaplinghat, L. Knox and Y.-S. Song Phys. Rev. Lett. 91, 241301 (2002), eprint astro-ph/0303344.