A dynamical analysis of the 14 Herculis planetary system

Precision radial velocity measurements of the Sun-like dwarf 14 Herculis published by Naef et al., Butler et al. and Wittenmyer, Endl & Cochran reveal a Jovian planet in a 1760-d orbit and a trend indicating the second distant object. On the grounds of dynamical considerations, we test a hypothesis that the trend can be explained by the presence of an additional giant planet. We derive dynamical limits to the orbital parameters of the putative outer Jovian companion in an orbit within ∼13 au. In this case, the mutual interactions between the Jovian planets are important for the long-term stability of the system. The best self-consistent and stable Newtonian fit to an edge-on configuration of Jovian planets has the outer planet in 9-au orbit with a moderate eccentricity of ∼0.2 and confined to a zone spanned by the low-order mean motion resonances 5:1 and 6:1. This solution lies in a shallow minimum of (χ²_ν)^1/2 and persists over a wide range of the system inclination. Other stable configurations within 1σ confidence interval of the best fit are possible for the semimajor axis of the outer planet in the range of (6,13) au and the eccentricity in the range of (0, 0.3). The orbital inclination cannot yet be determined but when it decreases, both planetary masses approach ∼10 m _J and for i ∼ 30° the hierarchy of the masses is reversed.     

GJ 581 update: Additional evidence for a Super‐Earth in the habitable zone

We present an analysis of the significantly expanded HARPS 2011 radial velocity data set for GJ 581 that was presented by Forveille et al. (2011). Our analysis reaches substantially different conclusions regarding the evidence for a Super‐Earth‐mass planet in the star's Habitable Zone. We were able to reproduce their reported χ²_ν and RMS values only after removing some outliers from their models and refitting the trimmed down RV set. A suite of 4000 N‐body simulations of their Keplerian model all resulted in unstable systems and revealed that their reported 3.6σ detection of e = 0.32 for the eccentricity of GJ 581e is manifestly incompatible with the system's dynamical stability. Furthermore, their Keplerian model, when integrated only over the time baseline of the observations, significantly increases the χ²_ν and demonstrates the need for including non‐Keplerian orbital precession when modeling this system. We find that a four‐planet model with all of the planets on circular or nearly circular orbits provides both an excellent self‐consistent fit to their RV data and also results in a very stable configuration. The periodogram of the residuals to a 4‐planet all‐circular‐orbit model reveals significant peaks that suggest one or more additional planets in this system. We conclude that the present 240‐point HARPS data set, when analyzed in its entirety, and modeled with fully self‐consistent stable orbits, by and of itself does offer significant support for a fifth signal in the data with a period near 32 days. This signal has a false alarm probability of <4% and is consistent with a planet of minimum mass 2.2 M_⊕, orbiting squarely in the star's habitable zone at 0.13 AU, where liquid water on planetary surfaces is a distinct possibility (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new analysis of the GJ581 extrasolar planetary system

We have done a new analysis of the available observations of the GJ581 exoplanetary system. Today this system is controversial due to choices that can be done in the orbital determination. The main ones are the occurrence of aliases and the additional bodies??the planets f and g??announced in Vogt et?al. (Astrophys J 723:954?C965, 2010). Any dynamical study of exoplanets requires the good knowledge of the orbital elements and the investigations involving the planet g are particularly interesting, since this body would lie in the habitable zone (HZ) of the star GJ581. This region, for this system, is very attractive of the dynamical point of view due to several resonances of two and three bodies present there. In this work, we investigate the conditions under which the planet g may exist. We stress the fact that the planet g is intimately related with the orbital elements of the planet d; more precisely, we conclude that it is not possible to disconnect its existence from the determination of the eccentricity of the planet d. Concerning the planet f, we have found one solution with period ??450?days, but we are judicious about any affirmation concerning this body because its signal is in the threshold of detection and the high period is in a spectral region where the occurrence of aliases is very common. Besides, we outline some dynamical features of the HZ with the dynamical map and point out the role played by some resonances laying there.     

Can GJ 876 host four planets in resonance?

Prior to the detection of its outermost Uranus-mass object, it had been suggested that GJ 876 could host an Earth-sized planet in a 15-day orbit. Observation, however, did not support this idea, but instead revealed evidence for the existence of a larger body in a ~125-day orbit, near a three-body resonance with the two giant planets of this system. In this paper, we present a detailed analysis of the dynamics of the four-planet system of GJ 876, and examine the possibility of the existence of other planetary objects interior to its outermost body. We have developed a numerical scheme that enables us to search the orbital parameter-space very effectively and, in a short time, identify regions where an object may be stable. We present details of this integration method and discuss its application to the GJ 876 four-planet system. The results of our initial analysis suggested possible stable orbits at regions exterior to the orbit of the outermost planet and also indicated that an island of stability may exist in and around the 15-day orbit. However, examining the long-term stability of an object in that region by direct integration revealed that the 15-day orbit becomes unstable and that the system of GJ 876 is most likely dynamically full. We present the results of our study and discuss their implications for the formation and final orbital architecture of this system.     

Global dynamics of the Gliese 876 planetary system

The Gliese 876 planetary system consists of two Jupiter-like planets having a nearly commensurate 2:1 orbital periods ratio. Because the semimajor axes of the planets are very small (of the order 0.1 au and 0.2 au, respectively), and the eccentricity of the inner companion is ≃0.3, the mutual perturbations are extremely large. However, many authors claim the long-term orbital stability of the system, at least over 500 Myr for initial conditions found by Rivera & Lissauer. Results of investigations of a migration of initially separated planets into the close 2:1 mean motion resonance lock from Lee & Peale also support the conclusion that the system should be stable for the lifetime of the parent star. Initial conditions of the system, found from non-linear N -body fits by Laughlin & Chambers and Rivera & Lissauer, to the radial velocity curve, formally allow for a variety of orbital configurations of the GJ 876 system, e.g. coplanar, with planetary inclinations in the range [≃30°, 90°], and with relative inclinations of orbital planes as high as 80°. Our work is devoted to the stability investigation of the systems originating from the fitted initial conditions. We study neighbourhoods of these initial states in the orbital parameter space. We found estimations of the 2:1 mean motion resonance width and dynamical limitations on the planetary masses. We also obtain a global representation of the domains of the orbital parameters space in which initial conditions leading to stable evolutions can be found. Our results can be useful in localization of the best, stable fits to the observational data. In our investigations we use the MEGNO technique (the Mean Exponential Growth factor of Nearby Orbits) invented by Cincotta & Simó. It allows us to distinguish efficiently and precisely between chaotic and regular behaviour of a planetary system.     

On the chaotic diffusion in multidimensional Hamiltonian systems

We present numerical evidence that diffusion in the herein studied multidimensional near-integrable Hamiltonian systems departs from a normal process, at least for realistic timescales. Therefore, the derivation of a diffusion coefficient from a linear fit on the variance evolution of the unperturbed integrals fails. We review some topics on diffusion in the Arnold Hamiltonian and yield numerical and theoretical arguments to show that in the examples we considered, a standard coefficient would not provide a good estimation of the speed of diffusion. However, numerical experiments concerning diffusion would provide reliable information about the stability of the motion within chaotic regions of the phase space. In this direction, we present an extension of previous results concerning the dynamical structure of the Laplace resonance in Gliese-876 planetary system considering variations of the orbital parameters accordingly to the error introduced by the radial velocity determination. We found that a slight variation of the eccentricity of planet c would destabilize the inner region of the resonance that, though chaotic, shows stable when adopting the best fit values for the parameters.     

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).     

A study on dynamic processes at late stages in the formation of planetary systems in gas and dust disks

The discovered exoplanetary systems have highly diverse dynamic properties, which differ from those of the Solar System. A single model including planet migration effects and their gravitational interaction is used to investigate the features of dynamic processes that lead to the formation of giant-planet systems with different orbital characteristics. It is shown for a system of four giant planets similar to the Solar System how Type I migration could lead to all the planets being captured into resonant configurations. The resonant motion can continue for a long period of time after the transition to Type II migration and after the dissipation of the gas-and-dust disk. The three-planet system of GJ 876 is used to investigate the migration of the planets inward the orbit of the most massive planet and their capture into low-order resonant configurations under the conditions of Type II migration. A system similar to the exoplanetary system of HD 102272 is used to study the capture into high-order resonances followed by an increase in the orbit’s eccentricity.     

Eccentricity evolution of giant planet orbits due to circumstellar disk torques

The extrasolar planets discovered to date possess unexpected orbital elements. Most orbit their host stars with larger eccentricities and smaller semi-major axes than similarly sized planets in our own Solar System do. It is generally agreed that the interaction between giant planets and circumstellar disks (Type II migration) drives these planets inward to small radii, but the effect of these same disks on orbital eccentricity, ?, is controversial. Several recent analytic calculations suggest that disk-planet interactions can excite eccentricity, while numerical studies generally produce eccentricity damping. This paper addresses this controversy using a quasi-analytic approach, drawing on several preceding analytic studies. This work refines the current treatment of eccentricity evolution by removing several approximations from the calculation of disk torques. We encounter neither uniform damping nor uniform excitation of orbital eccentricity, but rather a function d?/dt that varies in both sign and magnitude depending on eccentricity and other Solar System properties. Most significantly, we find that for every combination of disk and planet properties investigated herein, corotation torques produce negative values of d?/dt for some range in ? within the interval [0.1, 0.5]. If corotation torques are saturated, this region of eccentricity damping disappears, and excitation occurs on a short timescale of less than 0.08 Myr. Thus, our study does not produce eccentricity excitation on a timescale of a few Myr—we obtain either eccentricity excitation on a short time scale, or eccentricity damping on a longer time scale. Finally, we discuss the implications of this result for producing the observed range in extrasolar planet eccentricity.     

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.     

A Study on the Characteristics of the Structure of Vega's Debris Disk

The clumpy structure in the Vega's debris disk was reported at millimeter wavelengths previously, and attributed to the concentration of dust grains trapped in resonances with a potential high-eccentricity planet. However, current imaging at multi-wavelengths with higher sensitivity indicates that the Vega's debris disk has a smooth structure. But a planet orbiting Vega could not be neglected, and the present-day observations may place a severe constraint on the orbital parameters for the potential planet. Herein, we utilize the modi- fied MERCURY codes to numerically simulate the Vega system, which consists of a debris disk and a planet. In our simulations, the initial inner and outer boundaries of the debris disk are assumed to be 80 AU and 120 AU, respectively. The dust grains in the disk have the sizes from 10 μm to 100 μm, and the nearly coplanar orbits. From the outcomes, we show that the evolution of debris disk is consistent with recent observations, if there is no planet orbiting Vega. However, if Vega owns a planet with a high eccentricity (e.g., e = 0.6), the planet's semi- major axis cannot be larger than 60 AU, otherwise, an aggregation phenomenon will occur in the debris disk due to the existence of the postulated planet. In addition, the 2:1 mean motion resonances may play a significant role in forming the structure of debris disk.     

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.     

Radiation pressure and dust particle dynamics

The dynamics of small dust grains orbiting a planet are investigated when solar radiation pressure forces are added to the planet's gravitational central field. In the first part a set of differential equations is derived in a reference frame linked to the solar motion. The complete solution of these equations is given for particles lying in the planet's orbital plane, and we show that the orbital eccentricity may undergo considerable variation. At the same time the pericenter longitude librates or circulates according to initial conditions. With this result we establish a criterion for any orbiting particle (because of its highly eccentric orbit) to collide with its planet's atmosphere. The case of inclined orbit is studied through a numerical integration and allows us to draw conclusions related to the stability of the orbital plane. All solutions are periodic, with the period being independent of the initial conditions. This last point permits us to investigate the different time scales involved in that problem. Finally, the Poynting-Robertson drag is included, along with the radial radiation pressure forces, and the secular trend is considered. A coupling effect between the two components is ascertained, yielding a systematic behavior in the eccentricity and thus in the pericenter distance. Our solutions generalize the results of S. J. Peale (1966, J. Geophys. Res.71, 911–933) and J. A. Burns, P. Lamy, and S. Soter (1979, Icarus40, 1–48) by allowing eccentricities to be large (of order 1) and inclinations to be nonzero and by considering Poynting-Robertson drag.     

On the chaotic orbital dynamics of the planet in the system 16 Cyg

The chaotic orbital dynamics of the planet in the wide visual binary star system 16 Cyg is considered. The only planet in this system has a significant orbital eccentricity, e = 0.69. Previously, Holman et al. suggested the possibility of chaos in the orbital dynamics of the planet due to the proximity of 16 Cyg to the separatrix of the Lidov–Kozai resonance. We have calculated the Lyapunov characteristic exponents on the set of possible orbital parameters for the planet. In all cases, the dynamics of 16 Cyg is regular with a Lyapunov time of more than 30 000 yr. The dynamics is considered in detail for several possible models of the planetary orbit; the dependences of Lyapunov exponents on the time of their calculation and the time dependences of osculating orbital elements have been constructed. Phase space sections for the system dynamics near the Lidov–Kozai resonance have been constructed for all models. A chaotic behavior in the orbital motion of the planet in 16 Cyg is shown to be unlikely, because 16 Cyg in phase space is far from the separatrix of the Lidov–Kozai resonance at admissible orbital parameters, with the chaotic layer near the separatrix being very narrow.     

Five new and three improved mutual orbits of transneptunian binaries

We present three improved and five new mutual orbits of transneptunian binary systems (58534) Logos-Zoe, (66652) Borasisi-Pabu, (88611) Teharonhiawako-Sawiskera, (123509) 2000 WK₁₈₃, (149780) Altjira, 2001 QY₂₉₇, 2003 QW₁₁₁, and 2003 QY₉₀ based on Hubble Space Telescope and Keck II laser guide star adaptive optics observations. Combining the five new orbit solutions with 17 previously known orbits yields a sample of 22 mutual orbits for which the period P, semimajor axis a, and eccentricity e have been determined. These orbits have mutual periods ranging from 5 to over 800 days, semimajor axes ranging from 1600 to 37,000 km, eccentricities ranging from 0 to 0.8, and system masses ranging from 2 × 10¹⁷ to 2 × 10²² kg. Based on the relative brightnesses of primaries and secondaries, most of these systems consist of near equal-sized pairs, although a few of the most massive systems are more lopsided. The observed distribution of orbital properties suggests that the most loosely-bound transneptunian binary systems are only found on dynamically cold heliocentric orbits. Of the 22 known binary mutual orbits, orientation ambiguities are now resolved for 9, of which 7 are prograde and 2 are retrograde, consistent with a random distribution of orbital orientations, but not with models predicting a strong preference for retrograde orbits. To the extent that other perturbations are not dominant, the binary systems undergo Kozai oscillations of their eccentricities and inclinations with periods of the order of tens of thousands to millions of years, some with strikingly high amplitudes.     

The effect of orbital element variations on the mean seasonal daily insolation on Mars

In this paper we briefly study changes in the mean seasonal insolations on the planet Mars caused by significant large-scale variations in the following orbital elements: the eccentricity (e), the obliquity (ε) and the longitude of perihelion (λ _p ). Three orbital configurations have been investigated. In the first, the eccentricity equals successively 0, 0.075, and 0.15, whereas for the obliquity and the longitude of perihelion we took the present values which amount, respectively to 25° and 250°. In the second situation, ε=15, 25, and 35° for a circular orbit (e=0) and with λ _p =250°. In the last model we have sete=0.075 and ε=25° for λ _p =?90,0, and 90°. Although long-term periodic oscillations ofe (first case) and λ _p (third case) produce, respectively, very small or no variations in the average yearly insolation, fluctuations of the above mentioned planetary data strongly effect the mean summer and winter daily insolations. Indeed, the calculations reveal that between the two extreme values of the orbital elements used, the seasonal insolations exhibit a change in amplitude of about 15 to 20% difference over the entire latitude interval. Considering more particularly the second case it is found that the summertime insolation experiences a nearly similar variation as the mean annual daily insolation — i.e., a decrease of about 7% at the equator and a more than twofold increase at the poles. The corresponding mean winter daily insolation varies maximally by approximately 60% in the 60–80° latitude range.     

The Effect of Tidal Interaction with a Gas Disk on Formation of Terrestrial Planets

We have performed N-body simulation on final accretion stage of terrestrial planets, including the effect of damping of eccentricity and inclination caused by tidal interaction with a remnant gas disk. As a result of runway and oligarchic accretion, about 20 Mars-sized protoplanets would be formed in nearly circular orbits with orbital separation of several to ten Hill radius. The orbits of the protoplanets would be eventually destabilized by long-term mutual gravity and/or secular resonance of giant gaseous planets. The protoplanets would coalesce with each other to form terrestrial planets through the orbital crossing. Previous N-body simulations, however, showed that the final eccentricities of planets are around 0.1, which are about 10 times higher than the present eccentricities of Earth and Venus. The obtained high eccentricities are the remnant of orbital crossing. We included the effect of eccentricity damping caused by gravitational interaction with disk gas as a drag force (“gravitational drag”) and carried out N-body simulation of accretion of protoplanets. We start with 15 protoplanets with 0.2M⊕ and integrate the orbits for 10⁷ years, which is consistent with the observationally inferred disk lifetime (in some runs, we start with 30 protoplanets with 0.1M⊕). In most runs, the damping time scale, which is equivalent to the strength of the drag force, is kept constant throughout each run in order to clarify the effects of the damping. We found that the planets' final mass, spatial distribution, and eccentricities depend on the damping time scale. If the damping time scale for a 0.2M⊕ mass planet at 1 AU is longer than 10⁸ years, planets grow to Earth's size, but the final eccentricities are too high as in gas-free cases. If it is shorter than 10⁶ years, the eccentricities of the protoplanets cannot be pumped up, resulting in not enough orbital crossing to make Earth-sized planets. Small planets with low eccentricities are formed with small orbital separation. On the other hand, if it is between 10⁶ and 10⁸ years, which may correspond to a mostly depleted disk (0.01-0.1% of surface density of the minimum mass model), some protoplanets can grow to about the size of Earth and Venus, and the eccentricities of such surviving planets can be diminished within the disk lifetime. Furthermore, in innermost and outermost regions in the same system, we often find planets with smaller size and larger eccentricities too, which could be analogous to Mars and Mercury. This is partly because the gravitational drag is less effective for smaller mass planets, and partly due to the “edge effect,” which means the innermost and outermost planets tend to remain without collision. We also carried out several runs with time-dependent drag force according to depletion of a gas disk. In these runs, we used exponential decay model with e-folding time of 3×10⁶ years. The orbits of protoplanets are stablized by the eccentricity damping in the early time. When disk surface density decays to ?1% of the minimum mass disk model, the damping force is no longer strong enough to inhibit the increase of the eccentricity by distant perturbations among protoplanets so that the orbital crossing starts. In this disk decay model, a gas disk with 10⁻⁴-10⁻³ times the minimum mass model still remains after the orbital crossing and accretional events, which is enough to damp the eccentricities of the Earth-sized planets to the order of 0.01. Using these results, we discuss a possible scenario for the last stage of terrestrial planet formation.     

On the variations in the insolation at Mercury resulting from oscillations of the orbital eccentricity

This paper describes variations in the insolation on Mercury resulting from fluctuations of the orbital eccentricity (0.11≤e≤0.24) of the planet. Equations for the instantaneous and the daily insolation are briefly discussed and several numerical examples are given illustrating the sensitivity of the solar radiation to changes ine. Special attention is paid to the behavior of the solar radiation distribution curves near sunrise and sunset which at the warm pole of Mercury (longitudes ±90°) occur as the planet goes through perihelion. It has been found that for eccentricities larger than about 0.194 there exists two permanent thermal bulges on opposite sides of the Mercurian surface that alternately point to the Sun at every perihelion passage. The critical value ofe past which the Sun shortly sets after perihelion is near 0.213.