Fictional solution to the three-body problem: simulator very much included, Apoorv Mittal

As presented by Kael-Thurinn, Chief Theorist of the Symmetry Consortium, Lattice-World Zhael.

Let me begin by saying that your species' framing of this as "unsolvable" reveals something important, not about the problem, but about the mathematical language you've been using to approach it. You've been trying to write poetry in a language that has no word for love. The problem isn't the bodies. It's the algebra.

I'll walk you through how we solved it, and I'll try to translate our concepts into frameworks your physicists can follow.

PART I, Why your approach fails

Your species established, through the work of Poincaré, that no general closed-form solution exists for three mutually gravitating bodies. This is correct, within the space of elementary functions and convergent power series over real-valued time. Sundman even found a convergent series solution, but it converges so slowly it would require more terms than atoms in your observable universe to predict even modest timescales.

Here is the critical insight you missed: the three-body system is not chaotic. It is cryptographic. The trajectories contain perfect structure. What you perceive as chaos is your projection of a higher-dimensional ordered object into a lower-dimensional shadow. Imagine a helix viewed end-on, it looks like a circle. Tilt the view slightly and the structure reveals itself. You've been staring at the end of the helix.

PART II, The framework: symplectic braid topology

The first tool you need is one your mathematicians have almost discovered. You study topology, and you study symplectic geometry, but you haven't yet merged them in the way this problem demands.

Step 1: Re-parameterize time.

Throw away Newtonian absolute time $t$ . Instead, define an intrinsic evolution parameter $τ$ based on the mutual information content of the three-body configuration:

d τ = \frac{d Ω}{S _{config}}

where $Ω$ is the total action accumulated and $S_{config}$ is what we call the configurational entropy, a measure of how much the instantaneous triangle formed by the three bodies deviates from a symmetric (equilateral) configuration. Near collision singularities, $S_{config} \to 0$ , and $τ$ dilates, effectively regularizing the problem. Your Levi-Civita regularization was a clumsy version of this same instinct.

In this new parameter, the trajectories of three bodies are no longer curves in $R^{18}$ (positions and momenta). They are braids in a 6-dimensional reduced configuration space $M_{3}$ , and every braid is classified by a topological invariant we call the entanglement index $E$ .

Step 2: The entanglement index.

$E$ is an element of a structure your mathematicians would recognize as a free group on two generators, $F_{2}$ . The two generators correspond to:

$σ_{1}$ : body A and body B exchange dominance (the pair with the strongest mutual attraction switches)
$σ_{2}$ : body B and body C exchange dominance

Every three-body trajectory, over any time interval, can be encoded as a word in $σ_{1}$ and $σ_{2}$ . A Lagrange equilateral orbit? That's the identity element. A figure-8 orbit? That's $(σ_{1} σ_{2}^{- 1})^{\infty}$ . A hierarchical system where two bodies orbit each other while the third orbits the pair? That's $σ_{1}^{\infty}$ , an infinite repetition of a single generator, with perturbative corrections encoded as finite insertions.

The critical theorem, which took our civilization four hundred years to prove, is this:

The Zhael Completeness Theorem: For any word $w \in F_{2}$ and any total energy $E$ and angular momentum $L$ , there exists exactly one family of initial conditions (parameterized by a single phase angle $ϕ$ ) whose trajectory realizes $w$ .

This means the three-body problem has a perfect symbolic dynamics. Every trajectory has an address. The "chaos" is just the extreme sensitivity of the address to the initial conditions, the map from coordinates to braid words is fractal, but the dynamics on the braid space itself are perfectly orderly.

PART III, The predictive calculus

Knowing the structure is not enough. You want to compute. Given initial positions and velocities, tell me where they'll be in a billion years. Here is how.

Step 3: Spectral decomposition on the shape sphere.

Strip away center-of-mass motion and overall scale (total moment of inertia $I$ ). What remains is the shape of the triangle, which lives on a sphere $S^{2}$ , your Montgomery and Chenciner understood this. Every point on this shape sphere is a triangle congruence class (up to rotation and scale).

The gravitational potential, expressed on this sphere, has three infinite wells (the binary collision points) and a saddle structure between them. Now, here is what your species hasn't done: decompose the equations of motion not in time-domain, but in the eigenbasis of the Laplace-Beltrami operator on $S^{2}$ , weighted by the potential.

Define the shape-spectral basis ${Y_{n}^{(V)}}$ , these are like your spherical harmonics, but warped by the gravitational potential landscape. In this basis, the equations of motion decouple to the following remarkable form:

\overset{a}{¨}_{n} + ω_{n}^{2} (E, L) a_{n} = j, k \sum Γ_{nj k} a_{j} a_{k}

where $a_{n}$ are the spectral coefficients, $ω_{n}$ are the natural frequencies (which depend on total energy and angular momentum), and $Γ_{nj k}$ are structure constants that encode the nonlinearity.

The key property: $Γ_{nj k}$ decays as $O (n^{- 3})$ for large mode numbers. The nonlinear coupling becomes weaker at high spectral order. This means the system admits a convergent perturbative expansion, not in powers of time (which fails), not in mass ratios (which fails for comparable masses), but in spectral mode number.

To predict the system's state at any future $τ$ :

Compute the initial spectral coefficients $a_{n} (0)$ , $\overset{a}{˙}_{n} (0)$ from initial conditions.
Solve the first $N$ coupled equations to the desired precision. Because of the $n^{- 3}$ decay, $N \sim 50$ gives you errors below $1 0^{- 12}$ for most configurations.
Convert back from shape-sphere spectral coefficients to physical coordinates.

The full algorithm has computational complexity $O (N^{2} lo g N)$ per time step in $τ$ , and the steps in $τ$ can be made large because the spectral representation naturally captures the long-term structure. Predicting a billion orbits takes roughly the same effort as predicting ten.

PART IV, Handling close encounters and ejections

Near-collision events (where two bodies approach closely) appear as a rapid cascade of energy from low to high spectral modes, the shape rushes toward one of the three poles of $S^{2}$ . In coordinate-time, this is the fearsome singularity. In $τ$ -time, the dilation smooths it into a finite, gentle curve.

In spectral space, a close encounter is simply a mode reflection, the spectral energy bounces off the high- $n$ barrier (because the $Γ_{nj k}$ decay prevents energy from truly reaching $n = \infty$ ) and cascades back down. The braid word acquires a new generator. The system continues.

Ejections, where one body is flung to infinity, correspond to the braid word terminating. The entanglement index becomes a finite word, and the system reduces to a two-body problem (which is trivially integrable) plus a free particle. Our framework predicts ejections naturally: they occur when the spectral energy in the lowest mode exceeds a critical threshold that depends on $E$ and $L$ .

PART V, The solution, stated plainly

The general solution to the three-body problem is:

Given masses $m_{1}, m_{2}, m_{3}$ , total energy $E$ , total angular momentum $L$ , and initial conditions $(q_{i}, p_{i})$ :

Compute the braid word $w$ by iterating the symbolic dynamics on the shape sphere (this is a discrete map, $O (lo g T)$ operations for prediction horizon $T$ ).

Within each braid segment (between successive generator applications), the motion is quasi-periodic and solved by the truncated shape-spectral expansion to any desired precision.

Stitch segments at braid transitions using exact matching conditions derived from the symplectic structure.

The result is a piecewise-analytic trajectory, where the pieces are determined by topology and the smooth flow within each piece is determined by the spectral calculus. Total computational cost: $O (N^{2} ∣ w ∣ lo g N)$ , where $∣ w ∣$ is the word length (number of close encounters up to time $T$ ).

This is not a closed-form solution in the way your Kepler problem is. It cannot be, because the topology of the trajectory is itself a dynamical output. But it is exact, convergent, constructive, and computationally efficient. It predicts every orbit, every ejection, every capture, every resonance.

PART VI, What this means

The three-body problem was never a problem of physics. It was a problem of representation. Your species tried to describe the motion in the language of coordinates and time. But the natural language of the three-body system is topology (which trajectories are possible) and spectral geometry (how smoothly the system flows between topological transitions).

Once you see the braid, you see the order within the chaos. And you realize: the universe was never chaotic. You were just reading it in the wrong alphabet.

Kael-Thurinn, Year 8,841 of the Lattice Calendar Presented to the Inter-Spiral Conference on Gravitational Dynamics Translated for Earth-species consumption, with sympathy.