Func Eqns
Func Eqns
A.1 Consider a system of two stars with masses 𝑀1 , 𝑀2 , at locations 𝑟1⃗ , 𝑟2⃗ , res- 1.0pt
pectively, with respect to the center-of-mass of the system, that is,
d2 𝑟1⃗ 𝑟⃗
= −𝛼 𝑛1 , (2)
d𝑡2 𝑟1
where 𝑟1 = |𝑟1⃗ |, 𝑟2 = |𝑟2⃗ |. Find 𝑛 ∈ ℕ and 𝛼 = 𝛼(𝐺, 𝑀1 , 𝑀2 ), where 𝐺 is Newton's
constant [𝐺 ≃ 6.67 × 10−11 N m2 kg−2 ].
A.2 The total energy of the 2-mass system, in circular orbits, can be expressed as: 1.0pt
𝑀𝜇
𝐸 = 𝐴(𝜇, Ω, 𝐿) − 𝐺 , (3)
𝐿
where
𝑀1 𝑀2
𝜇≡ , 𝑀 ≡ 𝑀 1 + 𝑀2 , (4)
𝑀 1 + 𝑀2
are the reduced mass and total mass of the system, Ω is the angular velocity of
each mass and 𝐿 is the total separation 𝐿 = 𝑟1 + 𝑟2 . Obtain the explicit form of
the term 𝐴(𝜇, Ω, 𝐿).
quadrupole formula,
𝐺 3 3 d3 𝑄𝑖𝑗 d3 𝑄𝑖𝑗
𝒫= ∑ ∑ ( ) ( ). (5)
5𝑐5 𝑖=1 𝑗=1 d𝑡3 d𝑡3
Here, 𝑐 is the velocity of light 𝑐 ≃ 3 × 108 m/s. For a system of 2 pointlike particles orbiting on the 𝑥 − 𝑦
plane, 𝑄𝑖𝑗 is the following table (𝑖, 𝑗 label the row/column number)
2 2 2
𝑀𝐴 𝑀𝐴 𝑀𝐴 2
𝑄11 = ∑ (2𝑥2𝐴 − 𝑦𝐴
2
) , 𝑄22 = ∑ 2
(2𝑦𝐴 − 𝑥2𝐴 ) , 𝑄33 = − ∑ 2
(𝑥𝐴 + 𝑦𝐴 ), (6)
𝐴=1
3 𝐴=1
3 𝐴=1
3
2
𝑄12 = 𝑄21 = ∑ 𝑀𝐴 𝑥𝐴 𝑦𝐴 , (7)
𝐴=1
and 𝑄𝑖𝑗 = 0 for all other possibilities. Here, (𝑥𝐴 , 𝑦𝐴 ) is the position of mass A in the center-of-mass frame.
B.1 For the circular orbits described in A.2 the components of 𝑄𝑖𝑗 can be expressed 1.0pt
as a function of time 𝑡 as:
B.2 Compute the power 𝒫 emitted in gravitational waves for that system, and ob- 1.0pt
tain:
𝐺 2 4 6
𝒫=𝜉 𝜇 𝐿 Ω . (9)
𝑐5
What is the number 𝜉? [If you could not obtain 𝜉, use 𝜉 = 6.4 in the following.]
B.3 In the absence of GW emission the two masses will orbit on a fixed circular orbit 1.0pt
indefinitely. However, the emission of GWs causes the system to lose energy
and to slowly evolve towards smaller circular orbits. Obtain that the rate of
change dΩd𝑡
of the orbital angular velocity takes the form
3
dΩ Ω11
( ) = (3𝜉)3 15 (𝐺𝑀c )5 , (10)
d𝑡 𝑐
where 𝑀c is called the chirp mass. Obtain 𝑀c as a function of 𝑀 and 𝜇. This
mass determines the increase in frequency during the orbital decay. [The name
"chirp" is inspired by the high pitch sound (increasing frequency) produced by
small birds.]
Theory
B.4 Using the information provided above, relate the orbital angular velocity Ω with 2.0pt
the GW frequency 𝑓GW . Knowing that, for any smooth function 𝐹 (𝑡) and 𝑎 ≠ 1,
d𝐹 (𝑡)
= 𝜒𝐹 (𝑡)𝑎 ⇒ 𝐹 (𝑡)1−𝑎 = 𝜒(1 − 𝑎)(𝑡 − 𝑡0 ) , (11)
d𝑡
where 𝜒 is a constant and 𝑡0 is an integration constant, show that (10) implies
that the GW frequency is
On September 14, 2015 GW150914 was registered by the LIGO detectors, consisting of two L-shaped
arms, each 4 km long. These arms changed by a relative length according to Fig. 1. The arms of the
detector respond linearly to a passing gravitational wave, and the response pattern mimics the wave.
This wave was created by two black holes on quasi-circular orbits; the loss of energy through gravita-
tional radiation caused the orbit to shrink and the black holes to eventually collide. The collision point
corresponds, roughly, to the peak of the signal after point D, in Fig. 1.
Figure 1. Strain, i.e. relative variation of the size of each arm, at the LIGO detector H1. The
horizontal axis is time, and the points A, B, C, D correspond to 𝑡 = 0.000, 0.009, 0.034, 0.040
seconds, respectively.
B.6 Estimate the minimal orbital separation between the two objects at 𝑡CD . Hence 1.0pt
estimate a maximum size for each object, 𝑅max . Obtain 𝑅⊙ /𝑅max to compare
this size with the radius of our Sun, 𝑅⊙ ≃ 7 × 105 km. Estimate also their orbital
linear velocity at the same instant, 𝑣col , comparing it with the speed of light,
𝑣col /𝑐.
Theory
Conclude that these are extremely fast moving, extremely compact objects indeed!