Tuesday, November 15, 2016

mathematical physics - Can we infer the existence of periodic solutions to the three-body problem from numerical evidence?


I recently found out about the discovery of 13 beautiful periodic solutions to the three-body problem, described in the paper



Three Classes of Newtonian Three-Body Planar Periodic Orbits. Milovan Šuvakov and V. Dmitrašinović. Phys. Rev. Lett. 110 no. 11, 114301 (2013). arXiv:1303.0181.



I am particularly impressed by how elaborate the solutions are, and I'm struck by the tantalizing hint of an infinity of other distinct orbits given by the analogy with a free group. The solutions can be viewed in the Three-Body Galery, which has animations of the new orbits in real space and in something called the 'shape sphere', which is described in the paper.


enter image description here




I was aware already of the figure-of-eight solution, which is described nicely in




A new solution to the three body problem - and more. Bill Casselman. AMS Feature Column.



and which was discovered numerically by Christopher Moore (Phys. Rev. Lett. 70, 3675 (1993)). I understand that the figure-of-eight solution has been proven to actually exist as a solution of the ODE problem, in



A Remarkable Periodic Solution of the Three-Body Problem in the Case of Equal Masses. Alain Chenciner and Richard Montgomery. Ann. Math 152 no. 3 (2000), pp. 881-901.



There is also a large class of solutions called $N$-body choreographies by Carles Simó, in which a number of bodies - possibly more than three - all follow the same curve. Simó found a large class of them in 2000 (DOI/pdf), though this nice review (DOI) seems to imply that formal theorematic proof that they exist as periodic solutions of the ODE problem is still lacking.




So, this brings me to my actual question. For the numerical simulations, however well you do them, in the end you will only have a finite-precision approximation to a solution of a differential equation that is propagated for a finite time. Additionally, you might do a numerical stability analysis that strongly suggests (or rigorously proves?) that you are (or are not) in a stable orbit. However, this is quite a ways away from a rigorous existence theorem for a periodic orbit with that symmetry.



With this in mind, then, in what spirit are these simulations done? Is it purely a numerical approach, in the hope that good numerics do indicate existence, but with a rigorous proof left to the mathematicians through whatever other means they can manage? Or is there some overarching theorem that indicates the existence of a truly periodic solution after a given threshold? What tools are there for proving existence theorems of periodic solutions?



Answer



It seems like they were able to rigorously prove the existence of N-body choreographies by using interval Krawczyk method to show that a minimum exist to the variational problem solved in the subspace of the full phase space satisfying some symmetry conditions.


Following the links given I found this paper where they explain the method. It's not exactly a light reading material but on page 6 they say: "If all these conditions all fulfilled, then from Theorem 4.5 we are sure that in the set $Z \times \{c_0\}$ there is an initial condition for the choreography. Moreover, as the set Z is usually very small, the shape of the proved choreography is very similar to our first approximation."


It sounds like starting with "an initial guess", they are able to show that there exist an "exact solution" very close to this initial guess. And one can probably obtain a curve that is arbitrarily close to the actual solution by doing more and more precise calculation. But the existence of the choreography is established rigorously with the assist of their numerical method.


Note that in the beginning of the paper, they mention the solutions obtained by the usual numerical methods as "solutions produced in a non-rigorous numerical way."


No comments:

Post a Comment

classical mechanics - Moment of a force about a given axis (Torque) - Scalar or vectorial?

I am studying Statics and saw that: The moment of a force about a given axis (or Torque) is defined by the equation: $M_X = (\vec r \times \...