Lesson 15 of 15

Binary Inspiral Time

Binary Inspiral Time

How long does it take two compact objects to spiral together from an initial separation r0r_0? Integrating the orbital decay rate dr/dt=(64/5)G3m1m2(m1+m2)/(c5r3)|dr/dt| = (64/5) G^3 m_1 m_2 (m_1+m_2)/(c^5 r^3) from r0r_0 to 00 gives the Peters formula for a circular orbit:

Tmerge=5c5256G3r04m1m2(m1+m2)T_{\rm merge} = \frac{5\,c^5}{256\,G^3} \cdot \frac{r_0^4}{m_1 m_2 (m_1+m_2)}

This r04r_0^4 dependence is dramatic: halving the initial separation cuts the merger time by 16.

Separation as a Function of Time

Running the integral in reverse, the separation remaining at time tt before merger is:

r(t)=(r042565G3m1m2(m1+m2)c5t)1/4r(t) = \left(r_0^4 - \frac{256}{5} \frac{G^3 m_1 m_2 (m_1+m_2)}{c^5}\,t\right)^{1/4}

At t=0t=0 this gives r0r_0; at t=Tmerget = T_{\rm merge} the bracket vanishes and r=0r = 0.

Physical Timescales

Systemr0r_0TmergeT_{\rm merge}
Earth–Sun1.496×10111.496 \times 10^{11} m1023\sim 10^{23} years
NS binary at 0.01 AU1.496×1091.496 \times 10^{9} m5.8×108\sim 5.8 \times 10^{8} years
30+30 MM_\odot BH at 10810^8 m10810^8 m1.2\sim 1.2 years

Constants (define inside each function)

  • G=6.674×1011G = 6.674 \times 10^{-11} m³ kg⁻¹ s⁻²
  • c=299792458c = 299792458 m/s

Your Task

Implement the three functions below.

Python runtime loading...
Loading...
Click "Run" to execute your code.