Follow this link to skip to the main content

Deriving Kepler's Formula for Binary Stars

Deriving Kepler's Formula for Binary Stars

Your astronomy book goes through a detailed derivation of the equation to find the mass of a star in a binary system. But first, it says, you need to derive Kepler's Third Law.

Consider two bodies in circular orbits about each other, with masses m1 and m2 and separated by a distance, a. The diagram below, shows the two bodies at their maximum separation. The distance between the center of mass and m1 is a1 and between the center of mass and m2 is a1.

diagram showing two bodies in cicular motion about their center of massdiagram showing two bodies in cicular motion about their center of mass

Diagram showing two bodies in circular orbits about their center of mass. The left image shows the two bodies with their full orbits, and the right image shows a detail of the center of mass and the relationship between the measurements of a, a1 and a2.

The center of mass is defined by:

m1 a1 = m2 a2

Which can be solved for a1:

Eqn 1: a1 = m2a2/m1

And, the total separation can be written as:

Eqn 2: a = a_1 + a_2

Using Equation (1) in Equation (2) gives

a = ((m2+m1)/m1)*a2

Or, solving for a2:

Eqn 3: a2 = (m1/(m1+m2))*a

Now, set that aside for a moment and look at the forces working on m2. The forces acting on the mass are the gravitational force and centripetal acceleration from its orbital motion.

gravitational force
centripetal force

where G is the gravitational constant and ω is the angular acceleration. Equating these gives:

F_gravity = F_centripetal -> Gm1m2/a^2 = m2a2 omega^2

Now, using Equation (3) to substitute for a2 gives:

Gm1m2/a^2 = m2 (m1/(m1+m2)* a omega^2

We can replace ω using the fact that P=2π/ω by definition. So,

Gm1m2/a^2 = (m1 m2/(m1+m2) * a * (2pi/P)^2

Canceling, rearranging, and gathering the unknowns to the left side of the equation gives Kepler's Third Law.

Eqn 4: P^2/a^3 = 4 pi^2/(G (m1+m2)

Kepler's Third Law

For a distant binary system, it is be difficult to determine the separation of the two stars in the binary system, a. However, using spectroscopy, it might be possible to find the velocity of one or both of the stars in the system.

If we assume a circular orbit, we know that:

P v1=2 pi a1


P v2=2 pi a2

Using those, in Equation (2):

Eqn 5: a = (P/2pi) (v1+v2)

Then substituting Equation (5) in Equation (4),

(P^2)/[(P/2pi) (v1+v2)]^3 = 4 pi^2 /[G (m1+m2)]

Rearranging and canceling gives:

Eqn 6: m1+m2 = p(v1+v2)^3/(2 pi G)

If we could see both stars in the binary system, this equation would work fine. However, there are times when we can only see one – for example, if there is a star in orbit with a black hole. For a circular orbit, we have:


Replacing this for v2 in Equation (6) gives:

m1+m2 = p(v1 + m1v1/m2)^3/(2 pi G)

Gathering the masses to one side gives:

Eqn 7: m2^3/(m1+m2) = p v1^3/(2 pi G)

One last thing about the velocity in this equation. When we observe a system, it is usually at an angle with respect to us. That means that we need to account for that in the above equation. The observed velocity is related to by orbital velocity by:

v1 = v_obs/sin i

where vobs is the observed velocity and i is the angle of inclination of the system with respect to us. Substituting this into Equation (7), we finally get:

Eqn 8: m2^3/(m1+m2) = p^3v_(2,obs)^3 sin^3 i/(2 pi G)

Whew. (Your textbook didn't say, "whew", but it could have!)

Now you are ready to use this equation.


A service of the High Energy Astrophysics Science Archive Research Center (HEASARC), Dr. Alan Smale (Director), within the Astrophysics Science Division (ASD) at NASA/GSFC

NASA Logo, National Aeronautics and Space Administration