Comment by deskamess

9 hours 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?

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

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.

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.

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.

  • 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.

    • 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.

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.

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.

  • > 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.

That's exactly the case

  • 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).