30 S. Alexander, N. Yunes / <strong>Physics</strong> <strong>Reports</strong> 480 (2009) 1–55where α = κ, commas in index lists stand for partial differentiation, η is conformal time and ˜ɛ ijk = ˜ɛ ηijk . Variation of thelinearized action with respect to the metric perturbation yields the linearized field equations, namely [29]¯□h j i := 1 a 2 ˜ɛpk(j [( ϑ ,ηη − 2Hϑ ,η)hi)k,pη + ϑ ,η ¯□h i)k,p], (177)where ¯□ is the D’Alembertian operator associated with the background, namely¯□f = f ,ηη + 2Hf ,η − δ ij f ,ij , (178)with f some function of all coordinates and the conformal Hubble parameter H := a ,η /a. One could have, of course, obtainedthe same linearized field equations by perturbatively expanding the C-tensor.One can see from Eq. (177) that the evolution of GW perturbations is governed by second and third derivatives of theGW tensor. Jackiw and Pi [22] were the first to point out that for the canonical choice of ϑ the GW evolution is governedby the D’Alembertian of flat space only, if we neglect corrections due to the expansion history of the Universe (i.e. if thisvanishes, then the linearized <strong>modified</strong> field equations for the GW perturbation are satisfied to linear order). Such a resultimplies there are two linearly independent polarizations that propagate at the speed of light.Let us now concentrate on gravitational wave perturbations, for which one can make the ansatzh ij = A ija(η) exp [ −i ( φ(η) − κn k χ k)] , (179)where the amplitude A ij , the unit vector in the direction of wave propagation n k and the conformal wavenumber κ > 0 areall constant. It is convenient to decompose the amplitude into definite parity states, such asA ij = A R e R ij + A Le L ij(180)where the circular polarization tensors e R,Lijare given in terms of the linear ones e +,×ijby [21]e R = 1 (kl√ e+kl+ ie × )kl2(181)e L kl = 1 √2(e+kl− ie × kl). (182)These polarization tensors satisfy the conditionn i ɛ ijk e R,Lkl= iλ R,L(ejl) R,L, (183)where λ R = +1 and λ L = −1.With this decomposition, the linearized <strong>modified</strong> field equations [Eq. (177)] reduce to[iφ,ηη R,L + ( )φ,ηR,L 2] (+ H ,η + H 2 − κ 2 1 − λ )R,Lκϑ ,η= iλ R,Lκ ( ) (ϑ,ηη − 2Hϑ ,η φR,La 2 a 2,η − iH) . (184)Before attempting to solve this equation for arbitrary ϑ(η) and a(η), it is instructive to take the flat-space limit, that isa → 1, and thus, ȧ → 0. Assuming further that time derivatives of the CS scalar do not scale with H, one finds that theabove equation reduces to [29,33](i ¨φ R,L + ˙φ 2 R,L − k2) ( 1 − λ R,L k ˙ϑ ) = iλ R,L k ¨ϑ ˙φ R,L , (185)where k is the physical wavenumber 3-vector, t stands for cosmic time and overhead dots stand for partial differentiationwith respect to time. Let us further assume that the GW phase satisfies φ ,tt /φ 2 ,t≪ 1, which is the standard short-wavelengthapproximation, as well as ¨ϑ = ϑ 0 = const. Then the above equation can be solved to first order in ϑ to findφ(t) = φ 0 + kt + iλ R,Lϑ 0kt + O(ϑ) 2 , (186)2where φ 0 is a constant phase offset and the uncontrolled remainder O(ϑ) 2 stands for terms of the form k 2 ˙ϑ 2 or k ˙ϑ ¨ϑ. Theimaginary correction to the phase then implies an exponential enhancement/suppression effect of the GW amplitude, as thispropagates in CS <strong>modified</strong> gravity. Recall that here we are interested in the propagation of GWs, which is why the right-handside of Eq. (185) implicitly omits the stress–energy tensor. We shall see in the next section that if there is a stress–energytensor, then the CS correction depends both on the first and second derivatives of the CS scalar. Lastly, if we had not assumedthat ¨ϑ = const. then the solution would have becomeφ(t) = φ 0 + kt + iλ R,L2 k ˙ϑ(t) + O(ϑ) 2 , (187)which still exhibits the exponential suppression/enhancement effect.
S. Alexander, N. Yunes / <strong>Physics</strong> <strong>Reports</strong> 480 (2009) 1–55 31The exponentially growing modes could be associated with instabilities in the solution to gravitational wave propagation.One must be careful, however, to realize that the results above have been obtained within the approximation k 2 ˙ϑ 2 ≪ 1 ≫k ˙ϑ ¨ϑ. Thus, provided ˙ϑ is smaller than the age of the universe, then the instability time scale will not have enough time toset in. For larger values of ˙ϑ, the approximate solutions we presented above break down and one must account for higherorder corrections. A final caveat to keep in mind is that these results are derived within the non-dynamical formulation ofthe theory; GW solutions in the fully dynamical theory are only now being actively investigated.Let us now return to the field equation for the phase in an FRW background [Eq. (184)]. The solution to this equation isnow complicated by the fact that the scale factor also depends on conformal time, and thus, one cannot find a closed formsolution prior to specifying the evolution of a(η). Let us then choose a matter-dominated cosmological model, in whicha(η) = a 0 η 2 = a 0 /(1 + z), where a 0 is the value of the scale factor today and z is the redshift. It then follows that theconformal Hubble parameter is simply given by H = 2/η = 2(1 + z) 1/2 . With this choice, one can now compute the CScorrection to the accumulated phase as the plane-wave propagates from some initial conformal time to η, namely [33]∫ 1[ 1∆φ (R,L) = iλ R,L kH 04 ϑ ,ηη(η) − 1 ] dηη ϑ ,η(η)η + 4 O(ϑ)2 , (188)ηwhich one can check reduces to the flat space result of above in the right limit. Note that the exponentialenhancement/suppression effect now depends on an integrated measure of the evolution of the CS scalar and the scale factor.The CS correction to the GW amplitude derives from a modification to the evolution equations of the gravitationalperturbation, but it also leads to important observational consequences. One of these can be understood by considering aGW generated by a binary black hole system in the early inspiral phase. The GW produced by such a system can be describedas followsh R,L = √ 2 M ( )Mk0 (t) 2/3(1 + λR,L cos ι ) 2 [ ( )]exp −i Ψ (t) + ∆φR,L , (189)d L 2where d L = a 0 η (1 + z) is the luminosity distance to the binary’s center of mass, Ψ (t) is the GW phase described by GR,k 0 (t) is the instantaneous wave number of the gravitational wavefront passing the detector and M is the comoving chirpmass, which is a certain combination of the binary mass components. The inclination angle ι, the angle subtended by theorbital angular momentum and the observer’s line of sight, can be isolated ash R= 1 + cos ι [ ] 2k(t)1 − cos ι exp ζ = 1 + cos ῑ1 − cos ῑ , (190)h LH 0where we have defined∫ 1[ 1ζ := H 2 0η 4 ϑ ,ηη(η) − 1 ] dηη ϑ ,η(η)η . (191)4We see then that from an GW observational standpoint, the CS correction leads to an apparent inclination angle ῑ, whicheffectively modifies the actual inclination angle by a factor that depends on the integrated history of the CS correction:sinhcos ῑ =cosh(k 0 (t)ξ(z))H 0+ cosh( )k 0 (t)ξ(z)H 0+ sinh(k 0 (t)ξ(z))cosH 0ι( )k 0 (t)ξ(z)cosH 0ι∼ cos ι + k 0(t)ξ(z)H 0sin 2 ι + O ( ξ 2) . (192)We see then that the CS correction effectively introduces an apparent evolution of the inclination angle, which tracks thegravitational wave frequency.The interpretation of the CS correction as inducing an effective inclination angle should be interpreted with care. In GR, if agravitational wave propagates along the line of sight, such that the actual inclination angle is (eg. 0 or π), then the amplitudeis a maximum. In CS gravity, however, the amplitude can be either enhanced or suppressed, depending on whether the waveis right- or left-circularly polarized. When the CS effect suppresses the GW amplitude, one can think of this as an effectivemodification of the inclination angle away from the maximum. However, when the CS effect enhances the amplitude, thereis no real inclination angle that can mimic this effect (i.e. the effective angle would have to be imaginary).The evolution equation for the gravitational wave perturbation depends sensitively on the scale factor evolution[see, eg. Eq. (177)]. Alexander and Martin [29] have investigated gravitational wave solutions when the scale factorpresents an inflationary behavior. Suffice it to say in this section that Eq. (177) can be recast in the form of a parametricoscillator equation, with a non-trivial effective potential. In certain limits appropriate to inflation, one can solve thisdifferential equation in terms of Whittaker functions, which can be decomposed into products of trigonometric functionsand exponentials. In essence, the solutions present the same structure as that of a matter-dominated cosmology. Once thegravitational wave modes have been computed, one can proceed to calculate the power spectrum, but these results will bediscussed further in Section 8.1.
- Page 1 and 2: Physics Reports 480 (2009) 1-55Cont
- Page 3 and 4: S. Alexander, N. Yunes / Physics Re
- Page 5 and 6: S. Alexander, N. Yunes / Physics Re
- Page 7 and 8: S. Alexander, N. Yunes / Physics Re
- Page 9 and 10: S. Alexander, N. Yunes / Physics Re
- Page 11 and 12: S. Alexander, N. Yunes / Physics Re
- Page 13 and 14: S. Alexander, N. Yunes / Physics Re
- Page 15 and 16: S. Alexander, N. Yunes / Physics Re
- Page 17 and 18: S. Alexander, N. Yunes / Physics Re
- Page 19 and 20: S. Alexander, N. Yunes / Physics Re
- Page 21 and 22: S. Alexander, N. Yunes / Physics Re
- Page 23 and 24: S. Alexander, N. Yunes / Physics Re
- Page 25 and 26: S. Alexander, N. Yunes / Physics Re
- Page 27 and 28: S. Alexander, N. Yunes / Physics Re
- Page 29: S. Alexander, N. Yunes / Physics Re
- Page 33 and 34: S. Alexander, N. Yunes / Physics Re
- Page 35 and 36: S. Alexander, N. Yunes / Physics Re
- Page 37 and 38: S. Alexander, N. Yunes / Physics Re
- Page 39 and 40: S. Alexander, N. Yunes / Physics Re
- Page 41 and 42: S. Alexander, N. Yunes / Physics Re
- Page 43 and 44: S. Alexander, N. Yunes / Physics Re
- Page 45 and 46: S. Alexander, N. Yunes / Physics Re
- Page 47 and 48: S. Alexander, N. Yunes / Physics Re
- Page 49 and 50: S. Alexander, N. Yunes / Physics Re
- Page 51 and 52: S. Alexander, N. Yunes / Physics Re
- Page 53 and 54: S. Alexander, N. Yunes / Physics Re
- Page 55: S. Alexander, N. Yunes / Physics Re