Relativistic Precession of Elliptical Orbits of a Circumbinary Planet can be Cancelled

In publications presenting analytical results on the non-coplanar motion of a circumbinary planet it was shown that the unperturbed elliptical orbit of the planet undergoes simultaneously two kinds of the precession: the precession of the orbital plane and the precession of the orbit in its own plane. It is also well-known that there is also the relativistic precession of the planetary orbit in its own plane. In the present paper we study a combined effect of the all of the above precessions. For the general case, where the planetary orbit is not coplanar with the stars orbits, we analyzed the dependence of the critical inclination angle i_c, at which the precession of the planetary orbit in its own plane vanishes, on the angular momentum L of the planet. We showed that the larger the angular momentum, the smaller the critical inclination angle becomes. We presented the analytical result for i_c(L) and calculated the value of L, for which the critical inclination value becomes zero. For the particular case, where the planetary orbit is not coplanar with the stars orbits, we demonstrated analytically that at a certain value of the angular momentum of the planet, the elliptical orbit of the planet would become stationary: no precession. In other words, at this value of the angular momentum, the relativistic precession of the planetary orbit and its precession, caused by the fact that the planet revolves around a binary (rather than single) star, cancel each other out. This is a counterintuitive result.     

About Hypothesis of the Superconducting Origin of the Saturn’s Rings

Hypothesis of possible superconductivity of the iced matter of the rings of Saturn (based on the data of Voyager and Pioneer space missions) allow us to explain many phenomena which have not been adequately understood earlier. Introducing into planetary physics force of magnetic levitation of the superconducting iced particle of the rings, which interact with magnetosphere of the planet, becomes to be possible to explain origin, evolution, and dynamics of the rings; to show how the consequent precipitation of the rings’ matter upon the planet was concluded; how the rings began their rotation; how they were compressed by the magnetic field into the thin disc, and how this disc was fractured into hundreds of thousands of separated rings; why in the ring B do exist “spokes”; why magnetic field lines have distortion near by ring F; why there is a variable azimuth brightness of the ring A; why the rings reflected radio waves so efficiently; why exists strong electromagnetic radiation of the rings in the 20.4 kHz–40.2 MHz range and Saturnian kilometric radiation; why there is anomalous reflection of circularly polarized microwaves; why there are spectral anomalies of the thermal radiation of the rings; why the matter of the various rings does not mix but preserves its small-scale color differences; why there is an atmosphere of unknown origin nearby the rings of Saturn; why there are waves of density and bending waves within Saturn’s rings; why planetary rings in the solar system appear only after the Belt of Asteroids (and may be the Belt of Asteroids itself is a ring for the Sun); why our planet Earth has no rings of its own.     

Resonant Periodic Motion and the Stability of Extrasolar Planetary Systems

Families of nearly circular periodic orbits of the planetary type are studied, close to the 3/1 mean motion resonance of the two planets, considered both with finite masses. Large regions of instability appear, depending on the total mass of the planets and on the ratio of their masses.Also, families of resonant periodic orbits at the 2/1 resonance have been studied, for a planetary system where the total mass of the planets is the 4% of the mass of the sun. In particular, the effect of the ratio of the masses on the stability is studied. It is found that a planetary system at this resonance is unstable if the mass of the outer planet is smaller than the mass of the inner planet.Finally, an application has been made for the stability of the observed extrasolar planetary systems HD82943 and Gliese 876, trapped at the 2/1 resonance.     

On the stability of dust orbits in mean-motion resonances perturbed by from an interstellar wind

Circumstellar dust particles can be captured in a mean-motion resonance (MMR) with a planet and simultaneously be affected by non-gravitational effects. It is possible to describe the secular variations of a particle orbit in the MMR analytically using averaged resonant equations. We derive the averaged resonant equations from the equations of motion in near-canonical form. The secular variations of the particle orbit depending on the orientation of the orbit in space are taken into account. The averaged resonant equations can be derived/confirmed also from Lagrange’s planetary equations. We apply the derived theory to the case when the non-gravitational effects are the Poynting–Robertson effect, the radial stellar wind, and an interstellar wind. The analytical and numerical results obtained are in excellent agreement. We found that the types of orbits correspond to libration centers of the conservative problem. The averaged resonant equations can lead to a system of equations which holds for stationary points in a subset of resonant variables. Using this system we show analytically that for the considered non-gravitational effects, all stationary points should correspond to orbits which are stationary in interplanetary space after an averaging over a synodic period. In an exact resonance, the stationary orbits are stable. The stability is achieved by a periodic repetition of the evolution during the synodic period. Numerical solutions of this system show that there are no stationary orbits for either the exact or non-exact resonances.     

Computer simulation of planet formation in a binary star system: Terrestrial planets

A model of planetary formation in a binary system with a small relative mass of primary is computed on the assumption of a mass transfer from the less massive component to the more massive one with no mass and angular momentum carried away from the system under consideration. At the last stage of mass transfer the condensed Moon-like objects (planetoids) are ejected through the inner Lagrange point of the primary Roche lobe with the outflow of gaseous matter.The whole system is considered in the plane of binary star rotation. Newtonian equations of motion are integrated with the initial conditions for the planetoids referred to as the coordinates and velocity of the inner Lagrangian point at the moments of planetoid ejections, all the pairwise gravitational interactions being included in computations but without a gas-drag. The mass transfer ceases at the primary relative mass 10^–3 which corresponds to the present Sun-Jupiter system. The total mass of planetoids approximates that of the terrestrial planets. Those are formed through coagulation of the planetoids with the effective radius of capture cross-section as an input parameter in the computer simulation. When the minimum separation between the pair of bodies becomes less than this radius they coalesce into a single body with their masses and momenta summed. If the effective radius value is under a certain limit the computer simulation yields the planetary system like that of terrestrial planets of the present Sun system.Numerical computations reveal the division of the planetoids into 4 groups along their distances from the Sun. Further, each group forms a single planet or a planet and a less massive body at the nearest orbits. The parameters of simulated planet orbits are close to the present ones and the interplanetary spacings are in accord with the Titius-Bode law.     

Refined model for the evolution of distant satellite orbits

We consider a model that describes the evolution of distant satellite orbits and that refines the solution of the doubly averaged Hill problem. Generally speaking, such a refinement was performed previously by J. Kovalevsky and A.A. Orlov in terms of Zeipel’s method by constructing a solution of the third order with respect to the small parameter m, the ratio of the mean motions of the planet and the satellite. The analytical solution suggested here differs from the solutions obtained by these authors and is closest in form to the general solution of the doubly averaged problem (∼m ²). We have performed a qualitative analysis of the evolutionary equations and conditions for the intersection of satellite orbits with the surface of a spherical planet with a finite radius. Using the suggested solution, we have obtained improved analytical time dependences of the elements of evolving orbits for a number of distant satellites of giant planets compared to the solution of the doubly averaged Hill problem and, thus, achieved their better agreement with the results of our numerical integration of the rigorous equations of perturbed motion for satellites.     

On the migration of a system of protoplanets

The evolution of a system consisting of a protoplanetary disc with two embedded Jupiter-sized planets is studied numerically. The disc is assumed to be flat and non-self-gravitating; this is modelled by the planar (two-dimensional) Navier–Stokes equations. The mutual gravitational interaction of the planets and the star, and the gravitational torques of the disc acting on the planets and the central star are included. The planets have an initial mass of one Jupiter mass M _Jup each, and the radial distances from the star are one and two semimajor axes of Jupiter, respectively.
During the evolution a joint wide annular gap is created by the planets. Both planets increase their mass owing to accretion of gas from the disc: after about 2500 orbital periods of the inner planet it has reached a mass of 2.3  M _Jup, while the outer planet has reached a mass of 3.2  M _Jup. The net gravitational torques exerted by the disc on the planets result in an inward migration of the outer planet on time-scales comparable to the viscous evolution time of the disc. The semimajor axis of the inner planet remains constant as there is very little gas left in its vicinity to induce any migration. When the distance of close approach eventually becomes smaller than the mutual Hill radius, the eccentricities increase strongly and the system may become unstable.
If disc depletion occurs rapidly enough before the planets come too close to each other, a stable system similar to our own Solar system may remain. Otherwise the orbits may become unstable and produce systems like υ And.     

Kuiper Belt structure as a reflection of the migration process of a planet

One of the main particular features of the structure of the Kuiper Belt is that it contains clusters of objects of small orbital eccentricity and inclination (“cold population”). In order to solve the problem of the origin of the objects, we considered the process of the gravitational interaction of a comparatively small-mass planet with a planetesimal disk. We found that one particular property of the process is that the planet changes its direction of migration. The interaction with the planet results in the transportation of a considerable portion of planetesimals from the inner zone out to the Kuiper Belt. After such a transition of the objects, the planet returns to the inner regions of the planetesimal disk. Numerical simulations show that the reversible migration of a planet of a mass similar to that of the Earth can explain the main properties of the Kuiper Belt’s cold population orbit distribution.     

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.     

UV Capabilities to Probe the Formation of Planetary Systems: From the ISM to Planets

Planetary systems are angular momentum reservoirs generated during star formation. Solutions to three of the most important problems in contemporary astrophysics are needed to understand the entire process of planetary system formation: The physics of the ISM. Stars form from dense molecular clouds that contain ∼ 30% of the total interstellar medium (ISM) mass. The structure, properties and lifetimes of molecular clouds are determined by the overall dynamics and evolution of a very complex system – the ISM. Understanding the physics of the ISM is of prime importance not only for Galactic but also for extragalactic and cosmological studies. Most of the ISM volume (∼ 65%) is filled with diffuse gas at temperatures between 3000 and 300 000 K, representing about 50% of the ISM mass. The physics of accretion and outflow. Powerful outflows are known to regulate angular momentum transport during star formation, the so-called accretion–outflow engine. Elementary physical considerations show that, to be efficient, the acceleration region for the outflows must be located close to the star (within 1 AU) where the gravitational field is strong. According to recent numerical simulations, this is also the region where terrestrial planets could form after 1 Myr. One should keep in mind that today the only evidence for life in the Universe comes from a planet located in this inner disk region (at 1 AU) from its parent star. The temperature of the accretion–outflow engine is between 3000 and 10 ⁷ K. After 1 Myr, during the classical T Tauri stage, extinction is small and the engine becomes naked and can be observed at ultraviolet wavelengths. The physics of planet formation. Observations of volatiles released by dust, planetesimals and comets provide an extremely powerful tool for determining the relative abundances of the vaporizing species and for studying the photochemical and physical processes acting in the inner parts of young planetary systems. This region is illuminated by the strong UV radiation field produced by the star and the accretion–outflow engine. Absorption spectroscopy provides the most sensitive tool for determining the properties of the circumstellar gas as well as the characteristics of the atmospheres of the inner planets transiting the stellar disk. UV radiation also pumps the electronic transitions of the most abundant molecules (H ₂, CO, etc.) that are observed in the UV.Here we argue that access to the UV spectral range is essential for making progress in this field, since the resonance lines of the most abundant atoms and ions at temperatures between 3000 and 300 000 K, together with the electronic transitions of the most abundant molecules (H ₂, CO, OH, CS, S ₂, CO ₂ ⁺, C ₂, O ₂, O₃, etc.) are at UV wavelengths. A powerful UV-optical instrument would provide an efficient mean for measuring the abundance of ozone in the atmosphere of the thousands of transiting planets expected to be detected by the next space missions (GAIA, Corot, Kepler, etc.). Thus, a follow-up UV mission would be optimal for identifying Earth-like candidates.     

Gravitational interaction between a planet and an optically thin disc

We investigate the gravitational interaction between a planet and an optically thin protoplanetary disc, performing local three-dimensional hydrodynamical simulations. In the present study, we take account of radiative energy transfer in optically thin discs. Before the stage of planetary accretion, dust opacity is expected to decrease significantly because of grain growth and planetesimal formation. Thus, it would be reasonable to consider optically thin discs in the disc–planet interaction. Furthermore, we focus on small planets that can neither capture disc gas nor open a disc gap. The one-sided torque exerted on a planet by an optically thin disc is examined for various values of the disc optical thickness (<1). In optically thin discs, the temperature behind the density waves is lower than the unperturbed value because of radiative cooling. Heating due to shock dissipation is less effective than radiative cooling. Because of radiative cooling, the density distribution around the planet is not axisymmetric, which exerts an additional torque on the planet. The torque enhancement becomes maximum when the cooling time is comparable with the Keplerian period. The enhancement is significant for low-mass planets. For planets with  3 M_⊕  , the additional one-sided torque can be 40 per cent of the torque in the isothermal case. The radiative cooling is expected to change the differential torque and the migration speed of planets, too.     

Plasma interactions of exoplanets with their parent star and associated radio emissions

The relatively high contrast between planetary and solar low-frequency radio emissions suggests that the low-frequency radio range may be well adapted to the direct detection of exoplanets. We review the most significant properties of planetary radio emissions (auroral as well as satellite induced) and show that their primary engine is the interaction of a plasma flow with an obstacle in the presence of a strong magnetic field (of the flow or of the obstacle). Scaling laws have been derived from solar system planetary radio emissions that relate the emitted radio power to the power dissipated in the various corresponding flow–obstacle interactions. We generalize these scaling laws into a “radio-magnetic” scaling law that seems to relate output radio power to the magnetic energy flux convected on the obstacle, this obstacle being magnetized or unmagnetized. Extrapolating this scaling law to the case of exoplanets, we find that hot Jupiters may produce very intense radio emissions due to either magnetospheric interaction with a strong stellar wind or to unipolar interaction between the planet and a magnetic star (or strongly magnetized regions of the stellar surface). In the former case, similar to the magnetosphere–solar wind interactions in our solar system or to the Ganymede–Jupiter interaction, a hecto-decameter emission is expected in the vicinity of the planet with an intensity possibly 10³–10⁵ times that of Jupiter's low frequency radio emissions. In the latter case, which is a giant analogy of the Io–Jupiter system, emission in the decameter-to-meter wavelength range near the footprints of the star's magnetic field lines interacting with the planet may reach 10⁶ times that of Jupiter (unless some “saturation” mechanism occurs). The system of HD179949, where a hot spot has been tentatively detected in visible light near the sub-planetary point, is discussed in some details. Radio detectability is addressed with present and future low-frequency radiotelescopes. Finally, we discuss the interests of direct radio detection, among which access to exoplanetary magnetic field measurements and comparative magnetospheric physics.     

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.     

Effect of Planetary Eccentricity on Ballistic Capture in the Solar System

Ballistic capture of spacecraft and celestial bodies by planets of the solar system is studied considering the elliptic restricted three body model. A preferential region, due to the eccentricity of the planet and the Sun-gravity-gradient effect is found for the capture phenomenon. An analytical formula is derived which determines the limiting value of the satellite capture eccentricity e_c as a function of the pericenter distance x_p and planet’s true anomaly. The analytic values e_c are tested by a numerical propagator, which makes use of planetary ephemeris, and only a small difference with respect to numerical integration is found. It turns out that lower values of e_c occur when the planet anomaly is close to zero; that is, capture is easier when the planet is at its perihelion. This fact is confirmed by the capture of celestial bodies. It is shown that Jupiter comets are generally captured when Jupiter is in its perihelion region. Ballistic capture is also important in interplanetary missions. The propellant saved using the minimum ballistic capture eccentricity is evaluated for different missions and compared with respect to the case in which the insertion orbit is a parabola: a significant saving can be accomplished.     

CCD photometry inVRI bands of the galactic cluster NGC 2818

bdAbstract The open cluster NGC 2818 containing a planetary nebula has been observed inVRI bands using the CCD system at prime focus of the 2.3-metre Vainu Bappu Telescope. The study extending to starsV ∼ 21 magnitude establishes the distance modulus as(m-M) ₀ = 12.9 ±0.1 for the cluster. Based on the fitting of theoretical isochrones computed for solar metallicity, an age of 5(±1) × 10₈ years has been assigned to the cluster. Association of the planetary nebula with the cluster indicates that the progenitor mass of the planetary nebula on the main sequence is ≥2.5M_⊙ Based on observations obtained with the Vainu Bappu Telescope.     

The generalized non-conservative model of a 1-planet system revisited

We study the long-term dynamics of a planetary system composed of a star and a planet. Both bodies are considered as extended, non-spherical, rotating objects. There are no assumptions made on the relative angles between the orbital angular momentum and the spin vectors of the bodies. Thus, we analyze full, spatial model of the planetary system. Both objects are assumed to be deformed due to their own rotations, as well as due to the mutual tidal interactions. The general relativity corrections are considered in terms of the post-Newtonian approximation. Besides the conservative contributions to the perturbing forces, there are also taken into account non-conservative effects, i.e., the dissipation of the mechanical energy. This dissipation is a result of the tidal perturbation on the velocity field in the internal zones with non-zero turbulent viscosity (convective zones). Our main goal is to derive the equations of the orbital motion as well as the equations governing time-evolution of the spin vectors (angular velocities). We derive the Lagrangian equations of the second kind for systems which do not conserve the mechanical energy. Next, the equations of motion are averaged out over all fast angles with respect to time-scales characteristic for conservative perturbations. The final equations of motion are then used to study the dynamics of the non-conservative model over time scales of the order of the age of the star. We analyze the final state of the system as a function of the initial conditions. Equilibria states of the averaged system are finally discussed.     

Interaction of free-floating planets with a star–planet pair

The recent discovery of free-floating planets and their theoretical interpretation as celestial bodies, either condensed independently or ejected from parent stars in tight clusters, introduced an intriguing possibility. Namely, that some exoplanets are not condensed from the protoplanetary disk of their parent star. In this novel scenario a free-floating planet interacts with an already existing planetary system, created in a tight cluster, and is captured as a new planet. In the present work we study this interaction process by integrating trajectories of planet-sized bodies, which encounter a binary system consisting of a Jupiter-sized planet revolving around a Sun-like star. To simplify the problem we assume coplanar orbits for the bound and the free-floating planet and an initially parabolic orbit for the free-floating planet. By calculating the uncertainty exponent, a quantity that measures the dependence of the final state of the system on small changes of the initial conditions, we show that the interaction process is a fractal classical scattering. The uncertainty exponent is in the range (0.2–0.3) and is a decreasing function of time. In this way we see that the statistical approach we follow to tackle the problem is justified. The possible final outcomes of this interaction are only four, namely flyby, planet exchange, capture or disruption. We give the probability of each outcome as a function of the incoming planet’s mass. We find that the probability of exchange or capture (in prograde as well as retrograde orbits and for very long times) is non-negligible, a fact that might explain the possible future observations of planetary systems with orbits that are either retrograde (see e.g. Queloz et?al. Astron. Astrophys. 417, L1, 2010) or tight and highly eccentric.     

Holomorphy of coordinates with respect to eccentricities in the planetary three-body problem

The boundaries of the domains of holomorphy of the coordinates of unperturbed elliptic motion with respect to the eccentricities of planetary orbits are determined for the cases when any of the five anomalies of one of the planets-eccentric, true, tangential, or one of two mutual anomalies suggested by M.F. Subbotin—is used as an independent variable. The resulting equations are a generalization of the known equations for the boundaries of the domains of the holomorphy of coordinates for the cases when the time is the independent variable and determine the bisymmetric ovals, whose size and shape depend on the eccentricities and on the ratio of the planetary mean motions. The largest domains of holomorphy are obtained when the tangential anomaly or one of the Subbotin mutual anomalies is used. A function was found that conformally maps the domain of holomorphy to the unit disk. It was demonstrated that the application of any anomaly of the outer planet as the independent variable can result in a significant shrinking of the domain of the holomorphy of the coordinates of the inner planet, so that the analytic continuation of the initial power series with the center at the origin of the coordinates of a complex plane becomes impossible.