Arxiv:1910.01880V2 [Gr-Qc] 30 Oct 2019
Total Page:16
File Type:pdf, Size:1020Kb
2 II. INDUCED GRAVITATIONAL WAVES We start with deriving the expression for the energy density of induced GWs. The evolution equation for the induced GWs hij(η; ~x) is of the form 00 0 2 hij + 2Hhij − @ hij = Λij;klSkl; (1) where a dash denotes differentiation with respect to the conformal time η and Λij;kl is the projection tensor which extracts the transverse-traceless part of the source term Skl(η; ~x). A more explicit expression for Λij;klSkl will be given shortly below. We consider a universe filled with a perfect fluid whose equation of state parameter w is a constant. The scale factor is then given by a / η2=(1+3w) and hence H := a0=a = [2=(1 + 3w)]/η. The Fourier components of hij are defined by Z 3 X d k (A) ~ h (η; ~x) = h (η;~k)e (~k)eik·~x (A = +; ×): (2) ij (2π)3=2 A ij A (A) ~ The two polarization tensors eij (k) are defined as (+) ~ 1 h ~ ~ ~ ~ i e (k) = p ei(k)ej(k) − ei(k)ej(k) ; (3) ij 2 (×) ~ 1 h ~ ~ ~ ~ i e (k) = p ei(k)ej(k) + ei(k)ej(k) ; (4) ij 2 ~ ~ ~ where ei(k) and ej(k) are unit vectors orthogonal to each other and to k. It follows from the definition of the 0 i (A) (A) ~ (A ) ~ polarization tensors that k eij = 0 and eij (k)eij (k) = δAA0 . Using the polarization tensors one can write the right-hand side of Eq. (1) as Z 3 X d k (A) (A) Λ S (η; ~x) = e (~k)e (~k)S (η;~k); (5) ij;kl kl (2π)3=2 ij lm lm A where 3 Z d x ~ S (η;~k) := S (η; ~x)e−ik·~x (6) ij (2π)3=2 ij ~ (A) ~ ~ is the Fourier transform of Sij(η; ~x). We write SA(η; k) := eij (k)Sij(η; k). The formal solution to Eq. (1) in the Fourier domain is given by Z η ~ 1 ~ hA(η; k) = Gk(η; η¯)a(¯η)SA(¯η; k)d¯η; (7) a(η) 0 where the Green's function is expressed in terms of the Bessel functions as π G (η; η¯) = η1=2η¯1=2 [Y (kη)J (kη¯) − J (kη)Y (kη¯)] ; (8) k 2 ν ν ν ν with 3(1 − w) ν := : (9) 2(1 + 3w) The energy density of GWs is given by M 2 X 1 ρ (η) = Pl hh0 (η; ~x)h0 (η; ~x)i; (10) GW 2 a2 A A A where h· · · i denotes a spatial average and 3 Z d k ~ h (η; ~x) := h (η;~k)eik·~x: (11) A (2π)3=2 A 3 Since 3 3 0 Z d k d k ~ ~ 0 hh0 (η; ~x)h0 (η; ~x)i = h0 (η;~k)h∗ 0(η;~k0)hei(k−k )·~xi; (12) A A (2π)3=2 (2π)3=2 A A and 3 ~ ~ 0 1 Z ~ ~ 0 (2π) hei(k−k )·~xi = ei(k−k )·~xd3x = δ(3)(~k − ~k0); (13) V V where V is a volume whose size is much larger than the wavelengths of interest, we have M 2 1 X Z ρ (η) = Pl d3kjh0 (η;~k)j2: (14) GW 2 V a2 A A In the subhorizon regime, kη 1, the time derivative of hA is approximately given by Z η 0 ~ 1 ~ hA(η; k) ' @ηGk(η; η¯)a(¯η)SA(¯η; k)d¯η: (15) a(η) 0 We find the following approximate expression for the time derivative of the Green's function, π @ G (η; η¯) 'G (η; η¯) := (kη)1=2(kη¯)1=2 [Y (kη)J (kη¯) − J (kη)Y (kη¯)] : (16) η k k 2 ν−1 ν ν−1 ν More explicitly, for ν = 1=2 (w = 1=3) we have Gk(η; η¯) = cos[k(η − η¯)]; (17) and for ν = 3=2 (w = 0) we have sin[k(η − η¯)] G (η; η¯) = cos[k(η − η¯)] + : (18) k kη¯ It then follows that M 2 1 X Z Z η Z η ρ (η) = Pl d3k dη0 dη00a(η0)a(η00)G (η; η0)G (η; η00)S (η0;~k)S∗ (η00;~k): (19) GW 2 V a4(η) k k A A A 0 0 We are interested in the case where SA is of the form Z d3q S (η;~k) = e(A)(~k) qiqjA(~q)A(~k − ~q)F (~k; ~q; η); (20) A ij (2π)3=2 where A(~q) is a Gaussian random field and F (~k; ~q; η) is some function. In the actual calculation, A(~q) will be the primordial amplitude of scalar perturbations and F contains the information of their time evolution. The power spectrum P(k) is defined by 2π2 hA(~k)A∗(~q)i = δ(3)(~k − ~q)P(k): (21) k3 (Now h· · · i denotes an average over the whole distribution.) Using Wick's theorem, the two-point correlator of the source term SA can be written as Z 3 4 0 ~ ∗ 00 ~ 0 d q (A) ~ i j 2 8π (3) ~ ~ 0 ~ ~ 0 ~ 00 hSA(η ; k)SA(η ; k )i = (e (k)q q ) δ (k − k )P(q)P(jk − ~qj)F (k; ~q; η )F (k; ~q; η ): (22) (2π)3 ij q3j~k − ~qj3 Using this and [δ(3)(~k − ~k0)]2 = [V=(2π)3]δ(3)(~k − ~k0), one can write the energy density of GWs (19) as 2 Z 3 Z 3 Z η Z η 2MPl d k d q 0 00 0 00 0 00 ρGW(η) = 4 3 3 dη dη a(η )a(η )Gk(η; η )Gk(η; η ) a (η) (2π) (2π) 0 0 π4 × q4 sin4 θ P(q)P(j~k − ~qj)F (~k; ~q; η0)F (~k; ~q; η00); (23) q3j~k − ~qj3 4 where θ is the angle between ~q and ~k. The final step is to extract from the above expression the energy density parameter of GWs, ΩGW(η; k), defined by Z 2 2 ρGW(η) = 3MPlH ΩGW(η; k)d ln k: (24) In practice, F depends on ~k, ~q, and η through qη and j~k − ~qjη. Therefore, it is convenient to introduce dimensionless variables u := j~k − ~qj=k and v := q=k. Using these variables, the energy density parameter is expressed as Z 1 Z 1+v Z η Z η 1 1 0 00 0 00 0 00 ΩGW(η; k) = 2 4 dv du dη dη a(η )a(η )Gk(η; η )Gk(η; η ) 3H a (η) 0 j1−vj 0 0 " #2 k4 v2 1 + v2 − u2 2 × 1 − P(ku)P(kv)F (u; v; kη0)F (u; v; kη00) 4 u2 2v 2 2 Z 1 Z 1+v 2 " 2 2 2# k v 1 + v − u 2 = 2 dv du 2 1 − P(ku)P(kv)I (u; v; kη); (25) 12H 0 j1−vj u 2v where Z η k 0 0 0 0 I(u; v; kη) = dη a(η )Gk(η; η )F (u; v; kη ): (26) a(η) 0 The information of the initial conditions for scalar perturbations is encoded in the power spectrum P, while the time evolution of the perturbations determines the form of the integral I. The two distinct effects are thus separated. The gauge difference is essentially imprinted in I. This can be evaluated analytically for w = 0 and w = 1=3 [25] and numerically for the other values of w. III. ACTION APPROACH TO INDUCED GRAVITATIONAL WAVES To derive the gauge-ready form of the source term for induced GWs in a universe filled with an irrotational barotropic perfect fluid, it is convenient to employ the action-based approach, describing the fluid in terms of a shift-symmetric k-essence field. Our action is given by Z p M 2 S = d4x −g Pl R + P (X) ; (27) 2 µν where X := −g @µφ∂ν φ/2. The k-essence field is equivalent to a cosmological perfect fluid whose energy density and pressure are given respectively by ρ = 2XPX − P and p = P . (Here and hereafter we write @P=@X = PX .) Therefore, a w = p/ρ = const fluid can be mimicked by [26, 27] P / X(1+w)=2w: (28) Hereafter we will assume that w ≥ 0. Note that the w = 0 case appears to be singular, but a careful inspection shows that the limit in fact makes sense [23]. The metric in the 3 + 1 Arnowitt-Deser-Misner (ADM) form is given by 2 2 2 i i j j ds = −N dt + gij dx + N dt dx + N dt : (29) We only fix the spatial coordinate, while leaving the temporal gauge degree of freedom unfixed, as we are interested only in the latter gauge difference. We therefore write the ADM variables in terms of the scalar and tensor perturbations as 2 −2 h N = 1 + δN; Ni = @iχ, gij = a e (e )ij; (30) h where (e )ij = δij + hij + hikhkj=2 + ··· . The perturbed scalar field is written as φ = φ¯(t) + δφ. (31) 5 We will omit the bar from the background value if unnecessary. Using the temporal gauge degree of freedom one can eliminate one of χ, , and δφ. The background equations are given by 2 2 3MPlH = 2XPX − P; (32) 2 2 _ MPl 3H + 2H = −P; (33) d a3φP_ = 0: (34) dt X In order to derive the equations of motion for the perturbations, we substitute the metric (30) and the scalar field (31) to the action (27) and expand it to third order.