Planetans—oceanic planets

The development of principles, systems, and instruments enable the detection of exoplanets with 6–8 Earth masses or less. The launches of specialized satellites, such as CoRoT (2006) and Kepler (2009), into orbits around the Earth have enabled the discovery of new exoplanetary systems. These missions are searching for relatively low-mass planets by observing their transits over the disks of their parent stars. At the same time, supporting studies of exoplanets using ground-based facilities (that measure Keplerian components of radial velocities) are in progress. The properties of at least two objects discovered by different methods, Kepler-22 and GJ 1214b, suggested that there was another class of celestial bodies among the known types of extrasolar planets: planetans, or oceanic planets. The structure of Kepler-22 and GJ 1214b suggest that they can be these oceanic planets. In this paper, we consider to what extent this statement is valid. The consideration of exoplanet Gl 581g as an oceanic planet is more feasible. Some specific features of the physical nature of these unusual planets are presented.     

The road to Earth twins – Karl Schwarzschild Award Lecture 2010

A rich population of low‐mass planets orbiting solar‐type stars on tight orbits has been detected by Doppler spectroscopy. These planets have masses in the domain of super‐Earths and Neptune‐type objects, and periods less than 100 days. In numerous cases these planets are part of very compact multiplanetary systems. Up to seven planets have been discovered orbiting one single star. These low‐mass planets have been detected by the HARPS spectrograph around 30 % of solar‐type stars. This very high occurrence rate has been recently confirmed by the results of the Kepler planetary transit space mission. The large number of planets of this kind allows us to attempt a first characterization of their statistical properties, which in turn represent constraints to understand the formation process of these systems. The achieved progress in the sensitivity and stability of spectrographs have already led to the discovery of planets with masses as small as 1.5 M_⊕ (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

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.     

Weighing the cusp at the Galactic Centre

As stars close to the galactic centre have short orbital periods it has been possible to trace large fractions of their orbits in the recent years. Previously the data of the orbit of the star S2 have been fitted with Keplerian orbits corresponding to a massive black hole (MBH) with a mass of M_BH = 3–4 × 10⁶M_⊙ implying an insignificant cusp mass. However, it has also been shown that the central black hole resides in a ∼1″ diameter stellar cluster of a priori unknown mass. In a spherical potential which is neither Keplerian nor harmonic, orbits will precess resulting in inclined rosetta shaped trajectories on the sky. In this case, the assumption of non‐Keplerian orbits is a more physical approach. It is also the only approach through which cusp mass information can be obtained via stellar dynamics of the cusp members. This paper presents the first exemplary modelling efforts in this direction. Using positional and radial data of star S2, we find that there could exist an unobserved extended mass component of several 10⁵M_⊙ forming a so‐called ‘cusp’ centered on the black hole position. Considering only the fraction of the cusp mass Mequation/tex2gif-inf-4.gif within the apo‐center of the S2 orbit we find as an upper limit that Mequation/tex2gif-inf-6.gif/(M_BH + Mequation/tex2gif-inf-9.gif) ≤ 0.05. A large extended cusp mass, if present, is unlikely to be composed of sub‐solar mass constituents, but could be explained rather well by a cluster of high M/L stellar remnants, which we find to form a stable configuration. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Orbital structure of the GJ876 extrasolar planetary system based on the latest Keck and HARPS radial velocity data

We use full available array of radial velocity data, including recently published HARPS and Keck observatory sets, to characterize the orbital configuration of the planetary system orbiting GJ876. First, we propose and describe in detail a fast method to fit perturbed orbital configuration, based on the integration of the sensitivity equations inferred by the equations of the original N-body problem. Further, we find that it is unsatisfactory to treat the available radial velocity data for GJ876 in the traditional white noise model, because the actual noise appears autocorrelated (and demonstrates non-white frequency spectrum). The time scale of this correlation is about a few days, and the contribution of the correlated noise is about 2 m/s (i.e., similar to the level of internal errors in the Keck data). We propose a variation of the maximum-likelihood algorithm to estimate the orbital configuration of the system, taking into account the red noise effects. We show, in particular, that the non-zero orbital eccentricity of the innermost planet d, obtained in previous studies, is likely a result of misinterpreted red noise in the data. In addition to offsets in some orbital parameters, the red noise also makes the fit uncertainties systematically underestimated (while they are treated in the traditional white noise model). Also, we show that the orbital eccentricity of the outermost planet is actually ill-determined, although bounded by ~0.2. Finally, we investigate possible orbital non-coplanarity of the system, and limit the mutual inclination between the planets b and c orbits by 5°?C15°, depending on the angular position of the mutual orbital nodes.     

A model for accretion of the terrestrial planets

One possible origin of the terrestrial planets involves their formation by gravitational accretion of particles originally in Keplerian orbits about the sun. Some implications of this theory are considered. A formal expression for the rate of mass accretion by a planet is developed. The formal singularity of the gravitational collision cross-section for low relative velocities is shown to be without physical significance when the accreting bodies are in heliocentric orbits. The distribution of particle velocities relative to an accreting planet is considered; the mean velocity increases with time. The internal temperature of an accreting planet is shown to depend simply on the accretion rate. A simple and physically reasonable approximate expression for a planetary accretion rate is proposed.     

Cosmogony of the solar system

In Sections 1–6, we determine an approximate analytical model for the density and temperature distribution in the protoplanetary could. The rotation of the planets is discussed in Section 7 and we conclude that it cannot be determined from simple energy conservation laws.The velocity of the gas of the protoplanetary cloud is found to be smaller by about 5×10³ cm s^–1 in comparison to the Keplerian circular velocity. If the radius of the planetesimals is smaller than a certain limitr ₁, they move together with the gas. Their vertical and horizontal motion for this case is studied in Sections 8 and 9.As the planetesimals grow by accretion their radius becomes larger thanr ₁ and they move in Keplerian orbits. As long as their radius is betweenr ₁ and a certain limitr ₂ their gravitational interaction is negligible. In Section 10, we study the accretion for this case.In Section 11, we determine the change of the relative velocities due to close gravitational encounters. The principal equations governing the late stages of accretion are deduced in Section 12, In Section 13 there are obtained approximate analytical solutions.The effect of gas drag and of collisions is studied in Sections 14 and 15, respectively. Numerical results and conclusions concerning the last and principal stage of accretion are drawn in Section 16.     

Dynamical history of coplanar two-satellite systems

One of the possible early states of the Earth-Moon system was a system of several large satellites around the Earth. The dynamical evolution of coplanar three-body systems is studied; a planet (Earth) and two massive satellites (proto-moons) with geocentric orbits of slightly different radii. Such configurations may arise in multiple satellite systems receding from a planet due to tidal friction. The numerical integration of the equations of motion shows that initially circular Keplerian orbits are soon transformed into disturbed elliptic orbits which are intersecting. The life-time of such a coplanar system between two probable physical collisions of satellites is roughly from one day to one year for satellite systems with radii less than 20R⊕, and may reach 100 yr for three-dimensional systems. This time-scale is short in comparison with the duration of the removal of satellites due to tides raised on the planet, which is estimated as 10⁶–10⁸ yr for the same orbital dimensions. Therefore, the life-time of a system of several proto-moons is mainly determined by their tidal interactions with the Earth. For conditions which we have considered, the most probable result of the evolution was coalescence of satellites as the consequence of the collisions.     

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 habitability of the OGLE-2006-BLG-109L planetary system

We investigate the dynamics of putative Earth-mass planets in the habitable zone (HZ) of the extrasolar planetary system OGLE-2006-BLG-109L, a close analogue of the Solar system. Our work is inspired by the work of Malhotra & Minton. Using the linear Laplace–Lagrange theory, they identified a strong secular resonance that may excite large eccentricity of orbits in the HZ. However, due to uncertain or unconstrained orbital parameters, the subsystem of Jupiters may be found in a dynamically active region of the phase space spanned by low-order mean-motion resonances. To generalize this secular model, we construct a semi-analytical averaging method in terms of the restricted problem. The secular orbits of large planets are approximated by numerically averaged osculating elements. They are used to calculate the mean orbits of terrestrial planets by means of a high-order analytic secular theory developed in our previous works. We found regions in the parameter space of the problem in which stable, quasi-circular orbits in the HZ are permitted. The excitation of eccentricity in the HZ strongly depends on the apsidal angle of jovian orbits. For some combinations of that angle, eccentricities and semimajor axes consistent with the observations, a terrestrial planet may survive in low eccentric orbits. We also study the effect of post-Newtonian gravity correction on the innermost secular resonance.     

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.     

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.     

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.     

HD 12661行星系统构形的形成与稳定性

最近的多普勒观测表明恒星HD 12661周围存在两颗中等偏心率轨道上运行的行星,内行星的最小质量为2．3木星质量,轨道周期为263．6天;外行星的最小质量为1．57木星质量,轨道周期为1444．5天．该系统的稳定性要求两颗行星处在平运动轨道共振．用数值方法研究了该系统形成初期在恒星气体盘作用下的轨道迁移与稳定性,计算了行星在迁移中被平运动共振俘获的概率．发现这两颗行星目前很可能正处在11:2平运动共振边缘,且运动是混沌的,从而澄清了关于系统目前构形的不同说法,并且很可能在系统形成后行星迁移到目前构形时,气体盘几乎消失了．     

The dynamics of two massive planets on inclined orbits

The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.     

Generalized Averaging Principle and the Secular Evolution of Planet Crossing Orbits

Planet crossing orbits give rise to mathematical singularities that make it not possible to apply the classical averaging principle to study the qualitative evolution of Near Earth Asteroids (NEAs). Recently this principle has been generalized to deal with crossings in a mathematical model with the planets on circular coplanar orbits. More accuracy is needed to compute the averaged evolution of planet crossing orbits for different purposes: computing reliable crossing times for the averaged motion, writing more precise proper elements and frequencies for NEAs, etc. In this paper we present the generalization of the averaging principle using a model where the eccentricity and the inclination of the planets are taken into account.     

Young Jupiters are faint: new models of the early evolution of giant planets

Here we show preliminary calculations of the cooling and contraction of a 2 M_J planet. These calculations, which are being extended to 1–10 M_J, differ from other published “cooling tracks” in that they include a core accretion‐gas capture formation scenario, the leading theory for the formation of gas giant planets.We find that the initial post‐accretionary intrinsic luminosity of the planet is ∼3 times less than previously published models which use arbitrary initial conditions. These differences last a few tens of millions of years. Young giant planets are intrinsically fainter than has been previously appreciated. We also discuss how uncertainties in atmospheric chemistry and the duration of the formation time of giant planets lead to challenges in deriving planetary physical properties from comparison with tabulated model values. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Tidal Evolution of Low-mass Kepler Planets

The planets with a radius &lt; 4 R_⊕ observed by the Kepler mission exhibit a unique feature, and propose a challenge for current planetary formation models. The tidal effect between a planet and its host star plays an essential role in reconfiguring the final orbits of the short-period planets. In this work, based on various initial Rayleigh distributions of the orbital elements, the final semi-major axis distributions of the planets with a radius &lt; 4 R_⊕ after suffering tidal evolutions are investigated. Our simulations have qualitatively revealed some statistical properties: the semi-major axis and its peak value all increase with the increase of the initial semi-major axis and eccentricity. For the case that the initial mean semi-major axis is less than 0.1 au and the mean eccentricity is larger than 0.25, the results of numerical simulation are approximately consistent with the observation. In addition, the effects of other parameters, such as the tidal dissipation coefficient, stellar mass and planetary mass, etc., on the final semi-major axis distribution after tidal evolution are all relatively small. Based on the simulation results, we have tried to find some clues for the formation mechanism of low-mass planets. We speculate that these low-mass planets probably form in the far place of protoplanetary disk with a moderate eccentricity via the type I migration, and it is also possible to form in situ.     

Stability of inclined orbits of terrestrial planets in habitable zones

In long-term stability studies of terrestrial planets moving in the habitable zone (HZ) of a sun-like star, we distinguish four different configurations: (i) planets moving in binary star systems, (ii) the inner type (where the gas giant moves outside the HZ), (iii) the outer type (where the gas giant is closer to the star, than the HZ) and (iv) the Trojan type (where the gas giant moves in the HZ). Since earlier calculations indicated, that the stability of the motion in the HZ also depends on the inclination of the terrestrial planet orbits, we present a detailed numerical investigation to show correlations between the eccentricity, the mass and the distance of the giant planet for various inclinations of the terrestrial planets. The orbital stability of the HZ was examined for all four configurations stated above. While we could find hardly any stable orbits for the first three types for inclinations higher than 40^°, the Trojan planets can be stable up to an inclination of 60^°. Additionally, we could also find some stabilizing effects of the inclination for the first three types. As dynamical model we used the elliptic restricted three-body problem, which consists of two massive and one mass-less body. This allows an application to all detected and future extrasolar single planet systems.