Abstract:The three-body problem is arguably the oldest open question in astrophysics and has resisted a general analytic solution for centuries. Various implementations of perturbation theory provide solutions in portions of parameter space, but only where hierarchies of masses or separations exist. Numerical integrations1 show that bound, non-hierarchical triple systems of Newtonian point particles will almost2 always disintegrate into a single escaping star and a stable bound binary3,4, but the chaotic nature of the three-body problem5 prevents the derivation of tractable6 analytic formulae that deterministically map initial conditions to final outcomes. Chaos, however, also motivates the assumption of ergodicity7-9, implying that the distribution of outcomes is uniform across the accessible phase volume. Here we report a statistical solution to the non-hierarchical three-body problem that is derived using the ergodic hypothesis and that provides closed-form distributions of outcomes (for example, binary orbital elements) when given the conserved integrals of motion. We compare our outcome distributions to large ensembles of numerical three-body integrations and find good agreement, so long as we restrict ourselves to 'resonant' encounters10 (the roughly 50 per cent of scatterings that undergo chaotic evolution). In analysing our scattering experiments, we identify 'scrambles' (periods of time in which no pairwise binaries exist) as the key dynamical state that ergodicizes a non-hierarchical triple system. The generally super-thermal distributions of survivor binary eccentricity that we predict have notable applications to many astrophysical scenarios. For example, non-hierarchical triple systems produced dynamically in globular clusters are a primary formation channel for black-hole mergers11-13, but the rates and properties14,15 of the resulting gravitational waves depend on the distribution of post-disintegration eccentricities.

Statistics of N-Body Simulations. III. Unequal Masses

A comprehensive Nbody study of mass segregation in star clusters: Energy equipartition and escape

Rapid Mass segregation in small stellar clusters

A Comparison of Direct N-body Integration with Anisotropic Gaseous Models of Star Clusters

Signatures of Mass Segregation from Competitive Accretion and Monolithic Collapse

Gravitational collapse from cold uniform asymmetric initial conditions

Stellar, Brown Dwarf, and Multiple Star Properties from Hydrodynamical Simulations of Star Cluster Formation

Comparative study between N-body and Fokker-Planck simulations for rotating star clusters: I. Equal-mass system

Shortening of Core Collapse Time in Star Clusters with Two Mass Components

Testing lowered isothermal models with direct N-body simulations of globular clusters - II: Multimass models

Gravitational conundrum? Dynamical mass segregation versus disruption of binary stars in dense stellar systems

Spherical symmetry breaking in cold gravitational collapse of isolated systems

A statistical solution to the chaotic, non-hierarchical three-body problem

Evolution of the Stellar Mass Function in Multiple-Population Globular Clusters

Direct Collisional Simulation of 10000 Particles Past Core Collapse

The Lifetimes of Star Clusters Born with a Top-heavy IMF

Preparing the next gravitational million-body simulations: evolution of single and binary stars in nbody6++gpu, mocca, and mcluster

The Global Mass Functions of 35 Galactic globular clusters: II. Clues on the Initial Mass Function and Black Hole Retention Fraction

Million-body star cluster simulations: comparisons between Monte Carlo and direct N-body

Comparative Study Between N-body and Fokker–Planck Simulations for Rotating Star Clusters – II. Two-component Models