An atlas of periodic solutions to the three-body problem

368 pointsposted 5 days ago
by danielmorozoff

91 Comments

pixelpoet

3 days ago

I have something similar at https://gravitoy.xyz

With it I have discovered up to 11-dimensional choreographies, see https://lycium.github.io/hyperchoreography/ and code at https://github.com/lycium/hyperchoreography/

Exposition video: https://youtube.com/watch?v=sIfff10hYZA

Example rendered output from Gravitoy (not of a choreography though): https://www.youtube.com/watch?v=N3BwCoiwsGk

Looks like I need to update my catalogue to take into account the many different 2D choreographies from the references on this site!

emmelaich

3 days ago

That Pythagorean one with the ejection is cool. I wonder whether you could use it to accelerate a spacecraft.

SiempreViernes

2 days ago

Cool, are any of these particularly stable? Like, is there any choreography that could plausibly survive in practice?

pixelpoet

2 days ago

Most of them are extremely unstable, which is what the Morse number encodes - the number of directions you can roll down doing action minimisation. I don't think any N-body choreographies can exist in nature, because it requires point masses, exactly equal masses, Newtonian gravity etc.

BTW, I've wondered if it's possible to have N-body choreographies in general relativity; maybe some brave person can explore that :)

mr_mitm

3 days ago

It's well organized, the animations are smooth, and it looks beautiful... I'm not sure what to do with the information, but it's mesmerizing and fascinating. Great find!

ivanjermakov

3 days ago

> I'm not sure what to do with the information

Get inspired and spend a rest of your life looking for new 3bp solutions.

dice

3 days ago

Handy if you get really good at it because "the rest of your life" will be indicated by the unavoidable countdown timer that shows up in your vision!

ironSkillet

3 days ago

Something cool and interesting for its own sake. A rare find.

addaon

3 days ago

Pretty useful if you’re a Puppeteer, though.

doormatt

3 days ago

Klemperer rosette FTW.

dreamcompiler

3 days ago

I took graduate orbital mechanics from Roger Broucke. He was one of my best professors. Not only did I learn from him what orbital elements were, but he also taught me the Runge-Kutta numerical integration method.

I didn't learn until years later that he had discovered several of the periodic solutions to the three-body problems. You'll see his name on this page.

deskamess

3 days ago

I guess I misunderstood or mis-scoped the problem. Does the 3-body problem state 'the general case' has no solution, but that does not preclude some configurations from having a solution?

layer8

3 days ago

The three-body problem only states the problem to solve, it doesn’t itself state anything about the existence or non-existence of solutions. It has been proven that there is no general closed-form solution. And there are obvious solutions for trivial special cases, such as three equal masses in an equilateral triangle rotating around each other.

Further reading: https://en.wikipedia.org/wiki/Three-body_problem#Solutions

Sharlin

3 days ago

There is no closed-form solution for finding the roots of >4th degree polynomials in general, but that doesn’t preclude many families of >4th degree polynomials from having closed-form solutions. As a trivial example, x^5 - 1. The exact same thing with the three-body problem.

PaulHoule

3 days ago

In the case of a planet going around the sun we know that the planet travels in an ellipse, more generally a conic section including interstellar comets. Orbital periods and everything else can be computed straightforwardly with formulas.

In the three body problem you can always do a numerical integration (e.g. simulation) and this is valid for a certain amount of time but will not be valid forever because of: (1) chaotic motion which amplifies small errors exponentially over time and (2) celestial mechanics is symplectic which means it conserves certain geometric properties and most integrators are not symplectic and must have different long term dynamics. There are symplectic integrators but they don’t work as well overall as ordinary integrators.

We do not know, for instance, if the solar system is stable. In the short term the planets seem to be basically doing their own thing in their orbits with just minor perturbations. We know the orbits vary a bit over millions of years. We aren’t sure which side of the sun Pluto will be on in 30 million years. It’s very believable that the planets are doing the same thing in 4 billion years but we can’t rule out that the orbits could change in a big way or one could get ejected.

SideQuark

3 days ago

> In the case of a planet going around the sun we know that the planet travels in an ellipse

That's a Newtonian approximation. Under General Relativity it is not true - the possible solutions are time varying, and I'm not even sure there are periodic orbits, nor even stable orbits. All N-body patterns radiate gravitational energy till the system collapses. I'm pretty sure there are no general closed form solution, and I think it's true there is not even a single case where the path is a closed form solution.

PaulHoule

2 days ago

In general relativity it is closer to quasiperiodic.

Classical motion under an 1/r^2 field is really strange as a dynamical system because the periods of rotation, in-out and up-down motion are all the same which is why the motion closes as an ellipse. It is beautiful in a lot of ways but a godawful mess from the viewpoint of perturbation theory because these are always in resonance as opposed to only in resonance occasionally.

Ignoring gravitational radiation, in GR (or with secular perturbations from other planets) those periods no longer match up so the orbit goes in-and-out in not quite the same time it goes around and then you get a precession so it is still basically an ellipse but the angle of the ellipse changes and you'd see something like a spirograph if you draw it. It's a quasi-periodic orbit.

You can write down closed forms (often with special-function integrals) for things like 1/r^3 and it is the same story, the variables separate nicely because angular momentum is conserved.

kurthr

3 days ago

It's worth noting that the Sun-Earth-Moon "system" is much more complex than a three-body problem just due to momentum transfer from tides/bending. That's ignoring the effect of all of the other planets and asteroids or solar wind.

The 3BP is just the simplest chaotic system showing the limits of simple models and approximations.

PaulHoule

3 days ago

Well there are simpler chaotic Hamiltonians like Henon–Heiles.

There are a lot of things wrong with how we teach classical mechanics and one of them is that the two index problems that are used in the undergraduate course are the harmonic oscilator and celestial dynamics and these are the worst non-generic problems that there are. The harmonic oscillator is generic in quantum theory and is a good place to start doing perturbation theory from but the fact that the frequency doesn't change with amplitude makes it non-starter in classical perturbation theory.

In the 2-body case the periods to: (1) go around the sun, (2) go in and out towards the sun, and (3) go up and down out of the plane are all the same which again breaks perturbation theory... and of course classical perturbation theory is difficult and doesn't work that well even in cases where it does work. Contrast that to quantum mechanics where you can start doing simple calculation with perturbation theory, like to calculate the lifetime of an excited state, right away with paper and pencil.

Note that, in general, you can get chaos with as few as two position variables. In the case of the two body problem you have six position variables, but because of conservation of momentum the center of mass doesn't really matter, so you can pick a coordinate frame where the total angular momentum is zero and there are just three variables that matter which is the vector between the two bodies. There are a lot of conserved quantities here, especially angular momentum so the (1) and (2) and (3) motions all do the same thing coincidentally with the same period!

One you add a third body you are adding three more variables but not any more conserved quantities so it is a problem with a lot of dimensions.

Mapping out the orbits of the three body problem has gone pretty slowly because, compared to simpler Hamiltonians, we have to search for those orbits in a high dimensional space. You can find a periodic orbit numerically if you know where it is, but it takes a systematic approach to find many of them.

btilly

3 days ago

The general case always has a solution. At least until the point where two of the three bodies meet (which is a singularity). We can approximate that solution numerically.

The problem is that the solutions very strongly tend to be chaotic. Meaning that small differences in initial conditions, tend to grow exponentially with time. Which means that if you measure everything to 3 digits of precision, in finite time it will stop looking like the actual solution. Every additional digit of precision adds a similar finite time to how long the approximation is good for.

So when finally found, say, the 1953 BC conjunction described in https://en.wikipedia.org/wiki/Conjunction_%28astronomy%29?#N... - that was a very good stress test for our estimated planetary data. Because surprisingly small errors in modern data would have kept that conjunction from happening.

dreamcompiler

3 days ago

Yes. It's kind of like the halting problem: You cannot write a general computer program that will analyze the source code of any random other computer program and tell you if it will halt.

You can write a program that will analyze the code of a few specific other programs and tell you if they will halt. You just can't do it in general.

The 3-body problem is like that. Except it's much harder to find stable 3-body problems than computer programs that are predictable.

IAmBroom

3 days ago

> Except it's much harder to find stable 3-body problems than computer programs that are predictable.

Proving that statement is true might be harder than either of the other two issues.

mr_mitm

3 days ago

I'm pretty sure there is always a unique solution to the equations of motions (safe for some pathological edge cases perhaps). Classical mechanics is deterministic, after all. But for more than two bodies, there is in general no solution in closed form, and it's often chaotic, so not even computeable for arbitrary time frames.

The "about" info states that all of these are computed numerically.

coderenegade

3 days ago

There is a general solution to the three body problem due to Sundman, in the form of an infinite series. Unfortunately it converges so slowly that it's useless in practice.

Being unable to state the solution of a problem in closed form doesn't prevent solutions from existing.

incognito124

3 days ago

That's exactly the case

isolli

3 days ago

It's also not computable, as in chaotic. Small differences in initial positions will lead to unpredictably large differences in trajectory (with small and large having specific meanings to match the formal definition of a chaotic system).

cobbzilla

3 days ago

It’s my understanding that most 3-body orbits are unstable; minor random perturbations will set them adrift.

Are there any 3-body orbits with a natural resonance that maintains the shape of the orbits?

If so, how large of a disturbance can the most-stable 3-body orbit withstand?

chris_st

3 days ago

A few that I looked at are annotated, "It is linearly stable (largest Floquet multiplier 1.000), so a small nudge only makes it wobble." Doesn't define "small nudge", alas...

cobbzilla

3 days ago

further research indicates some do exist:

- equilateral triangle is the simple case

- some figure eight configurations

- some configurations that act like two nested binary systems

there might be others. fascinating stuff.

slwvx

3 days ago

One interesting scenario is where the Sun, Jupiter, and another object are the three bodies in question. I think we could put an object at the L4 or L5 point of Jupiter's orbit (see the Trojans) and have it be stable. How large could this object be? Would an Earth-sized object be stable there? If so, how close would a passing star have to be to disturb that stability?

threebodyorbits

3 days ago

Hi all, I just added a new feature to the site. You can now contribute compute via your browser to help find novel orbits. If you find a linearly stable orbit, you get to name it!

https://www.threebodyorbits.com/hunt

Also thank you so much for all the positive feedback, I truly appreciate it. If you have any additional suggestions please let me know

msuvakov

3 days ago

Really cool project! On a personal note, thanks for the reference! I'm one of the authors of the 2013 paper you cited, and it's awesome seeing how many new solutions have been found since then.

I have a quick question: how are you currently checking if a new candidate is genuinely novel rather than a rotated, time-shifted, or rescaled duplicate of an existing solution?

If you are not already computing it, topological classification (mapping orbits to free group words) would be a great way to quickly prune duplicates before heavy numerical checks, and it's also really useful metadata to display in the atlas itself. Also, showing shape sphere projections of the orbits would be a nice visual addition.

threebodyorbits

2 days ago

Thank you very much, I'm so excited that an original author of one of the three body problem papers that this atlas is based on found this.

After you sent your question I noticed an issue with the "hunt" settings. Orbit novelty was not properly assessed, and it just kept re-discovering already existing orbit families. I updated the search based on your suggestions and just submitted an update that now checks candidates like this:

- Does the motion really repeat? An independent more high precision periodicity check is performed to confirm that bodies truly return to their starting positions and velocities - Is it a different form of an existing orbit? Compares against reference orbits from the atlas, allowing for rotation, different starting points in time, rescaling, reflection, and relabeling of bodies -Is it just the same orbit traced several times? Checks for shorter repeated periods -Now topological information is taken into account as well. In the current version, I use reduced syzygy sequences that record which body is in the middle whenever the three bodies are lined up. This was easier to implement than the familiar a, b, A, B notation.

Also, new candidates here are not automatically declared as new discoveries. They will undergo additional validation before they are added to the atlas.

The shape-sphere suggestion sounds great - would you suggest adding this to the "hunt" page, or to every orbit in the atlas?

altairprime

3 days ago

How did you end up choosing the color scheme for this site?

threebodyorbits

2 days ago

This whole project was inspired by the Cixin Liu Three Body Problem books. I tried to capture the essence and feelings that I had when reading them in the color scheme and design of the website. I wanted a black background, representing the darkness, vastness, and mystery of space. For the orbits and other elements, I wanted the dominant color to represent danger and death, but also feel like a celebration of life and beauty. I thought blood-related colors best fit that purpose

jrflo

3 days ago

Quick, someone tell the Trisolarans!

hypersoar

3 days ago

The Trisolarians are supposedly at Alpha Centauri, a trinary star that isn't actually that chaotic. It has two sun-like stars orbiting each other quite close, and then a red dwarf, Proxima Centauri, orbiting that pair really really far away.

m4rtink

3 days ago

Given all the illogical insanity they are doing (like, trying to fight an interstellar war instead of just the simple stuff like moving to habitats, improving their bodies or even so,me stellar lifting) I don't think the will listen. ;-)

criemen

3 days ago

While I love the books as much as everyone else, the pedant in me has to point out that the Trisolareans have a four body problem at hand (3 suns and their planet), not a three body problem.

Rexxar

3 days ago

I think we can consider the planet(s) negligible from a mass perspective.

joquarky

3 days ago

Tell that to butterflies.

WA

3 days ago

What's the chance that either one of the four bodies crash into another?

inatreecrown2

3 days ago

Very cool visuals and site! Could I make a suggestion: You show the masses (1,1,1), but not the starting positions, which alter the course of events too.

MarkusQ

2 days ago

Any set of points (w. velocities) along the trajectories could be the "starting positions".

RALaBarge

3 days ago

This is an amazing looking website, I like it a lot.

MeteorMarc

3 days ago

I assume some of the solutions are stable against small perturbations, while others are not. That would be interesting to see.

QuesnayJr

3 days ago

If you go to the individual solutions, the text description tells you if it's stable. There's also a slider that allows you to perturb the orbit so you can see for yourself when you perturb it.

summa_tech

3 days ago

I think that's what "STABLE ONLY" clickable text filters by.

ur-whale

3 days ago

Oh, this is so nice, but man is it hugely frustrating not to be able to rotate the thing in 3D.

Or did I not find the controls?

hakuseki

3 days ago

I was surprised that I couldn't find any simple-looking solutions in this atlas. At first I was looking for Lagrange orbits, but maybe it makes sense to exclude them if zero-mass bodies aren't allowed. I think the equilateral triangle ought to be included though.

pletnes

3 days ago

They’re not stable except at L4 and L5, and they all assume oke body to be massless. Arguably they are 2-body orbits for that reason. Not sure but suspect that this atlas contains non-massless bodies.

moritzwarhier

3 days ago

Wow. This is really cool. Deterministic chaos is my absolute favorite in all the nerdy things there are to like in the abstract world.

edbaskerville

3 days ago

Deterministic chaos is cool. But this is even more special in a sense: for a problem where random initial conditions are almost always chaotic, this is a catalog of periodic orbits—these are all non-chaotic.

It would be cool to pair this with a numerical simulator that shows what happens when you perturb any of them.

EDIT: oh, it already does this, thanks other comments

IshKebab

3 days ago

I assume this is at least partially vibe coded, but this is the first good vibe coded website I've seen. Amazing work.

nautilus12

3 days ago

Is this assumed to be 2D? I was going to ask if there are any observed examples of 3 body equilibrium observed in nature.

raverbashing

3 days ago

Yes. Because 3 points are coplanar, so every "3D problem" with 3 objects can be turned into a 2D problem on the correspondent plane

(of course in real life your plane would keep changing, and probably some other complicated math I can't think right now)

petsfed

3 days ago

Put another way, while their positions are one set of 3 points, their momenta are another set of 3 points, and there is no requirement that 6 points will always lay on the same plane.

I wonder what phantom forces would appear when the reference frame changes in some complicated fashion. We get centrifugal "force" when we reconstruct F=dP/dt in a rotating reference frame, what would the 3-body "force" look like?

nautilus12

3 days ago

Oh makes perfect sense, do you have thoughts about real life examples? I did some research using AI and it said there were examples of restricted 3 body problems like the trojan asteroids, but no examples in real life similar to what is in this web app

raincole

3 days ago

Why would the plane keep changing? If there are only these three objects, won't the vectors of their gravitational pull to each other all be on this plane too?

mr_mitm

3 days ago

If you define the initial conditions such that their relative velocity is zero or parallel to the plane, yes. But that's not the case in general.

raverbashing

3 days ago

Not from an external point of view, as you might have a momentum component perpendicular to that plane

(but yes I think you might be right if we're centered on the CG)

twnettytwo

3 days ago

One of the necessary conditions here is that the three objects return to their exact initial position, and so does the centre of mass. Initial conditions with non-zero momentum must trivially be ruled out. But this doesn't stop them from having velocities perpendicular to the initial plane that cancel out perfectly, so this doesn't refute the assertion that the planes keep changing.

btilly

3 days ago

Sorry, but this is a word salad.

Most of the solutions always have non-zero momentum, including in the initial conditions. And https://numericaltank.sjtu.edu.cn/three-body/three-body.htm includes periodic solutions that move in all three dimensions.

drdeca

3 days ago

They presumably meant non-zero total momentum. If the total momentum were non-zero, then the center of mass will be moving in a straight line, and will not return to where it began, and therefore the orbit would not be periodic.

PaulHoule

3 days ago

Big picture the center of mass momentum is concerned and it does not matter if the system as a whole is moving up or down or to the right or the left. Like the Earth is basically orbiting the sun in an ellipse [1] so far as the sun is concerned and from the viewpoint of the solar system not care so much that it is moving around the galaxy unless we are interested that orbit being perturbed by other stars that we pass near over millions and milions of years.

[1] ignoring the parameters of that ellipse changing slightly and slowly thanks to the other planets

btilly

3 days ago

Oh. That makes sense. Yes. The total momentum has to be zero.

All of the things in the solution can constantly have momentum.

twnettytwo

3 days ago

Oh sorry, yeah - I wasn't fully awake yet when I wrote it. As a sibling comment mentioned, I was talking about the total momentum. Anyway, it was in response to GP's assertion that GGP may be correct if "we're centered on the CG", which doesn't make sense because the CG can't move.

crisf

3 days ago

Each orbit is a spinning top. You pull on a top from the side, it's spin axis precesses.

hammock

3 days ago

It’s an arbitrary plane, chosen at each moment just so you can flatten it

btilly

3 days ago

https://numericaltank.sjtu.edu.cn/three-body/three-body.htm shows that there are three dimensional solutions.

We have no observed examples in nature of three body equilibrium. But then again, all places we have looked are either influenced by the chaotic orbits around them of the Solar System, our surrounding galaxy, or nearby galaxies in a cluster.

There aren't a lot of orbiting three bodies without external gravitational influences disturbing them.

IAmBroom

3 days ago

> There aren't a lot of orbiting three bodies without external gravitational influences disturbing them.

This is of course a relative statement. Every object affects every other object, subject to the limitations of lightspeed propagation of gravity waves through expanding space.

But as you point out, we still haven't noticed any examples that are stable short-term.

btilly

3 days ago

That depends on what you mean by stable short-term.

The Sun, Jupiter, and any small asteroid at L4 of L5 is, by itself, linearly stable. Meaning that any small perturbation will only grow linearly.

But, of course, this system interacts with Saturn. The arrangement of those three objects is still remarkably stable. But the interaction with Saturn makes for a chaotic system again.

larodi

3 days ago

So proud to see the prof. who oversaw my masters thesis has their paper (and solutions) featured here. :D

sajithdilshan

3 days ago

Very nicely animated. Also didn’t know there can be so many stable solutions for 3 body problem

khalic

3 days ago

I'm going to spend so much time on this website, very good work

boringg

3 days ago

This is pretty cool to be able to see all the varieties.

kmitz

3 days ago

Fantastic UX, so smooth even on smartphone ! Congrats

RugnirViking

3 days ago

very ai but also pretty cool. wheres the data source? could I find my own periodic solution?

freakynit

3 days ago

This is beautiful. And fast!!!

sscaryterry

3 days ago

Holy fuck this is one of the coolest things I've seen.