Expansion of the perturbing function of a satellite restricted four-body problem

A satellite four-body problem is the problem of motion of an artificial satellite of a planet in a region of the space where perturbations due to the gravitational field of the planet are of the same order as perturbations due to influences of two perturbing bodies. In this paper an expansion of the perturbing function into a Fourier series in terms of angular Keplerian elements ( _j , _j ,M _j :j=0,1,2) (designations are standard) is obtained taking into account a sharp commensurability of the typen/ ₀=(p+q)/p (n is the mean motion of the artificial satellite and ₀ is the angular velocity of rotation of the planet,p andq are integers).The coefficients of the Fourier series are the functions of the positional Keplerian elements (a _j ,e _j ,i _j ;j=0, 1, 2) (designations are standard) and, in particular, are series in terms ofe _j that, generally speaking, can be written out to an accuracy ofe _j ¹⁹ .The expansion obtained can be used for the construction of a semianalytical theory of motion of resonant satellites on the basis of conditionally periodic solutions of the restricted four-body problem.     

The lunar orbit revisited,III

In this paper we present an investigation on the tidal evolution of a system of three bodies: the Earth, the Moon and the Sun. Equations are derived including dissipation in the planet caused by the tidal interaction between the planet and the satellite and between the planet and the sun. Dissipation within the Moon is included as well. The set of differential equations obtained is valid as long as the solar disturbances dominate the perturbations on the satellite's motion due to the oblateness of the planet, namelya/R _e greater than 15, and closer than that point equations derived in a preceding paper are used.The result shows the Moon was closer to the Earth in the past than now with an inclination to the ecliptic greater than today, whereas the obliquity was smaller. Toward the past, the inclination to the Earth's equator begins decreasing to 12° fora/R _e=12 and suddenly grows. During the first stage the results are weakly dependant on the magnitude of the dissipation within the satellite, whereas the distance of the closest approach and the prior history are strongly dependent on that dissipation. In particular, the crossing of the Roche limit can be avoided.     

Constructive analytical solution of the evolution hill problem

A new analytical solution of the system of differential equations describing secular perturbations and long-period solar perturbations of mean orbits of outer satellites of giant planets was obtained. As distinct from other solutions, the solution constructed using von Zeipel’s method approximately takes into account, in the secular part of the perturbing function, the totality of fourth order with respect to the small parameter m of the ratio of the mean motions of the primary planet and the satellite. This enables us to describe more accurately the evolution of satellite orbits with large apocentric distances, which in the course of evolution may exceed the halved radius of the Hill sphere of the planet with respect to the Sun. Among these are the orbits of the two outermost Neptunian satellites N10 (Psamathe) and N13 (Neso). For these satellites, the parameter m amounts to 0.152 and 0.165, respectively. Different from a purely analytical solution, the proposed solution requires preliminary calculations for each satellite. More precisely, in doing so, we need to construct some simple functions to approximate more complex ones. This is why we use the phrase “constructive analytical.” To illustrate the solution, we compare it with the results of the numerical integration of the strict motion equations of the satellites N10 and N13 over time intervals 5–15 thousand years.     

Eccentricities of Planets in Binary Systems

The most puzzling property of the extrasolar planets discovered by recent radial velocity surveys is their high orbital eccentricities, which are very difficult to explain within our current theoretical paradigm for planet formation. Current data reveal that at least 25% of these planets, including some with particularly high eccentricities, are orbiting a component of a binary star system. The presence of a distant companion can cause significant secular perturbations in the orbit of a planet. At high relative inclinations, large-amplitude, periodic eccentricity perturbations can occur. These are known as “Kozai cycles” and their amplitude is purely dependent on the relative orbital inclination. Assuming that every planet host star also has a (possibly unseen, e.g., substellar) distant companion, with reasonable distributions of orbital parameters and masses, we determine the resulting eccentricity distribution of planets and compare it to observations? We find that perturbations from a binary companion always appear to produce an excess of planets with both very high (?0.6) and very low (e ? 0.1) eccentricities. The paucity of near-circular orbits in the observed sample implies that at least one additional mechanism must be increasing eccentricities. On the other hand, the overproduction of very high eccentricities observed in our models could be combined with plausible circularization mechanisms (e.g., friction from residual gas) to create more planets with intermediate eccentricities (e? 0.1–0.6).     

Love numbers of the moon and of the terrestrial planets

In the IERS Standards (1989), for the Moon the adopted value of the tide Love number, k ₂, is equal to 0.0222. In this paper using the latest geodetic parameters of the Moon a group of internal structure models are constructed for this celestial body (see Table V), then the dependence of the Moon's core size on calculated value of k ₂ is explored. The obtained results indicate that the second degree Love number, k ₂ = 0.02664, of the lunar model 91–04 is near its observed value (0.027 ± 0.006). This implies that the Moon may possess an outer core of 660 km radius and of 300 kbar mean rigidity. With the same method the static Love numbers from degree 2 to 30 are computed for the terrestrial planets — Mercury, Venus, and Mars (see Table VII), and the influence of some parameters (such as the rigidity) of the outer core on low degree Love numbers is discussed. Finally, the likely range of the second degree Love numbers is determined for the terrestrial planets (see Table XI). It seems that if low degree Love numbers of a terrestrial planet can be detected in the future space explorations, there is some possibility to improve the planetary internal structure model. For example, as soon as space techniques yield an observed value of k ₂ > 0.10 for Mercury, there will be reason to anticipate that a partly melted iron core exists in this planet.     

Velocity distribution of meteoroids in the vicinity of planets and satellites

Collisions in the Solar System play an important role in its history. Impact processes depend essentially on the velocity distribution of meteoroids colliding with a chosen planet. According to Carleman's theorem it is sufficient to find the set of M _k = mathematical expectation of v ^k , v being the collisional velocity. We suppose that M _k for meteoroids of asteroidal nature differs slightly from that for asteroids themselves. So among all numbered minor planets we select those which may potentially collide with the chosen major planet. Then we calculate v at intersection points and count the average over all such points and all selected asteroids. The gravitation of a body-target may be taken into account or not. Numerical results are collected in four Tables.St.Petersburg University     

Numerical integration of dynamical systems with Lie series

The integration of the equations of motion in gravitational dynamical systems—either in our Solar System or for extra-solar planetary systems—being non integrable in the global case, is usually performed by means of numerical integration. Among the different numerical techniques available for solving ordinary differential equations, the numerical integration using Lie series has shown some advantages. In its original form (Hanslmeier and Dvorak, Astron Astrophys 132, 203 1984), it was limited to the N-body problem where only gravitational interactions are taken into account. We present in this paper a generalisation of the method by deriving an expression of the Lie terms when other major forces are considered. As a matter of fact, previous studies have been done but only for objects moving under gravitational attraction. If other perturbations are added, the Lie integrator has to be re-built. In the present work we consider two cases involving position and position-velocity dependent perturbations: relativistic acceleration in the framework of General Relativity and a simplified force for the Yarkovsky effect. A general iteration procedure is applied to derive the Lie series to any order and precision. We then give an application to the integration of the equation of motions for typical Near-Earth objects and planet Mercury.     

The effects of nebula surface density profile and giant‐planet eccentricities on planetary accretion in the inner solar system

Abstract— We describe results of 32 N‐body planetary accretion simulations that investigate the dependence of terrestrial‐planet formation on nebula surface density profile σ and evolution of the eccentricities of Jupiter and Saturn e_j,s. Two surface density profiles are examined: a decaying profile with σ ∝ 1/a, where a is orbital semi‐major axis, and a peaked profile in which σ increases for a < 2 AU and decreases for a > 2 AU. The peaked profiles are generated by models of coagulation in an initially hot nebula. Models with initial e_j,s = 0.05 (the current value) and 0.1 are considered. Simulations using the decaying profile with e_j,s = 0.1 produce systems most like the observed planets in terms of mass‐weighted mean a and the absence of a planet in the asteroid belt. Simulations with doubled σ produce planets roughly twice as massive as the nominal case. Most initial embryos are removed in each simulation via ejection from the solar system or collision with the Sun. The asteroid belt is almost entirely cleared on a timescale of 10–100 Ma that depends sensitively on e_j,s. Most initial mass with a < 2 AU survives, with the degree of mass loss increasing with a. Mass loss from the terrestrial region occurs on a timescale that is long compared to the mass loss time for the asteroid belt. Substantial radial mixing of material occurs in all simulations, but is greater in simulations with initital e_j,s = 0.05. The degree of mixing is equivalent to a feeding zone of half width 1.5 and 0.9 AU for an Earth mass planet at 1 AU for the cases e_j,s = 0.05 and 0.1, respectively. In simulations with e_j,s = 0.05, roughly one‐third and 5–10% of the mass contained in final terrestrial planets originated in the region a > 2.5 AU for the decaying and peaked profiles, respectively. In the case e_j,s = 0.1, the median mass accreted from a > 2.5 AU is zero for both profiles.     

Perturbations of a close-earth satellite due to sunlight diffusely reflected from the Earth

A semianalytical method has been developed to calculate the radiation-pressure perturbations of a close-Earth satellite due to sunlight reflected from the Earth. It is assumed that the satellite is spherically symmetric and that the solar radiation is reflected from the Earth according to Lambert's Law. To account for the increasing reflectivity of the Earth toward the poles, its albedo is assumed to have a latitudinal dependence given bya=a ₀ +a ₂ sin². The effect of the terminator on the perturbations has been neglected. The perturbations within a particular revolution are given analytically, while the long-range perturbations are obtained by accumulation.     

The 3:2 spin-orbit resonant motion of Mercury

Our purpose is to build a model of rotation for a rigid Mercury, involving the planetary perturbations and the non-spherical shape of the planet. The approach is purely analytical, based on Hamiltonian formalism; we start with a first-order basic averaged resonant potential (including J ₂ and C ₂₂, and the first powers of the eccentricity and the inclination of Mercury). With this kernel model, we identify the present equilibrium of Mercury; we introduce local canonical variables, describing the motion around this 3:2 resonance. We perform a canonical untangling transformation, to generate three sets of action-angle variables, and identify the three basic frequencies associated to this motion. We show how to reintroduce the short-periodic terms, lost in the averaging process, thanks to the Lie generator; we also comment about the damping effects and the planetary perturbations. At any point of the development, we use the model SONYR to compare and check our calculations.     

Frozen orbits for satellites close to an Earth-like planet

We say that a planet is Earth-like if the coefficient of the second order zonal harmonic dominates all other coefficients in the gravity field. This paper concerns the zonal problem for satellites around an Earth-like planet, all other perturbations excluded. The potential contains all zonal coefficientsJ ₂ throughJ ₉. The model problem is averaged over the mean anomaly by a Lie transformation to the second order; we produce the resulting Hamiltonian as a Fourier series in the argument of perigee whose coefficients are algebraic functions of the eccentricity — not truncated power series. We then proceed to a global exploration of the equilibria in the averaged problem. These singularities which aerospace engineers know by the name of frozen orbits are located by solving the equilibria equations in two ways, (1) analytically in the neighborhood of either the zero eccentricity or the critical inclination, and (2) numerically by a Newton-Raphson iteration applied to an approximate position read from the color map of the phase flow. The analytical solutions we supply in full to assist space engineers in designing survey missions. We pay special attention to the manner in which additional zonal coefficients affect the evolution of bifurcations we had traced earlier in the main problem (J ₂ only). In particular, we examine the manner in which the odd zonalJ ₃ breaks the discrete symmetry inherent to the even zonal problem. In the even case, we find that Vinti's problem (J ₄+J ₂ ² =0) presents a degeneracy in the form of non-isolated equilibria; we surmise that the degeneracy is a reflection of the fact that Vinti's problem is separable. By numerical continuation we have discovered three families of frozen orbits in the full zonal problem under consideration; (1) a family of stable equilibria starting from the equatorial plane and tending to the critical inclination; (2) an unstable family arising from the bifurcation at the critical inclination; (3) a stable family also arising from that bifurcation and terminating with a polar orbit. Except in the neighborhood of the critical inclination, orbits in the stable families have very small eccentricities, and are thus well suited for survey missions.     

An analytical method to compute comet cloud formation efficiency and its application

A quick analytical method is presented for calculating comet cloud formation efficiency in the case of a single planet or multiple-planet system for planets that are not too eccentric (e _p ≲ 0.3). A method to calculate the fraction of comets that stay under the control of each planet is also presented, as well as a way to determine the efficiency in different star cluster environments. The location of the planet(s) in mass-semi-major axis space to form a comet cloud is constrained based on the conditions developed by Tremaine (1993) together with estimates of the likelyhood of passing comets between planets; and, in the case of a single, eccentric planet, the additional constraint that it is, by itself, able to accelerate material to relative encounter velocity U ~ 0.4 within the age of the stellar system without sweeping up the majority of the material beforehand. For a single planet, it turns out the efficiency is mainly a function of planetary mass and semi-major axis of the planet and density of the stellar environment. The theory has been applied to some extrasolar systems and compared to numerical simulations for both these systems and the Solar System, as well as a diffusion scheme based on the energy kick distribution of Everhart (Astron J 73:1039–1052, 1968). The analytic results are in good agreement with the simulations.     

Delivery of icy planetesimals to inner planets in the Proxima Centauri planetary system

The estimates of the delivery of icy planetesimals from the feeding zone of Proxima Centauri c (with mass equal to 7m_E, m_E is the mass of the Earth) to inner planets b and d were made. They included the studies of the total mass of planetesimals in the feeding zone of planet c and the probabilities of collisions of such planetesimals with inner planets. This total mass could be about 10–15m_E. It was estimated based on studies of the ratio of the mass of planetesimals ejected into hyperbolic orbits to the mass of planetesimals collided with forming planet c. At integration of the motion of planetesimals, the gravitational influence of planets c and b and the star was taken into account. In most series of calculations, planetesimals collided with planets were excluded from integrations. Based on estimates of the mass of planetesimals ejected into hyperbolic orbits, it was concluded that during the growth of the mass of planet c the semi-major axis of its orbit could decrease by at least a factor of 1.5. Depending on possible gravitational scattering due to mutual encounters of planetesimals, the total mass of material delivered by planetesimals from the feeding zone of planet c to planet b was estimated to be between 0.002m_E and 0.015m_E. Probably, the amount of water delivered to Proxima Centauri b exceeded the mass of water in Earth's oceans. The amount of material delivered to planet d could be a little less than that delivered to planet b.     

Perturbation of a satellite orbiter by the oblateness of the primary planet

The perturbation of an orbiter around a large satellite of a giant planet (Jupiter, Saturn, Uranus or Neptune) produced by the oblateness of the planet is investigated. The perturbing force of theJ ₂-term (general case) and theJ ₄-term (special case of small eccentricity and inclination) is expanded in an appropriate form and the main term and the parallactic term are given explicitly. The variations of the orbital elements are derived using the stroboscopic method. An example shows that the perturbation of the orbit cannot be neglected.     

Search for circum‐planetary material and orbital period variations of short‐period Kepler exoplanet candidates

A unique short‐period (P = 0.65356(1) d) Mercury‐size Kepler exoplanet candidate KIC012557548b has been discovered recently by Rappaport et al. (2012). This object is a transiting disintegrating exoplanet with a circum‐planetary material–comet‐like tail. Close‐in exoplanets, like KIC012557548b, are subjected to the greatest planet‐star interactions. This interaction may have various forms. In certain cases it may cause formation of the comet‐like tail. Strong interaction with the host star, and/or presence of an additional planet may lead to variations in the orbital period of the planet. Our main aim is to search for comet‐like tails similar to KIC012557548b and for long‐term orbital period variations. We are curious about frequency of comet‐like tail formation among short‐period Kepler exoplanet candidates. We concentrate on a sample of 20 close‐in candidates with a period similar to KIC012557548b from the Kepler mission. We first improved the preliminary orbital periods and obtained the transit light curves. Subsequently we searched for the signatures of a circum‐planetary material in these light curves. For this purpose the final transit light curve of each planet was fitted with a theoretical light curve, and the residuals were examined for abnormalities. We then searched for possible long‐term changes of the orbital periods using the method of phase dispersion minimization. In 8 cases out of 20 we found some interesting peculiarities, but none of the exoplanet candidates showed signs of a comet‐like tail. It seems that the frequency of comet‐like tail formation among short‐period Kepler exoplanet candidates is very low. We searched for comet‐like tails based on the period criterion. Based on our results we can conclude that the short‐period criterion is not enough to cause comet‐like tail formation. This result is in agreement with the theory of the thermal wind and planet evaporation (Perez‐Becker & Chiang 2013). We also found 3 cases of candidates which showed some changes of the orbital period. Based on our results we can see that orbital period changes are not caused by comet‐like tail disintegration processes, but rather by possible massive outer companions. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Averaging on the motion of a fast revolving body. Application to the stability study of a planetary system

Exploring the global dynamics of a planetary system involves computing integrations for an entire subset of its parameter space. This becomes time-consuming in presence of a planet close to the central star, and in practice this planet will be very often omitted. We derive for this problem an averaged Hamiltonian and the associated equations of motion that allow us to include the average interaction of the fast planet. We demonstrate the application of these equations in the case of the μ Arae system where the ratio of the two fastest periods exceeds 30. In this case, the effect of the inner planet is limited because the planet’s mass is one order of magnitude below the other planetary masses. When the inner planet is massive, considering its averaged interaction with the rest of the system becomes even more crucial.     

The structure of Phaethon and detonation of its icy envelope

It is suggested that Phaethon — a hypothetical planet whose breakup gave rise to the asteroid belt — has a structure similar to that of Callisto, and thus consisted of a rocky core (40%; in mass) and an ice envelope (60%). Total breakup of the planet becomes possible in an explosition of the electrolysis products accumulated in the ice in the form of a solid solution if the planet mass M 0.5 M. Assuming M = 0.5 M we obtain 1750 km for the planet's radius with the envelope's thickness of 750km. Application of the hydrodynamic theory of detonation to the (2H₂ + O₂) solution in ice reveals that depending on the actual critical temperature which for conventional explosives lies in the range 700–900 K the minimum (2H₂ + O₂) concentration in ice required for its explosion is 13–18%;.     

Against all odds? Forming the planet of the HD 196885 binary

HD 196885 Ab is the most ??extreme?? planet-in-a-binary discovered to date, whose orbit places it at the limit for orbital stability. The presence of a planet in such a highly perturbed region poses a clear challenge to planet-formation scenarios. We investigate this issue by focusing on the planet-formation stage that is arguably the most sensitive to binary perturbations: the mutual accretion of kilometre-sized planetesimals. To this effect we numerically estimate the impact velocities dv amongst a population of circumprimary planetesimals. We find that most of the circumprimary disc is strongly hostile to planetesimal accretion, especially the region around 2.6 AU (the planet??s location) where binary perturbations induce planetesimal-shattering dv of more than 1 kms^?1. Possible solutions to the paradox of having a planet in such accretion-hostile regions are (1) that initial planetesimals were very big, at least 250 km (2) that the binary had an initial orbit at least twice the present one, and was later compacted due to early stellar encounters (3) that planetesimals did not grow by mutual impacts but by sweeping of dust (the ??snowball?? growth mode identified by Xie et al., in Astrophys J 724:1153, 2010b), or (4) that HD 196885 Ab was formed not by core-accretion but by the concurrent disc instability mechanism. All of these 4 scenarios remain however highly conjectural.     

An Upper Limit of Mass for a Stable Minor Planet in the Main Asteroid Belt

The stability of an imaginary planet located in the present main asteroid belt is studied with a 7-body model (Sun, Mars, Jupiter, Saturn, Uranus, Neptune and the imaginary planet). The fourth-order Hermite algorithm P(EC)³ is used, which has a very small secular energy error for the integration of periodic orbits with a constant time-step. The evolution of orbits is followed up to 10⁸ years. Our numerical results show that the low-order resonances with Jupiter can enhance the stability of the imaginary planet in some cases. The survival probability of the imaginary planet decreases with the planet mass. The upper limit of the imaginary planet's mass that can survive in the main belt is around 10²⁵ kg, i.e., about the Earth's mass.     

Mass and orbit estimation of Planet X via a family of comets

The crucial assumption of this paper is, that the observed clustering of aphelion distances of intermediate-period comets in the 70–90 AU range is due to the influence of a tenth planet, called Planet X. We contribute to the search for Planet X a new and extended evaluation of a family of comets assumed to be Planet X's family of comets.By averaging the aphelion distances of comets that belong to a transplutonic family of comets, we get Planet X's semi-major axis a _x = (83.0 ± 5.3) AU. The comets' orbits also yield the upper limit of the planet's orbital eccentricity e _x - 0.019. If this planet played an important part in sending quasi-periodic comet showers to the inner solar system, we can calculate its orbital inclination i _x = 46 .1 ± 3 .6. By distributing all planets' masses into the heliocentric, torus-like zones, in which they were formed, we get the density distribution of the primordial solar nebula. Extrapolating this distribution we find the mass of the planet M _x = (5.1_–2.4 ^+3.6 M _Earth. A few plausible assumptions (e.g. Uranus and Neptune perturbations being caused by Planet X) lead to Planet X's actual location with declination and eccliptic longitude being = 57 ± 17 and = 54 ± 34 , respectively (1989.5 position). In addition, we give Planet X's apparent brightness dependent on its unknown albedo. All those properties and predictions are more or less in agreement with earlier work on Planet X.