Gravitational Waves Theory, Sources and Detection
Kostas Kokkotas
September 28, 2016
Kostas Kokkotas Gravitational Waves Suggested Reading
I Books
I Gravitational Waves: Volume 1: Theory and Experiments by Michele Maggiore Oxford University Press (2007)
I Gravitation and Spacetime by Ohanian, Hans C. and Ruffini, Remo Cambridge University Preess (Sep 20, 2013)
I Gravitation by Charles W. Misner, Kip S. Thorne and John Archibald Wheeler (Sep 15, 1973) W.H. Freeman
I Review articles
I Gravitational wave astronomy F.F. Schutz, Class. Quantum Grav. 16 (1999) A131– A156
I Gravitational wave astronomy: in anticipation of first sources to be detected L P Grishchuk, V M Lipunov, K A Postnov, M E Prokhorov, B S Sathyaprakash, Physics- Uspekhi 44 (1) 1 -51 (2001)
I The basics of gravitational wave theory E.E. Flanagan and S.A. Hughes, New Journal of Physics 7 (2005) 204
I Gravitational Wave Astronomy K.D. Kokkotas Reviews in Modern Astrophysics, Vol 20, ”Cosmic Matter”, WILEY-VCH, Ed.S. Roeser (2008) arXiv:0809.1602 [astro-ph]
Kostas Kokkotas Gravitational Waves Tensors I
Tensor Transformations X ∂x ν X ∂x˜µ b˜ = b and˜ aµ = aν (1) µ ∂x˜µ ν ∂x ν ν ν
∂x˜α ∂x˜β ∂x˜α ∂x ν ∂x µ ∂x ν T˜ αβ = T µν , T˜ α = T µ & T˜ = T ∂x µ ∂x ν β ∂x µ ∂x˜β ν αβ ∂x˜α ∂x˜β µν Covariant Derivative
φ;λ = φ,λ (2) ρ Aλ;µ = Aλ,µ − ΓµλAρ (3) λµ λµ λ αµ µ λα T ;ν = T ,ν + Γαν T + Γαν T (4)
Parallel Transport
λ ν δaµ = Γ µν aλdx for covariant vectors (5) µ µ λ ν δa = −Γ λν a dx for contravariant vectors (6)
Kostas Kokkotas Gravitational Waves Tensors II
The distance ds of two points P(x µ) and P0(x µ + dx µ) is given by 2 2 2 ds2 = dx 1 + dx 2 + dx 3 (7) in two different coordinate frames , 2 µ ν α β ds =g ˜µν dx˜ dx˜ = gαβ dx dx . (8)
I Metric element for Minkowski spacetime ds2 = −c2dt2 + dx 2 + dy 2 + dz2 (9) ds2 = −c2dt2 + dr 2 + r 2dθ2 + r 2 sin2 θdφ2 (10)
I For a sphere with radius R : ds2 = R2 dθ2 + sin2 θdφ2 (11)
I The metric element of a torus with radii a and b ds2 = a2dφ2 + (b + a sin φ)2 dθ2 (12)
I The Schwarzschild metric: 2 GM 2 GM −1 ds2 = 1 − c2dt2 − 1 − dr 2 − r 2 dθ2 + sin2 θdφ2 c2 r c2 r
Kostas Kokkotas Gravitational Waves Tensors III
Curvature Tensor or Riemann Tensor
λ λ λ µ λ µ λ R βνσ = −Γ βν,σ + Γ βσ,ν − Γ βν Γ µσ + Γ βσΓ µν (13)
Christoffel Symbols
1 Γα = g αν (g + g − g ) (14) µρ 2 µν,ρ νρ,µ ρµ,ν
Geodesic Equations
duρ d 2x ρ dx µ dx ν + Γρ uµuν = 0 or + Γρ = 0 ds µν ds2 µν ds ds
Kostas Kokkotas Gravitational Waves The Ricci and Einstein Tensors
The contraction of the Riemann tensor leads to Ricci Tensor
λ λµ Rαβ = R αλβ = g Rλαµβ µ µ µ ν µ ν =Γ αβ,µ − Γαµ,β + Γαβ Γνµ − Γαν Γβµ (15)
which is symmetric Rαβ = Rβα. Further contraction leads to the Ricci or Curvature Scalar
α αβ αβ µν R = R α = g Rαβ = g g Rµανβ . (16) The following combination of Riemann and Ricci tensors is called Einstein Tensor 1 G = R − g R (17) µν µν 2 µν with the very important property:
µ µ 1 µ G ν;µ = R ν − δ ν R = 0 . (18) 2 ;µ
This results from the Bianchi Identity (how?)
λ R µ[νρ;σ] = 0 (19)
Kostas Kokkotas Gravitational Waves Flat & Empty Spacetimes
I When Rαβµν = 0 the spacetime is flat
I When Rµν = 0 the spacetime is empty
Prove that : λ λ λ κ a ;µ;ν − a ;ν;µ = −R κµν a
Kostas Kokkotas Gravitational Waves Einstein’s equations
Thus the typical form of Einstein’s equations is: 1 Rµν − g µν R + Λg µν = κT µν . (20) 2
8πG whereΛ= c2 ρv is the so called cosmological constant. They can also be written as: 1 R = −κ T − g T (21) µν µν 2 µν
Kostas Kokkotas Gravitational Waves GW: Linear Theory I
Weak gravitational fields can be represented by a slightly deformed Minkowski spacetime : 2 gµν ' ηµν + hµν + O(hµν ) , |hµν | 1 (22)
here hµν is a small metric perturbation. The indices will be raised and lowered by ηµν i.e. αβ αµ βν h = η η hµν (23) µν h = η hµν (24) g µν = ηµν − hµν (25)
and we will define the traceless (φµν ) tensor: 1 φ = h − η h . (26) µν µν 2 µν The Christoffel symbols & the Ricci tensor will become : 1 Γλ ≈ ηλρ (h + h − h ) (27) µν 2 ρν,µ µρ,ν µν,ρ 1 R ≈ Γα − Γα ≈ (hα + hα − h α − hα )(28) µν µν,α µα,ν 2 ν,µα µ,να µν,α α,µν µν αβ α,β R = η Rµν ≈ hαβ − hα,β (29)
Kostas Kokkotas Gravitational Waves Finally, Einstein tensor gets the form: 1 G (1) = (hα + hα − h α − hα ) − η h ,αβ − hα,β (30) µν 2 ν,µα µ,να µν,α α,µν µν αβ α,β Einstein’s equations reduce to (how?):
,α ,αβ ,α ,α −φµν ,α − ηµν φαβ + φµα,ν + φνα,µ = κTµν (31) Then by using the so called Hilbert (or Harmonic or De Donder) gauge similar α ,α 1 to Lorenz gauge( A ,α = Aα = 0) in EM
µα ,α φ ,α = φµα = 0 (32) we come to the following equation:
1 ∂2 φ ,α ≡ φ ≡ − − ∇2 φ = −κT (33) µν ,α µν c2 ∂t2 µν µν which is a simple wave equations describing ripples of spacetime propagating with the speed of light (why?). These ripples are called gravitational waves.
1The De Donder gauge is defined in a curved background by the condition µν √ ∂µ(g −g) = 0 Kostas Kokkotas Gravitational Waves GW: about Gauge conditions (*)
By careful choice of coordinates the linearized Einstein equations can be simplified. We can fix ηµν = diag(−1, 1, 1, 1) and make small changes in the coordinates that leave ηµν unchanged but induce small changes in hµν . For example lets consider a change of the form:
x 0µ = x µ + ξµ (34)
where ξµ are 4 small arbitrary functions of the same order as hµν . Then
∂x 0µ ∂x µ = δµ + ∂ ξµ and = δµ − ∂ ξµ ∂x ν ν ν ∂x 0ν ν ν Thus, the metric transforms as: ∂x ρ ∂x σ g 0 = g = δρ − ∂ ξρ (δσ − ∂ ξσ)(η + h ) µν ∂x 0µ ∂x 0ν ρσ µ µ ν ν ρσ ρσ 0 ≈ ηµν + hµν − ∂µξν − ∂ν ξµ = ηµν + hµν (35)
Then in the new coordinate system we get
0 hµν = hµν − ξµ,ν − ξν,µ (36)
This transformation is called gauge transformation.
Kostas Kokkotas Gravitational Waves GW: about Gauge conditions II (*)
This is analogous to the gauge transformation in Electromagnetism. If Aµ is a solution of the EM field equations then another solution that describes precisely the same physical situation is given by
(new) Aµ = Aµ − ψ,µ (37)
where ψ is any scalar field. (new)µ (new) ,µ Then the gauge condition A ,µ = A µ = 0 means that ,µ ,µ µ ψ ,µ = ψ,µ =ψ = A ,µ.
From (36) it is clear that if hµν is a solution to the linearised field equations then the same physical situation is also described by
(new) φµν = φµν −ξµ,ν − ξν,µ = φµν −Ξµν (38)
NOTE • This is a gauge transformation and not a coordinate one • We are still working on the same set of coordinates x µ and have defined a (new) new tensor φµν whose components in this basis are given by (38).
Kostas Kokkotas Gravitational Waves GW: about Gauge conditions III (*)
We can easily see that from (36) or (38) we can get
(new) µρ µρ µ φ ,ρ = φ ,ρ − ξ (39) Therefore, if we choose the function ξµ so that to satisfy
µ µρ ξ = φ ,ρ (40) we get the Hilbert gauge (new) µρ φ ,ρ = 0 (41) NOTE: This gauge condition is preserved by any further gauge transformation µ µ of the form (38) provided that the functions ξ satisfy ξ = 0 or equivalently µν Ξ = 0.
Kostas Kokkotas Gravitational Waves µν I The choice of the Hilbert gauge φ ,ν = 0, gives 4 conditions that reduces the 10 independent components of the symmetric tensor hµν to 6!
I Eqn (38) tells us that, from the 6 independent components of φµν which satisfy φµν = 0, we can subtract the functionsΞ µν , which depend on 4 independent arbitrary functions ξµ satisfying the same equation Ξµν = 0.
I This means that we can choose the functions ξµ so that as to impose 4 conditions on φµν . 0 I We can choose ξ such that the trace φ = 0. Note that if φ = 0 then φµν = hµν . 0i I The 3 functions ξi can be chosen so that φ = 0. ,0 ,i I Then the Hilbert condition for µ = 0 will be written φ00 + φ0i = 0. 0i ,0 But since we fixed φ = 0 we get φ00 = 0, i.e. φ00 is a constant in time.
I A time-independent part term φ00 corresponds to the static part of the grav. interactions i.e. to the Newtonian potential of the source.
I The GW itself is the time-dependent part and therefore as far as ,0 the GW concerns h00 = 0 means h00 = 0. In conclusion, we set
0µ i ij h = 0 , h i = 0 , h ,i = 0 (42)
Kostas Kokkotas Gravitational Waves GW: Properties
Equation (33) is the basis for computing the generation of GWs within the linearised theory. To study the propagation of GWs as well as the interaction with test masses (and therefore the GW detector) we are interested for the equations outside the source, i.e. where Tµν = 0. GWs are periodic changes of spacetime curvature and for weak gravitational fields far away from sources they described by a simple wave equations which admits a solution of the form:
α φµν = Aµν cos (kαx ) , (43)
where Aµν is a symmetric tensor called polarization tensor including information of the amplitude and the polarization properties of the GWs. kα ≡ (k0 = ω/c,~k) is the wave-vector. This solution satisfies Hilbert’s gauge condition, that is:
,ν ν α 0 = φµν = −Aµν k sin (kαx )
which lead to the orthogonality condition
ν Aµν k = 0 . (44)
Kostas Kokkotas Gravitational Waves From the wave equation (33) we get
,α α α α 0 = φµν ,α = −Aµν k kα cos (kαx ) ⇒ k kα = 0. (45) This relation suggests that the wave vector kα is null i.e. gravitational waves are propagating with the speed of light. But, (45) implies that ω2 = c2|~k|2 i.e. both group and phase velocity of GWs are equal to the speed of light. ∂ω λ ω vgroup = ∂k and vphase = T = k
Kostas Kokkotas Gravitational Waves GW: The Transverse - Traceless (TT) Gauge
Based on the gauge freedom which allows to choose ξµ we derived the following relations 0µ i ij h = 0 , h i = 0 , h ,i = 0 (46) which define the so-called Transverse - Traceless (TT) Gauge. Then for a GW propagating in the z direction i.e. it has a wave vector of the form kµ = (ω/c, 0, 0, −ω/c) where k0 = ω/c is the frequency of the wave that:
0 0 0 0 0 h h 0 hµν ≡ + × cos[ω(t − z/c)] (47) 0 h× −h+ 0 0 0 0 0
While h+ and h×, are the amplitudes of the gravitational waves in the two polarizations. The GWs described in this spacific gauge are Transverse and Traceless, and we TT will use the notation hµν .
Kostas Kokkotas Gravitational Waves GW: Effects...
We will study the effect of GWs on particles. A static or slowly moving particle has velocity vector uµ ≈ (1, 0, 0, 0) and one can assume that τ ≈ t. Then in linearized gravity the geodesic equation will be written as: duµ 1 = − (h + h − h ) uαuβ (48) dt 2 µα,β βµ,α αβ,µ leading to duµ 1 = − h − h . (49) dt µ0,0 2 00,µ
If we now use the T-T gauge (h00 = hµ0 = 0) we conclude that GWs do not affect isolated particles! If instead we consider a pair of test particles on the cartesian axis Ox being at distances x0 and −x0 from the origin and we assume a GW traveling in the z-direction then their distance will be given by the relation:
2 µ ν d` = gµν dx dx = ... 2 2 2 = −g11(dx) = (1 − h11)(2x0) = (1 − h+ cos ωt) (2x0) (50)
or approximatelly 1 ∆` ≈ 1 − h cos ωt (2x ) . (51) 2 + 0
Kostas Kokkotas Gravitational Waves GW: Effects...
In a similar way we can show for two particles on the Oy axis that: 1 ∆` ≈ 1 + h cos ωt (2y ) . (52) 2 + 0 In other wards the coordinate distance of two particles is varying periodically with the time
Figure: The effect of a travelling GW on a ring of particles
Kostas Kokkotas Gravitational Waves GW: Effects...
Figure: The effect of a travelling GW on a ring of particles
Kostas Kokkotas Gravitational Waves Geodesic deviation (*)
In a curved spacetime two geodesics that can be “parallel” initially will either converge or diverge depending on the local curvature. Consider two neighbouring geodesics G given by x α(σ) and G˜ given by˜x α(σ) where σ is an affine parameter. If ξα(σ) is a small vector connecting points of the two geodesics for the same values of σ i.e. x˜α(σ) = x α(σ) + ξα(σ) If we construct local geodesic coordinates about the point P, the Christoffel symbols will vanish but its derivatives will be non-zero there.
Figure: (Left) Two neighbouring geodesics. (Right) Converging geodesics on the surface of a sphere. Kostas Kokkotas Gravitational Waves In this coordinate system we will get
d 2x α d 2x˜α dx˜µ dx˜ν = 0 , + Γα = 0 (53) dσ2 dσ2 µν dσ dσ P Q But since ξα is small:
α α α λ α λ [Γ µν ]Q = [Γ µν ]P + [Γ µν,λ]P ξ =[Γ µν,λ]P ξ by subtracting the two equations in (53) we get (to 1st order, at P):
d 2ξα dx µ dx ν + Γα ξλ = 0 dσ2 µν,λ dσ dσ However, in our geodesic coordinates the 2nd order absolute (intrinsic) derivative of ξα at P is:
D2ξα d dξα d 2ξα dx µ dx λ = + Γα ξµx ν = + Γα ξν Dσ2 dσ dσ µν dσ2 µν,λ dσ dσ α where we have used the fact thatΓ µν (P) = 0. By combining the last two equations we get:
D2ξα dx µ dx λ +[Γ α − Γα ]ξν = 0 Dσ2 µλ,ν µν,λ dσ dσ which will give
Kostas Kokkotas Gravitational Waves D2ξα dx µ dx ν + Rα ξν = 0 (54) Dσ2 µνλ dσ dσ because the term in the square brackets is the Riemann tensor in local geodesic coordinates.
Kostas Kokkotas Gravitational Waves Tidal forces in a curved spacetime
Tidal forces deform the shape of bodies as they freely move in a gravitational field. Thus two nearby particles with trajectories x i (t) and˜x i (t) (in Cartesian coordinates) will be separated by a vector ξi = x i − x˜i d 2ξ ∂2Φ = − ξj (55) dt2 ∂x i ∂x j (why?) whereΦ is the Newtonian gravitational potential.
Figure: Tidal effects on a cloud of particles
Kostas Kokkotas Gravitational Waves Tidal effects can be also estimated in GR for two particles moving along timelike geodesics x µ(τ) and˜x µ(τ)( τ is the proper time of the 1st particle). The separation vector between the worldlines of the 2 particles is ξµ(τ) =x ˜µ − x µ: D2ξµ = Rµ uσuρξν ≡ Sµ ξν (56) Dτ 2 σρν ν µ σ σ where S ν is the so called tidal stress tensor and u = du /dτ. This is a fully covariant tensor equation and holds in any coordinate system.
Kostas Kokkotas Gravitational Waves Figure: The basis vectors of the instantaneous rest frame (IRF) at P. • ˆeα is a set of orthonormal basis vectors at P that define the IRF of the first particle (observer) withˆeα · ˆeβ = ηαβ . αˆ α α µ • ξ is a general connecting vector with ξ ≡ ˆe · ξ = (ˆe )µ ξ • ζ is the orthogonal connecting vector. For an observer sitting on the one of the particles it can be shown that in any orthonormal freely falling frame becomes: d 2ξµˆ = c2Rµˆ ξνˆ . (57) dτ 2 0ˆ0ˆˆν Newtonian limit (we will discuss the details later) ∂2Φ Rµˆ → (58) 0ˆ0ˆˆν ∂x i ∂x j
Kostas Kokkotas Gravitational Waves GW: Tidal forces
Riemann tensor is a ”measure” of spacetime’s curvature and in linearized gravity gets the form 1 R = (∂ h + ∂ h − ∂ h − ∂ h ) , (59) κλµν 2 νκ λµ λµ κν κµ λν λν κµ in the T-T gauge the Riemann tensor is considerably simplified
1 ∂2 RTT = − hTT, for j, k = 1, 2, 3. (60) j0k0 2 ∂t2 jk Actually, the Newtonian limit of the Riemann tensor is:
∂2U RTT ≈ , (61) j0k0 ∂x j ∂x k where U is the Newtonian potential. In other words the Riemann tensor has also a pure physical meaning i.e. it is a measure of the tidal gravitational acceleration. Then the distance between two nearby particles x µ(τ) will x µ(τ) + ξµ(τ) will be described by
2 k d ξ TT ≈ −Rk ξj . (62) dt2 0j0
Kostas Kokkotas Gravitational Waves The tidal force acting on a particle is (why?)
m d 2hTT f k ≈ −mRk ξj ≈ jk ξj (63) 0j0 2 dt2 where m is particle’s mass. This means that m m f x ≈ h ω2 cos[ω(t − z)]ξx , and f y ≈ − h ω2 cos[ω(t − z)]ξy . (64) 2 + 0 2 + 0 x y ~ ∂f ∂f ∇f = x + y = 0 . (65) ∂ξ0 ∂ξ0 Hence the divergence of the force ~f is zero, which tell us that the tidal force can be represented graphically by field lines.
Figure: The tidal field lines of force for a gravitational wave with polarization (+) (left panel) and (×) (right panel). The orientation of the field lines changes every half period producing the deformations as seen in Figure 1. Any point accelerates in the directions of the arrows, and the denser are the lines, the strongest is the acceleration. Since the acceleration is proportional to the distance from the center of mass, the force lines get denser as one moves away from the origin. For the polarization (×) the force lines undergo a 450 rotation.
Kostas Kokkotas Gravitational Waves GW: Properties
• GWs, once they are generated, propagate almost unimpeded. Indeed, they are even harder to stop than neutrinos! The only significant change they suffer as they propagate is the decrease in amplitude while they travel away from their source, and the redshift they feel (cosmological, gravitational or Doppler). • EM waves are fundamentally different, however, even though they share similar wave properties away from the source. • GWs are emitted by coherent bulk motions of matter (for example, by the implosion of the core of a star during a supernova explosion) or by coherent oscillations of spacetime curvature, and thus they serve as a probe of such phenomena. By contrast, • Cosmic EM waves are mainly the result of incoherent radiation by individual atoms or charged particles.
? As a consequence, from the cosmic electromagnetic radiation we mainly learn about the form of matter in various regions of the universe, especially about its temperature and density, or about the existence of magnetic fields.
Kostas Kokkotas Gravitational Waves GW: Properties
• Strong GWs are emitted from regions of spacetime where gravity is very strong and the velocities of the bulk motions of matter are near the speed of light. Since most of the time these areas are either surrounded by thick layers of matter that absorb EM radiation or they do not emit any at all (black holes), the only way to study these regions of the universe is via GWs.
Kostas Kokkotas Gravitational Waves GW: The energy of GWs
The fact that GWs carry energy and momentum is already clear from the discussion on the interaction with test masses. To get an explicit expression of the energy-momentum tensor of GWs we can follow two different routes one more geometrical and the other more field-theoretical A. According to GR, any form of energy contributes to the curvature of space-time, thus we can ask ”whether GWs are themselves a source of space-time curvature”. B. We can treat linearised gravity as any other classical filed theory and apply Noether’s theorem, the standard field-theoretical tool that answers this question.
Kostas Kokkotas Gravitational Waves GW: The energy of GWs
In order to include the contribution of the energy-momentum associated with the gravitational field itself one must modify the linearise Einstein’s equations to read 8πG G (1) = − (T + T ) (66) µν c4 µν µν where Tµν is the energy-momentum tensor of any matter present and Tµν is the energy-momentum tensor of the gravitational field itself. On the other hand Einstein’s equations may expand beyond first order to obtain 8πG G ≡ G (1) + G (2) + ··· = − T (67) µν µν µν c4 µν This suggest that, to a good approximation, we should make the identification
c4 c4 T ≡ G (2) = ... = G (2) (68) µν 8πG µν 8πG µν
Kostas Kokkotas Gravitational Waves GW: Energy
GWs carry energy. The stress-energy carried by GWs cannot be localized within a wavelength. Instead, one can say that a certain amount of stress-energy is contained in a region of the space which extends over several wavelengths. The stress-energy tensor can be written as: 1 h 1 i T µν = 2φαβ,µφ ,ν − φ,µφ,ν − ηµν φαβ,σφ − φ φ,σ (69) 4 αβ αβ,σ 2 ,σ which in the TT gauge of the linearized theory becomes (HOW?)
c4 T GW = ∂ hTT ∂ hTT . (70) µν 32πG µ ij ν ij where the angular brackets indicate averaging over several wavelengths. For the special case of a plane wave propagating in the z direction, the stress-energy tensor has only three non-zero components, which take the simple form
T GW T GW 1 c2 T GW = zz = − 0z = ω2 h2 + h2 , (71) 00 c2 c 32π G + ×
GW GW GW where T00 is the energy density, Tzz is the momentum flux and T0z the energy flow along the z direction per unit area and unit time.
Kostas Kokkotas Gravitational Waves Gravitational Waves: Nature of
? Electromagnetic radiation emitted by slowly varying charge distributions can be decomposed into a series of multipoles, where the amplitude of the 2`-pole (` = 0, 1, 2, ...) contains a small factor a`, with a equal to the ratio of the diameter of the source to the typical wavelength, namely, a number typically much smaller than 1.
From this point of view the strongest EM radiation would be expected for monopolar radiation( ` = 0), but this is completely absent, because the EM monopole moment is proportional to the total charge, which does not change with time (it is a conserved quantity).
Therefore, EM radiation consists only of ` ≥ 1 multipoles, the strongest being the electric dipole radiation (` = 1), followed by the weaker magnetic dipole & electric quadrupole radiation (` = 2).
Kostas Kokkotas Gravitational Waves Gravitational Waves: Nature of
For Gravitational Waves, it can be shown that mass conservation (which is equivalent to charge conservation in EM theory) will exclude monopole radiation.
Also, the rate of change of the mass dipole moment is proportional to the linear momentum of the system, which is a conserved quantity, and therefore there cannot be any mass dipole radiation in Einstein’s relativity theory.
The next strongest form of EM radiation is the magnetic dipole. For the case of gravity, the change of the “magnetic dipole” is proportional to the angular momentum of the system, which is also a conserved quantity and thus there is no dipolar grav. radiation of any sort. It follows that grav. radiation is of quadrupolar (or higher nature) and is directly linked to the quadrupole moment of the mass distribution.
Kostas Kokkotas Gravitational Waves Newtonian Gravity
Poisson equation Z ρ(~x) ∇2U(~x) = 4πGρ(~x) → U(~x) = −G d 3x~0 |~x − x~0| For a spherically symmetric mass distribution of radius R 1 Z R U(r) = − r 02ρ(r 0)dr 0 for r > R r 0 1 Z r Z R U(r) = − r 02ρ(r 0)dr 0 − r 0ρ(r 0)dr 0 for r < R r 0 r
Kostas Kokkotas Gravitational Waves For a non-spherical distribution the term1 /|~x − x~0| can be expanded as
k 0k k l 1 1 X x x 1 X X 0k 0l 02 l x x = + + 3x x −~r δk + ... |~x − x~0| r r 3 2 r 5 k k l
GM G X G X x k x l U(~x) = − − x k Dk − Qkl + ... r r 3 2 r 5 k kl Gravitational Multipoles
Z M = ρ(x~0)d 3x 0 Mass Z Dk = x 0k ρ(x~0)d 3x 0 Mass Dipole moment2
Z 2 kl 0k 0l ~0 l ~0 3 0 3 Q = 3x x − r δk ρ(x )d x Mass Quadrupole tensor
2If the center of mas is chosen to coincide with the origin of the coordinates then Dk = 0 (no mass dipole). 3If Qkl 6= 0 the potential will contain a term proportional to ∼ 1/r 3 and the gravitational force will deviate from the inverse square law by a term ∼ 1/r 4. Kostas Kokkotas Gravitational Waves GW: Generation
Einstein (1918) derived the quadrupole formula for gravitational radiation by solving the linearized form of his equations µν µν φ (t,~x) = −κT (t,~x) . (72) The solution is: κ Z T µν (t − |~x − ~x 0|,~x 0) µν ~ 3 0 φ (t, x) = − 0 d x , (73) 4π V |~x − ~x |
This solution suggests that φij is proportional to the second time derivative of the quadrupole moment of the source (WHY?): 2 G Z 1 φ = Q¨ TT (t − r/c) where QTT (x) = ρ x i x j − δij r 2 d 3x (74) ij r c4 ij ij 3
TT where, Qij is the quadrupole moment in the TT gauge, evaluated at the retarded time t − r/c 4. The energy radiated by the system per unit solid angle and unit time in the direction ns is d 2E − = r 2T 0s ns (75) dtdΩ 4This result is quite accurate for all sources, as long as the reduced wavelength λ˜ = λ/2π is much longer than the source size R. Kostas Kokkotas Gravitational Waves GW: Emission of Energy, Angular and Linear Momentum
Using the formulae (70) and (71) for the energy carried by GWs, one can derive the luminosity in GWs as a function of the third-order time derivative of the quadrupole moment tensor. This is the well-known quadrupole formula for the Energy emission dE 1 G ...... L = − = Q · Q (76) GW dt 5 c5 ij ij Angular momentum emission
dJ GW 2 X ... i = Q¨ · Q (77) dt 5 ijk j` `k jk`
Linear momentum emission GW ...... dPi 2 X 16 X = Q · Q + ijk Q · P `k (78) dt 63 jk jki 45 j` jk jk`
where Qijk : mass octupole moment Pij : current quadrupole moment
Kostas Kokkotas Gravitational Waves GW: Energy Flux
The energy flux has all the properties one would anticipate by analogy with electromagnetic waves: (a) it is conserved (the amplitude dies out as1 /r, the flux as1 /r 2), (b) it can be absorbed by detectors, and (c) it can generate curvature like any other energy source in Einstein’s formulation of relativity. Estimate the energy flux in GWs from the collapse of the core of a supernova to create a 10 M black hole at a distance of ∼15 Mpc from the earth (at the distance of the Virgo cluster of galaxies). An optimistic estimate of the amplitude of the GWs on Earth is of the order ofh ≈ 10−22 (at a frequency of about 1kHz). This corresponds to a flux of about 3 ergs/cm2 sec. This is an enormous amount of energy flux and is about ten orders of magnitude larger than the observed energy flux in electromagnetic waves! The basic difference is the duration of the two signals; GW signal will last a few milliseconds, whereas an EM signal lasts many days. This example provides us with a useful numerical formula for the energy flux:
f 2 h 2 ergs F = 3 , (79) 1kHz 10−22 cm2sec from which one can easily estimate the flux on Earth, given the amplitude (on Earth) and the frequency of the waves.
Kostas Kokkotas Gravitational Waves GW: Order of magnitude estimates
The quadrupole moment of a system is approximately equal to the mass M of the part of the system that moves, times the square of the size R of the system. This means that the 3rd-order time derivative of the quadrupole moment is ∂3Q MR2 MV2 E ij ∼ ∼ ∼ ns , (80) ∂t3 T 3 T T
where V is the mean velocity of the moving parts, Ens is the kinetic energy of the component of the source’s internal motion which is non-spherical, and T is the time scale for a mass to move from one side of the system to the other. The time scale (or period) is actually proportional to the inverse of the square root of the mean density of the system (why?) T ∼ pR3/GM. (81) This relation provides a rough estimate of the characteristic frequency of the system f = 2π/T . The luminosity of GWs of a given source is approximately G 4 M 5 G M 2 c5 R 2 V 6 L ∼ ∼ V6 ∼ Sch (82) GW c5 R c5 R G R c 2 where RSch = 2GM/c is the Schwarzschild radius of the source. It is obvious that the maximum value of the luminosity in GWs can be achieved if the source’s dimensions are of the order of its Schwarzschild radius and the typical velocities of the components of the system are of the order of the speed of light. Kostas Kokkotas Gravitational Waves GW: Order of magnitude estimates II
The above formula sets also an upper limit on the power emitted by a source, which for R ∼ RSch and V ∼ c is: 5 59 LGW ∼ c /G = 3.6 × 10 ergs/sec. (83) This is an immense amount of power, often called the luminosity of the universe. Using the above order-of-magnitude estimates, we can get a rough estimate of the amplitude of GWs at a distance r from the source: G E G εE h ∼ ns ∼ kin (84) c4 r c4 r
where εEkin (with 0 ≤ ε ≤ 1), is the fraction of kinetic energy of the source that is able to produce GWs. The factor ε is a measure of the asymmetry of the source and implies that only a time varying quadrupole moment will emit GWs. Another formula for the amplitude of GW relation can be derived from the flux formula (79). If, for example, we consider an event (perhaps a supernovae −4 explosion) at the Virgo cluster during which the energy equivalent of 10 M is released in GWs at a frequency of 1 kHz, and with signal duration of the order of 1 msec, the amplitude of the gravitational waves on Earth will be 1/2 −1 −1 −1/2 −22 EGW f τ r h ≈ 10 −4 . (85) 10 M 1kHz 1msec 15Mpc Kostas Kokkotas Gravitational Waves GW: Order of magnitude estimates III
For a detector with arm length of 4 km we are looking for changes in the arm length of the order of
∆` = h · ` = 10−22 · 4 km = 4 × 10−17cm!!!
These numbers shows why experimenters are trying so hard to build ultra-sensitive detectors and explains why all detection efforts till recently were not successful. Finally, it is useful to know the damping time, that is, the time it takes for a source to transform a fraction1 /e of its energy into gravitational radiation. One can obtain a rough estimate from the following formula
E 1 R 3 τ = kin ∼ R . (86) LGW c RSch For example, for a non-radially oscillating neutron star with a mass of roughly 1.4M and a radius of 12Km, the damping time will be of the order of ∼50msec. Also, by using formula (81), we get an estimate for the frequency of oscillation which is directly related to the frequency of the emitted gravitational waves, roughly 2kHz for the above case.
Kostas Kokkotas Gravitational Waves Example: Quadrupole Moment Tensor
We will calculate the mass quadrupole moment tensor of a homogeneous x2 y2 z2 0 0 0 triaxial ellipsoid a2 + b2 + c2 = 1. By setting x = x/a, y = y/a and z = z/a the volume integration over the ellipsoid reduces to that over the unit sphere ZZZ ZZZ 2 2 2 2 2 Q11 = ρ 3x − r dxdydz = ρ 2x − y − z dxdydz ZZZ = ρabc 2a2x 02 − b2y 02 − c2z02 dx 0dy 0dz0 ZZZ = ρabc 2a2 − b2 − c2 z02dx 0dy 0dz0
Z 2π Z π Z 1 = ρabc 2a2 − b2 − c2 r 4dr cos2 θ sin θdθdφ 0 0 0 m = 2a2 − b2 − c2 (87) 5 4 where m = 3 πabcρ is the mass of the ellipsoid. The other two non-vanishing components of the mass quadrupole tensor are 5: m m Q = −a2 + 2b2 − c2 and Q = −a2 − b2 + 2c2 22 5 33 5 5 By definition, the mass quad. moment tensor is traceless, Qjj = Q11 + Q22 + Q33. Kostas Kokkotas Gravitational Waves Example: Vibrating Quadrupole I
The distance of the masses from the center varies periodically as z = ±(b + a sin ωt). The quadrupole moment tensor for the pair of equal masses m is:
2 2 ! − 3 mb 0 0 ij 2 2 Q(0) ≡ 0 − 3 mb 0 (88) 4 2 0 0 3 mb Then the retarded value of the quadrupole tensor is: h 2a i Qij (t − r) ≈ 1 + sin ω(t − r) Q(0)ij (89) b The radiated gravitational field is: 2 G 2 G 2a φij = Q¨ ij (t − r) = ω2 sin ω(t − r)Q(0)ij (90) r c4 r c4 b
Kostas Kokkotas Gravitational Waves Example: Vibrating Quadrupole II
The energy radiated by the system per unit solid angle and unit time in the direction ns is d 2E − = r 2T 0s ns dtdΩ 1 κ 2 κ 2 2 = − Q = 2mabω3 cos ω(t − r) sin4 θ (91) 18 8π 8π
2 2 2 2 2 2 ...11 ...22 ...33 ...11 ...22 ...33 Q = Q + Q + Q − 2 Q n1 − 2 Q n2 − 2 Q n3
1 ...11 ...22 ...33 + Q n1n1 + Q n2n2 + Q n3n3 2
Kostas Kokkotas Gravitational Waves
Figure: The radiation pattern of emission of gravitational radiation by a quadrupole oscillating along the z-axis. Example: Vibrating Quadrupole III
The total emitted power is:
dE 1 G ...... 32 G 2 L = − = Q · Q = mabω3 cos ω(t − r) GW dt 5 c5 ij ij 15 c5 16G ≈ (mab)2 ω6 (92) 15c5 and the damping time of the oscillator, due to the emission of GWs is :
1 1 dE 16 G 2 4 1 2 2 = − = 5 mb ω where E = mω a (93) τrad E dt 15 c 2 The above formulae give an order of magnitude estimate for the GW emission of by any vibrating elastic body, provided that the vibrations are not spherical.
Kostas Kokkotas Gravitational Waves Example: Two-body collision I
We assume that a particle of mass m starts from infinity with zero velocity 1 2 GmM 2 ... 3/2 7/2 ( 2 mz˙ = |z| ,z ¨ = GM/z , z = (2GM) /|z| ) falls towards a massive body of mass M .
Radiated power dE 8 G ... − = m2 (3˙zz¨ + z z )2 (94) dt 15 c5 The energy during the plunge from z = ∞ to z = R
4 G m2(2GM)5/2 −∆E = (95) 105 c5 R7/2
If R = RSchw (M = 10M & m = 1M ) m −∆E =0 .019mc2 (96) M m −∆E =0 .0104mc2 → 2 × 1051erg (97) M Most radiation during the2 R → R phase ∆t ∼ R/ν ∼ R/c ∼ 30km/c ∼ 10−4sec → f ∼ 104Hz (98)
Kostas Kokkotas Gravitational Waves Example: Two-body collision III
Figure: Spectrum of GW emitted by a particle of mass m falling radially into a BH of mass M. The quantity dE/dω gives the amount of energy radiated Figure: The signal of a ringing per unit frequency interval. The curves black-hole. The signal can be produced marked L = 2, 3, 4 correspond to by a small body falling into a quadrupole, ... radiation. Note that black-hole. most of the radiation is emitted with frequency ω ∼ 0.3 − 0.5c3/GM.
Kostas Kokkotas Gravitational Waves GW: Binaries, an example I
If we assume that the two bodies m1 and m2 making up the binary lie in the x − y plane at distances a1 and a2 from the center of mass, their orbits are circular and rotating at angular frequencyΩ.
Then the only non-vanishing components of the quadrupole tensor are (why?): 1 Q = −Q = a2M + a2M cos2 Ωt = µa2 cos 2Ωt, (99) xx yy 1 1 2 2 2 1 Q = Q = µa2 sin 2Ωt, (100) xy yx 2
where a = a1 + a2, a1M1 = a2M2 = aµ. Here µ = M1M2/M is the reduced mass of the system and M = M1 + M2 its total mass.
Kostas Kokkotas Gravitational Waves GW: Binaries, an example II
The GW luminosity of the system is (we use Kepler’s third law,Ω 2 = GM/a3) (how?)
dE 1 G 1 2 LGW = − = (2Ω)2 a2µ sin2 2Ωt + sin2 2Ωt + 2 cos 2Ωt dt 5 c5 2 32 G 32 G 4 M3µ2 = µ2a4Ω6 = . (101) 5 c5 5 c5 a5 The total energy of the binary system can be written as (why?):
1 1 GM M 1 GµM E = M a2 + M a2 Ω2 − 1 2 = − (102) 2 1 1 2 2 2 a 2 a
Kostas Kokkotas Gravitational Waves GW: Binaries, an example III
As the gravitating system loses energy by emitting radiation, the distance between the two bodies shrinks at a rate dE 1 GµM da da 64 G 3 µM2 = ⇒ = − , (103) dt 2 a2 dt dt 5 c5 a3 and the orbital frequency increases accordingly (f˙/f = (3/2)˙a/a). If, the present separation of the two stars is a0, then the binary system will coalesce after a time 5 c5 a4 τ = 0 (104) 256 G 3 µM4 Finally, the amplitude of the GWs is (why?)
2/3 M µ f 2/3 15Mpc h = 5 × 10−22 . (105) 2.8M 0.7M 100Hz r
NOTE: if m1 = m2 = 30M then the amplitude will be higher by about 160 times from a binary with m1 = m2 = 1.4M ! In all these formulae we have assumed that the orbits are circular. In general, the orbits of the two bodies are approximately ellipses, but it has been shown that long before the coalescence of the two bodies, the orbits become circular, at least for long-lived binaries, due to gravitational radiation.
Kostas Kokkotas Gravitational Waves GW: Binaries, an example IV
The angular distribution of the radiated power, is given by
dP r 2c3 = h˙ 2 + h˙ (106) dΩ 16πG + x or dP 2Gµ2a2ω6 = g(θ) (107) dΩ πc5 2 1 + cos2 θ g(θ) = + cos2 θ (108) 2
Figure: The function g(θ) in polar coordinates
Kostas Kokkotas Gravitational Waves GW: Binaries, an example V
I The amplitude of the emitted GWs depends on the angle between the line of sight and the axis of angular momentum; formula (105) refers to an observer along the axis of the orbital angular momentum.
I The complete formula for the amplitude contains angular factors of order 1. The relative strength of the two polarizations depends on that angle as well.
I If 3 or more detectors observe the same signal it is possible to reconstruct the full waveform and deduce many details of the orbit of the binary system.
I As an example, we will provide some details of the well-studied pulsar PSR 1913+16 (the Hulse-Taylor pulsar), which is expected to coalesce after ∼ 3.5 × 108 years. The binary system is roughly 5kpc away from Earth, the masses of the two neutron stars are estimated to be ∼1.4M each, and the present period of the system is ∼7h and 45min. The predicted rate of period change is T˙ = −2.4 × 10−12sec/sec, while the corresponding observed value is in excellent agreement with the predictions, i.e., T˙ = (−2.30 ± 0.22) × 10−12sec/sec; finally the present amplitude of gravitational waves is of the order of h ∼ 10−23 at a frequency of ∼ 7 × 10−5Hz.
Kostas Kokkotas Gravitational Waves GW: Basic Formulae
For a gravitational wave binary with masses m1 and m2, in a circular orbit with gravitational wave frequency f = 2Ω =2 pGM/a3, then:
(m m )3/5 M = 1 2 ≡ µ3/5M2/5 chirp mass (109) c 1/5 (m1 + m2) 4G M G 2/3 h = c πf M scaling amplitude (110) o c5 D c3 c 3 8/3 ˙ 96 c f G f = 3 πf Mc chirp (111) 5 G Mc c The chirp indicates that as GWs are emitted, they carry energy away from the binary. The gravitational binding energy decreases, and the orbital frequency increases. The GW phase φ(t) evolves in time as
1 f˙ φ(t) = 2πf t + t2 + φ (112) 2 f 0
where F˙ is the chirp given above, and φ0 is the initial phase of the binary.
Kostas Kokkotas Gravitational Waves GW: Basic Formulae
A phenomenological form of the waveform then is given by