THE HYPERBOLIC RESTRICTED THREE-BODY PROBLEM AND APPLICATIONS TO CELESTIAL MECHANICS
HARINI CHANDRAMOULI
Contents 1. Motivation 2 2. The Kepler Problem 2 2.1. One-Body Central Force Problem 3 3. The Restricted 3-Body Problem 5 3.1. The Circular Restricted 3-Body Problem 6 3.2. The Sitnikov Problem 6 4. Hamiltonian Systems 8 4.1. Symplectic Transformations. 8 5. Regularizations 9 5.1. Levi-Civita Regularization. 9 5.2. Moser’s Regularization. 10 5.3. Belbruno’s Regularization. 14 6. The Scattering Problem 14 6.1. Computing Orbital Elements 14 6.2. Sorokovich 15 6.3. Current and Future Work 15 References 20
1 1. Motivation In our solar system, we have a few planets, but none as large as Jupiter. As Jupiter orbits our sun, it travels along an elliptic path in a plane at an inclination to a fixed reference plane. The magnitude of Jupiter’s size means that any effect that happens on Jupiter’s inclination or the eccentricity of its orbit will have a ripple effect on that of the other planets. It has been shown that the inclination and the eccentricity of a planet have the greatest effect on the climate of that planet. Inclination and eccentricity are two examples of what are called orbital elements, parameters which define the orbit of a planet around the sun.
Consider a star passing through our solar system “close enough” to disturb these elements of Jupiter. If a star passed within the heliosphere of our sun, and had large enough mass, what effect would it have on Jupiter? If we want to model such a situation we need to take into account three bodies: our Sun, the star passing through the system, and Jupiter. Here we have a three-body problem, a case of the famous n-body problem.
We want to study how Jupiter’s orbital elements change as time progresses from the past to the future. Assuming that our star doesn’t pass close enough to capture Jupiter, we hypothesize that the effect of the star could be felt while the star is nearby but disappear with the star. If these effects are large enough to change the climate on other planets, such as Earth, then we could try to see if any past climate data matches with such an event happening.
In this paper, we will first explore the Kepler problem and the restricted 3-body problem as mathematical ways to phrase the scattering problem. We will look at an example of a restricted 3-body problem to get a feel for how such a problem could be solved. We will also look at some properties of Hamiltonian systems, as we are working with one, and regularization methods as ways to deal with potential singularities. Lastly, we will utilize all this information to make some progress on our scattering problem.
2. The Kepler Problem Before the star appears, and after the star passes through the system, we have only two bodies: Jupiter and the Sun. Before studying what happens in the system with the star, we want to look into the system of celestial bodies. In classical mechanics, the two-body problem is to determine the position and speed of two bodies interacting with each other given their masses, initial positions, and initial velocities. The gravitational two-body problem is a special case of the two-body problem in which the two bodies interact by a central force, F , that varies in strength as the inverse square of the distance, r, between them.
Let m1 and m2 denote the masses of the two bodies, x1 = (x1, y1), x2 = (x2, y2) their po- p 2 2 sitions, and r = (x1 − x2) + (y1 − y2) the distance between them. By Newton’s law of universal gravitation, the force of attraction between the particles is given by Gm m F = 1 2 , r2 where G = 6.67408 × 10−11m3kg−1s−2 is the gravitation constant. The differential equations for the gravitational two-body problem are given by x − x m x¨ = F · 2 1 1 1 r (2.1) x − x m x¨ = F · 1 2 . 2 2 r 2 where it is assumed that x1(t0), x2(t0), x˙ 1(t0), x˙ 2(t0) are given at some initial time t0.
2.1. One-Body Central Force Problem. It is possible to reduce the gravitational two-body problem above to a one-body central force problem, also known as the Kepler problem. Let q = x2 − x1 and |q| = q. Then, divide the first of Eqs. (2.1) by m1 and the second by m2 and subtract the first from the second. This yields the equation q q (2.2) q¨ = −F · = −G(m + m ) . q 1 2 q3 Here we have what’s called the Kepler problem. Once we find out what q is, we can solve for both x1 and x2. The original two body gravitational problem has been reduced to an equivalent one body problem. The parameter q in the one body problem, called a relative coordinate, represents the distance between the “single body” and the central point. To solve (2.2), we let µ = G(m1 + m2) and take the resulting second order equation and turn it into a system of first order equations. We thus get the system ( q˙ = p (2.3) p˙ = −µq−3q.
Now to solve this, we need to consider both the conservation of angular momentum and conser- vation of energy. The techniques used to solve this problem can be found in many textbooks, including [12].
2.1.1. Conserved Quantities. Angular momentum is the rotational analog of linear momentum. It is defined for a point particle to be the vector q × p. So, we observe
q × p˙ = −µq−3 (q × q) = 0 d =⇒ (q × p) = q × p˙ + p × p = 0. dt Thus we have that q × p = c, a constant vector, and we observe that angular momentum is conserved. We will refer to c as the angular momentum.
In our original gravitational two-body problem (2.1), the gravitational interaction was an in- ternal conservative force and thus we can conclude that energy in (2.3) is constant as well. We can see this mathematically by first taking the second of equations (2.3) and taking the inner product on both sides with p. First, note that since q2 = q · q, we can see by taking the derivative that qq˙ = q · q˙. Thus, dq (2.4) p˙ · p = −µq−3 (q · p) = −µq−3(q · q˙) = −µq−3qq˙ = −µq−2 dt Now, if we integrate both sides of (2.4) and we find that 1 p2 = µq−1 + h, 2 where h is a constant.
Along with angular velocity, there is one other vector which stays constant – the vector rep- resenting the eccentric axis, denoted by e. We will derive it below. Recall that if r 6= 0, then 3 using q2 = q · q, qq˙ = q · q˙, and the vector identity (a × b) × c = (a · c)b − (b · c)a, we get d q qq˙ − qq˙ = dt q q2 q·q q·q˙ q q˙ − q q = q2 (q · q˙)q˙ − (q · q˙)q = q3 (q × q˙) × q = q3 c × q = . q3 Take both sides and multiply by −µ, d q −µ = c × (−µq−3q) dt q d q µ = p˙ × c. dt q If we integrate both sides with respect to t, we get that q (2.5) µ e + = p × c, q where e is the constant of integration.
2.1.2. Solution of the Kepler Problem. To solve the Kepler problem, we must consider two cases from above: when c = 0 and when c 6= 0. q When c = 0, from (2.5) we get that = −e. Thus, e lies on the line of motion and has length 1. q Now let us consider what happens when c 6= 0. We know that q · c = 0 and so this im- plies that e · c = 0, that is, e and c are perpendicular and so e lies in the plane of motion. Taking the inner product of both sides of (2.5) with q gets us (2.6) µ (e · q + q) = q · p × c = q × p · c = c · c c2 (2.7) =⇒ e · q + q = . µ c2 Suppose e = 0, then we have that q = and the motion is circular with uniform angular µ velocity. Now suppose e 6= 0. Let the fixed angle from the x-axis to e be denoted by ω. Then, let (r, θ) represent the polar coordinates of the particle from the positive x-axis. Now define f := θ − ω, which we will call the true anomaly. If we use e as the axis of coordinates, we can represent the particle’s location with the coordinates (r, f) instead. We know that e · r = er cos(f), and so applying this to the equation above, we get [12] c2/µ (2.8) r = . 1 + e cos(f) Equation (2.8) is the solution to the Kepler problem. Kepler’s first law tells us that the particle moves on a conic section of eccentricty e and thus we know that (2.8) must represent a conic section. We recall that 0 < e < 1 represents an ellipse, e = 1 represents a parabola, e > 1 represents a hyperbola, and e = 0 represents a circle. 4 3. The Restricted 3-Body Problem As the star gets closer to Jupiter and the Sun, we introduce a third body into the system. Now we want to work with the three-body problem. The three-body problem is the problem of determining the motions of the three bodies given an initial set of data that specifies the position, mass, and velocity of the bodies at some fixed time. If we say that the position of the masses mi is given by xi, then we get the following three coupled second-order differential equations:
x1 − x2 x1 − x3 x¨1 = −Gm2 3 − Gm3 3 |x1 − x2| |x1 − x3| x2 − x1 x2 − x3 x¨2 = −Gm1 3 − Gm3 3 |x2 − x1| |x2 − x3| x3 − x1 x3 − x2 x¨3 = −Gm1 3 − Gm2 3 . |x3 − x1| |x3 − x2| It has been shown that there is no general analytic solution for the three-body problem by alge- braic expressions and integrals, but we can classify types of solutions that exist. In this paper, we are concerned with solutions near infinity. What kinds of solutions can we see when the third body moves further and further away from the primaries? These solutions can be classified by the behavior of the distance between this third body and the center of mass of the primaries. If this rate is large enough, the distance will approach infinity while the velocity approaches a limit. If this limit is zero, we call the corresponding solution parabolic. If this limit is positive, we call such a solution hyperbolic. If this rate is not large enough to overcome the attraction of the third body to the primaries, the distance will start to decrease at some point and the third mass approaches the primaries once more. Such a solution is called elliptic. An escape orbit is a solution where the distance between the primaries and the third body is bounded in backward time but goes to infinity in forward time while a capture orbit is the exact opposite - the distance between the third body and the primaries is bounded in forward time and infinite in backward time. If the distance is infinite in both forward and backward time, we have a scattering solution.
In our problem, as described in the beginning, we have that the distance of the star to the primaries, Jupiter and the Sun, starts out infinite and it comes closer, lastly going off to infinity again. Our solution is a little different from a scattering solution, seeing as one of the bodies with mass is the one that comes closer and then goes off to infinity.
The restricted three-body problem is a special case of the three-body problem where one of the bodies is said to be so much smaller than than the other two that we say it has negligi- ble mass. So, we study the motion of this small particle moving under the influence of the gravitational attraction of the other two bodies. Thus, the equations become x − x x¨ = −Gm 1 2 1 2 3 |x1 − x2| x2 − x1 (3.1) x¨2 = −Gm1 , |x − x |3 2 1 x3 − x1 x3 − x2 x¨3 = −Gm1 3 − Gm2 3 |x3 − x1| |x3 − x2| where x3 is the particle with mass 0. We quickly notice that the first two equations are just the two-body problem from before, for which we know the solution. The two primaries, the bodies with mass m1 and m2, have a solution that is a conic section, as we learned from the previous section. The third body, with mass m3, simply moves in the gravitational field of the other two bodies. It does not exert influence on the paths of the primaries.
5 3.1. The Circular Restricted 3-Body Problem. We could solve the Kepler problem using classical mechanics but we need to be more clever with the restricted three body problem. While it has been shown that there are periodic solutions and solutions in special cases, a general so- lution does not exist. While the following section of this paper is not immediately relevant to the scattering problem, it provides insightful background into understanding and working with these systems.
The restricted three-body problem is often analyzed in a rotating coordinate system. In a coordinate system rotating with angular velocity c, the two primaries are kept stationary and the third body moves about in the gravitational field. So first, when looking at the solution to the system of equations (3.1), we first consider the case where the two primaries are moving in circular orbits about their center of mass. We can plug this solution in to the third equation of the system (3.1) and solve numerically for the motion of the smaller mass, m3, in the presence of the gravitational field created by m1 and m2 - generally an analytic solution is not possible. 3.1.1. Rotating Coordinate System. We will now move into the rotating coordinate system where m1 and m2 are stationary. To find our new equations, we will first consider xi ∈ C. We will also take the center of mass of the bodies m1 and m2 to be the origin, that is, let m1x1 + m2x2 = 0. The solution to the first two equations of system (3.1) are just circular orbits, and so they are given by iωt iωt x1 = αe and x2 = βe , where α, β ∈ R. We want to rewrite the motion of the third body, described by the third equation of system (3.1) in terms of our rotating coordinates. For this, we need to see how the iωt velocity and acceleration have changed in our new coordinate system. So, let x3 = γe , where now γ = γ(t) ∈ C. Then, we see that the velocity and acceleration change as follows: iωt iωt iωt x˙3 = γ · iωe +γe ˙ = e (γ ˙ + iωγ) iωt iωt iωt 2 x¨3 = e (¨γ + iωγ˙ ) + (γ ˙ + iωγ) · iωe = e γ¨ + 2iωγ˙ − ω γ . In the equation of the second derivative, 2iωγ˙ term represents the Coriolis effect1 and −ω2γ represents the centripetal acceleration2. Setting this second equation equal to the third equation in the system (3.1), we get iωt iωt iωt 2 e (γ − α) e (γ − β) e γ¨ + 2iωγ˙ − ω γ = −Gm1 − Gm2 , |γ − α|3 |γ − β|3 which implies