Prepared for submission to JCAP
Galaxy number counts to second order and their bispectrum
Enea Di Dio a , Ruth Durrer a , Giovanni Marozzi a , Francesco Montanari a
aUniversit´ede Gen`eve, D´epartement de Physique Th´eoriqueand CAP, 24 quai Ernest- Ansermet, CH-1211 Gen`eve 4, Switzerland E-mail: [email protected], [email protected], [email protected], [email protected]
Abstract. We determine the number counts to second order in cosmological perturbation theory in the Poisson gauge and allowing for anisotropic stress. The calculation is performed using an innovative approach based on the recently proposed ”geodesic light-cone” gauge. This allows us to determine the number counts in a purely geometric way, without using Einstein’s equation. The result is valid for general dark energy models and (most) modified gravity models. We then evaluate numerically some relevant contributions to the number counts bispectrum. In particular we consider the terms involving the density, redshift space distortion and lensing. arXiv:1407.0376v4 [astro-ph.CO] 8 Jun 2015 Contents
1 Introduction1
2 Galaxy Number Counts3
3 From the Geodesic Light-Cone to the Poisson gauge4
4 Evaluation of the number counts to second order8 4.1 Redshift density perturbation8 4.2 Volume perturbation 10 4.3 Number Counts 19 4.4 Number Counts: leading terms 23
5 The Bispectrum and its numerical evaluation 24 5.1 Density 25 5.2 Redshift space distortions 27 5.3 Lensing 28 5.4 Numerical Results 29
6 Conclusions 32
A Useful relations 34
1 Introduction
Cosmology has entered the precision era. So far, mainly via observations of the cosmic microwave background (CMB) most recently performed by the Planck satellite [1,2]. But also high precision observations of cosmic large scale structure (LSS), i.e. the distribution of matter in the Universe, are being performed, see e.g. [3,4], and are under construction or planning, see [5,6] 1 and references therein. The interpretation of observations of LSS is more complicated than the CMB. First of all we observe galaxies and not the density field, this is the biasing problem, and sec- ondly, fluctuations on small scales become large so that non-linearities become important. Nevertheless, if we can handle these problems, a detailed observation of LSS is very interest- ing because it contains much more information than the CMB. LSS is a three dimensional dataset, and if we can observe it from the scale of about 1h−1Mpc out to the Hubble scale, we have in principle ' (3000)3 = 2.7 × 1010 independent modes at our disposition. This has to be compared with the CMB with about 107 modes (including E-polarization). Furthermore, non-linearities in the CMB are very small and CMB fluctuations are very Gaussian. Hence the power spectrum captures nearly all the information and (reduced) higher order correlations are very small. For large scale structure, at least on intermediate to small scales, non-linearities become important. The analysis of the highly non-linear regime can probably only be done via full numerical simulations, but in the intermediate, weakly non-linear regime higher order perturbation theory can be applied. As soon as non-linearities
1See also: http://www.darkenergysurvey.org/science/ and http://www.euclid-ec.org
– 1 – are present, higher order correlations develop i.e. non-trivial three-point and reduced four- point functions are generated. These functions contain additional interesting information on gravitational clustering which can be very important for the Dark Energy problem or the contribution from massive neutrinos. The present work inscribes in this direction. We first compute a basic LSS observable, namely the galaxy number counts to second order in perturbation theory and then determine the induced bispectrum. Our approach is novel in several aspects. We describe the number counts as a function of direction and observed redshift. This has the big advantage that it is directly related to observables and that it is model indepen- dent contrary to the galaxy correlation function in real space or its power spectrum. The latter need the knowledge of a background Universe to convert measured redshifts and angles into distances and are therefore not well adapted, e.g., to estimate cosmological parameters. This method was introduced in [7] and implemented and tested numerically in [8,9]. Galaxy number counts have the advantage of being directly measured and their correlation func- 0 0 tion ξ(θ, z, z ) or the power spectrum in the form of C`(z, z ) can be directly compared to observations. The first order fluctuations of the number counts have been calculated in [7, 10–12]. In this work we go to second order. Several second order calculations have already been per- formed. Expressions for the redshift perturbations and cosmological distances were obtained in [13–15] and applications are discussed in [16–18]. Similar expressions are also obtained in [19, 20]. Results for the Newtonian density fluctuation were obtained in [21] and second order perturbation theory of lensing is discussed in [22]. Also expressions for the second order number counts have just appeared in [23–25]. However, since these calculations are very in- volved, it is useful to have results from independent groups, obtained in different ways, which can then be compared. Here we employ a novel method to calculate the second order num- ber counts. We use the geodesic light-cone gauge introduced in [26]. Actually, if we would know the density fluctuation exactly, our expression for the number counts would be fully non-perturbative. Then we perform a second order coordinate transformation to translate the result to the familiar Bardeen potentials, i.e. the metric perturbations in longitudinal or Poisson gauge. We also allow for anisotropic stresses. We however neglect second order vector and tensor contributions which can be easily obtained from the first order results by inserting the second order vector and tensor perturbations induced by the Bardeen poten- tials, which have been calculated in the past [27–29]. We write the final result fully in terms of the Bardeen potentials, the peculiar velocities and the density fluctuations. Of course, the final expressions remain complicated, but they can be readily used and compared to other results in the literature. Finally, we use our expressions to compute some of the most important new terms in the bispectrum, namely the contribution from redshift space distortion and the one from lensing, and we evaluate them numerically and compare them with the usual term from weakly non-linear Newtonian gravitational clustering. The remainder of this paper is organised as follows: in Section2 we lay out the basic quantities which we want to compute. In Section3 we describe the geodesic light-cone gauge and its transformation to the Poisson or longitudinal gauge. In Section4 we compute the second order number counts first in geodesic light-cone gauge and then we transform it to Poisson gauge. We also present a simplified formula containing only the dominant terms. In Section5 we derive and discuss some of the most relevant terms in the bispectrum for which we also plot numerical results using classgal [8]. In Section6 we conclude. Some useful
– 2 – relations are collected in AppendixA.
2 Galaxy Number Counts
What we really observe in galaxy surveys is the number N of galaxies in a redshift bin dz and solid angle dΩ as a function of the observed redshift z and of the observation direction defined by the unit vector −n (in this work n denotes the direction of propagation of the photon). Directly from these observable quantities we define the fluctuation of galaxy number counts as N (n, z) − hNi (z) ∆ (n, z) ≡ , (2.1) hNi (z) where h·i denotes the average over directions. If we neglect bias, the observed number of galaxies can be expressed in terms of the density ρ (n, z) and the volume V (n, z) at fixed redshift according to N (n, z) = ρ (n, z) V (n, z) . (2.2) To determine the galaxy number counts ∆ (n, z) up to second order in perturbation theory we need to expand the redshift density perturbation and the volume to second order. The volume and the density have to be defined in function of the observed redshift and direction of observation, then their expansion to second order can be given as ! δV (1) δV (2) V (n, z) = V¯ (z) + δV (1) (n, z) + δV (2) (n, z) = V¯ (z) 1 + + (2.3) V¯ V¯
(0) (1) (2) ρ (n, z) =ρ ¯ ηs 1 + δ + δ . (2.4) Here, z denotes the observed redshift, which is also perturbed, and we define the time coor- dinate of a fiducial background model of a source by a(η ) 1 + z = o , (2.5) s (0) a(ηs ) where ηo denotes the conformal time of the observer. We also note that the background (0) (0) conformal distance is rs = ηo − ηs . Hence to some extent, the perturbations of the volume and the density above come from the perturbations of the redshift. In the next section we shall present the second order perturbations of the redshift, obtained in [15], and Section3 is devoted to the calculation of the volume and the density fluctuations. Once these perturbations are calculated, inserted in the definition (2.1), they yield to second order in perturbation theory the following result " # δV (1) δV (1) δV (2) δV (1) δV (2) ∆ (n, z) = δ(1) + + δ(1) + δ(2) + − hδ(1) i − hδ(2)i − h i (2.6) V¯ V¯ V¯ V¯ V¯ Here we use the fact that the angular mean of first order perturbations vanishes. We then identify the first order and the second order contributions to galaxy number counts, respec- tively, δV (1) ∆(1) = δ(1) + , (2.7) V¯ δV (1) δV (2) δV (1) δV (2) ∆(2) = δ(1) + δ(2) + − hδ(1) i − hδ(2)i − h i . (2.8) V¯ V¯ V¯ V¯
– 3 – The first order contribution (2.7) has been already computed in [7, 10–12]. Here we are mainly interested in computing the second order term (2.8). For this purpose we adopt the so-called geodesic light-cone (GLC) coordinates defined in [26]. In terms of GLC coordinates we determine an exact solution which we then expand in a more conventional gauge, namely the Poisson gauge, to second order.
3 From the Geodesic Light-Cone to the Poisson gauge
Let us first introduce the geodesic light-cone (GLC) coordinates [26]. They consist of a timelike coordinate τ (which can be identified with the proper time of synchronous gauge and, therefore, describes a geodesic observer that is static in this gauge [30]), of a null coordinate w and of two angular coordinates θ˜a (a = 1, 2). We denote them with a tilde to distinguish these exact screen space angles from the background angles in the Friedmann metric which we introduce below. The line-element of the GLC metric takes the form
2 2 2 ˜a a ˜b b ds = Υ dw − 2Υdwdτ + γab(dθ − U dw)(dθ − U dw) , a, b = 1, 2 . (3.1)
It depends on six arbitrary functions (Υ, the 2-component vector U a and the symmetric 2×2 matrix γab). In matrix form:
0 −Υ ~0 −1 −Υ−1 −U b/Υ 2 2 µν −1 ~ gµν = −ΥΥ + U −Ub , g = −Υ 0 0 , (3.2) ~ T T a T ~ T ab 0 −Ua γab −(U ) /Υ 0 γ
ab 2 a b where γab and its inverse γ lower and raise the two-dimensional indices and U = γabU U . The 2-dimensional space described by the coordinates (θ˜1, θ˜2) is the exact screen space with metric γab. In an unperturbed Friedmann background, our coordinates correspond to cosmic time τ = t, background past light-cone w = η +r, where η is conformal time so that adη = dt, and the usual polar angles θ = θ˜1 and φ = θ˜2. The background metric components are Υ = a, ˜a ˜b 2 2 2 2 2 Ua = 0 and γabdθ dθ = a r (dθ + sin θdφ ). Here a denotes the cosmic scale factor. Like in synchronous gauge, also the GLC gauge has some residual gauge freedom (for µ details see [14]). The condition w = constant defines a null hypersurface (∂µw∂ w = 0), corresponding to the past light-cone of a given observer, hereafter chosen to be the µ geodesic observer which moves with 4-velocity uµ = −∂µτ. The geodesic equation for u ν is (∂ τ) ∇ν (∂µτ) = 0. Let us also recall that, in GLC gauge, the null geodesics connect- µ µν ing sources and observer are characterised simply by the tangent vector k = −ωg ∂νw = µw −1 µ −ωg = ωΥ δτ (where ω is an arbitrary normalization constant), meaning that photons travel at constant values of w and θ˜a. This renders the calculation of the redshift particularly simple in this gauge. We now give exact, non-perturbative, expressions for the redshift and the cosmological distances, like the luminosity and angular distances, in GLC gauge. As in the previous section, subscripts “o” and “s”, respectively, denote quantities evaluated at the observer and the source space-time positions. We consider a light ray emitted by a geodesic source with µ µν 4-velocity us = −(g ∂ντ)s lying on the past light-cone of the geodesic observer defined by w = wo and on the spatial hypersurface τ = τs. The light ray will be received by the geodesic
– 4 – observer at τ = τo > τs. The exact non-perturbative expression of the redshift zs associated with this light ray is then simply given by [26]
µ µ a (k uµ)s (∂ w∂µτ)s Υ(wo, τo, θe ) (1 + zs) = = = . (3.3) µ µ a (k uµ)o (∂ w∂µτ)o Υ(wo, τs, θe ) On the other hand, in [14] an exact expression for the so-called Jacobi Map [31] is derived in GLC gauge and the following non-perturbative solution for the luminosity distance dL and the angular diameter distance dA is obtained: √ 4 γ d2 = (1 + z )4d2 = (1 + z )4 s , (3.4) L s A s −1 ab 3/2 det uτ ∂τ γ γ o where γ denotes the determinant of the 2-dimensional matrix γab. We finally want to express our results in Poisson gauge. Neglecting vector and tensor contributions, the Poisson gauge (PG) metric [32, 33] (sometimes also denoted ’Newtonian gauge’ or ’longitudinal gauge’) takes the form
2 2 2 i j dsPG = a (η) −(1 + 2Φ)dη + (1 − 2Ψ)δijdx dx = a2(η) −(1 + 2Φ)dη2 + (1 − 2Ψ)(dr2 + r2d2Ω) (3.5) where the (generalized) Bardeen potentials Φ and Ψ are defined, up to second order, as follows: 1 1 Φ ≡ φ + φ(2) , Ψ ≡ ψ + ψ(2) , (3.6) 2 2 and we make no assumption on the anisotropic stress, so that Ψ and Φ may differ also at first order. To express the redshift and the luminosity distance given in Eqs.(3.3) and (3.4) in terms of standard PG variables we have to transform the GLC gauge quantities to quantities in PG. Following [15], and employing suitable boundary conditions: i) the transformation is non singular around the observer position at r = 0, and ii) the two-dimensional spatial sections r = const are locally parametrized at the observer position (for any time) by standard spherical coordinates (θ, φ) ≡ (θ1, θ2), the coordinate transformation to second order is given by [15]
τ = τ (0) + τ (1) + τ (2) Z η (0) 0 0 (1) a with τ = dη a(η ) = t − tin , τ = a(η)P (η, r, θ ) , ηin Z η 0 (2) 0 a(η ) h (2) 2 2 ab i 0 a τ = dη φ − φ + (∂rP ) + γ0 ∂aP ∂bP (η , r, θ ) , (3.7) ηin 2 w = w(0) + w(1) + w(2) (0) (1) a with w = η+ ≡ η + r , w = Q(η+, η−, θ ) , Z η− (2) 1 h (2) (2) 2 2 ab i a w = dx ψ + φ + 2(ψ − φ ) + 2(ψ + φ)∂+Q + γ0 ∂aQ ∂bQ (η+, x, θ ) , 4 ηo (3.8)
– 5 – θ˜a = θ˜a(0) + θ˜a(1) + θ˜a(2) Z η− ˜a(0) a ˜a(1) 1 h ab i a with θ = θ , θ = dx γ0 ∂bQ (η+, x, θ ) 2 ηo Z η− ˜a(2) h1 ac (2) ac (1) 1 dc (1) ˜a(1) 1 ˜a(1) θ = dx γ0 ∂cw + ψγ0 ∂cw + γ0 ∂cw ∂dθ + (ψ + φ)∂+θ ηo 2 2 2 a(1) (1) a(1)i a +(φ − ψ)∂−θ˜ − ∂+w ∂−θ˜ (η+, x, θ ) , (3.9)
ab −2 −2 −2 where (γ0 ) = diag(r , r sin θ), and ηin represents an early enough time when the perturbations (or better their integrands) were negligible. We have also introduced the zeroth-order light-cone variables η± = η ± r, with corresponding partial derivatives: ∂ 1 ∂η = ∂+ + ∂− , ∂r = ∂+ − ∂− , ∂± = = (∂η ± ∂r) . (3.10) ∂η± 2 Furthermore, we define
Z η 0 Z η− a 0 a(η ) 0 a a 1 a P (η, r, θ ) = dη φ(η , r, θ ) ,Q(η+, η−, θ ) = dx (ψ + φ)(η+, x, θ ) . ηin a(η) ηo 2 (3.11) Note that the integrals for τ go along the conformal time axis η at fixed r, while the integrals a for w and θ˜ are along the background light cone, η− at fixed η+. With this we can obtain the non-trivial entries of the GLC metric of Eq. (3.2) to second order in terms of the variables (η, r, θa):
1 1 2 1 Υ−1 = 1 + ∂ Q − ∂ P + (ψ − φ) + ∂ w(2) + ∂ τ (2) − φ(2) + 2φ2 − φ(φ + ψ) a(η) + r 2 η a − 2 1 3 −φ∂ Q + ∂ P φ − ψ − ∂ P ∂ Q − γab∂ P ∂ Q , (3.12) + r 2 2 r + 0 a b
1 1 1 U a = ∂ θ˜a(1) − γab∂ τ (1) + ∂ θ˜a(2) − γab∂ τ (2) − ∂ τ (1)∂ θ˜a(1) η a 0 b η a 0 b a r r 2 1 −φ∂ θ˜a(1) − ψγab∂ τ (1) − γcd∂ τ (1)∂ θ˜a(1) η a 0 b a 0 c d 1 1 + ∂ Q − ∂ P + (ψ − φ) −∂ θ˜a(1) + γab∂ τ (1) , (3.13) + r 2 η a 0 b ab −2 h ab (ac ˜b)(1) ab (2) 2 ˜a(1) ˜b(1) γ = a γ0 (1 + 2ψ) + 2γ0 ∂cθ + γ0 ψ + 4ψ − ∂ηθ ∂ηθ ˜a(1) ˜b(1) (ac ˜b)(1) cd ˜a(1) ˜b(1) (ac ˜b)(2)i +∂rθ ∂rθ + 4ψγ0 ∂cθ + γ0 ∂cθ ∂dθ + 2γ0 ∂cθ . (3.14)
Here X(a···Y ···b)··· denotes symmetrization and a and b. As a first example of the use of the coordinate transformation given above let us present the perturbed redshift up to second order in PG in the presence of anisotropic stress and in function of the observer’s angular coordinates [15]. This is easily done starting from the non-perturbative solution (3.3) and using Eq. (3.12) together with the fact that θ˜a actually are the standard angular coordinate at the observer position and are constant along the line of sight.
– 6 – Therefore, any perturbative quantity can be written as a function of the observer’s ˜a a angular coordinates, i.e. the direction of observation, by Taylor-expanding it around θ = θo . To this purpose, and considering general quantity with an angle-independent background value, it is enough to invert Eq. (3.9) to first order. Namely, we need the expansion
Z ηo Z ηo a a(0) a(1) a 0 ab 00 00 00 a θ = θ + θ = θo − dη γ0 ∂b dη (ψ + φ)(η , ηo − η , θo ) . (3.15) η η0 To obtain familiar expressions in terms also of peculiar velocities, we recall that τ plays the role of a gauge-invariant velocity potential (see [14]), so that the spatial components of the perturbed velocity vµ of the (geodesic) observer in polar coordinates in PG are given by
(2) (2) (vi) = (vr + vr , v⊥a + v⊥a) (1) (2) (2) (1) (2) (2) with vr = −∂rτ , vr = −∂rτ , v⊥a = −∂aτ , v⊥a = −∂aτ , (3.16) where τ (1) and τ (2) are the first- and second-order part of the coordinate transformation τ = τ(η, r, θa) given in Eq. (3.7). The unit vector nµ along the direction connecting the source to the observer can be then expanded, in polar coordinates and to first order (which is enough for our purpose),
1 nµ = 0, − (1 + ψ), 0, 0 , n = (0, −a(1 − ψ), 0, 0) . (3.17) a µ
Taking its scalar product with the spatial component of the perturbed velocity we obtain
Z η 0 (2) 1 0 a(η ) h (2) 2 2 ab i 0 a v·n = v||+v|| =∂rP +ψ∂rP + dη ∂r φ − φ + (∂rP ) + γ0 ∂aP ∂bP (η , r, θ ) . 2 ηin a(η) (3.18) For the velocity components which are perpendicular to the radial direction we have
a ab v⊥av⊥ = γ0 ∂aP ∂bP. (3.19)
Let us also define the following useful variables ψ + φ ψ − φ ψI = , ψA = , (3.20) 2 2 which define the isotropic and anisotropic part of the Bardeen potential; ψI is also the potential of the scalar part of the Weyl tensor [34]. In what follows we will use these variables to express our perturbed quantities. Furthermore, to simplify the notation, hereafter we only indicate the conformal time η as a the integration variable along the line of sight, instead of all the arguments like (η, ηo −η, θo ), 1 2 of course the result also depends on the direction fixed by (θo, θo). Finally, following [15], we obtain the observed redshift to second order and in function of the direction of observation
a(ηo) h (1) (2)i 1 + zs = 1 + δz + δz (3.21) a(ηs) with Z ηo (1) I A 0 I 0 δz = −v||s − ψs + ψs − 2 dη ∂η0 ψ η (3.22) ηs
– 7 – Z η0 (2) (2) 1 (2) 1 0 h (2) (2)i 0 1 2 1 I 2 δz = −v||s − φs − dη ∂η0 φ + ψ η + v||s + ψs 2 2 ηs 2 2 Z ηo Z ηo I I 0 I 0 1 a a 0 I 0 + −v||s − ψs −ψs − 2 dη ∂η0 ψ η + v⊥sv⊥a s − 2a v⊥s∂a dη ψ η ηs 2 ηs Z η0 Z ηo 0 I 0 I 0 I 0 00 I 00 +4 dη ψ η ∂η0 ψ η + ∂η0 ψ η dη ∂η00 ψ η 0 ηs η Z ηo Z ηo Z ηo I 0 00 2 I 00 ab 00 I 00 00 I 00 +ψ η dη ∂η00 ψ η − γ0 ∂a dη ψ η ∂b dη ∂η00 ψ η η0 η0 η0 Z ηo Z ηo I 0 ab 00 I 00 +2∂a v||s + ψs dη γ0 ∂b dη ψ η 0 ηs η Z ηo Z ηo Z ηo 0 I 0 00 ab 000 I 000 +4 dη ∂a ∂η0 ψ η dη γ0 ∂b dη ψ η 00 ηs ηs η Z ηo Z η0 A 3 A 2 I A 0 I A 0 A 0 I 0 −ψs v||s + (ψs ) − 3ψs ψs − 4 dη ∂η0 (ψ ψ ) η − 2ψs dη ∂η0 ψ η 2 ηs ηs Z ηo Z ηo A 0 ab 00 I 00 −2∂aψs dη γ0 ∂b dη ψ η . (3.23) 0 ηs η Let us stress that the above integrals are along the line of sight, f(η0) = f(η0, xi(η0)) and partial derivatives w.r.t. η0 are only derivatives with respect to the first dependence. Above and hereafter we neglect the terms from the gravitational potential and the ve- locity at the observer position. They cannot be calculated within perturbation theory. Note, however, that in addition to the monopole and the dipole terms, they now also contribute 2 2 to the quadrupole e.g. via terms of the form vko = (n · vo) . However, differences like C2(z2)−C2(z1) are meaningful and can be obtained from our formula, see [9] for a discussion of C0(z2) − C0(z1) for the power spectrum from linear perturbation theory. Furthermore, terms of the form XoYs which are a product of a term Xo at the observer position and a term Ys at the source position can contribute to all moments in the 4th order contributions to the power spectrum. Assuming that the correlation of Xo and Ys can be neglected, for 2 two sources at positions s1 and s2 they give correlators of the form hXo ihYs1 Ys2 i. We have checked that Xo is always either the velocity or the gravitational potential at the observer position. The first can be removed by going to the CMB frame while the second is of the order 10−10 times the power spectrum and hence negligible. Similar considerations apply for the (reduced) bispectrum.
4 Evaluation of the number counts to second order
4.1 Redshift density perturbation To second order the density perturbation at fixed PG coordinates can be written as a (1) a (2) a (1) (2) ρ (ηs, rs, θs ) =ρ ¯(ηs) + δρ (ηs, rs, θs ) + δρ (ηs, rs, θs ) =ρ ¯(ηs) 1 + δρ + δρ . (4.1) We want to evaluate this at fixed observed redshift and with respect to the observer’s an- a gular coordinates θo . Therefore, we define a fiducial background model with coordinates (0) (0) a (ηs , rs , θo ) for which the observed redshift and the past light-cone of our observer are given by a(η ) 1 + z = o , w = w = η = η(0) + r(0) . (4.2) s (0) o o s s a(ηs )
– 8 – We then expand conformal time and radial PG coordinates around the coordinates of the (0) (1) (2) (0) (1) (2) fiducial model as ηs = ηs + ηs + ηs and rs = rs + rs + rs , and obtain the terms of these expansions by perturbatively solving the following system of equations:
a(η ) a(η ) h i 1 + z = o = o 1 + δz(1) + δz(2) (4.3) s (0) a(ηs ) a(ηs) (0) (1) (2) w = ηo = w + w + w (4.4) where we use Eqs. (3.8), (3.22) and (3.23) for the perturbations appearing in Eqs. (4.3) and (4.4). We then obtain the perturbed conformal time and radial coordinate needed so that our measured redshift corresponds to the background redshift of our fiducial model,
(1) (1) δz ηs = (4.5) Hs 0 (2) 1 (2) (1→2) 1 Hs (1) 2 ηs = δz + δz − 1 + 2 (δz ) . (4.6) Hs 2 Hs and Z ηo (0) (0) (1) (1) 0 I 0 r = ηo − η , r = −η + 2 dη ψ η (4.7) s s s s (0) ηs (2) (2) (2) (1→2) rs = −ηs − ws − ws (4.8) 0 (0) where Hs = (a /a)(ηs ) is the comoving Hubble parameter of the fiducial background model, (0) evaluated at the background conformal time ηs , and all the quantities above are evaluated (0) (0) a (1→2) (1→2) at (ηs , rs , θo ) . The quantities δz and ws stand for the second order contribution (1) (1) coming from Taylor expanding δz and ws around the background fiducial model. They are given by
(1→2) (1) I A I I A δz = ηs ∂ηψs + ∂ηψs + Hsv||s + ∂rv||s + −2∂ηψs − ∂rψs + ∂rψs Z ηo Z ηo 0 2 I 0 0 I 0 −∂rv||s − 2 dη ∂η0 ψ (η ) 2 dη ψ (η ) (4.9) ηs ηs
Z ηo Z ηo Z ηo (1→2) I 0 I 0 0 I 0 I (1) 0 I 0 ws = −ψs − 2 dη ∂η0 ψ η 2 dη ψ (η ) + ψs 2ηs − 2 dη ψ (η ) ηs ηs ηs Z ηo Z ηo Z ηo 0 I 0 00 ab 000 I 000 +4 dη ∂aψ η dη γ0 ∂b dη ψ η , (4.10) 00 ηs ηs η
(0) (0) where to second order we are allowed to drop the superscript from ηs . Let us now expand the second order result (4.1) around the fiducial model, to obtain the second order density fluctuation at fixed observed redshift and direction of observation,
1 2 ρ (n, z ) =ρ ¯ + ∂ ρ¯ η(1) + η(2) + ∂2ρ¯ η(1) s η s s 2 η s (1) (1) (1) (1) (1) a(1) (1) (2) +δρ + ηs ∂ηδρ + rs ∂rδρ + θs ∂aδρ + δρ (4.11) =ρ ¯ 1 + δ(1) + δ(2) . (4.12)
– 9 – (0) (0) a Here all the quantities have to be evaluated at (ηs , rs , θo ). Finally, using Eqs. (4.5)-(4.7) dρ¯ ρ¯ and the relation dz = 3 1+z , we obtain the following perturbative result: at first order 1 δ(1) = η(1)∂ ρ¯ + δρ(1) ρ¯ s η (1) (1) = −3 δz + δρ (4.13) and at second order 1 1 2 δ(2) = η(2)∂ ρ¯ + η(1) ∂2ρ¯ + η(1)∂ δρ(1) + r(1)∂ δρ(1) + θa(1)∂ δρ(1) + δρ(2) ρ¯ s η 2 s η s η s r s a 2 (2) (1) 1 I A 1 (1) = −3δz + 6 δz − 3 ∂ηψs + ∂ηψs + v||s + ∂rv||s δz Hs Hs Z ηo Z ηo I I A 0 2 I 0 0 I 0 +6 2∂ηψs + ∂rψs − ∂rψs + ∂rv||s + 2 dη ∂η0 ψ η dη ψ η ηs ηs Z ηo 1 1 (1) (1) (1) (1) 0 I 0 + ∂η ρ¯ δρ − ∂rδρ δz + 2∂rδρ dη ψ η Hs ρ¯ ηs Z ηo Z ηo (1) 0 ab 00 I 00 (2) −2∂aδρ dη γ0 ∂b dη ψ η + δρ , (4.14) 0 ηs η where Eq. (4.13) agrees with the first order result of [7].
4.2 Volume perturbation In order to compute the contribution to the galaxy number counts due to the volume per- turbation, we start considering the 3-dimensional volume element seen by a source with 4-velocity uµ √ µ ν α β dV = −gµναβu dx dx dx . (4.15)
We then express this in terms of the observed quantities (z, θo, φo)
ν α β √ µ ∂x ∂x ∂x ∂ (θs, φs) dV = −gµναβu dzdθodφo ≡ v (z, θo, φo) dzdθodφo . (4.16) ∂z ∂θs ∂φs ∂ (θoφo) The volume perturbation is determined by δV v − v¯ δv = = . (4.17) V¯ v¯ v¯ Let us now express the above quantities in GLC coordinates, xµ = τ, w, θ˜a with a = 1, 2, to show how the volume perturbation can be easily given fully non-perturbatively ˜a ˜a a in GLC gauge. Considering that θ is constant along the geodesic of the light-cone θ = θo , defined through w = const, we have √ ∂τ dV = − −guw dzdθ dφ . (4.18) ∂z o o √ Using also uw = −gwτ = Υ−1 and −g = p|γ|Υ, we obtain
2 p dτ p dτ p Υs dV = |γ| − dzdθodφo , or v = |γ| − = |γ| , (4.19) dz dz Υo∂τ Υs
– 10 – where we made use of Eq. (3.3) to calculate dz/dτ. Note that Eq. (4.19) is again a non- perturbative expression for the volume element at the source in terms of the observed redshift and the observation angles in GLC gauge. To express the above quantity in Poisson gauge up to second order we have to calculate p dτ p |γ| and dz in this gauge. We first note that |γ| can be obtained, inverting Eq. (3.4), from the known result for the area distance given in [13–15]. In particular, neglecting the terms evaluated at the observer position, we have
2 (0) 2 2 p aors ¯(1) a ¯(1) a ¯(2) a |γ| = 2 sin θo 1 + 2δS (zs, θo ) + δS (zs, θo ) + 2δS (zs, θo ) , (4.20) (1 + zs) where we have set d (z , θa) d (z , θa) L s o = L s o = 1 + δ¯(1)(z , θa) + δ¯(2)(z , θa) , (4.21) (0) F LRW S s o S s o (1 + zs)aors dL (zs)
¯(1) a ¯(2) a and δS (zs, θo ) and δS (zs, θo ), the first and second order scalar perturbations of the lumi- nosity distance, are given in Sec. III of [15]. Finally, from the results of Sec. III of [15] and after some algebraic manipulations, we obtain the following useful perturbative result to first order (1) p ! |γ| Z ηo 1 (1) I A 4 0 I 0 = 2 1 − δz − 2 ψs + ψs + dη ψ η (0) (0) (0) (0) p|γ| Hsrs rs ηs
Z ηo 0 (0) 2 0 η − ηs I 0 − dη ∆2ψ η . (4.22) (0) (0) η − η0 rs ηs o
Here ∆2 denotes the angular Laplacian. To second order we have a more involved result (2) p " # |γ| 1 1 H0 1 2 2 = 2 1 − δz(2) + − 1 − s + 1 − δz(1) +Ξ +Ξ (0) H r H r H2 H r IS AS p|γ| s s s s s s s (4.23) where
Z ηo 0 Z ηo (2) 1 1 0 η − ηs h (2) (2)i 0 1 0 h (2) (2)i 0 ΞIS = −ψs − dη 0 ∆2 ψ + φ η + dη ψ + φ η 2 rs ηs ηo − η rs ηs 1 1 2 +2 1 − − v||s∂rv||s − v||s Hsrs Hs Z ηo Z ηo Z ηo 0 I 0 I 0 2 0 I 0 2 0 η − ηs I 0 0 + ψs − 2 dη ∂η ψ η − dη ψ η + dη 0 ∆2ψ η v||s ηs rs ηs rs ηs ηo − η Z ηo Z ηo I 0 I 0 0 I 0 1 + −ψs − 2 dη ∂η0 ψ η − 2Hs dη ψ η ∂rv||s ηs ηs Hs Z ηo Z ηo I I 0 2 I 0 0 I 0 + ∂rψs + 2∂ηψs + 2 dη ∂η0 ψ η −2 dη ψ η ηs ηs Z ηo Z ηo 0 Z ηo I 0 I 0 2 0 η − ηs I 0 2 0 I 0 0 − −ψs − 2 dη ∂η ψ η dη 0 ∆2ψ η − dη ψ η ηs rs ηs ηo − η rs ηs
– 11 – Z ηo Z ηo 2 0 I 0 2 I 1 I 1 I 1 0 I 0 + − dη ψ η + ψs + 2 ∂ηψs − ∂rψs − 2 dη ∆2ψ η 2v||s rs ηs Hsrs Hsrs Hs Hsrs ηs Z ηo Z ηo 2 Z ηo 2 ab 0 I 0 4 0 I 0 4 I 0 I 0 − ∂av||sγ0s dη ∂bψ (η ) + 2 dη ψ η − ψs dη ψ η Hs ηs rs ηs rs ηs Z ηo Z ηo I 0 I 0 I 2 2 I 1 I I 0 I 0 +8ψs dη ∂η0 ψ η + 4(ψs ) + ∂rψs − ∂ηψs −ψs − 2 dη ∂η0 ψ η ηs Hs Hsrs ηs Z ηo Z ηo Z ηo I 0 I 0 4 0 I 0 I 0 00 I 00 −4∂rψs dη ψ η + dη ψ η −ψ η − 2 dη ∂η00 ψ η 0 ηs rs ηs η Z ηo Z ηo ab 00 I 00 00 I 00 +γ0 ∂a dη ψ η ∂b dη ψ η η0 η0 Z ηo Z ηo 0 Z ηo 0 2 I 2 0 I 0 4 0 η − ηs I 0 1 0 η − ηs I 0 + ψs − dη ψ η dη 0 ∆2ψ η + 2 dη 0 ∆2ψ η rs ηs rs ηs ηo − η rs ηs ηo − η Z ηo Z ηo Z ηo 1 I 0 I 0 1 0 I 0 2 0 I 0 + −ψs − 2 dη ∂η0 ψ η − dη ψ η dη ∆2ψ η Hsrs ηs rs ηs rs ηs Z ηo Z ηo Z ηo Z ηo 0 1 I 0 0 1 00 I 00 0 I 0 −2 dη ψ (η ) + 2 dη dη ∂η00 ψ (η ) dη ∆2ψ (η ) 0 2 0 2 0 ηs (ηo − η ) ηs (ηo − η ) η ηs Z ηo Z ηo Z ηo Z ηo Z ηo I 0 ab 00 I 00 8 0 I 0 00 ab 000 I 000 +4∂aψs dη γ0 ∂b dη ψ η − dη ∂aψ η dη γ0 ∂b dη ψ η 0 00 ηs η rs ηs ηs η Z ηo Z ηo Z ηo Z ηo 0 I 0 0 1 ab 00 I 00 0 ab I 0 +2∂a dη ψ η 4 dη γ0 dη ∂bψ (η ) − 3 dη γ0 ∂bψ (η ) 0 0 ηs ηs (ηo − η ) η ηs Z ηo Z ηo 0 ab 00 I 00 −6 dη γ0 dη ∂b∂η00 ψ (η ) 0 ηs η Z ηo Z ηo Z ηo Z ηo 0 bc 00 I 00 0 ad 00 I 00 +2∂a dη γ0 ∂c dη ψ η ∂b dη γ0 ∂d dη ψ η 0 0 ηs η ηs η Z ηo Z ηo Z ηo 0 I 0 0 1 00 I 00 −4 dη ψ η dη − dη ∆2ψ η 0 3 0 ηs ηs (ηo − η ) η Z ηo Z ηo 0 2 1 1 I 0 00 I 00 2 1 0 η −ηs I 0 00 + ∆2ψ η + dη ∂η ∆2ψ η + 2 dη ∂θo ψ η 0 2 0 0 (ηo −η ) 2 η (sin θo) rs ηs ηo −η Z ηo Z ηo Z ηo 1 I 0 I 0 0 I 0 ab 0 I 0 +2∂a −ψs − 2 dη∂η0 ψ η − dη ψ (η ) γ0s dη ∂bψ (η ) Hs ηs ηs ηs Z ηo 0 Z ηo Z ηo 4 0 η − ηs I 0 00 ab 000 I 000 + dη ∂b ∆2 ψ (η ) dη γ0 ∂a dη ψ (η ) 0 00 rs ηs ηo − η ηs η Z ηo 0 Z ηo 2 0 η − ηs I 0 I 0 00 I 00 − dη ∆2 ψ η −ψ η − 2 dη ∂η00 ψ η 0 0 rs ηs ηo − η η Z ηo Z ηo ab 00 I 00 00 I 00 +γ0 ∂a dη ψ η ∂b dη ψ η η0 η0 Z ηo Z ηo 00 0 0 I 0 1 00 η − η I 00 −2 dη −2ψ η dη ∆2∂η00 ψ η 0 0 00 ηs ηo − η η ηo − η Z ηo Z ηo 00 0 ab 00 I 00 1 00 η − η I 00 +2γ0 ∂b dη ψ η 0 dη 00 ∂a∆2ψ η η0 ηo − η η0 ηo − η
– 12 – Z ηo Z ηo I 0 00 I 00 1 00 I 00 00 − −2ψ η − 2 dη ∂η ψ η 0 2 dη ∆2ψ η η0 (ηo − η ) η0 Z ηo Z ηo I 0 00 ab 000 I 000 −2∂aψ η dη γ0 ∂b dη ∂η000 ψ η η0 η00 Z ηo Z ηo Z ηo db 00 I 00 00 ac 000 I 000 +2∂a γ0 ∂b dη ψ η dη ∂d γ0 ∂c dη ψ η η0 η0 η00 Z ηo Z ηo ab I 0 00 I 00 00 I 00 +2γ0 ∂a ψ η + dη ∂η00 ψ η ∂b dη ψ η (4.24) η0 η0 and
A Z ηo 1 A ψs A 2 0 I 0 ΞAS = 2 1 − 2ψs v||s + ∂rv||s − ψs − dη ψ η Hsrs Hs rs ηs Z ηo Z ηo 0 0 I 0 2 0 η − ηs I 0 A2 I A 0 −2 dη ∂η ψ η + dη 0 ∆2ψ η − 2 ψs + 2ψs ψs ηs rs ηs ηo − η Z ηo A 1 A A 2 4 A 0 I 0 +2v||sψs + ∂ηψs − ∂rψs v||s − ∂rψs dη ψ η Hsrs Hs Hsrs ηs Z ηo Z ηo 0 A I 2 0 I 0 2 0 η − ηs I 0 +2ψs −2ψs − dη ψ η + dη 0 ∆2ψ η rs ηs rs ηs ηo − η Z ηo 1 A A 2 I A 0 I 0 − ∂ηψs − ∂rψs −(ψs − ψs ) − 2 dη∂η0 ψ η Hsrs Hs ηs Z ηo 1 I I 2 A 2 A 0 I 0 − ∂ηψs − ∂rψs ψs + 2 ψs dη ∆2ψ η Hsrs Hs Hsrs ηs Z ηo Z ηo Z ηo A 0 ab 00 I 00 2 A ab 0 I 0 +4∂aψs dη γ0 ∂b dη ψ η + ∂aψs γ0s ∂b dη ψ η 0 ηs η Hs ηs Z ηo Z ηo 0 8 0 I A 0 4 0 η − ηs I A 0 + dη ψ ψ η − dη 0 ∆2 ψ ψ η . (4.25) rs ηs rs ηs ηo − η
This second term ΞAS vanishes if no anisotropic stresses are present. Hereafter, we use the subscripts IS to denote terms in which the first order square contributions come entirely from ψI = (ψ + φ)/2 and AS for terms which come from the anisotropic stress, ∝ ψA = (ψ − φ)/2. dτ On the other hand, to compute dz to second order in perturbation theory in Poisson gauge we simply derive the coordinate transformation given in Eq.(3.7) with respect to the (0) (0) a redshift after expanding it around the fiducial background (ηs , rs , θo ). Proceeding step by step we first write (0) dτ dη(0) dτ 1 a ηs dτ = s = − . (4.26) (0) (0) dz dz dηs Hs a (ηo) dηs The second order expansion of τ around the fiducial model is given by
(η(1))2 τ(η , r , θa) = τ + (η(1) + η(2))∂ τ + r(1)∂ τ + θa(1)∂ τ + s ∂2τ , (4.27) s s s s s η s r s a 2 η
– 13 – (0) (0) a where now all the terms are evaluated at (ηs , rs , θo ). Then, using Eq.(3.7), and defining a (0) (1) (2) the perturbative expansion of τ(ηs, rs, θs ) = τs + τs + τs , we obtain the following results
(0) Z ηs (0) 0 0 τs = dη a(η ) (4.28) ηin
(0) Z ηs (1) (0) (1) 0 0 I A 0 (0) a τs = a(ηs )ηs + dη a(η )(ψ − ψ )(η , rs , θo ) (4.29) ηin
1 τ (2) = a(η(0))η(2) + H a(η(0))(η(1))2 + a(η(0))(ψI − ψA)η(1) + a(η(0))v r(1) s s s 2 s s s s s s s s ||s s (0) Z ηs 0 0 a(η ) h (2) I A 2 2 a i 0 (0) a + dη φ − (ψ − ψ ) + (v||) + v⊥av⊥ (η , rs , θo ) ηin 2 (0) Z ηs a(1) 0 0 I A 0 (0) a +θs ∂a dη a(η )(ψ − ψ )(η , rs , θo ) . (4.30) ηin
(0) We now derive the expressions above with respect to ηs . After a long but straightfor- ward calculation (see AppendixA for some useful relations), we obtain the following results
1 dτ (0) s = 1 (4.31) (0) (0) a(ηs ) dηs
1 dτ (1) H0 1 1 s = 1 − s δz(1) + ∂ ψI + ∂ v + ψI + ∂ ψA − ψA (4.32) (0) (0) 2 η s r ||s s η s s a(ηs ) dηs Hs Hs Hs
" # 1 dτ (2) H0 1 H0 3 H0 2 1 H00 2 s = 1 − s δz(2) + − s + s − s δz(1) (0) (0) 2 2 2 3 a(ηs ) dηs Hs 2 Hs 2 Hs 2 Hs 0 1 2 I 2 I I 1 I 2 3 Hs I + 2 ∂r ψs +∂η ψs −∂η∂rψs − ∂rψs + 1 − 2 ∂ηψs Hs Hs Hs 2 Hs 0 0 0 Hs 3 Hs 1 2 Hs I (1) − 2 v||s − 2 ∂rv||s − 2 ∂r v||s + 1 − 2 ψs δz Hs Hs Hs Hs Hs 0 1 2 A 2 A A 1 A 3 Hs A + 2 ∂η ψs −∂r ψs −∂η∂rψs + ∂rψs − 2 ∂ηψs Hs Hs Hs Hs 0 Hs A (1) − 1 − 2 ψs δz + ΠIS + ΠAS , (4.33) Hs
– 14 – where
0 Z ηo 0 1 (2) 1 (2) 1 (2) Hs 0 2 I 0 Hs I ΠIS = φs + ∂ηψs + ∂rv||s + 2 1 − 2 dη ∂η0 ψ η + 2 1 − 2 ∂ηψs 2 2Hs Hs Hs ηs Hs 0 0 Z ηo Hs I 1 I Hs 1 2 0 I 0 − 2 ∂rψs − ∂η∂rψs + 1 − 2 ∂rv||s − ∂r v||s −2 dη ψ η Hs Hs Hs Hs ηs Z ηo 1 1 I 0 I 0 1 I 2 2 I I + ∂rv||s + v||s + ∂ηψs −2 dη ∂η0 ψ η − (ψs ) + ψs ∂rv||s + ∂ηψs Hs Hs ηs 2 Hs 1 2 I 2 I 1 I 2 1 2 − (v||s) − ψs v||s + 2 ∂ηψs ∂rv||s + 2 ∂ηψs + 2 ∂rv||s 2 Hs Hs Hs Z ηo 1 4 a 2 a 4 a 0 I 0 1 a I − 1 − v⊥a sv⊥ s + av⊥ s∂av||s − av⊥ s∂a dη ψ η + av⊥ s∂aψs 2 Hsrs Hs Hsrs ηs Hs Z ηo 4 ab 0 I 0 1 1 I − γ0s ∂av||s∂b dη ψ η − v||s∂rv||s − v||s∂ηψs Hs ηs Hs Hs Z ηo Z ηo I 1 I 0 ab 00 I 00 +2∂a −ψs − ∂ηψs + ∂rv||s dη γ0 ∂b dη ψ η (4.34) 0 Hs ηs η and
0 Z ηo Z ηo 2 A Hs A 0 I 0 2 A 0 I 0 0 ΠAS = ∂η∂rψs − 2 2 ∂rψs dη ψ η − ∂ηψs dη ∂η ψ η Hs Hs ηs Hs ηs 2 A I 1 2 I A 1 A2 2 A I + ψs + ψs − v||s + 2 ∂ηψs + ∂rv||s ∂ηψs + 2 ∂ηψs + ψs ∂ηψs Hs Hs Hs Hs Hs Z ηo Z ηo A I 1 A2 A A 1 A 0 ab 00 I 00 +ψs ψs − ψs + ψs v||s + 2∂a ψs − ∂ηψs dηγ0 ∂b dη ψ η 0 2 Hs ηs η 1 a A + av⊥ s∂aψs . (4.35) Hs Finally, inserting the results given in Eqs. (4.22), (4.23), (4.32) and (4.33) in Eq. (4.19), we arrive at the following expressions for the first and second order volume perturbation ! (1) 0 Z ηo δV 2 Hs (1) I 4 0 I 0 = 3 − − δz − ψs + dη ψ η V¯ (0) H2 (0) (0) Hsrs s rs ηs
Z ηo 0 (0) 2 0 η − ηs I 0 1 I A 1 A − dη ∆2ψ η + ∂ηψs + ∂rv||s − 3ψs + ∂ηψs , (4.36) (0) (0) η − η0 H H rs ηs o s s
" δV (2) 2 H0 1 H0 3 H0 2 1 H00 1 H0 = 3 − − s δz(2) + − s + s − s − 1 − s ¯ 2 2 2 3 2 V Hsrs Hs 2 Hs 2 Hs 2 Hs Hsrs Hs # 2 0 2 1 1 Hs (1) + 1 − + 2 1 − 1 − 2 δz + ΛIS + ΛAS , (4.37) Hsrs Hsrs Hs
– 15 – where
Z ηo 0 (2) 1 (2) 1 (2) 1 (2) 1 1 0 η − ηs h (2) (2)i 0 ΛIS = −ψs + φs + ∂ηψs + ∂rv||s − dη 0 ∆2 ψ + φ η 2 2Hs Hs 2 rs ηs ηo − η Z ηo 1 0 h (2) (2)i 0 1 2 2 a + dη ψ + φ η + 2 1 − − v||s∂rv||s − (v||s) − v⊥a sv⊥ s rs ηs Hsrs Hs Z ηo Z ηo Z ηo 0 1 I 0 I 0 2 0 I 0 2 0 η − ηs I 0 0 + − ∂ηψs − 2 dη ∂η ψ η − dη ψ η + dη 0 ∆2ψ η v||s Hs ηs rs ηs rs ηs ηo − η Z ηo Z ηo I 0 I 0 0 I 0 1 + −2ψs − 4 dη ∂η0 ψ η − 2Hs dη ψ η ∂rv||s ηs ηs Hs Z ηo Z ηo Z ηo a 0 I 0 I I 0 2 I 0 0 I 0 +av⊥ s∂a dη ψ η + ∂rψs + 2∂ηψs + 2 dη ∂η0 ψ η −2 dη ψ η ηs ηs ηs Z ηo Z ηo 0 I 0 I 0 2 0 η − ηs I 0 0 − −ψs − 2 dη ∂η ψ η dη 0 ∆2ψ η ηs rs ηs ηo − η Z ηo 0 2 0 I 0 1 I I 3 a 2 a 1 Hs 2 − dη ψ η − ∂ηψs − ψs + v⊥a sv⊥ s + av⊥ s∂av||s + − + 2 (v||s) rs ηs Hs 2 Hs 2 Hs 0 Hs 1 1 h 2 2i 1 2 I 2 I I + −1 + 3 2 v||s∂rv||s + 2 v||s∂r v||s + ∂rv||s + − 2 ∂r ψs +∂η ψs −∂η∂rψs Hs Hs Hs Hs 0 0 Z ηo 1 I 3 Hs 2 1 I 4 Hs 0 I 0 − ∂rψs − 1 − 2 − ∂ηψs − 2 − 2 dη ψ η Hs Hs Hs 3 Hsrs rs Hs ηs 0 Z ηo 0 0 Z ηo Hs 2 0 η − ηs I 0 Hs 0 I 0 0 + 1 − 2 dη 0 ∆2ψ η − 2 1 − 2 dη ∂η ψ η Hs rs ηs ηo − η Hs ηs Z ηo 0 4 I 2 0 I 0 Hs 2 + ψs − 2 dη ∆2ψ η v||s + −2 1 − 2 − ∂rv||s Hsrs Hsrs ηs Hs Hsrs Z ηo 0 Z ηo 2 2 0 I 0 2 Hs 2 2 0 I 0 0 + ∂r v||s dη ψ η + − 1 − 3 2 ∂rv||s + 2 ∂r v||s dη ∂η ψ η Hs ηs Hs Hs Hs ηs Z ηo 0 0 2 1 0 η − ηs I 0 2 I 1 1 2 Hs I − ∂rv||s dη 0 ∆2ψ η + 2 ∂ηψs ∂rv||s + ∂r v||s + 3 2 ∂rv||s ψs Hs rs ηs ηo − η Hs Hs Hs Hs Z ηo Z ηo 2 a 0 I 0 1 a I 6 ab 0 I 0 − av⊥ s∂a dη ψ η + av⊥ s∂aψs − γ0s ∂av||s∂b dη ψ η Hsrs ηs Hs Hs ηs Z ηo 2 0 0 Z ηo 4 0 I 0 Hs I Hs 0 I 0 0 + 2 dη ψ η + 2 1 − 2 ψs + 4 1 − 2 dη ∂η ψ η rs ηs Hs Hs ηs 0 Z ηo 0 0 2 I 1 Hs 0 2 I 0 Hs I Hs I − ∂ηψs + 2 1 − 2 dη ∂η0 ψ η + 2 1 − 2 ∂ηψs + 2 − 2 ∂rψs Hs rs Hs ηs Hs Hs Z ηo 0 0 1 I 0 I 0 3 Hs 2 1 I 1 I Hs I − ∂η∂rψs −2 dη ψ η + 1− 2 − ∂ηψs + ∂rψs − 5− 2 ψs Hs ηs Hs Hs 3 Hsrs Hs Hs Z ηo 0 1 2 I 2 I I 0 I 0 Hs I 0 + 2 ∂r ψs +∂η ψs −∂η∂rψs −2 dη ∂η ψ η + 4 − 2 2 ψs Hs ηs Hs 0 Z ηo Z ηo Hs 0 I 0 8 0 I 0 0 +4 1 − 2 dη ∂η ψ η − dη ψ η Hs ηs rs ηs Z ηo 0 Z ηo 0 2 2 I 1 0 η − ηs I 0 1 0 η − ηs I 0 − ∂ηψs dη 0 ∆2ψ η + 2 dη 0 ∆2ψ η Hs rs ηs ηo − η rs ηs ηo − η
– 16 – Z ηo Z ηo Z ηo 1 I 0 I 0 1 0 I 0 2 0 I 0 + −ψs − 2 dη ∂η0 ψ η − dη ψ η dη ∆2ψ η Hs∆η ηs ∆η ηs rs ηs Z ηo Z ηo Z ηo Z ηo 0 1 I 0 0 1 00 I 00 0 I 0 −2 dη ψ (η ) + 2 dη dη ∂η00 ψ (η ) dη ∆2ψ (η ) 0 2 0 2 0 ηs (ηo − η ) ηs (ηo − η ) η ηs 0 5 Hs I 2 1 I 2 1 2 I 2 I I + − 2 ψs + 2 ∂ηψs + − 2 ∂r ψs +∂η ψs −∂η∂rψs 2 Hs Hs Hs 0 1 Hs 2 I 1 I I + −2 + 3 2 + ∂ηψs − ∂rψs ψs Hs Hs Hsrs Hs Z ηo Z ηo I 0 ab 00 I 00 −2∂aψs dη γ0 ∂b dη ψ η 0 ηs η Z ηo Z ηo Z ηo 2 I 0 ab 00 I 00 a 0 I 0 − ∂a ∂ηψs + ∂rv||s dη γ0 ∂b dη ψ η − 2av⊥ s∂a dη ψ η 0 Hs ηs η ηs Z ηo Z ηo 4 0 I 0 I 0 00 I 00 + dη ψ η −ψ η − 2 dη ∂η00 ψ η 0 rs ηs η Z ηo Z ηo Z ηo Z ηo ab 00 I 00 00 I 00 I 0 ab 00 I 00 +γ0 ∂a dη ψ η ∂b dη ψ η +4∂aψs dη γ0 ∂b dη ψ η 0 0 0 η η ηs η Z ηo Z ηo Z ηo 8 0 I 0 00 ab 000 I 000 − dη ∂aψ η dη γ0 ∂b dη ψ η 00 rs ηs ηs η Z ηo Z ηo Z ηo Z ηo 0 I 0 0 1 ab 00 I 00 0 ab I 0 +2∂a dη ψ η 4 dη γ0 dη ∂bψ (η ) − 3 dη γ0 ∂bψ (η ) 0 0 ηs ηs (ηo − η ) η ηs Z ηo Z ηo 0 ab 00 I 00 −6 dη γ0 dη ∂b∂η00 ψ (η ) 0 ηs η Z ηo Z ηo Z ηo Z ηo 0 bc 00 I 00 0 ad 00 I 00 +2∂a dη γ0 ∂c dη ψ η ∂b dη γ0 ∂d dη ψ η 0 0 ηs η ηs η Z ηo Z ηo Z ηo 0 I 0 0 1 00 I 00 1 1 I 0 −4 dη ψ η dη − dη ∆2ψ η + ∆2ψ η 0 3 0 0 2 ηs ηs (ηo − η ) η (ηo − η ) 2 Z ηo Z ηo 0 2 00 I 00 2 1 0 η − ηs I 0 00 + dη ∂η ∆2ψ η + 2 dη ∂θo ψ η 0 0 η (sin θo) rs ηs ηo − η Z ηo Z ηo Z ηo 1 I 0 I 0 0 I 0 ab 0 I 0 +2∂a −ψs − 2 dη∂η0 ψ η − dη ψ (η ) γ0s dη ∂bψ (η ) Hs ηs ηs ηs Z ηo 0 Z ηo Z ηo 4 0 η − ηs I 0 00 ab 000 I 000 + dη ∂b ∆2 ψ (η ) dη γ0 ∂a dη ψ (η ) 0 00 rs ηs ηo − η ηs η Z ηo 0 Z ηo 2 0 η − ηs I 0 I 0 00 I 00 − dη ∆2 ψ η −ψ η − 2 dη ∂η00 ψ η 0 0 rs ηs ηo − η η Z ηo Z ηo ab 00 I 00 00 I 00 +γ0 ∂a dη ψ η ∂b dη ψ η η0 η0 Z ηo Z ηo 00 0 0 I 0 1 00 η − η I 00 −2 dη −2ψ η dη ∆2∂η00 ψ η 0 0 00 ηs ηo − η η ηo − η Z ηo Z ηo 00 0 ab 00 I 00 1 00 η − η I 00 +2γ0 ∂b dη ψ η 0 dη 00 ∂a∆2ψ η η0 ηo − η η0 ηo − η
– 17 – Z ηo Z ηo I 0 00 I 00 1 00 I 00 00 − −2ψ η − 2 dη ∂η ψ η 0 2 dη ∆2ψ η η0 (ηo − η ) η0 Z ηo Z ηo I 0 00 ab 000 I 000 −2∂aψ η dη γ0 ∂b dη ∂η000 ψ η η0 η00 Z ηo Z ηo Z ηo db 00 I 00 00 ac 000 I 000 +2∂a γ0 ∂b dη ψ η dη ∂d γ0 ∂c dη ψ η η0 η0 η00 Z ηo Z ηo ab I 0 00 I 00 00 I 00 +2γ0 ∂a ψ η + dη ∂η00 ψ η ∂b dη ψ η (4.38) η0 η0 and A 1 A v||s A ψs ΛAS = 2 1 − 3ψs v||s − ∂ηψs + 2 ∂rv||s Hsrs Hs Hs Z ηo Z ηo Z ηo 0 A 2 0 I 0 0 I 0 2 0 η − ηs I 0 0 −ψs − dη ψ η − 4 dη ∂η ψ η + dη 0 ∆2ψ η rs ηs ηs rs ηs ηo − η Z ηo 1 A I A 0 I 0 1 A I A2 I A + ∂ηψs − ψs − ψs − 2 dη ∂η0 ψ η + ψs ∂ηψs − 3 ψs + 4ψs ψs Hs ηs Hs 0 1 2 A 2 A A Hs 2 1 A 3 A + − 2 ∂η ψs −∂r ψs −∂η∂rψs + −1+3 2 + ∂ηψs − ∂rψs Hs Hs Hsrs Hs Hs 0 0 Hs A 1 2 Hs 1 A 2 A + 6−4 2 ψs v||s + − 2 ∂r v||s − 2 + 3 2 ∂rv||s ψs + 2 ∂ηψs ∂rv||s Hs Hs Hs Hs Hs 0 Z ηo 0 1 a A 2 A Hs A 1 0 η − ηs I 0 + av⊥s∂aψs + − ∂ηψs + 4 + 2 2 ψs dη 0 ∆2ψ η Hs Hs Hs rs ηs ηo − η Z ηo 0 2 A 0 I 0 2 A Hs A 1 + 2 ψs dη ∆2ψ η + − ∂ηψs + 2 1 + 2 ψs Hsrs ηs Hs Hs rs 0 Z ηo 1 A Hs 2 A 0 I 0 − ∂η∂rψs + 2 + ∂rψs −2 dη ψ η Hs Hs Hsrs ηs 0 1 2 A 2 A A Hs 2 1 A 3 A + 2 ∂η ψs −∂r ψs −∂η∂rψs + 1 − 3 2 − ∂ηψs + ∂rψs Hs Hs Hsrs Hs Hs 0 Z ηo 0 Hs A 0 I 0 3 Hs A2 1 A2 0 −3 1 − 2 ψs −2 dη ∂η ψ η − − 3 2 ψs + 2 ∂ηψs Hs ηs 2 Hs Hs 2 A I 1 2 A 2 A A 1 2 I 2 I I + 2 ∂ηψs ∂ηψs + 2 ∂η ψs −∂r ψs −∂η∂rψs + 2 ∂r ψs +∂η ψs −∂η∂rψs Hs Hs Hs 0 0 0 Hs 2 1 A Hs 2 1 I Hs I 3 A − 3 2 + ∂ηψs + 2 − 3 2 − ∂ηψs − 1 + 2 ψs + ∂rψs Hs Hsrs Hs Hs Hsrs Hs Hs Hs 0 1 I A 1 2 A 2 A A Hs 2 1 A + ∂rψs ψs + − 2 ∂η ψs −∂r ψs −∂η∂rψs + 3 2 + ∂ηψs Hs Hs Hs Hsrs Hs Z ηo Z ηo 3 A I A 0 ab 00 I 00 − ∂rψs ψs + 6 ∂aψs dη γ0 ∂b dη ψ η 0 Hs ηs η Z ηo Z ηo Z ηo 2 A ab 0 I 0 2 A 0 ab 00 I 00 + ∂aψs γ0s ∂b dη ψ η − ∂a ∂ηψs dη γ0 ∂b dη ψ η 0 Hs ηs Hs ηs η Z ηo Z ηo 0 8 0 I A 0 4 0 η − ηs I A 0 + dη ψ ψ η − dη 0 ∆2 ψ ψ η . (4.39) rs ηs rs ηs ηo − η
– 18 – 4.3 Number Counts We have now all we need to evaluate the galaxy number counts at observed redshift and in the direction of observation. Using the relations given in Eqs. (2.7) and (2.8) and inserting the results obtained for the density, Eqs. (4.13) and (4.14), and for the volume fluctuation, Eqs. (4.36)-(4.39), at fixed observed redshift and observer direction, we obtain to first order ! 0 Z ηo Z ηo (1) 2 Hs I A 0 I 0 I 4 0 I 0 ∆ = + v||s + ψs − ψs + 2 dη ∂η0 ψ η − ψs + dη ψ η (0) H2 (0) (0) (0) Hsrs s ηs rs ηs
Z ηo 0 (0) 2 0 η − ηs I 0 1 I A 1 A (1) − dη ∆2ψ η + ∂ηψs + ∂rv||s − 3ψs + ∂ηψs + δρ , (0) (0) η − η0 H H rs ηs 0 s s (4.40) which agrees with [7], where it was derived directly in Poisson gauge. To second order we obtain of course a much more involved expression: ∆(2) = Σ − hΣi (4.41) where Σ = ΣIS + ΣAS (4.42) and 0 Z η0 2 Hs (2) 1 (2) 1 0 h (2) 0 (2) 0i 1 2 0 ΣIS = − − 2 −v||s − φs − dη ∂η φ η + ψ η + v||s Hsrs Hs 2 2 ηs 2 Z ηo 1 I 2 I I 0 I 0 1 a + ψs + −v||s − ψs −ψs − 2 dη ∂η0 ψ η + v⊥sv⊥a s 2 ηs 2 Z ηo Z η0 Z ηo a 0 I 0 0 I 0 I 0 I 0 00 I 00 −2a v⊥s∂a dη ψ η + 4 dη ψ η ∂η0 ψ η + ∂η0 ψ η dη ∂η00 ψ η 0 ηs ηs η Z ηo Z ηo Z ηo I 0 00 2 I 00 ab 00 I 00 00 I 00 +ψ η dη ∂η00 ψ η − γ0 ∂a dη ψ η ∂b dη ∂η00 ψ η η0 η0 η0 Z ηo Z ηo I 0 ab 00 I 00 +2∂a v||s + ψs dη γ0 ∂b dη ψ η 0 ηs η Z ηo Z ηo Z ηo 0 I 0 00 ab 000 I 000 +4 dη ∂a ∂η0 ψ η dη γ0 ∂b dη ψ η 00 ηs ηs η " 0 0 2 00 0 # 1 Hs 3 Hs 1 Hs 1 Hs 1 h 2 I 2 I + 2 + 2 − 3 + 1 + 3 2 + v||s + (ψs ) + 2ψs v||s 2 Hs 2 Hs 2 Hs Hsrs Hs Hsrs # Z ηo Z ηo 2 I 0 I 0 0 I 0 (2) 1 (2) 1 (2) +4 v||s + ψs dη ∂η0 ψ η + 4 dη ∂η0 ψ η − ψs + φs + ∂ηψs ηs ηs 2 2Hs Z ηo 0 Z ηo 1 (2) 1 1 0 η − ηs h (2) (2)i 0 1 0 h (2) (2)i 0 + ∂rv||s − dη 0 ∆2 ψ + φ η + dη ψ + φ η Hs 2 rs ηs ηo − η rs ηs Z ηo 1 2 2 a 1 I 0 I 0 +2 1 − − v||s∂rv||s − (v||s) − v⊥a sv⊥ s + − ∂ηψs − 2 dη ∂η0 ψ η Hsrs Hs Hs ηs Z ηo Z ηo 0 Z ηo 2 0 I 0 2 0 η − ηs I 0 I 0 I 0 0 − dη ψ η + dη 0 ∆2ψ η v||s + −2ψs − 4 dη ∂η ψ η rs ηs rs ηs ηo − η ηs Z ηo Z ηo 0 I 0 1 a 0 I 0 I I −2Hs dη ψ η ∂rv||s + av⊥ s∂a dη ψ η + ∂rψs + 2∂ηψs ηs Hs ηs
– 19 – Z ηo Z ηo 0 2 I 0 0 I 0 +2 dη ∂η0 ψ η −2 dη ψ η ηs ηs Z ηo Z ηo 0 Z ηo I 0 I 0 2 0 η − ηs I 0 2 0 I 0 0 − −ψs − 2 dη ∂η ψ η dη 0 ∆2ψ η − dη ψ η ηs rs ηs ηo − η rs ηs 0 1 I I 3 a 2 a 5 Hs 2 − ∂ηψs − ψs + v⊥a sv⊥ s + av⊥ s∂av||s + + 2 (v||s) Hs 2 Hs 2 Hs 0 Hs 1 1 h 2 2i 1 2 I 2 I I + 5 + 3 2 v||s∂rv||s + 2 v||s∂r v||s + ∂rv||s + − 2 ∂r ψs +∂η ψs −∂η∂rψs Hs Hs Hs Hs 0 0 Z ηo 1 I 3 Hs 2 1 I 4 Hs 0 I 0 − ∂rψs − −1 − 2 − ∂ηψs − −1 − 2 dη ψ η Hs Hs Hs 3 Hsrs rs Hs ηs 0 Z ηo 0 0 Z ηo Hs 2 0 η − ηs I 0 Hs 0 I 0 0 + −2 − 2 dη 0 ∆2ψ η − 2 −2 − 2 dη ∂η ψ η Hs rs ηs ηo − η Hs ηs Z ηo 0 4 I 2 0 I 0 Hs 2 + ψs − 2 dη ∆2ψ η v||s + 2 2 + 2 + ∂rv||s Hsrs Hsrs ηs Hs Hsrs Z ηo 0 Z ηo 2 2 0 I 0 2 Hs 2 2 0 I 0 0 + ∂r v||s dη ψ η + 5 + 3 2 ∂rv||s + 2 ∂r v||s dη ∂η ψ η Hs ηs Hs Hs Hs ηs Z ηo 0 2 1 0 η − ηs I 0 2 I 1 1 2 − ∂rv||s dη 0 ∆2ψ η + 2 ∂ηψs ∂rv||s + ∂r v||s Hs rs ηs ηo − η Hs Hs Hs 0 Z ηo Hs I 2 a 0 I 0 1 a I + 6 + 3 2 ∂rv||s ψs − av⊥ s∂a dη ψ η + av⊥ s∂aψs Hs Hsrs ηs Hs Z ηo Z ηo 2 0 6 ab 0 I 0 4 0 I 0 Hs I − γ0s ∂av||s∂b dη ψ η + 2 dη ψ η + 2 −2 − 2 ψs Hs ηs rs ηs Hs 0 Z ηo 0 Z ηo Hs 0 I 0 2 I 1 Hs 0 2 I 0 0 +4 −2 − 2 dη ∂η ψ η − ∂ηψs + 2 −2 − 2 dη ∂η0 ψ η Hs ηs Hs rs Hs ηs 0 0 Z ηo Hs I Hs I 1 I 0 I 0 +2 −2 − 2 ∂ηψs + −1 − 2 ∂rψs − ∂η∂rψs −2 dη ψ η Hs Hs Hs ηs 0 0 3 Hs 2 1 I 1 I Hs I + −1 − 2 − ∂ηψs + ∂rψs − 2 − 2 ψs Hs Hs 3 Hsrs Hs Hs Z ηo 0 1 2 I 2 I I 0 I 0 Hs I 0 + 2 ∂r ψs +∂η ψs −∂η∂rψs −2 dη ∂η ψ η + −2 − 2 2 ψs Hs ηs Hs 0 Z ηo Z ηo Hs 0 I 0 8 0 I 0 0 +4 −2 − 2 dη ∂η ψ η − dη ψ η Hs ηs rs ηs Z ηo 0 Z ηo 0 2 2 I 1 0 η − ηs I 0 1 0 η − ηs I 0 − ∂ηψs dη 0 ∆2ψ η + 2 dη 0 ∆2ψ η Hs rs ηs ηo − η rs ηs ηo − η Z ηo Z ηo Z ηo 1 I 0 I 0 1 0 I 0 2 0 I 0 + −ψs − 2 dη ∂η0 ψ η − dη ψ η dη ∆2ψ η Hs∆η ηs ∆η ηs rs ηs Z ηo Z ηo Z ηo Z ηo 0 1 I 0 0 1 00 I 00 0 I 0 −2 dη ψ (η ) + 2 dη dη ∂η00 ψ (η ) dη ∆2ψ (η ) 0 2 0 2 0 ηs (ηo − η ) ηs (ηo − η ) η ηs 0 1 Hs I 2 1 I 2 1 2 I 2 I I + − − 2 ψs + 2 ∂ηψs + − 2 ∂r ψs +∂η ψs −∂η∂rψs 2 Hs Hs Hs 0 1 Hs 2 I 1 I I + 4 + 3 2 + ∂ηψs − ∂rψs ψs Hs Hs Hsrs Hs
– 20 – Z ηo Z ηo I 1 I 0 ab 00 I 00 +2∂a −ψs − ∂ηψs + ∂rv||s dη γ0 ∂b dη ψ η 0 Hs ηs η Z ηo Z ηo Z ηo a 0 I 0 4 0 I 0 I 0 00 I 00 −2av⊥ s∂a dη ψ η + dη ψ η −ψ η − 2 dη ∂η00 ψ η 0 ηs rs ηs η Z ηo Z ηo Z ηo Z ηo ab 00 I 00 00 I 00 I 0 ab 00 I 00 +γ0 ∂a dη ψ η ∂b dη ψ η +4∂aψs dη γ0 ∂b dη ψ η 0 0 0 η η ηs η Z ηo Z ηo Z ηo 8 0 I 0 00 ab 000 I 000 − dη ∂aψ η dη γ0 ∂b dη ψ η 00 rs ηs ηs η Z ηo Z ηo Z ηo Z ηo 0 I 0 0 1 ab 00 I 00 0 ab I 0 +2∂a dη ψ η 4 dη γ0 dη ∂bψ (η ) − 3 dη γ0 ∂bψ (η ) 0 0 ηs ηs (ηo − η ) η ηs Z ηo Z ηo 0 ab 00 I 00 −6 dη γ0 dη ∂b∂η00 ψ (η ) 0 ηs η Z ηo Z ηo Z ηo Z ηo 0 bc 00 I 00 0 ad 00 I 00 +2∂a dη γ0 ∂c dη ψ η ∂b dη γ0 ∂d dη ψ η 0 0 ηs η ηs η Z ηo Z ηo Z ηo 0 I 0 0 1 00 I 00 1 1 I 0 −4 dη ψ η dη − dη ∆2ψ η + ∆2ψ η 0 3 0 0 2 ηs ηs (ηo − η ) η (ηo − η ) 2 Z ηo Z ηo 0 2 00 I 00 2 1 0 η − ηs I 0 00 + dη ∂η ∆2ψ η + 2 dη ∂θo ψ η 0 0 η (sin θo) rs ηs ηo − η Z ηo Z ηo Z ηo 1 I 0 I 0 0 I 0 ab 0 I 0 +2∂a −ψs − 2 dη∂η0 ψ η − dη ψ (η ) γ0s dη ∂bψ (η ) Hs ηs ηs ηs Z ηo 0 Z ηo Z ηo 4 0 η − ηs I 0 00 ab 000 I 000 + dη ∂b ∆2 ψ (η ) dη γ0 ∂a dη ψ (η ) 0 00 rs ηs ηo − η ηs η Z ηo 0 Z ηo 2 0 η − ηs I 0 I 0 00 I 00 − dη ∆2 ψ η −ψ η − 2 dη ∂η00 ψ η 0 0 rs ηs ηo − η η Z ηo Z ηo ab 00 I 00 00 I 00 +γ0 ∂a dη ψ η ∂b dη ψ η η0 η0 Z ηo Z ηo 00 0 0 I 0 1 00 η − η I 00 −2 dη −2ψ η dη ∆2∂η00 ψ η 0 0 00 ηs ηo − η η ηo − η Z ηo Z ηo 00 0 ab 00 I 00 1 00 η − η I 00 +2γ0 ∂b dη ψ η 0 dη 00 ∂a∆2ψ η η0 ηo − η η0 ηo − η Z ηo Z ηo I 0 00 I 00 1 00 I 00 00 − −2ψ η − 2 dη ∂η ψ η 0 2 dη ∆2ψ η η0 (ηo − η ) η0 Z ηo Z ηo I 0 00 ab 000 I 000 −2∂aψ η dη γ0 ∂b dη ∂η000 ψ η η0 η00 Z ηo Z ηo Z ηo db 00 I 00 00 ac 000 I 000 +2∂a γ0 ∂b dη ψ η dη ∂d γ0 ∂c dη ψ η η0 η0 η00 Z ηo Z ηo ab I 0 00 I 00 00 I 00 +2γ0 ∂a ψ η + dη ∂η00 ψ η ∂b dη ψ η η0 η0 0 Z ηo 2 Hs I 0 I 0 1 I 0 + + 2 v||s + ψs + 2 dη ∂η ψ η − 3v||s + ∂rv||s − 4ψs Hsrs Hs ηs Hs
– 21 – Z ηo Z ηo Z ηo 0 0 I 0 4 0 I 0 2 0 η − ηs I 0 1 I (1) 0 −6 dη ∂η ψ η + dη ψ η − dη 0 ∆2ψ η + ∂ηψs δρ ηs rs ηs rs ηs ηo − η Hs Z ηo 1 (1) (1) 1 I 0 I 0 + ∂η ρ¯ δρ − ∂rδρ −v||s − ψs − 2 dη ∂η0 ψ η ρ¯ Hs ηs Z ηo Z ηo Z ηo (1) 0 I 0 (1) 0 ab 00 I 00 (2) +2∂rδρ dη ψ η − 2∂aδρ dη γ0 ∂b dη ψ η + δρ . (4.43) 0 ηs ηs η The contributions proportional to ψA which vanish for ΛCDM are given by
0 Z ηo 2 Hs A 3 A 2 I A 0 I A 0 0 ΣAS = − − 2 −ψs v||s + (ψs ) − 3ψs ψs − 4 dη ∂η (ψ ψ ) η Hsrs Hs 2 ηs " Z η0 Z ηo Z ηo 0 0 2 A 0 I 0 A 0 ab 00 I 00 1 Hs 3 Hs −2ψs dη ∂η0 ψ η − 2∂aψs dη γ0 ∂b dη ψ η + + 0 2 2 ηs ηs η 2 Hs 2 Hs 00 0 Z ηo 1 Hs 1 Hs 1 A2 a I 0 I 0 0 − 3 + 1 + 3 2 + ψs + 2ψs −v||s − ψs − 2 dη ∂η ψ η 2 Hs Hsrs Hs Hsrs ηs A 1 A 1 A ψs +2 1 − 3ψs v||s − v||s∂ηψs + 2 ∂rv||s Hsrs Hs Hs Z ηo Z ηo Z ηo 0 A 2 0 I 0 0 I 0 2 0 η − ηs I 0 0 −ψs − dη ψ η − 4 dη ∂η ψ η + dη 0 ∆2ψ η rs ηs ηs rs ηs ηo − η Z ηo 1 A I A 0 I 0 1 A I A2 I A + ∂ηψs − ψs − ψs − 2 dη ψ η + ψs ∂ηψs − 3 ψs + 4ψs ψs Hs ηs Hs 0 1 2 A 2 A A Hs 2 1 A 3 A + − 2 ∂η ψs −∂r ψs −∂η∂rψs + 5 + 3 2 + ∂ηψs − ∂rψs Hs Hs Hsrs Hs Hs 0 0 Hs A 1 2 Hs 1 A 2 A − 6 + 4 2 ψs v||s + − 2 ∂r v||s − 8 + 3 2 ∂rv||s ψs + 2 ∂ηψs ∂rv||s Hs Hs Hs Hs Hs 0 Z ηo 0 1 a A 2 A Hs A 1 0 η − ηs I 0 + av⊥s∂aψs + − ∂ηψs + 10 + 2 2 ψs dη 0 ∆2ψ η Hs Hs Hs rs ηs ηo − η Z ηo 0 2 A 0 I 0 2 A Hs A 1 + 2 ψs dη ∆2ψ η + − ∂ηψs + 2 4 + 2 ψs Hsrs ηs Hs Hs rs 0 Z ηo 1 A Hs 2 A 0 I 0 − ∂η∂rψs + 3 + 2 + ∂rψs −2 dη ψ η Hs Hs Hsrs ηs 0 1 2 A 2 A A Hs 2 1 A 3 A + 2 ∂η ψs −∂r ψs −∂η∂rψs − 5 + 3 2 + ∂ηψs + ∂rψs Hs Hs Hsrs Hs Hs 0 Z ηo 0 Hs A 0 I 0 15 Hs A2 1 A2 0 + 3 2 + 2 ψs −2 dη ∂η ψ η + + 3 2 ψs + 2 ∂ηψs Hs ηs 2 Hs Hs 2 A I 1 2 A 2 A A 1 2 I 2 I I + 2 ∂ηψs ∂ηψs + 2 ∂η ψs −∂r ψs −∂η∂rψs + 2 ∂r ψs +∂η ψs −∂η∂rψs Hs Hs Hs 0 0 0 Hs 2 1 A Hs 2 1 I Hs I + 2−3 2 − ∂ηψs + 4−3 2 − ∂ηψs − 7+ 2 ψs Hs Hsrs Hs Hs Hsrs Hs Hs 3 A 1 I A 1 2 A 2 A A + ∂rψs + ∂rψs ψs + − 2 ∂η ψs −∂r ψs −∂η∂rψs Hs Hs Hs 0 Hs 2 1 A 3 A I + 6 + 3 2 + ∂ηψs − ∂rψs ψs Hs Hsrs Hs Hs
– 22 – Z ηo Z ηo Z ηo A 0 ab 00 I 00 2 A ab 0 I 0 +6 ∂aψs dη γ0 ∂b dη ψ η + ∂aψs γ0s ∂b dη ψ η 0 ηs η Hs ηs Z ηo Z ηo Z ηo 2 A 0 ab 00 I 00 8 0 I A 0 − ∂a ∂ηψs dη γ0 ∂b dη ψ η + dη ψ ψ η 0 Hs ηs η rs ηs Z ηo 0 0 4 0 η − ηs I A 0 2 Hs A 1 A (1) − dη 0 ∆2 ψ ψ η + − − 2 ψs + ∂ηψs δρ rs ηs ηo − η Hsrs Hs Hs 1 (1) (1) 1 A + ∂η ρ¯ δρ − ∂rδρ ψs . (4.44) ρ¯ Hs 4.4 Number Counts: leading terms Let us finally give a more concise expression for the second order number counts in Eqs. (4.43) and (4.44). We will neglect all the subleading terms, keeping only the leading (potentially observable) terms in the large redshift and small angle limit, where the pure Doppler and potential terms can be neglected. This means we keep only the terms with four spatial derivatives of metric perturbations (the velocity has one spatial derivative at first order, see Eq. (3.18)) and terms with the density and two spatial derivatives. We also keep the density, the redshift space distortions and the lensing at second order. Restricting to these dominant terms we obtain, after some algebraic manipulation, the following leading order result for Σ:
Z ηo 0 Leading 1 (2) 1 1 0 η − ηs (2) 0 (2) 0 (2) Σ ' ∂rv||s − dη 0 ∆2 ψ η + φ η + δρ Hs 2 rs ηs ηo − η Z ηo 0 1 h 2 2i 2 1 0 η − ηs I 0 + 2 v||s∂r v||s + ∂rv||s − ∂rv||s dη 0 ∆2ψ η Hs Hs rs ηs ηo − η Z ηo Z ηo Z ηo 0 2 2 0 ab 00 I 00 1 0 η − ηs I 0 − ∂a ∂rv||s dη γ0 ∂b dη ψ η + 2 dη ∆2ψ η 0 0 Hs ηs η rs ηs ηo − η Z ηo 0 Z ηo Z ηo 4 0 η − ηs I 0 00 ab 000 I 000 + dη ∂b ∆2ψ (η ) dη γ0 ∂a dη ψ (η ) 0 00 rs ηs ηo − η ηs η Z ηo Z ηo 1 ab 00 I 00 00 I 00 +∆2 − γ0 ∂a dη ψ η ∂b dη ψ η 2 η0 η0 Z ηo Z ηo Z ηo 00 0 0 ab 00 I 00 1 00 η − η I 00 −4 dη γ0 ∂b dη ψ η dη ∂a∆2ψ η 0 0 0 00 ηs η ηo − η η ηo − η η 0 1 2 Z o η − ηs 1 + ∂ v − dη0 ∆ ψI η0 δ(1) + v ∂ δ(1) r ||s 0 2 ρ ||s r ρ Hs rs ηs ηo − η Hs Z ηo Z ηo (1) 0 ab 00 I 00 −2∂aδρ dη γ0 ∂b dη ψ η . (4.45) 0 ηs η All these terms come from density fluctuations, redshift space distortions and lensing-like terms. The ΛCDM expression corresponding to Eq. (4.45) can be easily obtained substituting ψI with ψ (in this limit ψI = ψ and ψA = 0). The terms indicated in bold face are the once which we evaluate numerically in Section5 for ΛCDM. The other terms, which are really new physical effects (like e.g. second order lensing), can also contribute significantly as they have four transverse derivatives We shall study them in detail in a forthcoming publication [35].
– 23 – 5 The Bispectrum and its numerical evaluation
In this section we compute two new contributions to the bispectrum of the number count fluctuations and evaluate them numerically for some configurations, namely the contribution from redshift space distortions and the one from lensing. Since we consider the truly observed galaxy number count, and not fluctuations on some unobservable spatial hypersurface, the 3-point function is in general a function of 3 directions2 and 3 redshifts,
B (n1, n2, n3, z1, z2, z3) = h∆ (n1, z1) ∆ (n2, z2) ∆ (n3, z3)i , (5.1)
The leading non-vanishing terms are in the form
(2) (1) (1) (1) (1) h∆ (n1, z1) ∆ (n2, z2) ∆ (n3, z3)i = hΣ(n1, z1) ∆ (n2, z2) ∆ (n3, z3)i (1) (1) −hΣ(n1, z1)ih∆ (n2, z2) ∆ (n3, z3)i (5.2) plus the permutations with respect to the galaxy positions (ni, zi). Expanding the direction dependence of ∆ in spherical harmonics, X ∆ (n, z) = a`m(z)Y`m(n) , `,m we can write X B (n , n , n , z , z , z ) = Bm1m2m3 (z , z , z )Y (n )Y (n )Y (n ) , 1 2 3 1 2 3 `1`2`3 1 2 3 `1m1 1 `2m2 2 `3m3 3 `1,`2,`3,m1,m2,m3 (5.3) where
Bm1m2m3 (z , z , z ) = ha (z )a (z )a (z )i `1`2`3 1 2 3 `1m1 1 `2m2 2 `3m3 3 Z ∗ ∗ ∗ = dΩ1dΩ2dΩ3B (n1, n2, n3, z1, z2, z3) Y`1m1 (n1) Y`2m2 (n2) Y`3m3 (n3) . (5.4)
Statistical isotropy requires that
Bm1m2m3 (z , z , z ) = Gm1,m2,m3 b (z , z , z ) , (5.5) `1`2`3 1 2 3 `1,`2,`3 `1,`2,`3 1 2 3 where Gm1,m2,m3 is the Gaunt integral given by `1,`2,`3 Z Gm1,m2,m3 = dΩ Y (n) Y (n) Y (n) (5.6) `1,`2,`3 `1m1 `2m2 `3m3 r ` ` ` ` ` ` (2` + 1)(2` + 1)(2` + 1) = 1 2 3 1 2 3 1 2 3 . (5.7) 0 0 0 m1 m2 m3 4π On the second line we have expressed the Gaunt integral in terms of the Wigner 3j symbols, see e.g. [36]. This integral vanishes if `3 6∈ [|`1 − `2| , `1 + `2], if the sum `1 + `2 + `3 is odd, or if m1 + m2 + m3 6= 0. The quantity b`1,`2,`3 (z1, z2, z3) is the reduced bispectrum of the number counts. For Gaussian initial conditions it vanishes in linear perturbation theory
2We remember that n is the photon direction, so from the source to the observer. Hence we observer the galaxy in direction −n. But, for the sake of simplicity, we prefer to use n to indicate the galaxy position.
– 24 – (2) (1) (1) and its first contributions are terms of the form h∆ (n1, z1) ∆ (n2, z2) ∆ (n3, z3) + (1) (2) (1) (1) (1) (2) ∆ (n1, z1) ∆ (n2, z2) ∆ (n3, z3) + ∆ (n1, z1) ∆ (n2, z2) ∆ (n3, z3)i. Here we want to determine the presumably largest contributions to this expression. In subsequent work [35] we shall study it in more detail. We write the perturbations in terms of Fourier modes and we express them in terms of the power spectrum of the primordial curvature perturbation Rin (k). In this section, we assume simple adiabatic Gaussian initial perturbations from inflation which are given by the curvature power spectrum,