Equations of motion of elliptic restricted problem of three bodies with variable mass

The differential equations of motion of the elliptic restricted problem of three bodies with decreasing mass are derived. The mass of the infinitesimal body varies with time. We have applied Jeans' law and the space-time transformation of Meshcherskii. In this problem the space-time transformation is applicable only in the special case whenn=1,k=0,q=1/2. We have applied Nechvile's transformation for the elliptic problem. We find that the equations of motion of our problem differ from that of constant mass only by a small perturbing force.     

A note on the elliptic restricted three-body problem

This work considers periodic solutions, arc-solutions (solutions with consecutive collisions) and double collision orbits of the plane elliptic restricted problem of three bodies for =0 when the eccentricity of the primaries,e _p , varies from 0 to 1. Characteristic curves of these three kinds of solutions are given.     

The elliptic restricted problem at the 3 : 1 resonance

Four 3 : 1 resonant families of periodic orbits of the planar elliptic restricted three-body problem, in the Sun-Jupiter-asteroid system, have been computed. These families bifurcate from known families of the circular problem, which are also presented. Two of them, I _c , II _c bifurcate from the unstable region of the family of periodic orbits of the first kind (circular orbits of the asteroid) and are unstable and the other two, I _e , II _e , from the stable resonant 3 : 1 family of periodic orbits of the second kind (elliptic orbits of the asteroid). One of them is stable and the other is unstable. All the families of periodic orbits of the circular and the elliptic problem are compared with the corresponding fixed points of the averaged model used by several authors. The coincidence is good for the fixed points of the circular averaged model and the two families of the fixed points of the elliptic model corresponding to the families I _c , II _c , but is poor for the families I _e , II _e . A simple correction term to the averaged Hamiltonian of the elliptic model is proposed in this latter case, which makes the coincidence good. This, in fact, is equivalent to the construction of a new dynamical system, very close to the original one, which is simple and whose phase space has all the basic features of the elliptic restricted three-body problem.     

Out-of-plane motion about libration points: Nonlinearity and eccentricity effects

Out-of-plane motion about libration points is studied within the framework of the elliptic restricted three-body problem. Nonlinear motion in the circular restricted problem is given to third order in the out-of-plane amplitudeA _z by Jacobi elliptic functions. Linear motion in the elliptic problem is studied using Mathieu's and Hill's equations. Additional terms needed for a complete third-order theory are found using Lindsted's method. This theory is constructed for the case of collinear libration points; for the case of triangular points, a third-order nonlinear solution is given separately in terms of Jacobi elliptic functions.     

On the triangular libration points in photogravitational restricted three-body problem with variable mass

This paper investigates the triangular libration points in the photogravitational restricted three-body problem of variable mass, in which both the attracting bodies are radiating as well and the infinitesimal body vary its mass with time according to Jeans’ law. Firstly, applying the space-time transformation of Meshcherskii in the special case when q=1/2, k=0, n=1, the differential equations of motion of the problem are given. Secondly, in analogy to corresponding problem with constant mass, the positions of analogous triangular libration points are obtained, and the fact that these triangular libration points cease to be classical ones when α≠0, but turn to classical L ₄ and L ₅ naturally when α=0 is pointed out. Lastly, introducing the space-time inverse transformation of Meshcherskii, the linear stability of triangular libration points is tested when α>0. It is seen that the motion around the triangular libration points become unstable in general when the problem with constant mass evolves into the problem with decreasing mass.     

On the dynamics in the asteroids belt. Part II: Detailed study of the main resonances

The general theory exposed in the first part of this paper is applied to the following resonances with Jupiter's motion : 3/2, 2/1, 5/2, 3/1, 7/2, 4/1; these are the most relevant resonances for the asteroids. The whole analysis is performed in the framework of the spatial problem of three bodies, both in the circular and in the elliptic case. The results are also compared with the observed distribution of the asteroids.     

Bifurcations of plane with three-dimensional periodic orbits in the elliptic restricted problem

We show that the procedure employed in the circular restricted problem, of tracing families of three-dimensional periodic orbits from vertical self-resonant orbits belonging to plane families, can also be applied in the elliptic problem. A method of determining series of vertical bifurcation orbits in the planar elliptic restricted problem is described, and one such series consisting of vertical-critical orbits (a _v=+1) is given for the entire range (0,1/2) of the mass parameter . The initial segments of the families of three-dimensional orbits which bifurcate from two of the orbits belonging to this series are also given.     

Steady state obliquity of a rigid body in the spin–orbit resonant problem: application to Mercury

In this paper, we consider the elliptic collinear solutions of the classical n-body problem, where the n bodies always stay on a straight line, and each of them moves on its own elliptic orbit with the same eccentricity. Such a motion is called an elliptic Euler–Moulton collinear solution. Here we prove that the corresponding linearized Hamiltonian system at such an elliptic Euler–Moulton collinear solution of n-bodies splits into \((n-1)\) independent linear Hamiltonian systems, the first one is the linearized Hamiltonian system of the Kepler 2-body problem at Kepler elliptic orbit, and each of the other \((n-2)\) systems is the essential part of the linearized Hamiltonian system at an elliptic Euler collinear solution of a 3-body problem whose mass parameter is modified. Then the linear stability of such a solution in the n-body problem is reduced to those of the corresponding elliptic Euler collinear solutions of the 3-body problems, which for example then can be further understood using numerical results of Martínez et al. on 3-body Euler solutions in 2004–2006. As an example, we carry out the detailed derivation of the linear stability for an elliptic Euler–Moulton solution of the 4-body problem with two small masses in the middle.     

Recent progress in the theory and application of symplectic integrators

In this paper various aspect of symplectic integrators are reviewed. Symplectic integrators are numerical integration methods for Hamiltonian systems which are designed to conserve the symplectic structure exactly as the original flow. There are explicit symplectic schemes for systems of the formH=T(p)+V(q), and implicit schemes for general Hamiltonian systems. As a general property, symplectic integrators conserve the energy quite well and therefore an artificial damping (excitation) caused by the accumulation of the local truncation error cannot occur. Symplectic integrators have been applied to the Kepler problem, the motion of minor bodies in the solar system and the long-term evolution of outer planets.     

Equations of motion of the restricted problem of three bodies with variable mass

In this paper we have deduced the differential equations of motion of the restricted problem of three bodies with decreasing mass, under the assumption that the mass of the satellite varies with respect to time. We have applied Jeans law and the space time transformation contrast to the transformation of Meshcherskii. The space time transformation is applicable only in the special casen=1,k=0,q=1/2. The equations of motion of our problem differ from the equations of motion of the restricted three body problem with constant mass only by small perturbing forces.     

New stacked central configurations for the planar 5-body problem

A stacked central configuration in the n-body problem is one that has a proper subset of the n-bodies forming a central configuration. In this paper we study the case where three bodies with masses m ₁, m ₂, m ₃ (bodies 1, 2, 3) form an equilateral central configuration, and the other two with masses m ₄, m ₅ are symmetric with respect to the mediatrix of the segment joining 1 and 2, and they are above the triangle generated by {1, 2, 3}. We show the existence and non-existence of this kind of stacked central configurations for the planar 5-body problem.     

On a global expansion of the disturbing function in the planar elliptic restricted three-body problem

Starting with a simple Taylor-based expansion of the inverse of the distance between two bodies, we are able to obtain a series expansion of the disturbing function of the three-body problem (planar elliptic case) which is valid for all points of the phase space outside the immediate vicinity of the collision points. In particular, the expansion is valid for very high values of the eccentricity of the perturbed body. Furthermore, in the case of an interior mean-motion resonant configuration, the above-mentioned expression is easily averaged with respect to the synodic period, yielding once again a global expansion of (R) valid for very high eccentricities.Comparisons between these results and the numerically computed exact function are presented for various resonances and values of the eccentricity. Maximum errors are determined in each case and their origin is established. Lastly, we discuss the applicability of the present expansion to practical problems.     

Newtonian and relativistic periodic orbits around two fixed black holes

We compare families of simple periodic orbits of test particles in the Newtonian and relativistic problems of two fixed centers (black holes). The Newtonian problem is integrable, while the relativistic problem is highly non-integrable.The orbits are calculated on the meridian plane through the fixed centersM ₁ (atz=+1) andM ₂ (atz=–1) for energies smaller than the escape energyE=1. We use prolate spheroidal coordinates (, , =const) and also the variables =cosh and =–cos . The orbits are inside a curve of zero velocity (CZV). The Newtonian orbits are also limited by an ellipse and a hyperbola, or by two eillipses. There are 3 main types of periodic orbits (1) elliptic type (around both centers), (2) hyperbolic-type, and (3) resonant-type.The elliptic type orbits are stable in the Newtonian case and both stable and unstable in the relativistic case. From the stable orbits bifurcate double period orbits both symmetric and asymmetric with respect to thez-axis. There are also higher order bifurcations. The hyperbolic-type orbits are unstable. The Newtonian resonant orbits are defined by the ratiot _µ/t =n/m of oscillations along and during one period, and they are all marginally unstable. The corresponding relativistic orbits are stable, or unstable. The main families are figure eight orbits aroundM ₁, or aroundM ₂ (3/1 orbits); gamma, or inverse gamma orbits (4/2); higher resonant families 5/1,7/1,...,8/2,12/2,...;, more complicated orbits, like 5/3, and bifurcations from the above orbits. Satellite orbits aroundM ₁, orM ₂, and their bifurcations (e.g. double period) exist in the relativistic case but not in the Newtonian case. The characteristics of the various families are quite different in the Newtonian and the relativistic cases. The sizes of the orbits and their stabilities are also quite different in general. In the Appendix we study the various types of straight line orbits and prove that some subcases introduced by Charlier (1902) are impossible.     

A method for stability determinations in the elliptic restricted problem

In the ordinary restricted problem of three bodies, the first-order stability of planar periodic orbits may be determined by means of their characteristic exponents, as derived from the condition of a vanishing determinant for the coefficients of an infinite system of homogenous linear equations associated with the exponential series solutionu, v representing any initially small oscillations about the periodic solutionx, y. In the elliptic restricted problem, periodic solutions are possible only for periods which are equal to, or integral multiples of, the periodP of the elliptic motion of the two primary masses. It is shown that the infinite determinant approach to the determination of the characteristic exponents can be extended to the treatment of superposed free oscillations in the elliptic problem, and that in generaltwo exponents appear in any complete solutionu, v for eachone existing in the corresponding ordinary restricted problem. The value of each exponent depends on a series proceeding in even powers of the eccentricitye of the relative orbit of the two primaries, in addition to its basic dependence on the mass ratio . For stable periodic orbits, the oscillation frequenciesn ₁ (,e ²),n ₂ (,e ²) associated with these two exponents tend, withe0, to certain limiting valuesn ₁ (),n ₂(), which differ from each other by the amount of the frequencyN=2/P of the orbital motion of the primaries. One of the two frequencies, sayn ₁(), is identical with the frequency of the corresponding oscillations in the ordinary restricted problem, while the second one gives rise to oscillations only in the elliptic restricted problem, withe0.The method will be described in more detail, together with its application to two families of small periodie librations about the equilateral points of the elliptic restricted problem (E. Rabe: Two new Classes of Periodic Trojan Librations in the Elliptic Restricted Problem and their Stabilities) in theProceedings of the Symposium on Periodic Orbits, Stability and Resonances, held at the University of São Paulo, Brasil, 4–12 September, 1969.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 17–23, 1969.     

Stability of periodic resonance-orbits in the elliptic restricted 3-body problem

Results of families of periodic orbits in the elliptic restricted problem are shown. They are calculated for the mass ratios =0.5 and =0.1 for the primary bodies and for different values of the eccentricity of the orbit of the primaries which is the second parameter. The case =0.5 is also a good model for planetary orbits in binaries. Finally we show detailed stability diagrams and give results according to the stability classification of Contopoulos.     

Linear stability of the triangular equilibrium points in the photogravitational elliptic restricted problem,II

The linear stability of the triangular equilibrium points in the photogravitational elliptic restricted three-body problem is examined and the stability regions are determined in the space of the parameters of mass, eccentricity, and radiation pressure, in the case of equal radiation factors of the two primaries. The full range of values of the common radiation factor is explored, from the gravitational caseq ₁ =q ₂ =q = 1 down to the critical value ofq = 1/8 at which the triangular equilibria disappear by coalescing on the rotating axis of the primaries. It is found that radiation pressure exerts a significant influence on the stability regions. For certain intervals of radiation values these regions become qualitatively different from the gravitational case as well as the solar system case considered in Paper I. There exist values of the common radiation factor, in the range considered, for which the triangular equilibrium points are stable for the entire range of mass distribution among the primaries and for large eccentricities of their orbits.     

2D calculations of the collapse of magnetized gas cloud

All the families of planar symmetric simple-periodic orbits of the photogravitational restricted plane circular three-body problem, are determined numerically in the case when the primaries are of equal mass and radiate with equal radiation factors (q ₁=q₂=q). We obtain a global view of the possible patterns of periodic three-body motion while the full range of values of the common radiation factor is explored, from the gravitational case (q=1) down to near the critical value at which the triangular equilibria disappear by coalescing with the inner equilibrium pointL ₁ on the rotating axis of the primaries. It is found that for large deviations of its value from the gravitational case the radiation factorq can have a strong effect on the structure of the families.     

A self-consistent world model

The field equations of the generalized field theory constructed by Mikhail and Wanas have been applied to a well-established geometrical structure given earlier by H. P. Robertson in connection with the cosmological problem. A unique solution, representing a specified expanding Universe (withq ₀=0, ₀=0.75,k=–1) has been obtained. The model obtained has been compared with cosmological observations and with FRW-models of relativistic cosmology. It has been shown that the suggested model is free of particle horizons. The existence of singularities has been discussed.The two cases, when the associated Riemannian-space has a definite or indefinite metric are considered. The case of indefinite metric with signature (+ – – –) is found to be characterized byk=–1, while the case of +ve definite metric is characterized byk=+1. Apart from that difference, the two cases give rise to the same cosmological parameters. It has been shown that energy conditions are satisfied by the material contents in both cases.     

Particle Dynamics in a Maxwell’s Ring-Type Configuration with a Radiating Central Primary

Maxwell’s ring-type configuration (i.e. an N-body model where the ν = Ν − 1 bodies have equal masses and are located at the vertices of a regular ν-gon while the N-th body with a different mass is located at the center of mass of the system) has attracted special attention during the last 15 years and many aspects of it have been studied by considering Newtonian and post-Newtonian potentials (Mioc and Stavinschi 1998, 1999), homographic solutions (Arribas et al. 2007) and relative equilibrium solutions (Elmabsout 1996), etc. An equally interesting problem, known as the ring problem of (N + 1) bodies, deals with the dynamics of a small body in the combined force field produced by such a configuration. This is the problem we are dealing with in the present paper and our aim is to investigate the variations in the dynamics of the small body in the case that the central primary is also a radiating source and therefore acts on the particle with both gravitation and radiation. Based on the general outlines of Radzievskii’s model, we study the permitted and the existing trapping regions of the particle, its equilibrium locations and their parametric variations as well as the existence of focal points in the zero-velocity diagrams. The distribution of the characteristic curves of families of planar symmetric periodic orbits and their stability for various values of the radiation coefficient of the central body is additionally investigated.     

Photo—gravitational effects in the two—body problem

The effect of the radiation-pressure in the problem of the two bodies is analysed. It is found that such an effect is the same as that of reducing the stellar mass, m, by the factor p (0 ≤ p ≤ 1). The case p < 0 leads to stellar instability and, consequently, it is not considered here. The problem of the mass determination from spectroscopic observations is reviewed and an iterative method is recommended.