N-Body Gravitational Interaction Estimator

The N-body problem asks how N point masses move under their mutual gravitational attraction. For N = 2, an exact analytical solution exists (Kepler orbits). For N = 3 or more, no general closed-form solution exists and the system is typically chaotic. This calculator computes the instantaneous forces and accelerations at a snapshot in time, not the full orbital evolution.

Gravitational force follows an inverse-square law: F = Gm1m2/r^2. Doubling the separation reduces the force by a factor of four. Tripling it reduces the force by a factor of nine. This rapid falloff explains why the Moon's tidal effect on Earth is far stronger than the Sun's, even though the Sun is far more massive.

Gravitational potential energy U = -Gm1m2/r is defined to be zero when the bodies are at infinite separation. Because gravity is attractive, energy is released as the bodies approach each other, so U decreases below zero. A more negative U means the bodies are more tightly bound. The total mechanical energy of a bound orbit is negative.

Escape velocity is the minimum speed a small test mass would need at the given separation to escape the gravitational pull of the larger body to infinity. It equals sqrt(2GM/r). For Earth's gravity at 1 AU from the Sun, the escape velocity from the Sun is about 42 km/s, which is why spacecraft leaving the solar system need a large velocity kick.

In a three-body system, tiny differences in initial conditions grow exponentially over time, making long-term prediction impossible without numerical integration. Even small rounding errors in the starting positions lead to completely different trajectories after a few orbital periods. This sensitive dependence on initial conditions is the hallmark of deterministic chaos.

The gravitational force between Earth (5.972 x 10^24 kg) and the Moon (7.342 x 10^22 kg) at their average separation of 384,400 km is about 1.98 x 10^20 N. This force is what keeps the Moon in orbit and drives Earth's ocean tides. It accelerates the Moon toward Earth at about 2.70 mm/s^2.

Solar masses work best for stars. Jupiter masses (1 Mjup = 1.898 x 10^27 kg) are convenient for gas giants. Earth masses (1 Mearth = 5.972 x 10^24 kg) suit rocky planets and large moons. For sub-lunar bodies, switch to kilograms directly. The unit dropdown converts automatically to SI internally.

No. This tool computes instantaneous force magnitudes and accelerations at a single snapshot. It does not integrate the equations of motion over time, so it cannot tell you whether an orbit is stable, circular, or elliptical. For full trajectory simulation use REBOUND, N-body Gravity Simulator, or similar codes.

Newton's second law gives a = F/m. Both bodies experience the same force magnitude, but the lighter body accelerates far more. For the Earth-Sun pair, the Sun accelerates at 1.78 x 10^-8 m/s^2 while Earth accelerates at 5.93 x 10^-3 m/s^2. The ratio equals the inverse mass ratio: a_Earth/a_Sun = Msun/Mearth = 333,000.

Force F = Gm1m2/r^2 depends on both masses. Acceleration a = GM/r^2 (the gravitational field) depends only on the source mass M and distance r. A feather and a cannonball experience the same gravitational acceleration in free fall even though the cannonball experiences a far larger force, because it also has far more inertia.

The three-body mode sums all three pairwise potential energies: U_total = U12 + U23 + U13, where each U_ij = -Gmi*mj/rij. This total is the energy that would need to be supplied to move all three bodies to infinite mutual separation. It does not include kinetic energy, so it is not the total mechanical energy of the system.