Wolfram Data Repository

Immediate Computable Access to Curated Contributed Data

Collisionless Periodic Orbits in the Free-fall Three-body Problem

Source Notebook

The 315 Collisionless Periodic Orbits in the Free-fall Three-body Problem

Originator: X. Li and S. Liao

The three-body problem seeks solutions for the orbits of three masses under the mutual influence of gravity. Periodic solutions are very hard to find.

(316 elements)

Examples

Basic Examples

Extract velocities:

In[1]:=
ResourceData[ "Collisionless Periodic Orbits in the Free-fall Three-body Problem"]["F1(1,1,1)"]["Position"]
Out[1]=

Maximum value of period:

In[2]:=
Max[AssociationMap[ ResourceData[ "Collisionless Periodic Orbits in the Free-fall Three-body Problem"][#]["Period"] &, Keys[ResourceData[ "Collisionless Periodic Orbits in the Free-fall Three-body Problem"]]]]
Out[2]=

All initial velocities:

In[3]:=
Graphics[Point[ Values[AssociationMap[ ResourceData[ "Collisionless Periodic Orbits in the Free-fall Three-body Problem"][#1]["Position"] &, Keys[ResourceData[ "Collisionless Periodic Orbits in the Free-fall Three-body Problem"]]]]], AspectRatio -> 1, Frame -> True, FrameTicks -> True]
Out[3]=

Use NBodySimulation to solve an orbit:

In[4]:=
Module[{n = 1, view = "real space"}, With[{parameters = ResourceData[ "Collisionless Periodic Orbits in the Free-fall Three-body Problem"][ Keys[ResourceData[ "Collisionless Periodic Orbits in the Free-fall Three-body Problem"]][[n]]]}, Module[{data = NBodySimulation[ "InverseSquare", {<|"Mass" -> parameters["Mass"][[1]], "Position" -> {-.5, 0}, "Velocity" -> {0., 0.}|>, <|"Mass" -> parameters["Mass"][[2]], "Position" -> {.5, 0}, "Velocity" -> {0., 0.}|>, <|"Mass" -> parameters["Mass"][[3]], "Position" -> parameters["Position"], "Velocity" -> {0., 0.}|>}, parameters["Period"]]}, ParametricPlot[ Evaluate[data[All, "Position", t]], {t, 0, parameters["Period"]}]]]]
Out[4]=

Plot in a "shape sphere" using Jacobi three-body coordinates:

In[5]:=
data = Module[{n = 45}, With[{parameters = ResourceData[ "Collisionless Periodic Orbits in the Free-fall Three-body Problem"][ Keys[ResourceData[ "Collisionless Periodic Orbits in the Free-fall Three-body Problem"]][[n]]]}, Module[{data = NBodySimulation[ "InverseSquare", {<|"Mass" -> parameters["Mass"][[1]], "Position" -> {-.5, 0}, "Velocity" -> {0., 0.}|>, <|"Mass" -> parameters["Mass"][[2]], "Position" -> {.5, 0}, "Velocity" -> {0., 0.}|>, <|"Mass" -> parameters["Mass"][[3]], "Position" -> parameters["Position"], "Velocity" -> {0., 0.}|>}, parameters["Period"]]}, data]]]
Out[5]=
In[6]:=
Show[ParametricPlot3D[{(((x1 - x2) (x1 + x2 - 2 x3))/Sqrt[ 3] + ((y1 - y2) (y1 + y2 - 2 y3))/Sqrt[3])/( 1/2 (x1 - x2)^2 + 1/6 (x1 + x2 - 2 x3)^2 + 1/2 (y1 - y2)^2 + 1/6 (y1 + y2 - 2 y3)^2), (-(1/2) (x1 - x2)^2 + 1/6 (x1 + x2 - 2 x3)^2 - 1/2 (y1 - y2)^2 + 1/6 (y1 + y2 - 2 y3)^2)/( 1/2 (x1 - x2)^2 + 1/6 (x1 + x2 - 2 x3)^2 + 1/2 (y1 - y2)^2 + 1/6 (y1 + y2 - 2 y3)^2), ( 2 (-(((x1 + x2 - 2 x3) (y1 - y2))/( 2 Sqrt[3])) + ((x1 - x2) (y1 + y2 - 2 y3))/(2 Sqrt[3])))/( 1/2 (x1 - x2)^2 + 1/6 (x1 + x2 - 2 x3)^2 + 1/2 (y1 - y2)^2 + 1/6 (y1 + y2 - 2 y3)^2)} /. Thread[{x1, y1, x2, y2, x3, y3} -> Flatten[data[All, "Position", t]]], {t, 0, 32.2}, BoxRatios -> 1, PlotRange -> All, PlotPoints -> 120, AspectRatio -> 1, PlotStyle -> Yellow, MaxRecursion -> 9, Axes -> False, ImageSize -> 445, Boxed -> False, SphericalRegion -> True, Background -> Black], Graphics3D[{Red, PointSize[.03], Point[{{-(Sqrt[3]/2), -(1/2), 0}, {Sqrt[3]/2, -(1/2), 0}, {0, 1, 0}}], Opacity[.35], Green, Sphere[]}]]
Out[6]=

Plot orbits from different families:

In[7]:=
Table[With[{parameters = ResourceData[ "Collisionless Periodic Orbits in the Free-fall Three-body Problem"][n]}, Module[{data = NBodySimulation[ "InverseSquare", {<|"Mass" -> parameters["Mass"][[1]], "Position" -> {-.5, 0}, "Velocity" -> {0., 0.}|>, <|"Mass" -> parameters["Mass"][[2]], "Position" -> {.5, 0}, "Velocity" -> {0., 0.}|>, <|"Mass" -> parameters["Mass"][[3]], "Position" -> parameters["Position"], "Velocity" -> {0., 0.}|>}, parameters["Period"]]}, ParametricPlot[ Evaluate[data[All, "Position", t]], {t, 0, parameters["Period"]}]]], {n, {"F1(1,1,1)", "F1(1,0.8,0.8)", "F1(1,0.8,0.6)", "F1(1,0.8,0.4)", "F1(1,0.8,0.2)", "F1(0.6,0.8,1)", "F1(1,0.412,0.267)", "F1(0.393,0.288,0.662)", "F1(1,0.1,0.01)"}}]
Out[7]=

Enrique Zeleny, "Collisionless Periodic Orbits in the Free-fall Three-body Problem" from the Wolfram Data Repository (2021)  

Data Resource History

  • Date Created:

Source Metadata

  • Source: http://numericaltank.sjtu.edu.cn/free-fall-3b/free-fall-3b-movies.htm
  • Citation:
    • X. Li and S. Liao, "Collisionless periodic orbits in the free-fall three-body problem", 2018. arXiv:1805.07980

Data Downloads

Publisher Information

  • Prepared for the Wolfram Data Repository By: Enrique Zeleny
  • Publisher of Record: Enrique Zeleny