Medical, Pharma, Engineering, Science, Technology and Business

**Kryukov N and Oks E ^{*}**

Physics Department, Allison Lab, Auburn University, Auburn, USA

- *Corresponding Author:
- Oks E

Physics Department, 206 Allison Lab

Auburn University, Auburn, AL 36849, USA

**Tel:**3348444362

**Fax:**3348444613

**E-mail:**[email protected]

**Received Date:** March 23, 2017;** Accepted Date:** April 22, 2017;** Published Date:** April 27, 2017

**Citation:** Kryukov N, Oks E (2017) Conic-Helical Motion in the Three-Body Problem: Star-Planet-Moon Systems and Relativistic Effects in Binary-Star-Planet Systems. J Astrophys Aerospace Technol 5:144. doi: 10.4172/2329-6542.1000144

**Copyright:** © 2017 Kryukov N, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Astrophysics & Aerospace Technology

Previously it was shown analytically that it is possible for a planet around a binary star to have stable or metastable conic-helical orbits, whose axis of symmetry coincides with the interstellar axis. That study was performed in frames of the nonrelativistic classical mechanics. In the present paper, first, we extend that study to star-planet-moon systemsalso in frames of the nonrelativistic classical mechanics. We complement analytical results by simulations showing that the moon can have practically stable conic-helical orbits around the planet, the average plane of the orbits being perpendicular to the axis connecting the planet and the star. Second, we extend that study to the relativistic classical mechanics. We show that relativistic effects can become significant in conic-helical orbits of a planet around a binary star for the situations where the mass of the planet is relatively small (such planets are so-called planetoids). Again, we complement analytical results by simulations showing that the planet can have relatively stable conic-helical orbits around the lighter star, the average plane of the orbits being perpendicular to the interstellar axis.

Celestial mechanics; Planets and satellites: General; Binaries: General; Planetary systems

The subject of the study is the following variation of the general three-body problem in astrophysics. There are two relatively heavy bodies rotating around their barycentre, one of the bodies being significantly heavier than the other one. In addition, there is a third, much lighter body whose orbit is relatively close to the lighter of the two heavy bodies, the orbit being not in the plane of the rotation of the two heavy bodies, so that the motion of the entire three-body system is really three-dimensional. One particular kind of such systems are star-planet- moon systems. Another particular kind of such systems are binary-star-planet systems (i.e., a planet around a binary star). For brevity, we call such 3-dimensional 3-body systems “3D3BS”.

Conic-helical orbits in 3D3BS-the astrophysical analogue of 1-electron Rydberg quasi molecules studied previously by Oks [1,2] are of interest for several reasons. First, the possibility of stable (or metastable) conic-helical orbits in 3D3BS is a fundamental problem in its own right, as a relatively new chapter of the centuries-old classical three-body problem. Second, one of its applications-namely, to binarystar- planet systems-is especially significant for the search of an extraterrestrial life. Indeed, according to estimates [3-5], approximately 50% of binary stars could support habitable terrestrial planets within stable orbital ranges.

In paper by Oks [6], referred hereafter as paper I, it was shown that stable (or metastable) conic-helical orbits are possible for 3D3BS. Previously the scope of papers on the binary-star-planet systems was typically limited to the motion of these three bodies in the same plane [3-12].

The result in paper I was obtained by applying the standard general analytical method for systems that can be separated into rapid and slow subsystems. The method was applicable because there are ranges of parameters, specified in paper I, where the primary frequency Ω of the conic-helical motion of the planet in binary-star-planet systems is much greater than the Kepler frequency ω of the star’s revolution around their barycentre. Even with the allowance for the star’s rotation, the trajectory of the planet is still conic-helical. The plane of the quasicircular planetary orbit undergoes relatively small oscillations along the rotating interstellar axis and the radius of the planetary orbit also undergoes relatively small oscillations. The last, but not least: in paper I, there were also obtained positive results concerning the transitability (and thus, detectability) of such planets.

In the present paper, first, we extend that study to star-planetmoon systems-also in frames of the nonrelativistic classical mechanics. We complement analytical results, obtained by the standard accurate approximate method based on the separation of rapid and slow subsystems, by exact simulations showing that the moon can have practically stable conic helical orbits around the planet, the average plane of the orbits being perpendicular to the axis connecting the planet and the star. Second, we extend that study to the relativistic classical mechanics. We show that relativistic effects can become significant in conic-helical orbits of a planet around a binary star for the situations where the mass of the planet is relatively small (such planets are socalled planetoids). Again, we complement analytical results, obtained by the standard accurate approximate method based on the separation of rapid and slow subsystems, by exact simulations showing that the planet can have relatively stable conic helical orbits around the lighter star, the average plane of the orbits being perpendicular to the interstellar axis.

The paper is organized as follows. In Section 2 we apply those nonrelativistic results to star-planet-moon systems and demonstrate the stability of the moon in conic-helical orbits. In Section 3 we present the relativistic study of OBSS and of the stability of planetoids in conic helical orbits. In Section 4 we summarize the results into conclusions.

We consider, as an example, an Earth-like planet around a Sun-like star, so that the ratio of the masses μ/μ' of the planet and the star is 3 × 10^{–6}. The separation between the planet and the star is R=1 AU. The planet has a moon.

**Figure 1** shows the ratio Ω/ω of the frequency Ω of the moon revolution around the planet-star axis to the frequency ω of the planet rotation around the star versus the scaled projection w=z/R of the average plane of the moon’s orbit on the planet-star axis. We use the cylindrical coordinates with the z-axis along the planet-star axis in the direction from the planet to the star. It is seen that in the range of 0<w<10^{-4} presented in **Figure 1**, the ratio Ω/ω >> 1, so that the separation of rapid and slow subsystems is justified.

**Figure 2** presents the scaled amplitude δw_{0} of the oscillations of the moon’s orbit along the planet-star axis versus the scaled coordinate w=z/R. **Figure 3** similarly shows the scaled amplitude δv0 of the oscillations of the scaled radius v=ρ/R of the moon’s orbit in the plane perpendicular to the planet-star axis versus the scaled coordinate w. It is seen that these oscillations are small.

In addition to the above analytical results, obtained by the standard accurate approximate method based on the separation of rapid and slow subsystems, we also performed exact simulations of the moon’s motion, as the moon revolves around the planet-star axis while this axis rotates with the Kepler frequency corresponding to the star-planet twobody problem. **Figure 4** shows the dependence of m, the scaled (by R) y-coordinate in the Cartesian righthanded xyz-system with the z-axis being the planet-star axis, on the scaled time τ defined in Equation (A.4) from Appendix, for the initial value of the scaled projection of the average plane of the moon’s orbit on the planet-star axis w_{0}=10^{-4}, i.e., 150000 km, which is of the same order of magnitude as the distance between the Earth and the Earth’s Moon.

**Figure 5** similarly presents the dependence of n, the scaled x-coordinate in the Cartesian right-handed xyz-system with the z-axis being the planet-star axis, on the scaled time τ defined in Equation (A.4), for the initial value of the scaled projection of the average plane of the moon’s orbit on the planet-star axis w_{0}=10^{-4}.

**Figure 6** shows the dependence of w, the scaled planet-star axial coordinate of the moon (co-ordinate in the above-mentioned Cartesian system), on the scaled time τ defined in Equation (A.4) for the initial value of the scaled projection of the average plane of the moon’s orbit on the planet-star axis w_{0}=10^{–4}.

It is clearly seen from **Figures 4-7** that the conic-helical orbit of the moon in the star-planet-moon system is stable. Thus, under the influence of a relatively distant star, the plane of the moon’s orbit could orient itself to be perpendicular to the planet-star axis. The trajectory of the moon becomes conic-helical.

The relativistic force acting on the planet (scaled by its mass, i.e., the force divided by the mass of the planet) is given by the formula

(1)

where α V are the acceleration and the velocity of the planet, respectively, in the inertial reference frame, and The additional acceleration given by Equation (A.1) from Appendix can be substituted into (1) as the last factor, with V as the time derivative of the radius-vector in the rotating reference frame.

(2)

where the z-axis is along the interstellar axis in the direction from the lighter star of the mass μ to the heavier star of the mass μ'. After the calculation, we obtain the following values for the components of the force:

(3)

(here β =V/c). Since ω << Ω, the dominant term is the z-projection of the force (given by the third equation in (3)). By analogy with the non-relativistic case from paper I, we derive the small oscillations of the scaled axial w and radial v coordinates about their equilibrium values, due to the dominant force term (compare to Equation (A.3)):

(4)

where is the equilibrium value of ρ/R and the tilde above ω and Ω means the scaling by multiplying by (R^{3}/Z)^{1/2}. We can now find the conditions, under which the amplitude of the oscillations is small while the condition ω << Ω stays valid. In the relativistic case, the Hamiltonian, divided by the mass m of the planet, h=H/m for the 3D3BS case is given by

(5)

where G is the gravitational constant. In a circular state, the components of the linear momentum in the cylindrical coordinates p_{z}=p_{ρ}=0 and |p_{φ}|/m=const=L. Using the scaling

(6)

we write the scaled energy for the circular state:

(7)

For the relativistic motion one has:

(8)

Using the scaled radius v=ρ/R and the first formula in (6), we can find the speed of the planet in the relativistic case in the units of the speed of light:

(9)

p=v^{2}, v=ρ/R

From the first and the third formulas in (6), we get â„“^{2}=α/r^{2}, where α=(Z/(cL))^{2}. Thus,

(10)

The equilibrium values for p and r can be obtained from differentiating (7) with respect to w and v. The first differentiation yields the same relation as in the non-relativistic case, so the equilibrium value of p is the squared right-hand side of the second line in (4). The second differentiation (with respect to v), with the later substitution of the equilibrium value of v (or p), yields the equilibrium value of â„“, which is related to r by â„“^{2}=α/r^{2} , as mentioned above. Using the substitution [13]

(11)

which significantly simplifies the formulas in the two-Coulombcenter problem, we find the speed of the planet in the circular state (in the units of the speed of light):

(12)

From (8) with the substituted value of V=Ωρ, by using the first relation in (6), the last relation in (9), and the formula â„“^{2}=α/r^{2}, we find

(13)

and its scaled counterpart

(14)

As given in (5), and measuring now the mass of the star μ in the units of the mass of the Sun (in distinction to the nonrelativistic case in Section 2, where the star masses were measured in units of the mass of the planet), we get

(15)

where Substituting (15) into (14), we obtain

(16)

where the quantity (which is one half of the Schwarzschild radius of the Sun and is approximately equal to 147700 cm). As defined above, α^{1/2}=Z/(cL), and substituting the scaling relation for L from (6), â„“^{2}=α/r^{2} and (15), we have

(17)

We substitute (17) into (12), and then substitute the resulting equation for β and the solution for α into (16), obtaining the equation for the scaled frequency of the revolution of the planet that only depends on the ratio R/(μs) and γ (or w-see (11)). Finally, from (4), and because the scaled Kepler frequency is á¿¶=(1 + b)^{1/2} the amplitude of the small oscillations about the equilibrium on the w-axis and v-axis are

(18)

For deriving more explicit expressions for the amplitude of these small oscillations in the relativistic case, we use Equation (16) for the scaled frequency of the revolution of the planet, and Equation (12) for β, and Equation (A.4) for the scaled eigenfrequencies ω_{+} and ω_{–}, and the equilibrium value for p=v_{0}^{2} (v_{0} being given in Equation (4)), and the following equilibrium value of r:

(19)

As a result, we arrive to the following final formula for the amplitudes of the small oscillations of the planet on the w-axis and the v-axis in the relativistic case:

(20)

In turns out that the amplitude in the relativistic case is the same as in the non-relativistic case.

We checked that the ratio of the frequencies of the revolution of the planet Ω and the Kepler frequency ω of the revolution of the stars about their barycentre is much greater than 1. By using the previously found values of Ω and ω and substituting the equilibrium values of p, r, α and β, we derive the formula for the ratio k of the frequencies depending on the axial coordinate for the given interstellar distance R and the star masses μ and μ':

(21)

From Equation (21) it can be found out that k in the relativistic case is equal to the nonrelativistic k_{NR} multiplied by (1-β^{2})^{1/2}. **Figure 8** shows the ratio k/k_{NR} versus the scaled radius of the orbit v, for the masses of the stars μ=1 and μ'=100 (in the units of the mass of the Sun) and the interstellar distance R=100 a.u. It is seen that the relativistic effects become significant when v ~ 10^{–9} or smaller.

**Figure 8:** The ratio of the relativistic and non-relativistic k (each k being the ratio of the frequencies of the revolution of the planet Ω and the Kepler frequency ω of the revolution of the stars about their barycenter) versus the scaled radius of the planetary orbit v=ρ/R, for the masses of the stars μ=1 and μ'=100 (in the units of the mass of the Sun) and the interstellar distance R=100 a.u.

The range of v=(1 ÷ 2) × 10^{-9} for R=100 a.u. corresponds to the range of the radius of the planetary orbit ρ=(15 ÷ 30) km. Since the radius of the planet should be smaller than ρ, in this case it should be a planetoid.

**Figure 9** shows the dependence k(v) for the example where the mass of the lighter star μ=1 and the heavier star μ'=100 (in the units of the mass of the Sun) and for the interstellar distance R=100 a.u. (wide binary system). It is seen that the ratio of the frequencies is greater than 10^{12} for v < 2 × 10^{–9}, which means that such a system would be stable for a very long time.

**Figure 9:** The ratio of the revolution frequency of the planet and the Kepler frequency of the rotation of the stars versus the scaled orbit radius v=ρ/R in the relativistic case, for the example where the mass of the lighter star μ=1 and the heavier star μ'=100 (in the units of the mass of the Sun) and the interstellar distance R=100 a.u.

In **Figure 10** we present the plot of the scaled amplitude δ*w*_{0} of the axial oscillations of the planetary orbit from Equation 20 versus the scaled radial coordinate v for the mass of the lighter star μ=1 and the heavier star μ'=100, in the units of the mass of the Sun. It is seen that the amplitude of the oscillations is really much smaller-about quadrillion (!) times smaller-than the interstellar distance.

In addition to the above analytical results, obtained by the standard accurate approximate method based on the separation of rapid and slow subsystems, we also performed exact simulations of the planetary motion, as the planet revolves around the interstellar axis while this axis rotates with the Kepler frequency corresponding to the binary star twobody problem. **Figure 11** shows the dependence of m, the scaled (by R) y-coordinate of the planet in the Cartesian xyz-system with z-axis being the interstellar axis, on the scaled time τ defined in Equation (A.4) from Appendix, for the initial value of the scaled projection of the average plane of the planetary orbit on the interstellar axis w_{0}=10^{-25}, i.e., practically coinciding with the position of the lighter star in the binary system.

**Figure 12** similarly presents the dependence of n, the scaled (by R) x-coordinate of the planet in the above-mentioned Cartesian xyzsystem, on the scaled time τ defined in Equation (A.4), for the same initial conditions.

**Figure 13 **shows the dependence of w, the scaled z-coordinate of the planet on scaled time τ defined in Equation (A.4), for the same initial conditions. We remind that that the z-axis is along the interstellar axis.

It is seen from **Figures 11-13** that the conic-helical orbit of the planet in the binary-star-planet system is relatively stable. Thus, under the influence of the heavier star in the binary system, the plane of the planetary orbit, being at the lighter star, could orient itself to be perpendicular to the interstellar axis. The trajectory of the planet becomes conic-helical.

We extended the results of paper I to star-planet-moon systems-in frames of the nonrelativistic classical mechanics. We complemented analytical results, obtained by the standard accurate approximate method based on the separation of rapid and slow subsystems, by exact simulations. The exact simulations demonstrated that the moon can have practically stable conic-helical orbits around the planet, the average plane of the orbits being perpendicular to the axis connecting the planet and the star. It is counterintuitive that under the influence of a relatively distant star, the plane of the moon’s orbit could orient itself to be perpendicular to the planet-star axis and the orbit becomes conic-helical.

We also extended the nonrelativistic study from paper I to the relativistic classical mechanics (and also corrected some printing errors from paper I). We showed that relativistic effects in binarystar- planet systems can become significant for the situations where the mass of the planet is relatively small (such planets are so-called planetoids). In this case, the conic helical trajectory of the planetoid encloses the star of the lighter mass, the average plane of the planetary orbit practically coinciding with the position of the star of the lighter mass. It is counterintuitive that a distant heavier star can cause such a modification of the planetoid orbit around the lighter star. We showed that in this case the frequency of the rapid subsystem (i.e., the frequency of the planetoid revolution around the lighter star) exceeds the frequency of the slow subsystem (i.e., the Kepler frequency of the star’s rotation around their barycentre) by more than a trillion times. This means that the standard analytical method of separating rapid and slow subsystems, which we used, is well justified and that this configuration will remain stable for a very long time [14,15].

Finally, we complemented analytical results, obtained by the standard accurate approximate method based on the separation of rapid and slow subsystems, by exact simulations. The simulations confirmed that the planet can have relatively stable conic-helical orbits around the lighter star, the average plane of the orbits being perpendicular to the interstellar axis.

**Appendix: corrections of misprints/errors in some formulas from paper I and new results based on the corrected formulas**

In this Appendix we present corrections to three formulas from paper I. These are three formulas that are used in the main text of the present paper. Because of misprints/errors in these formulas in paper I, it is impossible to simply refer in the main text of the present paper to the corresponding formulas from paper I, which is why it is necessary to present here the corrected formulas numbered below as (A.1), (A.2), and (A.3).

The analysis in paper I was performed in a rotating frame. The frame was rotating at the same frequency as the Kepler frequency ω of the stars rotation. In this frame, there is an additional force (see, e.g., Landau & Lifshitz 1960; Goldstein 1980): [14,15]

(a.1)

The order of magnitude of the velocity v in Equation (A.1) is related to the average radius ρ of the planetary orbit as v ~ Ωρ. Here Ω is the main frequency of the planetary motion. Therefore, Equation (A.1) can be approximated as

(a.2)

where is the unit vector along the z-axis, which is the interstellar axis. The motion in the zρ-plane corresponded to a 2D-oscillator driven by the force given by Equation (A.2). In the scaled coordinates w = z/R, *v*=ρ/R, the solution for the oscillation amplitudes is

(a.3)

where,

(a.4)

are the scaled (dimensionless) time, as well as the scaled frequencies ω and Ω, respectively. The similarly scaled eigenfrequencies ω_{+} and ω_{-} of the above two-dimensional oscillator are

(a.5)

The equilibrium value of the scaled radial coordinate in Equations (A.3, A.5) was defined in Equation (4).

Based on the corrected formula (A.3) for the small oscillations of this driven two-dimensional oscillator, below we present new illustrations demonstrating the ranges of the stability of the conic-helical orbits for various parameters of this nonrelativistic 3-body problem. **Figure 14** shows the scaled amplitude of the oscillations of the planetary orbit along the interstellar axis versus the scaled projection w = z/R of the average plane of the planetary orbit on the interstellar axis, for three values of the ratio b=μ'/μ of the stellar masses. It is seen that δw_{0} << 1 (i.e., the amplitude of the oscillations is much smaller than the interstellar distance) for b greater or of the order of 10. For the validity condition Ω >> ω of this result to be satisfied with a large margin of “safety”, the average plane of the planetary orbit should be very close to the star of the smaller mass (and the closer it is in the range of w < 0.05, the smaller is the amplitude δw_{0}, as seen from **Figure 14**). As long as the primary frequency Ω of the planet revolution about the interstellar axis exceeds by many orders of magnitude the Kepler frequency ω of stars rotation, the conic-helical planetary orbit would remain stable for a very long time. (Rigorously speaking, the planet is in a metastable state, which is an analogy with atomic/molecular systems: they have metastable states living by many orders of magnitude longer than other states of the system.)

**Figure 14:** The scaled amplitude δw_{0} of the oscillations of the planetary orbit along the interstellar axis versus the scaled projection w=z/R of the average plane of the planetary orbit on the interstellar axis, for three values of the ratio b of the stellar masses: b=100 (solid line), b=30 (dashed line), and b=10 (dotted line).

**Figure 15** shows the scaled amplitude δ*v*=2v_{0}(w,b)ω_{s}f_{p}(sin 2α)/|ω_{+}^{2}-ω_{-}^{2}|of the oscillations of the radius of the planetary orbit versus the scaled projection w = z/R of the average plane of the planetary orbit on the interstellar axis, for three values of the ratio b=μ'/μ of the stellar masses. It is seen that for the range of b greater or of the order of 10 (required for having δw_{0} << 1), we have also δ*v *<< 1.

**Figure 15:** The scaled amplitude δv_{0} of the oscillations of the radius of the planetary orbit versus the scaled projection w=z/R of the average plane of the planetary orbit on the interstellar axis, for three values of the ratio b of the stellar masses: b=100 (solid line), b=30 (dashed line), and b=10 (dotted line).

- Oks E (2000) Are crossings of energy levels and charge exchange really quantum phenomena? Phys Rev Lett 85: 2084-2087.
- Oks E (2000) A purely classical description of crossings of energy levels and of spectroscopic signatures of charge exchange. J Phys B Atom Mol Opt Phys 33: 3319-3330.
- Quintana EV, Lissauer JJ (2010) Terrestrial planet formation in binary star systems. In: Naghighipour N (Ed.) Planets in Binary Star Systems 265-284 Springer Dordrecht.
- Fatuzzo M, Adams FC, Gauvin R, Proszkow EM (2006) A statistical stability analysis of earth-like planetary orbits in binary systems. PASP 118: 1510-1527.
- David EM, Quintana EV, Fatuzzo M, Adams FC (2003) Dynamical stability of earth-like planetary orbits in binary systems. PASP 115: 825-836
- Oks E (2015) Stable conic-helical orbits of planets around binary stars: Analytical results. Ap J 804: 106.
- Kaib NA, Raymond SN, Duncan M (2013) Planetary system disruption by Galactic perturbations to wide binary stars. Nature 493: 381-384.
- Desidera S, Barbieri M (2007) Properties of planets in binary systems. The role of binary separation. Astron Astrophys 462: 345-353.
- Quintana EV, Adams FC, Lissauer JJ, Chambers JE (2007) Terrestrial planer formation around individual stars within binary star systems Ap J 660: 807-822.
- Turrini D, Barbieri M, Marzari F, Thebault P, Tricarico P(2005) Planetary formation and orbital stability in binary star systems. Mem S A It Suppl 6: 172-177.
- Holman MJ, Wiegert PA (1999) Long term stability of planets in binary systems. Astron J 117: 621-628.
- Innanen KA, Zheng JQ, Mikkola S, Valtonen MJ (1997) The Kozai mechanism and the stability of planetary orbits in binary star systems. Astron J 113: 1915-1919.
- Kryukov N, Oks E (2013) Applications of the classical two-Coulomb-center systems to atomic/molecular physics. Intern Rev Atom Mol Phys 4: 121-177.
- Landau LD, Lifshitz EM (1960) Mechanics. Pergamon Oxford.
- Goldstein H (1980) Classical mechanics. Addison-Wesley Reading MA.

Select your language of interest to view the total content in your interested language

- Aerodynamic Flow Control
- Aerodynamics
- Aeroelasticity
- Aerospace Dynamics
- Aerospace Engineering techniques
- Air Safety
- Aircraft
- Aircraft Flight Mechanics
- Astrodynamics
- Aviation
- Avionics
- Carbon Nanotubes
- Computational Fluid Dynamics
- Cosmoparticle Physics
- Detectors and Optical Sensors
- Flight Dynamics
- Fusion Plasmas
- Luminescence
- Nuclear Physics
- Parallel Programming
- Quantum Gravity
- Rocketry
- Soil Physics
- Solid State Plasmas
- Space
- Space Plasmas
- Stellar Spectroscopy
- Suprnova Remnants
- Surface Physics
- Unmanned-Vehicles

- 5th International Conference and Exhibition on
**Mechanical**&**Aerospace Engineering**

October 02-04, 2017 Las Vegas USA - 2nd Internetional Conference on
**Astrophysics**and**Particle Physics**

November 13-15, 2017 San Antonio, Texas, USA

- Total views:
**155** - [From(publication date):

June-2017 - Sep 20, 2017] - Breakdown by view type
- HTML page views :
**129** - PDF downloads :
**26**

Peer Reviewed Journals

International Conferences 2017-18