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.     

Analytical solution for the three-dimensional motion of a circumbinary planet around a binary star

By using the method of separating rapid and slow subsystem, we obtain an analytical solution for a stable three-dimensional motion of a circumbinary planet around a binary star. We show that the motion of the planet is more complicated than it was obtained for this situation analytically by Farago and Laskar (2010). Namely, in addition to the precession of the orbital plane of the planet around the angular momentum of the binary (found by Farago and Laskar (2010)), there is simultaneously the precession of the orbital plane of the planet within the orbital plane. We show that the frequency of this additional precession is different from the frequency of the precession of the orbital plane around the angular momentum of the binary. We demonstrate that this problem is mathematically equivalent both to the problem of the motion of a satellite around an oblate planet and to the problem of a hydrogen Rydberg atom in the field of a high-frequency linearly-polarized laser radiation, thus discovering yet another connection between astrophysics and atomic physics. We point out that all of the above physical systems have a higher than geometrical symmetry, which is a counterintuitive result. In particular, it is manifested by the fact that, while the elliptical orbit of the circumbinary planet (around a binary star) or of the satellite (around an oblate planet) or of the Rydberg electron (in the laser field) undergoes simultaneously two types of the precession, the shape of the orbit does not change. The fact that a system, consisting of a circumbinary planet around a binary star, possesses the hidden symmetry should be of a general physical interest. Our analytical results could be used for benchmarking future simulations.     

Precession of the orbital nodes of Jupiter and Saturn triggered by the mutual perturbation: A model of two rings

The problem of the precession of the orbital planes of Jupiter and Saturn under the influence of mutual gravitational perturbations was formulated and solved using a simple dynamical model. Using the Gauss method, the planetary orbits are modeled by material circular rings, intersecting along the diameter at a small angle α. The planet masses, semimajor axes and inclination angles of orbits correspond to the rings. What is new is that each ring has an angular momentum equal to the orbital angular momentum of the planet. Contrary to popular belief, it was proved that the orbital resonance 5: 2 does not preclude the use of the ring model. Moreover, the period of averaging of the disturbing force (T ≈ 1332 yr) proves to be appreciably greater than a conventionally used period (≈900 yr). The mutual potential energy of rings and the torque of gravitational forces between the rings were calculated. We compiled and solved the system of differential equations for the spatial motion of rings. It was established that a perturbing torque causes the precession and simultaneous rotation of the orbital planes of Jupiter and Saturn. Moreover, the opposite orbit nodes on the Laplace plane coincide and perform a secular movement in retrograde direction with the same velocity of 25.6″/yr and the period T _J = T _S ≈ 50687 yr. These results are close to those obtained in the general theory (25.93″/yr), which confirms the adequacy of the developed model. It was found that the vectors of the angular velocity of orbital rings move counterclockwise over circular cones and describe circles on the celestial sphere with radii β₁ ≈ 0.8403504° (Saturn) and β₂ ≈ 0.3409296° (Jupiter) around the point which is located at an angular distance of 1.647607° from the ecliptic pole.     

Resonant and non-resonant gravity-gradient perturbations of a tumbling tri-axial satellite

Gravity-gradient perturbations of the attitude motion of a tumbling tri-axial satellite are investigated. The satellite center of mass is considered to be in an elliptical orbit about a spherical planet and to be tumbling at a frequency much greater than orbital rate. In determining the unperturbed (free) motion of the satellite, a canonical form for the solution of the torque-free motion of a rigid body is obtained. By casting the gravity-gradient perturbing torque in terms of a perturbing Hamiltonian, the long-term changes in the rotational motion are derived. In particular, far from resonance, there are no long-period changes in the magnitude of the rotational angular momentum and rotational energy, and the rotational angular momentum vector precesses abound the orbital angular momentum vector.At resonance, a low-order commensurability exists between the polhode frequency and tumbling frequency. Near resonance, there may be small long-period fluctuations in the rotational energy and angular momentum magnitude. Moreover, the precession of the rotational angular momentum vector about the orbital angular momentum vector now contains substantial long-period contributions superimposed on the non-resonant precession rate. By averaging certain long-period elliptic functions, the mean value near resonance for the precession of the rotational angular momentum vector is obtained in terms of initial conditions.     

Polar orbits around binary stars

Oks proposes the existence of a new class of stable planetary orbits around binary stars, in the shape of a helix on a conical surface whose axis of symmetry coincides with the interstellar axis, and rotates with the same orbital frequency as the binary pair. We show that this claim relies on the inappropriate use of an effective potential that is only applicable when the stars are held motionless. We also present numerical evidence that the only planetary orbits whose planes are initially orthogonal to the interstellar axis that remain stable on the time scale of the stellar orbit are ordinary polar orbits around one of the stars, and that the perturbations due to the binary companion do not rotate the plane of the orbit to maintain a fixed relationship with the axis.     

On the stability of circumbinary planetary systems

The dynamics of circumbinary planetary systems (the systems in which the planets orbit a central binary) with a small binary mass ratio discovered to date is considered. The domains of chaotic motion have been revealed in the “pericentric distance–eccentricity” plane of initial conditions for the planetary orbits through numerical experiments. Based on an analytical criterion for the chaoticity of planetary orbits in binary star systems, we have constructed theoretical curves that describe the global boundary of the chaotic zone around the central binary for each of the systems. In addition, based on Mardling’s theory describing the separate resonance “teeth” (corresponding to integer resonances between the orbital periods of a planet and the binary), we have constructed the local boundaries of chaos. Both theoretical models are shown to describe adequately the boundaries of chaos on the numerically constructed stability diagrams, suggesting that these theories are efficient in providing analytical criteria for the chaoticity of planetary orbits.     

On the relation between resonance and instability in planetary systems

The factors which affect the linear stability of a periodic planetary orbit in the plane are studied. It is proved that planetary systems with two planets describing nearly circular orbits in the same direction are linearly stable and no perturbation exists which destroys the stability, unless a resonance of the form 1/3, 3/5, 5/7, ... among the orbits of the planets occurs. This latter resonant case is always unstable. Retrograde motion is always linearly stable. Planetary systems with three or more planets in nearly circular orbits in the same direction are proved to be unstable, in the sense that a Hamiltonian perturbation always exists which destroys the stability. The generation of instability in the case of three or more planets is not only due to the existence of resonances, as in the case of two planets, but also to the nonexistence of integrals of motion, apart from the energy and angular momentum integrals. It is also proved that planetary systems with nearly elliptic orbits of the planets are unstable.     

Obliquity variations of terrestrial planets in habitable zones

We have investigated obliquity variations of possible terrestrial planets in habitable zones (HZs) perturbed by a giant planet(s) in extrasolar planetary systems. All the extrasolar planets so far discovered are inferred to be jovian-type gas giants. However, terrestrial planets could also exist in extrasolar planetary systems. In order for life, in particular for land-based life, to evolve and survive on a possible terrestrial planet in an HZ, small obliquity variations of the planet may be required in addition to its orbital stability, because large obliquity variations would cause significant climate change. It is known that large obliquity variations are caused by spin-orbit resonances where the precession frequency of the planet's spin nearly coincides with one of the precession frequencies of the ascending node of the planet's orbit. Using analytical expressions, we evaluated the obliquity variations of terrestrial planets with prograde spins in HZs. We found that the obliquity of terrestrial planets suffers large variations when the giant planet's orbit is separated by several Hill radii from an edge of the HZ, in which the orbits of the terrestrial planets in the HZ are marginally stable. Applying these results to the known extrasolar planetary systems, we found that about half of these systems can have terrestrial planets with small obliquity variations (smaller than 10°) over their entire HZs. However, the systems with both small obliquity variations and stable orbits in their HZs are only 1/5 of known systems. Most such systems are comprised of short-period giant planets. If additional planets are found in the known planetary systems, they generally tend to enhance the obliquity variations. On the other hand, if a large/close satellite exists, it significantly enhances the precession rate of the spin axis of a terrestrial planet and is likely to reduce the obliquity variations of the planet. Moreover, if a terrestrial planet is in a retrograde spin state, the spin-orbit resonance does not occur. Retrograde spin, or a large/close satellite might be essential for land-based life to survive on a terrestrial planet in an HZ.     

On possible circumbinary configurations of the planetary systems of α Centauri and EZ Aquarii

Possible configurations of the planetary systems of the binary stars α Cen A–BandEZAqr A–C are analyzed. The P-type orbits—circumbinary ones, i.e., the orbits around both stars of the binary, are studied. The choice of these systems is dictated by the fact that α Cen is closest to us in the Galaxy, while EZ Aqr is the closest system whose circumbinary planets, as it turns out, may reside in the “habitability zone.” The analysis has been performed within the framework of the planar restricted three-body problem. The stability diagrams of circumbinary motion have been constructed: on representative sets of initial data (in the pericentric distance–eccentricity plane), we have computed the Lyapunov spectra of planetary motion and identified the domains of regular and chaotic motion through their statistical analysis. Based on present views of the dynamics and architecture of circumbinary planetary systems, we have determined the most probable planetary orbits to be at the centers of the main resonance cells, at the boundary of the dynamical chaos domain around the parent binary star, which allows the semimajor axes of the orbits to be predicted. In the case of EZ Aqr, the orbit of the circumbinary planet is near the habitability zone and, given that the boundary of this zone is uncertain, may belong to it.     

Gravitational loss-cone instability in stellar systems with retrograde orbit precession

We study spherical and disc clusters in a near-Keplerian potential of galactic centres or massive black holes. In such a potential orbit precession is commonly retrograde, that is, the direction of the orbit precession is opposite to the orbital motion. It is assumed that stellar systems consist of nearly-radial orbits. We show that if there is a loss-cone at low angular momentum (e.g. due to consumption of stars by a black hole), an instability similar to loss-cone instability in plasma may occur. The gravitational loss-cone instability is expected to enhance black hole feeding rates. For spherical systems, the instability is possible for the number of spherical harmonics   l ≥ 3  . If there is some amount of counter-rotating stars in flattened systems, they generally exhibit the instability independent of azimuthal number m . The results are compared with those obtained recently by Tremaine for distribution functions monotonically increasing with angular momentum.
The analysis is based on simple characteristic equations describing small perturbations in a disc or a sphere of stellar orbits highly elongated in radius. These characteristic equations are derived from the linearized Vlasov equations (combining the collisionless Boltzmann kinetic equation and the Poisson equation), using the action-angle variables. We use two techniques for analysing the characteristic equations: the first one is based on preliminary finding of neutral modes, and the second one employs a counterpart of the plasma Penrose–Nyquist criterion for disc and spherical gravitational systems.     

The 1/1 resonance in extrasolar planetary systems

We study orbits of planetary systems with two planets, for planar motion, at the 1/1 resonance. This means that the semimajor axes of the two planets are almost equal, but the eccentricities and the position of each planet on its orbit, at a certain epoch, take different values. We consider the general case of different planetary masses and, as a special case, we consider equal planetary masses. We start with the exact resonance, which we define as the 1/1 resonant periodic motion, in a rotating frame, and study the topology of the phase space and the long term evolution of the system in the vicinity of the exact resonance, by rotating the orbit of the outer planet, which implies that the resonance and the eccentricities are not affected, but the symmetry is destroyed. There exist, for each mass ratio of the planets, two families of symmetric periodic orbits, which differ in phase only. One is stable and the other is unstable. In the stable family the planetary orbits are in antialignment and in the unstable family the planetary orbits are in alignment. Along the stable resonant family there is a smooth transition from planetary orbits of the two planets, revolving around the Sun in eccentric orbits, to a close binary of the two planets, whose center of mass revolves around the Sun. Along the unstable family we start with a collinear Euler–Moulton central configuration solution and end to a planetary system where one planet has a circular orbit and the other a Keplerian rectilinear orbit, with unit eccentricity. It is conjectured that due to a migration process it could be possible to start with a 1/1 resonant periodic orbit of the planetary type and end up to a satellite-type orbit, or vice versa, moving along the stable family of periodic orbits.     

Secular and tidal evolution of circumbinary systems

We investigate the secular dynamics of three-body circumbinary systems under the effect of tides. We use the octupolar non-restricted approximation for the orbital interactions, general relativity corrections, the quadrupolar approximation for the spins, and the viscous linear model for tides. We derive the averaged equations of motion in a simplified vectorial formalism, which is suitable to model the long-term evolution of a wide variety of circumbinary systems in very eccentric and inclined orbits. In particular, this vectorial approach can be used to derive constraints for tidal migration, capture in Cassini states, and stellar spin–orbit misalignment. We show that circumbinary planets with initial arbitrary orbital inclination can become coplanar through a secular resonance between the precession of the orbit and the precession of the spin of one of the stars. We also show that circumbinary systems for which the pericenter of the inner orbit is initially in libration present chaotic motion for the spins and for the eccentricity of the outer orbit. Because our model is valid for the non-restricted problem, it can also be applied to any three-body hierarchical system such as star–planet–satellite systems and triple stellar systems.     

The structure of motion in a 4-component galaxy mass model

We use a composite galaxy model consisting of a disk-halo, bulge, nucleus and dark-halo components in order to investigate the motion of stars in ther-z plane. It is observed that high angular momentum stars move in regular orbits. The majority of orbits are box orbits. There are also banana-like orbits. For a given value of energy, only a fraction of the low angular momentum stars — those going near the nucleus — show chaotic motion while the rest move in regular orbits. Again one observes the above two kinds of orbits. In addition to the above one can also see orbits with the characteristics of the 2/3 and 3/4 resonance. It is also shown that, in the absence of the bulge component, the area of chaotic motion in the surface of section increases, significantly. This suggests that a larger number of low angular momentum stars are in chaotic orbits in galaxies with massive nuclei and no bulge components.     

Modulation of circumstellar extinction in a young binary system with a low-mass companion in a noncoplanar orbit

We consider a model of cyclic brightness variations in a young star with a low-mass (q = M ₂/M ₁ ≤ 0.1) companion that accretes matter from the remnants of a protostellar cloud (circumbinary disk). We assume that the orbit of the companion is circular and that its plane does not coincide with the disk plane. We have computed grids of hydrodynamic models for such a binary by the SPH method based on which we have investigated the circumstellar extinction variations produced by the streams of matter and density waves excited in the circumbinary disk by the orbital motion of the companion. We show that, depending on the inclination and orientation of the binary’s line of nodes relative to the observer, the brightness of the primary component can undergo various (in shape and depth) oscillations with a period equal to the orbital one. In contrast to the models with coplanar circular orbits, the accretion rate onto the components of a binary with a noncoplanar orbit depends on the orbital phase. The results of our computations can be used to study the cyclic activity of UX Ori stars and young eclipsing binaries with anomalously long eclipses.     

Eccentric Extrasolar Planets: The Jumping Jupiter Model

Most extrasolar planets discovered to date are more massive than Jupiter, in surprisingly small orbits (semimajor axes less than 3 AU). Many of these have significant orbital eccentricities. Such orbits may be the product of dynamical interactions in multiplanet systems. We examine outcomes of such evolution in systems of three Jupiter-mass planets around a solar-mass star by integration of their orbits in three dimensions. Such systems are unstable for a broad range of initial conditions, with mutual perturbations leading to crossing orbits and close encounters. The time scale for instability to develop depends on the initial orbital spacing; some configurations become chaotic after delays exceeding 10⁸ y. The most common outcome of gravitational scattering by close encounters is hyperbolic ejection of one planet. Of the two survivors, one is moved closer to the star and the other is left in a distant orbit; for systems with equal-mass planets, there is no correlation between initial and final orbital positions. Both survivors may have significant eccentricities, and the mutual inclination of their orbits can be large. The inner survivor's semimajor axis is usually about half that of the innermost starting orbit. Gravitational scattering alone cannot produce the observed excess of “hot Jupiters” in close circular orbits. However, those scattered planets with large eccentricities and small periastron distances may become circularized if tidal dissipation is effective. Most stars with a massive planet in an eccentric orbit should have at least one additional planet of comparable mass in a more distant orbit.     

Retrograde precession of stellar orbits and gravitational loss-cone instability in galactic centers

We consider disk and spherical subsystems of stars with nearly radial orbits under conditions when the well-known radial orbit instability is not possible. This requires that the precession of stellar orbits be retrograde, i.e., in the direction opposite to the orbital rotation of stars. We show that an instability that is an analogue of the loss-cone instability known in plasma physics can then develop in the presence of a “loss cone” in the angular momentum distribution of stars, which ensures a deficit or even absence of stars with low angular momenta. Examples of systems with a loss cone are the centers of galaxies or star clusters with massive black holes. The instability can produce a flux of stars onto the galactic center, i.e., it can serve as a mechanism of fueling the nuclear activity of galaxies. Mathematically, the problem is reduced to analyzing simple characteristic equations that describe small perturbations in a disk and a sphere of radially highly elongated stellar orbits. In turn, these characteristics equations are derived through a number of successive simplifications of the general linearized Vlasov equations (i.e., the system that includes the collisionless Boltzmann kinetic equation and the Poisson equation) in action—angle variables.     

Analytical direct planetary perturbations on the orbital motion of satellites

An analytical expansion of the disturbing function arising from direct planetary perturbations on the motion of satellites is derived. As a Fourier series, it allows the investigation of the secular effects of these direct perturbations, as well as of every argument present in the perturbation. In particular, we construct an analytical model describing the evection resonance between the longitude of pericenter of the satellite orbit and the longitude of a planet, and study briefly its dynamic. The expansion developed in this paper is valid in the case of planar and circular planetary orbits, but not limited in eccentricity or inclination of the satellite orbit.     

Planetary dynamics in the system α Centauri: The stability diagrams

The stability of the motion of a hypothetical planet in the binary system ?? Cen A?CB has been investigated. The analysis has been performed within the framework of a planar (restricted and full) three-body problem for the case of prograde orbits. Based on a representative set of initial data, we have obtained the Lyapunov spectra of the motion of a triple system with a single planet. Chaotic domains have been identified in the pericenter distance-eccentricity plane of initial conditions for the planet through a statistical analysis of the data obtained. We have studied the correspondence of these chaotic domains to the domains of initial conditions that lead to the planet??s encounter with one of the binary??s stars or to the escape of the planet from the system. We show that the stability criterion based on the maximum Lyapunov exponent gives a more clear-cut boundary of the instability domains than does the encounterescape criterion at the same integration time. The typical Lyapunov time of chaotic motion is ??500 yr for unstable outer orbits and ??60 yr for unstable inner ones. The domain of chaos expands significantly as the initial orbital eccentricity of the planet increases. The chaos-order boundary has a fractal structure due to the presence of orbital resonances.     

Nodal and periastron precession of inclined orbits in the field of a rotating black hole

The inclination of low-eccentricity orbits is shown to significantly affect orbital parameters, in particular, the Keplerian, nodal precession, and periastron rotation frequencies, which are interpreted in terms of observable quantities. For the nodal precession and periastron rotation frequencies of low-eccentricity orbits in a Kerr field, we derive a Taylor expansion in terms of the Kerr parameter at arbitrary orbital inclinations to the black-hole spin axis and at arbitrary radial coordinates. The particle radius, energy, and angular momentum in the marginally stable circular orbits are calculated as functions of the Kerr parameter j and parameter s in the form of Taylor expansions in terms of j to within O[j ⁶]. By analyzing our numerical results, we give compact approximation formulas for the nodal precession frequency of the marginally stable circular orbits at various s in the entire range of the Kerr parameter.