Binary orbits as mass measurements

The relative orbit measures the separation of the components

A binary star system contains two stars bound through their mutual gravity. In a sufficiently isolated Newtonian two-body description, each component moves about their common centre of mass. The relative position of one component with respect to the other also follows a Keplerian orbit.

Let a₁ and a₂ be the physical semimajor axes of the components about that centre of mass. The relative semimajor axis is a = a₁ + a₂, and the barycentric relation is M₁a₁ = M₂a₂. The more massive component has the smaller barycentric semimajor axis.

For positive constant masses and a bound nonradial inverse-square orbit, the period P satisfies M₁ + M₂ = 4π²a³/(GP²). The a in this expression belongs to the physical relative orbit. One component’s semimajor axis or a projected separation measured in a single image cannot be substituted without further information.

This model neglects material exchange, appreciable extended-body effects, perturbations and relativistic corrections. It is useful when those effects are sufficiently small for the required accuracy. Interacting systems need a more complete physical account.

Different observations reveal different projections

Resolved positions over time can trace an orbit on the sky. An angular scale requires a distance to become a physical length, and a projected orbit requires a reconstruction of the spatial orientation. A single angular snapshot supplies neither the complete orbit nor its period.

Spectroscopy can reveal changing Doppler shifts of one or both components. When suitable spectral features from both stars can be followed, a double-lined spectroscopic orbit provides two velocity curves. When only one component is measurable, a single-lined orbit provides less information.

Write γ for the system’s line-of-sight velocity offset and K₁, K₂ for the positive orbital radial-velocity semi-amplitudes. Each K is half the fitted orbital maximum-minus-minimum velocity; it excludes γ. It is not the largest raw observed velocity when the entire system is also moving.

Define inclination i as the angle between the orbital plane and the sky plane: i = 0 is face-on and i = 90° is edge-on. For a circular orbit, K_j = 2πa_j sin i/P. Projection suppresses both component amplitudes when the orbit is seen closer to face-on.

Build a velocity table with all conventions stated

Choose arbitrary positive length and time units L₀ and T₀, and define M₀ = L₀³/(GT₀²). In an original isolated circular Newtonian example, set M₁ = 2M₀, M₂ = M₀ and relative semimajor axis a = 3L₀. These are stipulated inputs, not measured stars.

The barycentric axes are a₁ = L₀ and a₂ = 2L₀. The period relation gives P = 6πT₀ and angular orbital frequency 2π/P = 1/(3T₀). Adopt i = 30°, so sin i = 1/2.

The projected semi-amplitudes are then K₁ = L₀/(6T₀) and K₂ = L₀/(3T₀). For the velocity table, additionally set γ = L₀/T₀, positive velocity receding, and phase φ = 2πt/P with v_r,1 = γ + K₁ sin φ and v_r,2 = γ − K₂ sin φ.

The curves vary in opposite orbital directions about the common γ. Their nonzero offset is translation of the system in this teaching description; it does not increase the two orbital semi-amplitudes. At phase zero the radial components of orbital velocity vanish even though the circular orbital speeds do not.

The two amplitudes share the same projection

For this common two-body orbit, K₁/K₂ = a₁/a₂ = M₂/M₁. The example gives (1/6)/(1/3) = 1/2, so the mass ratio is recovered even before an independent inclination is supplied. It does not follow that both absolute masses have been measured.

In the circular model, the amplitudes and period give a sin i = P(K₁ + K₂)/(2π). Here that projected relative axis is 1.5L₀. Combining it with the period gives (M₁ + M₂)sin³i = 4π²(a sin i)³/(GP²) = 3M₀/8.

With the independently supplied i = 30°, divide by sin³i = 1/8 to recover total mass 3M₀. If instead the same projected observations were interpreted with i = 90°, the inferred total would be only 3M₀/8. The actual projection must be constrained rather than selected for convenience.

When only component 1 is measurable, its single-lined mass function is f = PK₁³/(2πG) = M₂³sin³i/(M₁ + M₂)² for a circular orbit. The supplied example gives f = M₀/72. This is one combination of masses and inclination, not an independent measurement of each mass.

An orbit fit needs physical and observational checks

Eccentric orbits require the appropriate velocity model. Their amplitude relation includes the eccentricity factor, and the circular mass-function expression acquires (1 − e²)³ᐟ² on its observational side. A sinusoid should not be forced onto every velocity history.

Astrometry can constrain orientation when sufficient positions, distance information and a valid model are available. Eclipse geometry adds another useful constraint, but the presence of eclipses does not establish an exactly 90° inclination. The next lesson examines what the combined light actually measures.

Calibration drift, stellar activity and changing or blended line profiles can affect measured velocities. Additional companions or a changing orbit can also undermine a simple two-body interpretation. Residual structure and independent observations help test whether the adopted model is adequate.

A good dynamical mass result records the orbit convention, data coverage, distance or inclination constraints, uncertainties and physical model. The strength of the method comes from connecting measured motion to tested dynamics while retaining the information that remains projected or unconstrained.

Worked examples

Try the reasoning

1. May a₁ alone replace a = a₁ + a₂ in the total-mass period formula?

   Show the answer and reasoning

   No.

   That formula uses the physical relative semimajor axis, not just one component’s barycentric orbit.

2. Does a nonzero systemic velocity γ have to be included in K?

   Show the answer and reasoning

   No.

   K is the semi-amplitude of the orbital variation about the systemic offset.

3. Does one single-lined circular mass function determine both component masses without further constraints?

   Show the answer and reasoning

   No.

   It supplies a combination of masses and inclination rather than separate unique values.

4. Does detecting an eclipse prove i = 90° exactly?

   Show the answer and reasoning

   No.

   Eclipses constrain orientation through a geometry model; finite stellar sizes allow a range of inclinations.

What to carry forward

• The relative physical orbit and period constrain a total mass under stated dynamics.
• Two component velocity amplitudes can constrain a mass ratio while leaving an inclination-dependent absolute scale.
• Astrometry, spectroscopy and eclipse modelling supply complementary evidence rather than interchangeable measurements.