Orbital Precession Calculator

Orbital precession in GR refers to the slow rotation of the orbit's orientation in space. The point of closest approach (periapsis) advances slightly each orbit due to spacetime curvature caused by the central mass. Over many orbits this traces out a rosette pattern rather than a closed ellipse. The GR precession rate per orbit is delta_phi = 6*pi*G*M / (c^2*a*(1-e^2)) radians.

Mercury has the fastest orbital period in the solar system (88 days), the highest eccentricity (0.206) among the inner planets, and orbits closest to the Sun where spacetime curvature is strongest. These three factors combine to give it the largest GR precession of any solar-system planet. The calculated GR value is 42.99 arcsec/century, matching the anomalous residual observed since the 1800s.

Newtonian gravity predicts closed elliptical orbits for a pure two-body problem with no precession. In the solar system, the gravitational pulls of other planets cause Newtonian apsidal precession (about 5,557 arcsec/century for Mercury). GR adds an additional curvature-driven precession on top of this. The GR component was the missing 43 arcsec/century that classical mechanics could not explain.

An arcsecond is 1/3600 of a degree, or about 1/206,265 of a radian. 43 arcseconds per century means Mercury's perihelion shifts by 43/3600 degrees (0.012 degrees) per 100 years, a tiny but measurable angle. The historical importance is that it was observed long before GR was developed, creating a known anomaly that any correct theory of gravity had to explain. Einstein's calculation of exactly this value in 1915 was a landmark triumph for GR.

The precession formula contains (1 - e^2) in the denominator. For a circular orbit (e = 0), this factor is 1 and gives the baseline precession. As eccentricity increases toward 1, (1 - e^2) approaches 0, making the denominator smaller and the precession much larger. PSR B1913+16 with e = 0.617 has (1 - e^2) = 0.619, boosting its precession by a factor of 1.62 compared to a circular orbit with the same semi-major axis.

By the mid-1800s, astronomers had measured a 43 arcsec/century residual in Mercury's precession that Newtonian gravity could not account for even after including all known planetary perturbations. Various ad-hoc explanations were proposed including an undiscovered inner planet (Vulcan). In November 1915, Einstein showed that his newly completed general theory of relativity predicted exactly 43 arcsec/century for Mercury, which he called the most exciting result of his life.

PSR B1913+16 (Hulse-Taylor pulsar) is a double neutron star system discovered in 1974. Each neutron star has a mass near 1.4 solar masses, giving a total mass of 2.83 Msun. The orbital period is only 7.75 hours (0.3233 days) and the eccentricity is 0.617. This combination of high mass, tiny orbit, and high eccentricity produces a GR precession of 4.2267 degrees per year (421.98 deg/century), measured via pulsar timing with extraordinary precision.

Yes, the Binary mode accepts any two masses and computes the total mass M = m1 + m2 for use in the formula. The formula remains accurate as long as the orbit is quasi-Keplerian: orbital velocity much less than c and orbital separation much greater than the Schwarzschild radius. For extreme mass-ratio inspirals or final merger phases (v/c > 0.1), higher-order post-Newtonian or numerical relativity methods are needed.

The periapsis is the point of closest approach in an orbit, called perihelion for orbits around the Sun and perigee for Earth orbits. In Newtonian gravity a two-body orbit is a closed ellipse and the periapsis returns to the same point each orbit. In GR, the additional gravitational effect of spacetime curvature acts like a small extra central force that causes the ellipse to rotate slightly each orbit, advancing the periapsis in the direction of orbital motion.

No. The formula delta_phi = 6*pi*G*M/(c^2*a*(1-e^2)) is derived for the Schwarzschild metric (non-rotating, spherically symmetric mass). For Kerr black holes, the precession depends on the spin parameter a* and the orbital plane orientation relative to the spin axis, requiring the full Kerr geodesic equations. The Schwarzschild formula gives a useful first approximation when a* is small.

Three factors drive the enormous precession of binary pulsars compared to planetary orbits. First, neutron stars are about 300,000 times more massive than the Earth. Second, their orbital separations are only a few million kilometers, compared to tens of millions of km for inner planets. Third, high eccentricities amplify the rate. Combined, these factors push PSR B1913+16 to 421.98 degrees per century versus Mercury's 0.0119 degrees per century.

The orbital period T = 2*pi*sqrt(a^3/(G*M_total)) is the full revolution period computed from Kepler's third law using the total mass M = m1 + m2. For PSR B1913+16 this gives 0.3233 days (7.75 hours). This period is also measurable directly from pulsar timing, providing an independent cross-check on the orbital parameters. The number of orbits per century (T_century / T_orbit) then sets how rapidly the annual precession accumulates.