Planet mass stability margin of two-planetary systems

The Lagrange stability of the Sun-Jupiter-Saturn and 47 UMa two-planetary systems at a time scale of 10⁶ yr was studied using the method of averaging. When the masses of Jupiter and Saturn increase by 19 times, these planets can closely converge. The study of Lagrange stability in the case of successive mass increase allows for the obtainment of upper estimates of possible masses of extrasolar planets. Close approachs in the 47 UMa system are possible when minimal masses increase by 38 times. approachs are revealed when analyzing osculating elements; in averaged elements, approachs are absent. Resonant properties of six extrasolar two-planetary systems where the outer planet is less massive than the inner one are studied. The values of semi-major axes of planet orbits in HD 82943 and HD 73526 systems lie in a narrow resonant zone; in 47 UMa, μ Ara and HD 108874 systems lie in a wide resonant zone. In the HD 12661, the system resonances of a lower order were not revealed.     

On the possibility of earth-type habitable planets around 47 UMa

We investigate whether Earth-type habitable planets can in principle exist in the planetary system of 47 UMa. The system of 47 UMa consists of two Jupiter-size planets beyond the outer edge of the stellar habitable zone, and thus resembles our own Solar System most closely compared to all exosolar planetary systems discovered so far. Our study of habitability deliberately follows an Earth-based view according to the concept of Franck and colleagues, which assumes the long-term possibility of photosynthetic biomass production under geodynamic conditions. Consequently, a broad variety of climatological, biogeochemical, and geodynamical processes involved in the generation of photosynthesis-driven life conditions is taken into account. The stellar luminosity and the age of the star/planet system are of fundamental importance for planetary habitability. Our study considers different types of planetary continental growth models and takes into account a careful assessment of the stellar parameters. In the event of successful formation and orbital stability, two subjects of intense research, we find that Earth-type habitable planets around 47 UMa are in principle possible! The likelihood of those planets is increased if assumed that 47 UMa is relatively young (?6 Gyr) and has a relatively small stellar luminosity as permitted by the observational range of those parameters.     

A planetary system with an escaping Mars

The chaotic behaviour of the motion of the planets in our Solar System is well established. In this work to model a hypothetical extrasolar planetary system our Solar System was modified in such a way that we replaced the Earth by a more massive planet and let the other planets and all the orbital elements unchanged. The major result of former numerical experiments with a modified Solar System was the appearance of a chaotic window at κ _E ∈ (4, 6), where the dynamical state of the system was highly chaotic and even the body with the smallest mass escaped in some cases. On the contrary for very large values of the mass of the Earth, even greater than that of Jupiter regular dynamical behaviour was observed. In this paper the investigations are extended to the complete Solar System and showed, that this chaotic window does still exist. Tests in different ‘Solar Systems’ clarified that including only Jupiter and Saturn with their actual masses together with a more ‘massive’ Earth (4 < κ _E < 6) perturbs the orbit of Mars so that it can even be ejected from the system. Using the results of the Laplace‐Lagrange secular theory we found secular resonances acting between the motions of the nodes of Mars, Jupiter and Saturn. These secular resonances give rise to strong chaos, which is the cause of the appearance of the instability window. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hydrogen Eos at Megabar Pressures and the Search for Jupiter’s Core

The interior structure of Jupiter serves as a benchmark for an entire astrophysical class of liquid–metallic hydrogen-rich objects with masses ranging from ~0.1M _J to ~80M _J (1M _J = Jupiter mass = 1.9e30 g), comprising hydrogen-rich giant planets (mass < 13M _J) and brown dwarfs (mass > 13M _J but ~ < 80M _J), the so-called substellar objects (SSOs). Formation of giant planets may involve nucleated collapse of nebular gas onto a solid, dense core of mass ~0.04M _J rather than a stellar-like gravitational instability. Thus, detection of a primordial core in Jupiter is a prime objective for understanding the mode of origin of extrasolar giant planets and other SSOs. A basic method for core detection makes use of direct modeling of Jupiter’s external gravitational potential terms in response to rotational and tidal perturbations, and is highly sensitive to the thermodynamics of hydrogen at multi-megabar pressures. The present-day core masses of Jupiter and Saturn may be larger than their primordial core masses due to sedimentation of elements heavier than hydrogen. We show that there is a significant contribution of such sedimented mass to Saturn’s core mass. The sedimentation contribution to Jupiter’s core mass will be smaller and could be zero.     

Core formation in giant gaseous protoplanets

Sedimentation rates of silicate grains in gas giant protoplanets formed by disk instability are calculated for protoplanetary masses between 1 MSaturn to 10 MJupiter. Giant protoplanets with masses of 5 MJupiter or larger are found to be too hot for grain sedimentation to form a silicate core. Smaller protoplanets are cold enough to allow grain settling and core formation. Grain sedimentation and core formation occur in the low mass protoplanets because of their slow contraction rate and low internal temperature. It is predicted that massive giant planets will not have cores, while smaller planets will have small rocky cores whose masses depend on the planetary mass, the amount of solids within the body, and the disk environment. The protoplanets are found to be too hot to allow the existence of icy grains, and therefore the cores are predicted not to contain any ices. It is suggested that the atmospheres of low mass giant planets are depleted in refractory elements compared with the atmospheres of more massive planets. These predictions provide a test of the disk instability model of gas giant planet formation. The core masses of Jupiter and Saturn were found to be ∼0.25 M_⊕ and ∼0.5 M_⊕, respectively. The core masses of Jupiter and Saturn can be substantially larger if planetesimal accretion is included. The final core mass will depend on planetesimal size, the time at which planetesimals are formed, and the size distribution of the material added to the protoplanet. Jupiter's core mass can vary from 2 to 12 M_⊕. Saturn's core mass is found to be ∼8 M_⊕.     

Planets in the asteroid belt

Abstract— The main asteroid belt has lost >99.9% of its solid mass since the time at which the planets were forming, according to models for the protoplanetary nebula. Here we show that the primordial asteroid belt could have been cleared efficiently if much of the original mass accreted to form planetsized bodies, which were capable of perturbing one another into unstable orbits. We provide results from 25 N‐body integrations of up to 200 planets in the asteroid belt, with individual masses in the range 0.017–0.33 Earth masses. In the simulations, these bodies undergo repeated close encounters which scatter one another into unstable resonances with the giant planets, leading to collision with the Sun or ejection from the solar system. In response, the giant planets' orbits migrate radially and become more circular. This reduces the size of the main‐belt resonances and the clearing rate, although clearing continues. If ～3 Earth masses of material was removed from the belt this way, Jupiter and Saturn would initially have had orbital eccentricities almost twice their current values. Such orbits would have made Jupiter and Saturn 10–100x more effective at clearing material from the belt than they are on their current orbits. The time required to remove 90% of the initial mass from the belt depends sensitively on the giant planets' orbits, and weakly on the masses of the asteroidal planets. 18 of the 25 simulations end with no planets left in the belt, and the clearing takes up to several hundred million years. Typically, the last one or two asteroidal planets are removed by interactions with planets in the terrestrial region     

A comparison of the interiors of Jupiter and Saturn

Interior models of Jupiter and Saturn are calculated and compared in the framework of the three-layer assumption, which rely on the perception that both planets consist of three globally homogeneous regions: a dense core, a metallic hydrogen envelope, and a molecular hydrogen envelope. Within this framework, constraints on the core mass and abundance of heavy elements (i.e. elements other than hydrogen and helium) are given by accounting for uncertainties on the measured gravitational moments, surface temperature, surface helium abundance, and on the inferred protosolar helium abundance, equations of state, temperature profile and solid/differential interior rotation. Results obtained solely from static models matching the measured gravitational fields indicate that the mass of Jupiter’s dense core is less than 14 M_⊕ (Earth masses), but that models with no core are possible given the current uncertainties on the hydrogen–helium equation of state. Similarly, Saturn’s core mass is less than 22 M_⊕ but no lower limit can be inferred. The total mass of heavy elements (including that in the core) is constrained to lie between 11 and 42 M_⊕ in Jupiter, and between 19 and 31 M_⊕ in Saturn. The enrichment in heavy elements of their molecular envelopes is 1–6.5, and 0.5–12 times the solar value, respectively. Additional constraints from evolution models accounting for the progressive differentiation of helium (Hubbard WB, Guillot T, Marley MS, Burrows A, Lunine JI, Saumon D, 1999. Comparative evolution of Jupiter and Saturn. Planet. Space Sci. 47, 1175–1182) are used to obtain tighter, albeit less robust, constraints. The resulting core masses are then expected to be in the range 0–10 M_⊕, and 6–17 M_⊕ for Jupiter and Saturn, respectively. Furthermore, it is shown that Saturn’s atmospheric helium mass mixing ratio, as derived from Voyager, Y=0.06±0.05, is probably too low. Static and evolution models favor a value of Y=0.11−0.25. Using, Y=0.16±0.05, Saturn’s molecular region is found to be enriched in heavy elements by 3.5 to 10 times the solar value, in relatively good agreement with the measured methane abundance. Finally, in all cases, the gravitational moment J₆ of models matching all the constraints are found to lie between 0.35 and 0.38×10⁻⁴ for Jupiter, and between 0.90 and 0.98×10⁻⁴ for Saturn, assuming solid rotation. For comparison, the uncertainties on the measured J₆ are about 10 times larger. More accurate measurements of J₆ (as expected from the Cassini orbiter for Saturn) will therefore permit to test the validity of interior models calculations and the magnitude of differential rotation in the planetary interior.     

Habitable zones of host stars during the post-MS phase

A star will become brighter and brighter with stellar evolution, and the distance of its habitable zone will become larger and larger. Some planets outside the habitable zone of a host star during the main sequence phase may enter the habitable zone of the host star during other evolutionary phases. A terrestrial planet within the habitable zone of its host star is generally thought to be suitable for the existence of life. Furthermore, a rocky moon around a giant planet may be also suitable for life to survive, provided that the planet–moon system is within the habitable zone of its host star. Using Eggleton’s code and the boundary flux of the habitable zone, we calculate the habitable zone of our Solar system after the main sequence phase. It is found that Mars’ orbit and Jupiter’s orbit will enter the habitable zone of the Solar system during the subgiant branch phase and the red giant branch phase, respectively. And the orbit of Saturn will enter the habitable zone of Solar during the He-burning phase for about 137 million years. Life is unlikely at any time on Saturn, as it is a giant gaseous planet. However, Titan, the rocky moon of Saturn, may be suitable for biological evolution and become another Earth during that time. For low-mass stars, there are similar habitable zones during the He-burning phase as our Solar, because there are similar core masses and luminosities for these stars during that phase.     

Computer studies of the evolution of planetary and satellite systems

This paper describes two computer experiments carried out with a CDC-Cyber 74 program (Barricelliet al., 1979) for computer simulation of a large number of objects in orbit about a central body or primary. The first experiment was started with 125 planets of which the two largest ones had coplanar orbits and masses comparable to those of Jupiter and Saturn, respectively. Their semimajor axes and eccentricities were, however, much larger. The smaller planets had a distribution promoting the formation of an axial meeting area. The experiment gives information relevant to the question of focusing of planetary orbits into a common plane and to the question of the formation and stability of an axial meeting area. Together with the next experiment, it also gives information about the development of commensurabilities (or resonances) with the largest planets.The second experiment started with 55 planets none of them with a mass greater than about 20% of Jupiter's but several of them with orbits close to a common plane. The aim of the experiment was to investigate whether successive captures followed by planetary fusion (Barricelli, 1972a) could lead to the formation of major planets comparable to Jupiter and Saturn, and in similar orbits. Also this experiment gives information relevant to the commensurability problem.     

Models and theoretical spectrum for the free oscillations of Saturn

We constructed new models of Saturn with an allowance made for a helium mass fraction of ～0.18–0.25 in its atmosphere. Our modeling shows that the composition of Saturn differs markedly from the solar composition; more specifically, during its formation, the planet was ～11–15 planetary masses short of the hydrogen-helium component. Saturn, as well as the other giant planets, must have been formed according to Schmidt’s scenario, through the formation of embryonic nuclei, rather than according to Laplace’s scenario. The masses of the embryonic nuclei themselves lie within the range (3.5–8) M_⊕. We calculated a theoretical free-oscillation spectrum for our models of Saturn, each of which fits all of the available observational data. The results of our calculations are presented graphically and in tables. Of particular interest are the diagnostic potentialities of the discontinuity gravitational modes related to density jumps in the molecular envelope of Saturn and at the interface between its molecular and metallic envelopes. When observational data appear, our results can be used both to identify the observed modes and to improve the models themselves. We discuss some of the cosmogonical aspects associated with the formation of the giant planets.     

On the Stability of the Terrestrial Planets as Models for Exosolar Planetary Systems

All results, achieved up to now, show the long term stability of our planetary system, although, especially the inner solar system is chaotic, due to some specific secular resonances. We study, by means of numerical integrations, the dynamical evolution of the planetary system where we concentrate on the stability of motion of the terrestrial planets Venus, Earth and Mars. Our model consists of a simplified planetary system with the inner planets Venus, Earth and Mars as well as Jupiter and Saturn. A mass factor was introduced to uniformly change the masses of the terrestrial planets; Jupiter and Saturn were involved in the system with their actual masses. We integrated the equations of motion with a Lie-integration method for a time interval of 10⁷ years. It turned out that when 220 < < 245 and > 250 the system became unstable due to the strong interactions between the planets. We discuss the model planetary systems for small mass-factors 0.5 10 and large ones 160 270 with the aid of several different numerical tools. These results can be applied to recently discovered exoplanetary systems, which configuration is comparable to our own.     

On the stability of Earth-like planets in multi-planet systems

We present a continuation of our numerical study on planetary systems with similar characteristics to the Solar System. This time we examine the influence of three giant planets on the motion of terrestrial-like planets in the habitable zone (HZ). Using the Jupiter–Saturn–Uranus configuration we create similar fictitious systems by varying Saturn’s semi-major axis from 8 to 11 AU and increasing its mass by factors of 2–30. The analysis of the different systems shows the following interesting results: (i) Using the masses of the Solar System for the three giant planets, our study indicates a maximum eccentricity (max-e) of nearly 0.3 for a test-planet placed at the position of Venus. Such a high eccentricity was already found in our previous study of Jupiter–Saturn systems. Perturbations associated with the secular frequency g ₅ are again responsible for this high eccentricity. (ii) An increase of the Saturn-mass causes stronger perturbations around the position of the Earth and in the outer HZ. The latter is certainly due to gravitational interaction between Saturn and Uranus. (iii) The Saturn-mass increased by a factor 5 or higher indicates high eccentricities for a test-planet placed at the position of Mars. So that a crossing of the Earth’ orbit might occur in some cases. Furthermore, we present the maximum eccentricity of a test-planet placed in the Earth’ orbit for all positions (from 8 to 11 AU) and masses (increased up to a factor of 30) of Saturn. It can be seen that already a double-mass Saturn moving in its actual orbit causes an increase of the eccentricity up to 0.2 of a test-planet placed at Earth’s position. A more massive Saturn orbiting the Sun outside the 5:2 mean motion resonance (a _S  ≥9.7 AU) increases the eccentricity of a test-planet up to 0.4.     

Formation of terrestrial planets in disks with different surface density profiles

We present the results of an extensive study of the final stage of terrestrial planet formation in disks with different surface density profiles and for different orbital configurations of Jupiter and Saturn. We carried out simulations in the context of the classical model with disk surface densities proportional to \({r^{-0.5}}, {r^{-1}}\) and \({r^{-1.5}}\), and also using partially depleted, non-uniform disks as in the recent model of Mars formation by Izidoro et al. (Astrophys J 782:31, 2014). The purpose of our study is to determine how the final assembly of planets and their physical properties are affected by the total mass of the disk and its radial profile. Because as a result of the interactions of giant planets with the protoplanetary disk, secular resonances will also play important roles in the orbital assembly and properties of the final terrestrial planets, we will study the effect of these resonances as well. In that respect, we divide this study into two parts. When using a partially depleted disk (Part 1), we are particularly interested in examining the effect of secular resonances on the formation of Mars and orbital stability of terrestrial planets. When using the disk in the classical model (Part 2), our goal is to determine trends that may exist between the disk surface density profile and the final properties of terrestrial planets. In the context of the depleted disk model, results of our study show that in general, the \(\nu _5\) resonance does not have a significant effect on the dynamics of planetesimals and planetary embryos, and the final orbits of terrestrial planets. However, \(\nu _6\) and \(\nu _{16}\) resonances play important roles in clearing their affecting areas. While these resonances do not alter the orbits of Mars and other terrestrial planets, they strongly deplete the region of the asteroid belt ensuring that no additional mass will be scattered into the accretion zone of Mars so that it can maintain its mass and orbital stability. In the context of the classical model, the effects of these resonances are stronger in disks with less steep surface density profiles. Our results indicate that when considering the classical model (Part 2), the final planetary systems do not seem to show a trend between the disk surface density profile and the mean number of the final planets, their masses, time of formation, and distances to the central star. Some small correlations were observed where, for instance, in disks with steeper surface density profiles, the final planets were drier, or their water contents decreased when Saturn was added to the simulations. However, in general, the final orbital and physical properties of terrestrial planets seem to vary from one system to another and depend on the mass of the disk, the spatial distribution of protoplanetary bodies (i.e., disk surface density profile), and the initial orbital configuration of giant planets. We present results of our simulations and discuss their implications for the formation of Mars and other terrestrial planets, as well as the physical properties of these objects such as their masses and water contents.     

Formation of the giant planets

Observational constraints on interior models of the giant planets indicate that these planets were all much hotter when they formed and they all have rock and/or ice cores of ten to thirty earth masses. These cores are probably soluble in the envelopes above, especially in Jupiter and Saturn, and are therefore likely to be primordial. They persist despite the continual upward mixing by thermally driven convection throughout the age of the solar system, because of the inefficiency of double-diffusive convection. Thus, these planets most probably formed by the hydrodynamic collapse of a gaseous envelope onto a core rather than by direct instability of the gaseous solar nebula. Recent calculations by Mizuno (1980, Prog. Theor. Phys.64, 544) show that this formation mechanism may explain the similarity of giant planet core masses. Problems remain however, and no current model is entirely satisfactory in explaining the properties of the giant planets and simultaneously satisfying the terrestrial planet constraints. Satellite systematics and protoplanetary disk nebulae are also discussed and related to formation conditions.     

Models of Jupiter and Saturn after Galileo mission

New models of Jupiter are based on observational data provided by the Galileo spaceprobe, which considerably improved previously existing estimates of the helium abundance in the atmosphere of Jupiter. These data yield for Jupiter’s atmosphere 20% of the solar oxygen abundance and do not agree with the results of the analysis of the collision of comet Shoemaker-Levy 9 with Jupiter (10 times the solar value). Therefore, both the models of Jupiter with water-depleted and water-enriched atmosphere are considered. By analogy with Jupiter, trial models of Saturn with a water-depleted external envelope are also developed. The molecular-metallic phase transition pressure of hydrogen P_m was taken to be 1.5, 2 and 3 Mbar. Since Saturn’s internal molecular envelope is noticeably enriched in the IR-component (its weight concentration, 0.25–0.30, being by a factor of 3–4 higher than in Jupiter), the phase transition pressure in Saturn can be lower than in Jupiter. In the constructed models, the IR-core masses are 3–3.5 M_⊕ for Jupiter and 3–5.5 M_⊕ for Saturn. Jupiter’s and Saturn’s IR-cores can be considered embryos onto which the accretion of the gas occurred during the formation of the planets. The mass of the hydrogen–helium component dispersed in the zone of planetary formation constitutes ≈2–5 planetary masses for Jupiter and ≈11–14 planetary masses for Saturn.     

Some dynamical aspects of the accretion of Uranus and Neptune: The exchange of orbital angular momentum with planetesimals

The final stage of the accretion of Uranus and Neptune is numerically investigated. The four Jovian planets are considered with Jupiter and Saturn assumed to have reached their present sizes, whereas Uranus and Neptune are taken with initial masses 0.2 of their present ones. Allowance is made for the orbital variation of the Jovian planets due to the exchange of angular momentum with interacting bodies (“planetesimals”). Two possible effects that may have contributed to the accretion of Uranus and Neptune are incorporated in our model: (1) an enlarged cross section for accretion of incoming planetesimals due to the presence of extended gaseous envelopes and/or circumplanetary swarms of bodies; and (2) intermediate protoplanets in mid-range orbits between the Jovian planets. Significant radial displacements are found for Uranus and Neptune during their accretion and scattering of planetesimals. The orbital angular momentum budgets of Neptune, Uranus, and Saturn turn out to be positive; i.e., they on average gain orbital angular momentum in their interactions with planetesimals and hence they are displaced outwardly. Instead, Jupiter as the main ejector of bodies loses orbital angular momentum so it moves sunward. The gravitational stirring of planetesimals caused by the introduction of intermediate protoplanets has the effect that additional solid matter is injected into the accretion zones of Uranus and Neptune. For moderate enlargements of the radius of the accretion cross section (2–4 times), the accretion time scale of Uranus and Neptune are found to be a few 10⁸ years and the initial amount of solid material required to form them of a few times their present masses. Given the crucial role played by the size of the accretion cross section, questions as to when Uranus and Neptune acquired their gaseous envelopes, when the envelopes collapsed onto the solid cores, and how massive they were are essential in defining the efficiency and time scale of accretion of the two outer Jovian planets.     

The influence of the great inequality on the secular disturbing function of the planetary system

This paper derives the contributionF ₂ ^* by the great inequality to the secular disturbing function of the principal planets. Andoyer's expansion of the planetary disturbing function and von Zeipel's method of eliminating the periodic terms is employed; thereby, the corrected secular disturbing function for the planetary system is derived. An earlier solution suggested by Hill is based on Leverrier's equations for the variation of elements of Jupiter and Saturn and on the semi-empirical adjustment of the coefficients in the secular disturbing function. Nowadays there are several modern methods of eliminating periodic terms from the Hamiltonian and deriving a purely secular disturbing function. Von Zeipel's method is especially suitable. The conclusion is drawn that the canonicity of the equations for the secular variation of the heliocentric elements can be preserved if there be retained, in the secular disturbing function, terms only of the second and fourth order relative to the eccentricity and inclinations.The Krylov-Bogolubov method is suggested for eliminating periodic terms, if it is desired to include the secular perturbations of the fifth and higher order in the heliocentric elements. The additional part of the secular disturbing functionF ₂ ^* derived in this paper can be included in existing theories of the secular effects of principal planets. A better approach would be to preserve the homogeneity of the theory and rederive all the secular perturbations of principal planets using Andoyer's symbolism, including the part produced by the great inequality.     

Symmetric horseshoe periodic orbits in the general planar three-body problem

We present some families of horseshoe periodic orbits in the general planar three-body problem for the case of two equal masses. The considered system is a symmetric version of the one formed by Saturn, Janus and Epimetheus. We use a mass ratio equal to 35×10⁻⁵, corresponding to 10⁵ times the Saturn-Janus mass parameter of the restricted case; for this mass ratio the satellites have a significantly bigger influence on the planet than in the classical Saturn, Janus and Epimetheus system. To obtain periodic orbits, we search those horseshoe orbits passing through two reversible configurations. A particular kind of periodic orbits where the minor bodies follow the same path is discussed.     

First Ground-Based Photometry and Light Curve Analysis of the Recently Discovered Contact Binary Hx Uma

Photoelectric UBV light curves of the recently discovered eclipsing binary HX UMa were obtained and studied to determine the preliminary physical parameters of the system for the first time. The observations were taken at the TüBİTAK¹ – Turkish National Observatory (TUG) on three nights in April 2003. A simultaneous analysis of the light and radial velocity curves yields a typical A-type contact binary with a high degree of overcontact. The influence of the close visual companion to the total light of the system was taken into account during the analysis. ¹TüBİTAK: The Scientific and Technical Research Council of Turkey.