Are redshift-space distortions actually a probe of growth of structure?
We present an impact of coupling between dark matter and a scalar field, which might be responsible for dark energy, on measurements of redshift-space distortions. We point out that, in the presence of conformal and/or disformal coupling, linearized continuity and Euler equations for total matter fluid significantly deviate from the standard ones even in the sub-horizon scales. In such a case, a peculiar velocity of total matter field is determined not only by a logarithmic time derivative of its density perturbation but also by density perturbations for both dark matter and baryon, leading to a large modification of the physical interpretation of observed data obtained by measurements of redshift-space distortions. We reformulate galaxy two-point correlation function in the redshift space based on the modified continuity and Euler equations. We conclude from the resultant formula that the true value of the linear growth rate of large-scale structure cannot be necessarily constrained by single-redshift measurements of the redshift-space distortions, unless one observes the actual time-evolution of structure.
The current cosmological observations, such as type Ia supernovae Perlmutter et al. (1997); Riess et al. (1998) and cosmic microwave background Ade et al. (2016), indicate the presence of dark matter and dark energy, which have not been identified yet. The existence of dark matter is also well established by astrophysical observations, which indicate dark matter as a non-luminous and pressure-less fluid with small dispersion velocity Zwicky (1933); Babcock (1939); Kahn and Woltjer (1959); Clowe et al. (2004). The dark energy is responsible for explaining the late-time accelerated expansion of the Universe, and numerous attempts to identify it have been intensively proposed in many literatures. One such candidate is to introduce a scalar degree of freedom as a new contribution to energy-momentum tensor or modification in a gravitational sector (see for reviews e.g. Tsujikawa (2010); Clifton et al. (2012)).
When the ordinary matter, baryon, directly couples with such a scalar degree of freedom, it induces the fifth force. While the fifth force between baryonic matter is tightly constrained by the solar-system experiments Will (2006), this is not true for the dark force that is active only between dark matter since the solar-system experiments do not probe such an interaction. Then, the natural arena for probing such interactions is cosmology. When additional interaction only between the cold dark matter (CDM) is present, the growth rate of the CDM density perturbations would be generically different from that of the baryon density perturbations. We then expect that observing the growth of the CDM density perturbation provides us with rich information about such an interaction.
In the standard treatments, the linear growth rate of large-scale structure is mainly obtained by observing galaxy peculiar velocity field through measurements of redshift-space distortions (RSDs) in galaxy survey. On large scales, where the linear perturbation theory is valid, the galaxy peculiar velocity field is considered to be identical to the velocity field of the total matter. Based on the continuity equation, the matter velocities should be given by the logarithmic time derivative of the density field, that is, the linear growth rate, . Galaxy maps produced by estimating distances from observed radial velocities, which include components from both the Hubble flow and peculiar velocities driven by the clustering of matter, show an anisotropic galaxy distribution. Due to such an effect, the galaxy power spectrum on large scales is known to be enhanced by the factor (named “Kaiser formula”), where with being the linear galaxy bias factor and is the cosine of the angle between the line of sight and the Fourier momentum Kaiser (1987). Although there is a degeneracy between the growth rate and the linear bias factor, this degeneracy can be in principle broken by using e.g., higher-order statistics Scoccimarro et al. (1999) and cross-correlations between other observables Hashimoto et al. (2016) by which the linear bias factor alone can be constrained. Hence, it is widely believed that measurements of RSDs even at single redshift allow direct constraints on the growth rate. Moreover, several attempts show that the relation between the peculiar velocity and the growth rate for each species, which is based on the continuity equation, is valid even for the wide range of cosmological scenarios including modified theories of gravity (see e.g., Gleyzes et al. (2016)). However, as we will show below, this relation is not necessarily correct in more general situations. In this Letter, we would like to address how the above Kaiser formula is modified when the CDM couples with the scalar field.
Let us consider the following invertible metric transformation Bekenstein (1993),
where is the original frame metric, and and are respectively called conformal and disformal factors, which are functions of the scalar field and its kinetic term . Here and hereafter, is a generic scalar field, and we do not specify it though the case with being responsible for dark energy is the most interesting. The action is given by
where represents a Lagrangian for scalar field and a total matter action. For simplicity, we consider the canonical scalar field: and assume the scalar field does not modify the gravitational sector, i.e., the absence of kinetic braiding Deffayet et al. (2010). As for the matter sector, we assume that the baryon is minimally coupled for simplicity while the CDM couples with the scalar field through the barred metric defined in (1). The total matter action is thus given by
where and represent the actions for baryon and CDM, respectively. Due to the non-minimal coupling between the dark matter and the scalar field, baryonic matter and dark matter do not move in the same way.
The variation with respect to the metric leads to the Einstein equations as usual,
Here and hereafter, and . The superscript represents , or for baryon, dark matter and total matter, respectively. The combination of the energy-momentum tensor for total matter and the scalar sector is conserved as The energy-momentum conservation for baryon also takes the familiar form, . On the other hand, the energy-momentum tensors for the scalar field and dark matter no longer satisfy the conservation law individually, and it rather takes the following form, .
The scalar equation is given by
where , which characterizes the coupling between CDM and the scalar field, is defined as
where , for . By the use of Eq. (5), the energy-momentum conservation for CDM and total matter can be recast as
We work on a spatially flat FLRW metric in Newtonian gauge,
and define the background and perturbations of the energy-momentum tensor for the baryon, the dark matter and the total matter as
where . We also split the scalar field as . The background part of the Einstein equation gives
The background equation of motion for (5) yields
and the energy-momentum conservation for baryon and CDM lead to the background equations
where , , and are evaluated at background fields.
In deriving perturbed equations, we use the quasi-static approximation,
which is applicable when the wavelength of perturbations is well inside
the sound Horizon of the scalar field, ,
where is the sound speed of the scalar field.
Although the sound speed of the scalar field generally differs from unity in our setup Kimura et al. (2017), we assume for simplicity.
Then we can neglect time-derivative terms of perturbations
while keeping spatial derivative terms
The continuity and Euler equations for baryon are the standard form;
while those for CDM get modified as follows
where the scalar field perturbation is determined by
In the quasi-static limit, the most relevant terms in can be extracted as
Physical meanings of these functions are as follows: characterizes the strength of the modification of the CDM velocity field from the disformal coupling, and represents the contribution from the scalar field perturbation in . We emphasize that the term gives the non-vanishing contributions only if the conformal and/or disformal factors depend on the kinetic term.
After eliminating the scalar field perturbations by using (26), we can rewrite the modified continuity equation in terms of the CDM density contrast and velocity field as
In the minimal coupling case (), all time dependent coefficients are zero, . When conformal and disformal factors depend only on , we have , , and . Thus the continuity equation is the same as the one in the minimal coupling case. One can verify this property even for a wider class of scalar-tensor theories Gleyzes et al. (2016). When at least one of and depend on , there arises a new contribution of in the continuity equation, and the CDM velocity can significantly differ from the standard case. An important implication from these equations is that the continuity equation for the total matter fluctuations, (13) and (14), is given by
which differs from the standard form by the presence of the non-minimal coupling. We also found that even when the contribution is negligible the standard form of the continuity equation cannot be reproduced due to the second term of the right-hand-side in the second equation, which originates from the deviation of the background energy density from the standard matter (see Eq. (19)). Therefore, we conclude that there are two possibilities to break the standard relation of the continuity equation for the total matter field: One comes from the term in the CDM continuity equation, which appears only when the coupling depends on the kinetic term, and the other corresponds to the deviation of the background dynamics from the standard one characterized by .
Combining all the perturbed equations to eliminate velocities as usual, we obtain two coupled second-order differential equation for the baryon and CDM density contrasts. Since the evolution equations for the baryon and CDM density contrasts are independent of the wavenumber , one can decompose the density contrasts into the (normalized) -independent linear growth factors and initial density contrasts for the baryon, CDM and total matter as follows,
Here, we have chosen the initial time to be much after the time of CMB decoupling () but much before the effect of the dark interaction becomes important () and assumed that the baryon density contrast has caught up with the CDM density contrast by the initial time. We also define the growth rate for each species, , as the logarithmic derivative of the linear growth, that is . Rewriting the continuity equations for the baryon, CDM and total matter in terms of the growth factors, we obtain the suggestive form of velocity potentials
Although one can easily see , from Eq. (30) the effective linear growth rate of the CDM, , can significantly differ from the standard one due to the contribution as
Moreover, by the use of Eq. (32), the effective growth rate of the total matter fluctuations can be written as
with . Although is naturally given by the growth-factor-weighted average of the effective growth rates for CDM and baryon, it does not in general coincide with . As discussed above, its deviation is due to the non-trivial terms in the CDM continuity equation and the background dynamics that produce the second and third terms in (36).
Modified interpretation of Kaiser formula.
In the above investigation, we found that the effective growth rate inferred from the peculiar velocities no longer coincides with the actual growth rate , namely measurements of the peculiar velocity field do not necessarily provide the growth rate of clustering directly. Our example vividly demonstrates that the standard dictionary translating the RSDs measurements into the growth rate is not universal and fails for some classes of theories beyond the CDM model. To see the impact of the breaking of the relation between the peculiar velocities and the actual growth rate, we now focus on the modification of the Kaiser formula as the simplest and most important observable effect of RSDs. The generalization to other observables related to the peculiar velocities is straightforward.
In addition to the Hubble expansion, the peculiar velocities of the galaxies relative to the Hubble expansion distort the distribution of galaxies in the 3-dimensional redshift space, and such effects must carefully be taken into account when comparing galaxy two-point correlation function with theoretical predictions Kaiser (1987). The mapping of the observed redshift position from the real space position is given by
where is a line-of-sight component of the peculiar velocity of a galaxy and is a unit vector of line-of-sight. In Eq. (37), we have assumed the plane-parallel approximation, so that the line-of-sight is taken as a fixed direction, . Recalling that the number of galaxies in the infinitesimal volume of both spaces is invariant, the overdensities in the redshift space and the real space are related through
The galaxy density contrast in the real space, , is related to the total matter density contrast given by Eq. (13), through the standard linear bias model on large scales. The peculiar velocity fields of the galaxies, , on large scales are expected to be related to the CDM and baryon fluid velocities, and the explicit relation is determined by imposing the reasonable physical condition, e.g., momentum conservation law for each galaxy Gleyzes et al. (2016). For simplicity, here, we assume that the peculiar velocity fields of galaxies on large scales are given by the total matter fluid velocities, , given in Eq. (14), as in the standard case: .
Therefore, the resultant galaxy power spectrum in redshift space is given by
where is the real-space galaxy power spectrum, is the power spectrum for the total matter density contrast, and
This is a generalization of the Kaiser formula. In fact, in the minimal coupling case, we have and , and hence Eq. (39) is reduced to the standard Kaiser formula. However, we found that in the presence of the coupling between the CDM and the scalar sector we have no longer the relation as we have discussed, and it means that the RSDs are not trustable probes of growth of structure. It is notable that the RSDs cannot provide the true value of the growth rate even in the simple case where the conformal and disformal factors depend only on . Since in this case, the deviation from the standard formula is proportional to , this effect is suppressed when the background evolution of the dark matter is almost same as the one of the baryon. On the other hand, there is a wider room for sizable modification of the standard Kaiser formula in our general setup; even when either or the baryonic contamination is negligibly small, can differ from by . To see this clearly, let us expand the formula (40) in terms of the baryon-CDM ratio to neglect the ambiguity from the baryon contribution. The leading term gives . This immediately shows the single-redshift RSDs measurements can not give a constraint on the linear growth rate unless the contributions from the couplings is fixed by using other observables. This fact demonstrates that one has to keep this new effect in mind when testing beyond CDM theories by the RSDs measurements. Even if the growth index is obtained from RSDs in future galaxy survey, it is still possible that the true theory is different from the standard CDM model. One way to obtain the actual growth rate of large-scale structure is to directly observe the time-evolution of structure by e.g. multiple redshift observations of galaxy power spectrum. In fact, we have a strong degeneracy between the growth of large scale structure and the redshift-dependence of the linear bias. Thus, to measure by multiple redshift observations, we need to fix the bias for each redshift by using other observations, i.e., cross-correlation between the clustering of galaxies and weak lensing (see, e.g., Hashimoto et al. (2016)). After evaluating the actual growth rate, one can compare the actual and effective growth rates to constrain the couplings between the CDM and scalar field.
We have shown that the additional interaction mediated by the scalar field that operates only between dark matter through conformal and disformal couplings changes the continuity and Euler equations for cosmological perturbations in a non-trivial manner and investigated its impact on RSDs measurements in galaxy survey. We found that the effects of such modifications appear even at sub-horizon scales in the presence of and -dependence of the conformal and/or disformal couplings. The effective linear growth rate, which is inferred from measurements of the peculiar velocities of the distributed galaxies, no longer corresponds to the logarithmic time derivative of the density perturbation and is rather characterized by both the density perturbations and their derivatives for each species in general situation. In other words, the information of the coupling is encoded in the peculiar velocity fields and the true value of the growth rate of large-scale structure cannot necessarily be constrained by the single-redshift RSDs measurements. It can be extracted by using multiple power spectra of the galaxy distribution at different redshift. This fact will play a vital role of measuring the linear growth rate by the RSDs measurement, and it will provide us a rich information of dark matter and dark energy.
We thank Toshifumi Noumi, Masamune Oguri, Jiro Soda, Masahiro Takada, and Kazuhiro Yamamoto for many useful comments and discussions. This work was supported in part by JSPS Grant-in-Aid for Scientific Research Nos. JP17K14304 (D.Y.), JP17K14276 (R.K.), JP25287054 (R.K. & M.Y.), JP26610062 (M.Y.), JP15K17632 (T.S.) and JP15K17659 (S.Y.), MEXT KAKENHI Grant-in-Aid for Scientific Research on Innovative Areas ”Cosmic Acceleration” Nos. JP15H05888 (M.Y.) and JP16H01103 (S.Y.), MEXT KAKENHI Grant Numbers JP17H06359 (T.S.).
While this paper was being completed, Ref. A Borges and Wands (2017) appeared, in which the redshift space distortions in the context of interacting dark matter and vacuum energy are discussed.
- preprint: RUP-17-19
- Note that the pressureless feature of the CDM is robust at least at first order of perturbations even if we take other definitions of energy momentum tensor such as and .
- We also neglect the mass of the scalar field, which could be crucial when the mass of the scalar field is large enough, i.e., . The full analysis including large mass effects will be reported in the full paper Kimura et al. (2017).
- S. Perlmutter et al. (Supernova Cosmology Project), Bull. Am. Astron. Soc. 29, 1351 (1997), eprint astro-ph/9812473.
- A. G. Riess et al. (Supernova Search Team), Astron. J. 116, 1009 (1998), eprint astro-ph/9805201.
- P. A. R. Ade et al. (Planck), Astron. Astrophys. 594, A13 (2016), eprint 1502.01589.
- F. Zwicky, Helvetica Physica Acta, 6, 110 (1933).
- H. W. Babcock, Lick Observatory Bulletin 19, 41 (1939).
- F. D. Kahn and L. Woltjer, ApJ 130, 705 (1959).
- D. Clowe, A. Gonzalez, and M. Markevitch, Astrophys. J. 604, 596 (2004), eprint astro-ph/0312273.
- S. Tsujikawa, Lect. Notes Phys. 800, 99 (2010), eprint 1101.0191.
- T. Clifton, P. G. Ferreira, A. Padilla, and C. Skordis, Phys. Rept. 513, 1 (2012), eprint 1106.2476.
- C. M. Will, Living Rev. Rel. 9, 3 (2006), eprint gr-qc/0510072.
- N. Kaiser, Mon Not R Astron Soc 227, 1 (1987).
- R. Scoccimarro, H. M. P. Couchman, and J. A. Frieman, Astrophys. J. 517, 531 (1999), eprint astro-ph/9808305.
- I. Hashimoto, A. Taruya, T. Matsubara, T. Namikawa, and S. Yokoyama, Phys. Rev. D93, 103537 (2016), eprint 1512.08352.
- J. Gleyzes, D. Langlois, M. Mancarella, and F. Vernizzi, JCAP 1602, 056 (2016), eprint 1509.02191.
- J. D. Bekenstein, Phys. Rev. D48, 3641 (1993), eprint gr-qc/9211017.
- C. Deffayet, O. Pujolas, I. Sawicki, and A. Vikman, JCAP 1010, 026 (2010), eprint 1008.0048.
- R. Kimura, T. Suyama, M. Yamaguchi, D. Yamauchi, and S. Yokoyama, in preparation (2017).
- H. A Borges and D. Wands, arXiv:1709.08933 (2017).