There are ten,not eight planets

In 1946, E. Sevin postulated the global vibrations of the Sun with a period P ₀ = 1/9 day and a “wavelength” L ₀ = c × P ₀ = 19.24 AU and predicted the tenth planet at a mean distance of 4.0 × L ₀ ≈ 77.0 AU from the Sun (c is the speed of light). The global vibrations of the Sun, precisely with the period of 1/9 day, were actually detected in 1974. Recently, the largest Kuiper Bell object 2003 UB₃₁₃, or Eris, with an orbital semimajor axis ≈ 3.5 × L ₀ ≈ 67.5 AU was discovered. We adduce arguments for the status of Eris as our tenth planet: (i) the object is larger and farther from the Sun than Pluto and (ii) the semimajor axis of Eris agrees well with the sequence of planetary distances that follows from the resonance spectrum of the Solar system dimensions (with the scale L ₀ and for all 11 orbits, including those of Pluto, Eris, and the asteroid belt). We point to a mistake of the Prague (2006) IAU Assembly, which excluded Pluto from the family of planets by introducing a new, highly controversial class of objects—“dwarf planets.”     

Orbital stability of systems of closely-spaced planets

An investigation of the stability of systems of 1 M_⊕ (Earth-mass) bodies orbiting a Sun-like star has been conducted for virtual times reaching 10 billion years. For the majority of the tests, a symplectic integrator with a fixed timestep of between 1 and 10 days was employed; however, smaller timesteps and a Bulirsch-Stoer integrator were also selectively utilized to increase confidence in the results. In most cases, the planets were started on initially coplanar, circular orbits, and the longitudinal initial positions of neighboring planets were widely separated. The ratio of the semimajor axes of consecutive planets in each system was approximately uniform (so the spacing between consecutive planets increased slowly in terms of distance from the star). The stability time for a system was taken to be the time at which the orbits of two or more planets crossed. Our results show that, for a given class of system (e.g., three 1 M_⊕ planets), orbit crossing times vary with planetary spacing approximately as a power law over a wide range of separation in semimajor axis. Chaos tests indicate that deviations from this power law persist for changed initial longitudes and also for small but non-trivial changes in orbital spacing. We find that the stability time increases more rapidly at large initial orbital separations than the power-law dependence predicted from moderate initial orbital separations. Systems of five planets are less stable than systems of three planets for a specified semimajor axis spacing. Furthermore, systems of less massive planets can be packed more closely, being about as stable as 1 M_⊕ planets when the radial separation between planets is scaled using the mutual Hill radius. Finally, systems with retrograde planets can be packed substantially more closely than prograde systems with equal numbers of planets.     

Close encounters of nearly parabolic comets and planets

An overview is given of close encounters of nearly parabolic comets (NPCs; with periods of P > 200 years and perihelion distances of q > 0.1 AU; the number of the comets is N = 1041) with planets. The minimum distances Δ_min between the cometary and planetary orbits are calculated to select comets whose Δ_min are less than the radius of the planet’s sphere of influence. Close encounters of these comets with planets are identified by numerical integration of the comets’ equations of motion over an interval of ±50 years from the time of passing the perihelion. Close encounters of NPCs with Jupiter in 1663–2011 are reported for seven comets. An encounter with Saturn is reported for comet 2004 F2 (in 2001).     

Relationship of comets and planets

We analyze our earlier data on the numerical integration of the equations of motion for 274 short-period comets (with the period P<200 yr) on a time interval of 6000 yr. As many as 54 comets had no close approaches to planets, 13 comets passed through the Saturnian sphere of action, and one comet passed through the Uranian sphere of action. The orbital elements of these 68 comets changed by no more than ±3 percent in a space of 6000 yr. As many as 206 comets passed close to Jupiter. We confirm Everhart’s conclusion that Jupiter can capture long-period comets with q = 4–6 AU and i < 9° into short-period orbits. We show that nearly parabolic comets cross the solar system mainly in the zone of terrestrial planets. No relationship of nearly parabolic comets and terrestrial planets was found for the epoch of the latest apparition of comets. Guliev’s conjecture about two trans-Plutonian planets is based on the illusory excess of cometary nodes at large heliocentric distances. The existence of cometary nodes at the solar system periphery turns out to be a solely geometrical effect.     

Relativistic effects and solar oblateness from radar observations of planets and spacecraft

We used more than 250 000 high-precision American and Russian radar observations of the inner planets and spacecraft obtained in the period 1961–2003 to test the relativistic parameters and to estimate the solar oblateness. Our analysis of the observations was based on the EPM ephemerides of the Institute of Applied Astronomy, Russian Academy of Sciences, constructed by the simultaneous numerical integration of the equations of motion for the nine major planets, the Sun, and the Moon in the post-Newtonian approximation. The gravitational noise introduced by asteroids into the orbits of the inner planets was reduced significantly by including 301 large asteroids and the perturbations from the massive ring of small asteroids in the simultaneous integration of the equations of motion. Since the post-Newtonian parameters and the solar oblateness produce various secular and periodic effects in the orbital elements of all planets, these were estimated from the simultaneous solution: the post-Newtonian parameters are β = 1.0000 ± 0.0001 and γ = 0.9999 ± 0.0002, the gravitational quadrupole moment of the Sun is J₂ = (1.9 ± 0.3) × 10^?7, and the variation of the gravitational constant is ?/G = (?2 ± 5) × 10^?14 yr^?1. The results obtained show a remarkable correspondence of the planetary motions and the propagation of light to General Relativity and narrow significantly the range of possible values for alternative theories of gravitation.     

Predictions on the core mass of Jupiter and of giant planets in general

More than 80 giant planets are known by mass and radius. Their interior structure in terms of core mass, number of layers, and composition however is still poorly known. An overview is presented about the core mass M _core and envelope mass of metals M _Z in Jupiter as predicted by various equations of state. It is argued that the uncertainty about the true H/He EOS in a pressure regime where the gravitational moments J ₂ and J ₄ are most sensitive, i.e. between 0.5 and 4 Mbar, is in part responsible for the broad range \(M_{\mathit{core}}=0{-}18\:M_{\oplus }\), \(M_{Z}=0{-}38\:M_{\oplus }\), and \(M_{\mathit{core}}+M_{Z}=14{-}38\:M_{\oplus }\) currently offered for Jupiter. We then compare the Jupiter models obtained when we only match J ₂ with the range of solutions for the exoplanet \(\mathrm{GJ}\:436\mathrm{b}\), when we match an assumed tidal Love number k ₂ value.     

Stability of Trojan planets in multi-planetary systems

Today there are more than 340 extra-solar planets in about 270 extra-solar systems confirmed. Besides the observed planets there exists also the possibility of a Trojan planet moving in the same orbit as the Jupiter-like planet. In our investigation we take also into account the habitability of a Trojan planet and whether such a terrestrial planet stays in the habitable zone. Its stability was investigated for multi-planetary systems, where one of the detected giant planets moves partly or completely in the habitable zone. By using numerical computations, we studied the orbital behaviour up to 10⁷ years and determined the size of the stable regions around the Lagrangian equilibrium points for different dynamical models for fictitious Trojans. We also examined the interaction of the Trojan planets with a second or third giant planet, by varying its semimajor axis and eccentricity. We have found two systems (HD 155358 and HD 69830) that can host habitable Trojan planets. Another aim of this work was to determine the size of the stable region around the Lagrangian equilibrium points in the restricted three body problem for small mass ratios μ of the primaries μ ≤ 0.001 (e.g. Neptune mass of the secondary and smaller masses). We established a simple relation for the size depending on μ and the eccentricity.     

Effect of the darkforce on the extra-anomalous apsidal precession of solar planets

We demonstrate that, while the proposed Gravitational Dark-force Theory (of Nyambuya (New Astron. 67:1, 2019b) here-in Paper II) predicts an extra-anomalous apsidal precession for Solar planets due to the gravitational dark-force on the orbits of these planets, the predicted extra-anomalous apsidal precession is so small—so much that—it can not account for the observed extra-anomalous apsidal precession of Solar planets. This null result is important in that it informs us that whatever may be the cause of the extra-anomalous apsidal precession, it is not the proposed gravitational dark-force.     

Secular tidal changes in lunar orbit and Earth rotation

Small tidal forces in the Earth–Moon system cause detectable changes in the orbit. Tidal energy dissipation causes secular rates in the lunar mean motion n, semimajor axis a, and eccentricity e. Terrestrial dissipation causes most of the tidal change in n and a, but lunar dissipation decreases eccentricity rate. Terrestrial tidal dissipation also slows the rotation of the Earth and increases obliquity. A tidal acceleration model is used for integration of the lunar orbit. Analysis of lunar laser ranging (LLR) data provides two or three terrestrial and two lunar dissipation parameters. Additional parameters come from geophysical knowledge of terrestrial tides. When those parameters are converted to secular rates for orbit elements, one obtains dn/dt = \(-25.97\pm 0.05 ''/\)cent\(^{2}\), da/dt = 38.30 ± 0.08 mm/year, and di/dt = ?0.5 ± 0.1 \(\upmu \)as/year. Solving for two terrestrial time delays and an extra de/dt from unspecified causes gives \(\sim \) \(3\times 10^{-12}\)/year for the latter; solving for three LLR tidal time delays without the extra de/dt gives a larger phase lag of the N2 tide so that total de/dt = \((1.50 \pm 0.10)\times 10^{-11}\)/year. For total dn/dt, there is \(\le \)1 % difference between geophysical models of average tidal dissipation in oceans and solid Earth and LLR results, and most of that difference comes from diurnal tides. The geophysical model predicts that tidal deceleration of Earth rotation is \(-1316 ''\)/cent\(^{2}\) or 87.5 s/cent\(^{2}\) for UT1-AT, a 2.395 ms/cent increase in the length of day, and an obliquity rate of 9 \(\upmu \)as/year. For evolution during past times of slow recession, the eccentricity rate can be negative.     

Long-period comet flux in the planetary region: Dynamical evolution from the Oort cloud

This study analyzes the evolution of 2 × 10⁵ orbits with initial parameters corresponding to the orbits of comets of the Oort cloud under the action of planetary, galactic, and stellar perturbations over 2 × 10⁹ years. The dynamical evolution of comets of the outer (orbital semimajor axes a > 10⁴ AU) and inner (5 × 10³ < a (AU) < 10⁴) parts of the comet cloud is analyzed separately. The estimates of the flux of “new” and long-period comets for all perihelion distances q in the planetary region are reported. The flux of comets with a > 10⁴ AU in the interval 15 AU < q < 31 AU is several times higher than the flux of comets in the region q < 15 AU. We point out the increased concentration of the perihelia of orbits of comets from the outer cloud, which have passed several times through the planetary system, in the Saturn-Uranus region. The maxima in the distribution of the perihelia of the orbits of comets of the inner Oort cloud are located in the Uranus-Neptune region. “New” comets moving in orbits with a < 2 × 10⁴ AU and arriving at the outside of the planetary system (q > 25 AU) subsequently have a greater number of returns to the region q < 35 AU. The perihelia of the orbits of these comets gradually drift toward the interior of the Solar System and accumulate beyond the orbit of Saturn. The distribution of the perihelia of long-period comets beyond the orbit of Saturn exhibits a peak. We discuss the problem of replenishing the outer Oort cloud by comets from the inner part and their subsequent dynamical evolution. The annual rate of passages of comets of the inner cloud, which replenish the outer cloud, in the region q < 1 AU in orbits with a > 10⁴ AU (～ 5.0 × 10^?14 yr^?1) is one order of magnitude lower than the rate of passage of comets from the outer Oort cloud (～ 9.1 × 10^?13 yr^?1).     

Evolutionary studies of the young star clusters: NGC 1960, NGC 2453 and NGC 2384

We report the analysis of the young star clusters NGC 1960, NGC 2453 and NGC 2384 observed in the J (1.12 μm), H (1.65 μm) and K′ (2.2 μm) bands. Estimates of reddening, distance and age as E(B?V)=0.25, d=1380 pc and t=31.6 to 125 Myr for NGC 1960, E(B?V)=0.47, d=3311 pc and t=40 to 200 Myr for NGC 2453 and E(B?V)=0.25, d=3162 pc and t=55 to 125 Myr for NGC 2384 have been obtained. Also, we have extended the color–magnitude diagrams of these clusters to the fainter end and thus extended the luminosity functions to fainter magnitudes. The evolution of the main sequence and luminosity functions of these clusters have been compared with themselves as well as Lyngå 2 and NGC 1582.     

Variability of the spectrum and brightness of RY Tau

The results of spectroscopic observations of the star RY Tau in the ultraviolet based on IUE data and in the visual spectral range obtained at ShAO are presented. Despite significant brightness variability in 1983–1984, the Mg II λ2800 Å emission doublet showed no synchronous variation with the UBV photometric data. Periodic variability of the Mg II λ2800 Å emission intensity with a period of 23 days has been detected for the first time. The periodicity is also observed for a group of such lines as CIV λ1450 Å, He II λ1640 Å, and S II λ1756 Å. The equivalent widths and shifts of the individual components of the Hα, H + H _ε , and CaII K lines also vary with the period found. The observed variability of the emission spectrum can be explained by the existence of a companion in the system in an orbit with a semimajor axis of about 0.13 AU.     

Astrometric Study of Two Multiple Stars: ADS 2668 and ADS 8236

Based on CCD observations with the Pulkovo 26-inch refractor in 2003–2018, we have obtained the orbit of the visual double star ADS 2668 AB (P = 947 yr, a = 2.9″, e = 0.41, ω = 246°, Ω = 131°, i = 114°, T = 1456 yr) for the first time by the apparent motion parameter (AMP) method, which is consistent with the inner orbit of ADS 2668 Aa-Ab, and improved the orbit of ADS 8236 AB (P = 1996 yr, a = 4.69″, e = 0.39, ω = 201°, Ω = 166°, i = 110°, T = 1246 yr). The inner orbit of the photocenter of ADS 8236 with a period of 4.627 yr has been calculated from the residuals. This orbit of ADS 8236 Ba–Bb supplements the spectroscopic orbit by the elements specifying the orbital plane (i and Ω). In both cases, the planes of the inner and outer orbits are noncoplanar. The presence of an additional companion in the system ADS 2668 is discussed.     

Capture of comets from the Oort cloud into Halley-type and Jupiter-family orbits

This paper analyzes the capture of comets into Halley-type and Jupiter-family orbits from the nearparabolic flux of the Oort cloud. Two types of capture into Halley-type orbits are found. The first type is the evolution of near-parabolic orbits into short-period orbits (with heliocentric orbital periods P < 200 years) as a result of close encounters with giant planets. This process is followed by a very slow drift of cometary orbits into the inner part of the Solar System. Only those comets may pass from short-period orbits into Halley-type and Jupiter-family orbits, which move in orbits with perihelion distances q < 13 au. In the second type of capture, the perihelion distances of cometary orbits become rather small (< 1.5 au) during the first stage of dynamic evolution under the action of perturbations from the Galaxy, and then their semimajor axes decrease as a result of diffusion. The capture takes place, on average, in 500 revolutions of the comet about the Sun, whereas in the first case, the comet is captured, on average, after 12500 revolutions. The region of initial orbital perihelion distances q > 4 au is found to be at least as important a source of Halley-type comets as the region of perihelion distances q < 4 au. More than half of the Halley-type comets are captured from the nearly parabolic flux with q > 4 au. The analysis of the dynamic evolution of objects moving in short-period orbits shows that the distribution of Centaurs orbits agrees well with the observed distribution corrected for observational selection effects. Hence, the hypothesis associating the origin of Centaurs with the Edgeworth-Kuiper belt and the trans-Neptunian region exclusively should be rejected.     

On the MOND external field effect in the solar system

In the framework of the MOdified Newtonian Dynamics (MOND), the internal dynamics of a gravitating system s embedded in a larger one S is affected by the external background field E of S even if it is constant and uniform, thus implying a violation of the Strong Equivalence Principle: it is the so-called External Field Effect (EFE). In the case of the solar system, E would be A _cen≈10^?10 m?s^?2 because of its motion through the Milky Way: it is orders of magnitude smaller than the main Newtonian monopole terms for the planets. We address here the following questions in a purely phenomenological manner: are the Sun’s planets affected by an EFE as large as 10^?10 m?s^?2? Can it be assumed that its effect is negligible for them because of its relatively small size? Does E induce vanishing net orbital effects because of its constancy over typical solar system’s planetary orbital periods? It turns out that a constant and uniform acceleration, treated perturbatively, does induce non-vanishing long-period orbital effects on the longitude of the pericenter ? of a test particle. In the case of the inner planets of the solar system and with E≈10^?10 m?s^?2, they are 4–6 orders of magnitude larger than the present-day upper bounds on the non-standard perihelion precessions \(\Delta\dot{\varpi}\) recently obtained with by E.V. Pitjeva with the EPM ephemerides in the Solar System Barycentric frame. The upper limits on the components of E are E _x ≤1×10^?15 m?s^?2, E _y ≤2×10^?16 m?s^?2, E _z ≤3×10^?14 m?s^?2. This result is in agreement with the violation of the Strong Equivalence Principle by MOND. Our analysis also holds for any other exotic modification of the current laws of gravity yielding a constant and uniform extra-acceleration. If and when other corrections \(\Delta\dot{\varpi}\) to the usual perihelion precessions will be independently estimated with different ephemerides it will be possible to repeat such a test.     

Dynamical lifetime survey of geostationary transfer orbits

In this paper, we study the long-term dynamical evolution of highly elliptical orbits in the medium-Earth orbit region around the Earth. The real population consists primarily of Geosynchronous Transfer Orbits (GTOs), launched at specific inclinations, Molniya-type satellites and related debris. We performed a suite of long-term numerical integrations (up to 200 years) within a realistic dynamical model, aimed primarily at recording the dynamical lifetime of such orbits (defined as the time needed for atmospheric reentry) and understanding its dependence on initial conditions and other parameters, such as the area-to-mass ratio (A / m). Our results are presented in the form of 2-D lifetime maps, for different values of inclination, A / m, and drag coefficient. We find that the majority of small debris (\(>70\%\), depending on the inclination) can naturally reenter within 25–90 years, but these numbers are significantly less optimistic for large debris (e.g., upper stages), with the notable exception of those launched from high latitude (Baikonur). We estimate the reentry probability and mean dynamical lifetime for different classes of GTOs and we find that both quantities depend primarily and strongly on initial perigee altitude. Atmospheric drag and higher A / m values extend the reentry zones, especially at low inclinations. For high inclinations, this dependence is weakened, as the primary mechanisms leading to reentry are overlapping lunisolar resonances. This study forms part of the EC-funded (H2020) “ReDSHIFT” project.     

Numerical periodic orbits of charged grains around magnetic planets

We apply a numerical searching method to investigate three-dimensional periodic orbits of charged dust particles in planetary magnetospheres. A classic generalized Stormer model of magnetic planets along with the parameters of Saturn is employed. More periodic orbits are found, besides the already known circular periodic orbits in or parallel to the equatorial plane. We divide all these orbits into six categories based on their appearances. By calculating the characteristic multipliers of the orbits, we investigate the stabilities of these periodic orbits.     

Inverse problem of central configurations and singular curve in the collinear 4-body problem

In this paper, we consider the inverse problem of central configurations of n-body problem. For a given \({q=(q_1, q_2, \ldots, q_n)\in ({\bf R}^d)^n}\), let S(q) be the admissible set of masses denoted \({ S(q)=\{ m=(m_1,m_2, \ldots, m_n)| m_i \in {\bf R}^+, q}\) is a central configuration for m}. For a given \({m\in S(q)}\), let S _m (q) be the permutational admissible set about m = (m ₁, m ₂, . . . , m _n ) denoted

$S_m(q)=\{m^\prime | m^\prime\in S(q),m^\prime \not=m \, {\rm and} \, m^\prime\,{\rm is\, a\, permutation\, of }\, m \}.$

The main discovery in this paper is the existence of a singular curve \({\bar{\Gamma}_{31}}\) on which S _m (q) is a nonempty set for some m in the collinear four-body problem. \({\bar{\Gamma}_{31}}\) is explicitly constructed by a polynomial in two variables. We proved:

1. (1)
   If \({m\in S(q)}\), then either ^# S _m (q) = 0 or ^# S _m (q) = 1.
    
2. (2)
   ^#S _m (q) = 1 only in the following cases:

   1. (i)
      If s = t, then S _m (q) = {(m ₄, m ₃, m ₂, m ₁)}.
       
   2. (ii)
      If \({(s,t)\in \bar{\Gamma}_{31}\setminus \{(\bar{s},\bar{s})\}}\), then either S _m (q) = {(m ₂, m ₄, m ₁, m ₃)} or S _m (q) = {(m ₃, m ₁, m ₄, m ₂)}.
       
    

     

The matter in the Big-Bang model is dust and not any arbitrary perfect fluid!

We consider a spherically symmetric general relativistic perfect fluid in its comoving frame. It is found that, by integrating the local energy momentum conservation equation, a general form of g ₀₀ can be obtained. During this study, we get a cue that an adiabatically evolving uniform density isolated sphere having ρ(r,t)=ρ ₀(t), should comprise “dust” having p ₀(t)=0; as recently suggested by Durgapal and Fuloria (J. Mod. Phys. 1:143, 2010) In fact, we offer here an independent proof to this effect. But much more importantly, we find that for the homogeneous and isotropic Friedmann-Robertson-Walker (FRW) metric having p(r,t)=p ₀(t) and ρ(r,t)=ρ ₀(t), \(g_{00} = e^{-2p_{0}/(p_{0} +\rho_{0})}\). But in general relativity (GR), one can choose an arbitrary t→t _?=f(t) without any loss of generality, and thus set g ₀₀(t _?)=1. And since pressure is a scalar, this implies that p ₀(t _?)=p ₀(t)=0 in the Big-Bang model based on the FRW metric. This result gets confirmed by the fact the homogeneous dust metric having p(r,t)=p ₀(t)=0 and ρ(r,t)=ρ ₀(t) and the FRW metric are exactly identical. In other words, both the cases correspond to the same Einstein tensor \(G^{a}_{b}\) because they intrinsically have the same energy momentum tensor \(T^{a}_{b}=\operatorname {diag}[\rho_{0}(t), 0,0, 0]\).