“Horseshoe” Structures in the Debris Disks of Planet-Hosting Binary Stars

The formation of a planetary system from the protoplanetary disk leads to destruction of the latter; however, a debris disk can remain in the form of asteroids and cometary material. The motion of planets can cause the formation of coorbital structures from the debris disk matter. Previous calculations have shown that such a ring-like structure is more stable if there is a binary star in the center of the system, as opposed to a single star. To analyze the properties of the coorbital structure, we have calculated a grid of models of binary star systems with a circumbinary planet moving in a planetesimal disk. The calculations are performed considering circular orbits of the stars and the planet; the mass and position of the planet, as well as the mass ratio of the stars, are varied. The analysis of the models shows that the width of the coorbital ring and its stability significantly depend on the initial parameters of the problem. Additionally, the empirical dependences of the width of the coorbital structure on the parameters of the system have been obtained, and the parameters of the models with the most stable coorbital structures have been determined. The results of the present study can be used for the search of planets around binary stars with debris disks.     

系外行星系统掩星周期变化研究进展

在掩星法发现的系外行星系统中,如果存在其他未知的伴星绕同一颗恒星运动,掩星行星由于受到伴星引力的影响,运动轨道将发生变化,轨道周期不再是常数,而是变化的。利用这种变化探测掩星系统中的其他行星,已成为一种新的方法。主要介绍了未知行星与掩星行星之间的引力作用引起的掩星周期变化效应,以及掩星周期变化法探测系外行星的理论和研究进展状况,最后简要讨论了几种影响掩星周期变化的其他因素:共轨行星、卫星、潮汐效应、相对论效应及恒星的引力四极矩等。     

On the dynamics of the exoplanetary system Kepler-413

This study is concerned with the stability of motion of the circumbinary exoplanet Kepler-413b. The analysis is performed within the framework of a flat restricted three-body problem. The stability diagram is plotted in the plane of initial conditions “pericentric distance—eccentricity” using mass calculations of Lyapunov exponents. According to the diagram, the Kepler-413b planet is located in a stable resonance cell, confined by the mean-motion resonances 6: 1 and 7: 1 with a central binary star, which agrees with the conclusions of Kostov et al. (2014) based on calculations of the MEGNO parameter. It is shown that the value of the critical semimajor axis acquired from the empirical formula of Holman and Wiegert (1999) almost coincides with the value obtained directly from the stability diagram; at low and moderate eccentricities of the planetary orbit, the position of the calculated boundary of the chaos zone is in close agreement with the boundary predicted by Shevchenko’s theory (2015). If the planet were in the instability zone, its characteristic Lyapunov time would be only ~1 year. In accordance with the conclusions of Kostov et al. (2014), it has been shown that the planet Kepler-413b is outside the habitability zone of the system.     

Debris disks: seeing dust, thinking of planetesimals and planets

Debris disks are optically thin, almost gas-free dusty disks observed arounda significant fraction of main-sequence stars older than about 10 Myr. Since the circumstellar dust is short-lived, the very existence of these disks is considered as evi-dence that dust-producing planetesimals are still present in mature systems, in whichplanets have formed – or failed to form – a long time ago. It is inferred that theseplanetesimals orbit their host stars at asteroid to Kuiper-belt distances and continuallysupply ...     

Kepler-9b、c的近2:1平运动共振分析

对多行星系统中行星周期比的统计发现,行星周期比在简单整数比2:1、3:2的右侧边缘处有明显聚集,而在其紧邻的左边有明显空缺.针对这一现象有各种不同的动力学解释. Kepler-9系统中已发现的3个行星中,行星b、c周期比约为2.03,是接近2:1共振的一个典型例子.利用关于偏心率的二阶哈密顿方程,针对只考虑长期作用和加入共振摄动两种情况,通过研究当前状态下系统在能量等高线图与相空间截面图中的位置,讨论了两行星可能的近共振状态.     

On the evolution of giant protoplanets forming in circumbinary discs

We present the results of hydrodynamic simulations of Jovian mass protoplanets that form in circumbinary discs. The simulations follow the orbital evolution of the binary plus protoplanet system acting under their mutual gravitational forces, and forces exerted by the viscous circumbinary disc. The evolution involves the clearing of the inner circumbinary disc initially, so that the binary plus protoplanet system orbits within a low density cavity. Continued interaction between disc and protoplanet causes inward migration of the planet towards the inner binary. Subsequent evolution can take three distinct paths: (i) the protoplanet enters the 4 : 1 mean motion resonance with the binary, but is gravitationally scattered through a close encounter with the secondary star; (ii) the protoplanet enters the 4 : 1 mean motion resonance, the resonance breaks, and the planet remains in a stable orbit just outside the resonance; (iii) when the binary has initial eccentricity   e _bin≥ 0.2  , the disc becomes eccentric, leading to a stalling of the planet migration, and the formation of a stable circumbinary planet.
These results have implications for a number of issues in the study of extrasolar planets. The ejection of protoplanets in close binary systems provides a source of 'free-floating planets', which have been discovered recently. The formation of a large, tidally truncated cavity may provide an observational signature of circumbinary planets during formation. The existence of protoplanets orbiting stably just outside a mean motion resonance (4 : 1) in the simulations indicate that such sites may harbour planets in binary star systems, and these could potentially be observed. Finally, the formation of stable circumbinary planets in eccentric binary systems indicates that circumbinary planets may not be uncommon.     

Delivery of icy planetesimals to inner planets in the Proxima Centauri planetary system

The estimates of the delivery of icy planetesimals from the feeding zone of Proxima Centauri c (with mass equal to 7m_E, m_E is the mass of the Earth) to inner planets b and d were made. They included the studies of the total mass of planetesimals in the feeding zone of planet c and the probabilities of collisions of such planetesimals with inner planets. This total mass could be about 10–15m_E. It was estimated based on studies of the ratio of the mass of planetesimals ejected into hyperbolic orbits to the mass of planetesimals collided with forming planet c. At integration of the motion of planetesimals, the gravitational influence of planets c and b and the star was taken into account. In most series of calculations, planetesimals collided with planets were excluded from integrations. Based on estimates of the mass of planetesimals ejected into hyperbolic orbits, it was concluded that during the growth of the mass of planet c the semi-major axis of its orbit could decrease by at least a factor of 1.5. Depending on possible gravitational scattering due to mutual encounters of planetesimals, the total mass of material delivered by planetesimals from the feeding zone of planet c to planet b was estimated to be between 0.002m_E and 0.015m_E. Probably, the amount of water delivered to Proxima Centauri b exceeded the mass of water in Earth's oceans. The amount of material delivered to planet d could be a little less than that delivered to planet b.     

Metallicity, debris discs and planets

We investigate the populations of main-sequence stars within 25 pc that have debris discs and/or giant planets detected by Doppler shift. The metallicity distribution of the debris sample is a very close match to that of stars in general, but differs with >99 per cent confidence from the giant planet sample, which favours stars of above average metallicity. This result is not due to differences in age of the two samples. The formation of debris-generating planetesimals at tens of au thus appears independent of the metal fraction of the primordial disc, in contrast to the growth and migration history of giant planets within a few au. The data generally fit a core accumulation model, with outer planetesimals forming eventually even from a disc low in solids, while inner planets require fast core growth for gas to still be present to make an atmosphere.     

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.     

Forming Jupiter, Saturn, Uranus and Neptune in few million years by core accretion

Giant planet formation process is still not completely understood. The current most accepted paradigm, the core instability model, explains several observed properties of the Solar System’s giant planets but, to date, has faced difficulties to account for a formation time shorter than the observational estimates of protoplanetary disks’ lifetimes, especially for the cases of Uranus and Neptune. In the context of this model, and considering a recently proposed primordial Solar System orbital structure, we performed numerical calculations of giant planet formation. Our results show that if accreted planetesimals follow a size distribution in which most of the mass lies in 30-100 m sized bodies, Jupiter, Saturn, Uranus and Neptune may have formed according to the nucleated instability scenario. The formation of each planet occurs within the time constraints and they end up with core masses in good agreement with present estimations.     

Envelope instability in giant planet formation

We compute the growth of isolated gaseous giant planets for several values of the density of the protoplanetary disk, several distances from the central star and two values for the (fixed) radii of accreted planetesimals. Calculations were performed in the frame of the core instability mechanism and the solids accretion rate adopted is that corresponding to the oligarchic growth regime. We find that for massive disks and/or for protoplanets far from the star and/or for large planetesimals, the planetary growth occurs smoothly. However, notably, there are some cases for which we find an envelope instability in which the planet exchanges gas with the surrounding protoplanetary nebula. The timescale of this instability shows that it is associated with the process of planetesimals accretion. The presence of this instability makes it more difficult the formation of gaseous giant planets.     

High-resolution simulations of the final assembly of Earth-like planets I. Terrestrial accretion and dynamics

The final stage in the formation of terrestrial planets consists of the accumulation of ∼1000-km “planetary embryos” and a swarm of billions of 1-10 km “planetesimals.” During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ∼0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: (1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. (2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. (3) If planetary embryos form relatively slowly, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. (4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.     

The Leonard Award Address Presented 1998 July 27, Dublin,Ireland: On the difficulties of making Earth-like planets

Abstract— Here I discuss the series of events that led to the formation and evolution of our planet to examine why the Earth is unique in the solar system. A multitude of factors are involved: These begin with the initial size and angular momentum of the fragment that separated from a molecular cloud; such random factors are crucial in determining whether a planetary system or a double star develops from the resulting nebula. Another requirement is that there must be an adequate concentration of heavy elements to provide the 2% “rock” and “ice” components of the original nebula. An essential step in forming rocky planets in the inner nebula is the loss of gas and depletion of volatile elements, due to early solar activity that is linked to the mass of the central star. The lifetime of the gaseous nebula controls the formation of gas giants. In our system, fine timing was needed to form the gas giant, Jupiter, before the gas in the nebula was depleted. Although Uranus and Neptune eventually formed cores large enough to capture gas, they missed out and ended as ice giants. The early formation of Jupiter is responsible for the existence of the asteroid belt (and our supply of meteorites) and the small size of Mars, whereas the gas giant now acts as a gravitational shield for the terrestrial planets. The Earth and the other inner planets accreted long after the giant planets, from volatile-depleted planetesimals that were probably already differentiated into metallic cores and silicate mantles in a gas-free, inner nebula. The accumulation of the Earth from such planetesimals was essentially a stochastic process, accounting for the differences among the four rocky inner planets—including the startling contrast between those two apparent twins, Earth and Venus. Impact history and accretion of a few more or less planetesimals were apparently crucial. The origin of the Moon by a single massive impact with a body larger than Mars accounts for the obliquity (and its stability) and spin of the Earth, in addition to explaining the angular momentum, orbital characteristics, and unique composition of the Moon. Plate tectonics (unique among the terrestrial planets) led to the development of the continental crust on the Earth, an essential platform for the evolution of Homo sapiens. Random major impacts have punctuated the geological record, accentuating the directionless course of evolution. Thus a massive asteroidal impact terminated the Cretaceous Period, resulted in the extinction of at least 70% of species living at that time, and led to the rise of mammals. This sequence of events that resulted in the formation and evolution of our planet were thus unique within our system. The individual nature of the eight planets is repeated among the 60-odd satellites—no two appear identical. This survey of our solar system raises the question whether the random sequence of events that led to the formation of the Earth are likely to be repeated in detail elsewhere. Preliminary evidence from the “new planets” is not reassuring. The discovery of other planetary systems has removed the previous belief that they would consist of a central star surrounded by an inner zone of rocky planets and an outer zone of giant planets beyond a few astronomical units (AU). Jupiter-sized bodies in close orbits around other stars probably formed in a similar manner to our giant planets at several astronomical units from their parent star and, subsequently, migrated inwards becoming stranded in close but stable orbits as “hot Jupiters”, when the nebula gas was depleted. Such events would prevent the formation of terrestrial-type planets in such systems.     

A matched filter method for ground-based sub-noise detection of terrestrial extrasolar planets in eclipsing binaries: application to CM Draconis

The photometric detection of extrasolar planets by transits in eclipsing binary systems can be significantly improved by cross-correlating the observational light curves with synthetic models of possible planetary transit features, essentially a matched filter approach. We demonstrate the utility and application of this transit detection algorithm for ground-based detections of terrestrial-sized (Earth-to-Neptune radii) extrasolar planets in the dwarf M-star eclipsing binary system CM Draconis. Preliminary photometric observational data of this system demonstrate that the observational noise is well characterized as white and Gaussian at the observational time steps required for precision photometric measurements. Depending on planet formation scenarios, terrestrial-sized planets may form quite close to this low-luminosity system. We demonstrate, for example, that planets as small as 1.4 Earth radii with periods on the order of a few months in the CM Draconis system could be detected at the 99.9% confidence level in less than a year using 1-m class telescopes from the ground. This result contradicts commonly held assumptions limiting present ground-based efforts to, at best, detections of gas giant planets after several years of observation. This method can be readily extended to a number of other larger star systems with the utilization of larger telescopes and longer observing times. Its extension to spacecraft observations should also allow the determination of the presence of terrestrial-sized planets in nearly 100 other known eclipsing binary systems.     

Terrestrial planet formation surrounding close binary stars

Most stars reside in binary/multiple star systems; however, previous models of planet formation have studied growth of bodies orbiting an isolated single star. Disk material has been observed around both components of some young close binary star systems. Additionally, it has been shown that if planets form at the right places within such disks, they can remain dynamically stable for very long times. Herein, we numerically simulate the late stages of terrestrial planet growth in circumbinary disks around ‘close’ binary star systems with stellar separations 0.05 AU?aB?0.4 AU and binary eccentricities 0?eB?0.8. In each simulation, the sum of the masses of the two stars is 1 M_⊙, and giant planets are included. The initial disk of planetary embryos is the same as that used for simulating the late stages of terrestrial planet formation within our Solar System by Chambers [Chambers, J.E., 2001. Icarus 152, 205-224], and around each individual component of the α Centauri AB binary star system by Quintana et al. [Quintana, E.V., Lissauer, J.J., Chambers, J.E., Duncan, M.J., 2002. Astrophys. J. 576, 982-996]. Multiple simulations are performed for each binary star system under study, and our results are statistically compared to a set of planet formation simulations in the Sun-Jupiter-Saturn system that begin with essentially the same initial disk of protoplanets. The planetary systems formed around binaries with apastron distances QB≡aB(1+eB)?0.2 AU are very similar to those around single stars, whereas those with larger maximum separations tend to be sparcer, with fewer planets, especially interior to 1 AU. We also provide formulae that can be used to scale results of planetary accretion simulations to various systems with different total stellar mass, disk sizes, and planetesimal masses and densities.     

The formation and habitability of terrestrial planets in the presence of close-in giant planets

‘Hot jupiters,’ giant planets with orbits very close to their parent stars, are thought to form farther away and migrate inward via interactions with a massive gas disk. If a giant planet forms and migrates quickly, the planetesimal population has time to re-generate in the lifetime of the disk and terrestrial planets may form [P.J. Armitage, A reduced efficiency of terrestrial planet formation following giant planet migration, Astrophys. J. 582 (2003) L47-L50]. We present results of simulations of terrestrial planet formation in the presence of hot/warm jupiters, broadly defined as having orbital radii ?0.5 AU. We show that terrestrial planets similar to those in the Solar System can form around stars with hot/warm jupiters, and can have water contents equal to or higher than the Earth's. For small orbital radii of hot jupiters (e.g., 0.15, 0.25 AU) potentially habitable planets can form, but for semi-major axes of 0.5 AU or greater their formation is suppressed. We show that the presence of an outer giant planet such as Jupiter does not enhance the water content of the terrestrial planets, but rather decreases their formation and water delivery timescales. We speculate that asteroid belts may exist interior to the terrestrial planets in systems with close-in giant planets.     

Simulations of planet migration driven by planetesimal scattering

Evidence has mounted for some time that planet migration is an important part of the formation of planetary systems, both in the Solar System [Malhotra, R., 1993. Nature 365, 819-821] and in extrasolar systems [Mayor, M., Queloz, D., 1995. Nature 378, 355-359; Lin, D.N.C., Bodenheimer, P., Richardson, D.C., 1996. Nature 380, 606-607]. One mechanism that produces migration (the change in a planet's semi-major axis a over time) is the scattering of comet- and asteroid-size bodies called planetesimals [Fernandez, J.A., Ip, W.-H., 1984. Icarus 58, 109-120]. Significant angular momentum exchange can occur between the planets and the planetesimals during local scattering, enough to cause a rapid, self-sustained migration of the planet [Ida, S., Bryden, G., Lin, D.N.C., Tanaka, H., 2000. Astrophys. J. 534, 428-445]. This migration has been studied for the particular case of the four outer planets of the Solar System (as in Gomes et al. [Gomes, R.S., Morbidelli, A., Levison, H.F., 2004. Icarus 170, 492-507]), but is not well understood in general. We have used the Miranda [McNeil, D., Duncan, M., Levison, H.F., 2005. Astron. J. 130, 2884-2899] computer simulation code to perform a broad parameter-space survey of the physical variables that determine the migration of a single planet in a planetesimal disk. Migration is found to be predominantly inwards, and the migration rate is found to be independent of planet mass for low-mass planets in relatively high-mass disks. Indeed, a simple scaling relation from Ida et al. [Ida, S., Bryden, G., Lin, D.N.C., Tanaka, H., 2000. Astrophys. J. 534, 428-445] matches well with the dependencies of the migration rate:

(1)     

Limits on the Orbits of Possible Eccentric and Inclined Moons of Extrasolar Planets Orbiting Single Stars

Limits are placed on the range of orbits and masses of possible moons orbiting extrasolar planets which orbit single central stars. The Roche limiting radius determines how close the moon can approach the planet before tidal disruption occurs; while the Hill stability of the star–planet–moon system determines stable orbits of the moon around the planet. Here the full three-body Hill stability is derived for a system with the binary composed of the planet and moon moving on an inclined, elliptical orbit relative the central star. The approximation derived here in Eq. (17) assumes the binary mass is very small compared with the mass of the star and has not previously been applied to this problem and gives the criterion against disruption and component exchange in a closed form. This criterion was applied to transiting extrasolar planetary systems discovered since the last estimation of the critical separations (Donnison in Mon Not R Astron Soc 406:1918, 2010a) for a variety of planet/moon ratios including binary planets, with the moon moving on a circular orbit. The effects of eccentricity and inclination of the binary on the stability of the orbit of a moon is discussed and applied to the transiting extrasolar planets, assuming the same planet/moon ratios but with the moon moving with a variety of eccentricities and inclinations. For the non-zero values of the eccentricity of the moon, the critical separation distance decreased as the eccentricity increased in value. Similarly the critical separation decreased as the inclination increased. In both cases the changes though very small were significant.     

Reducing the probability of capture into resonance

A migrating planet can capture planetesimals into mean motion resonances. However, resonant trapping can be prevented when the drift or migration rate is sufficiently high. Using a simple Hamiltonian system for first- and second-order resonances, we explore how the capture probability depends on the order of the resonance, drift rate and initial particle eccentricity. We present scaling factors as a function of the planet mass and resonance strength to estimate the planetary migration rate above which the capture probability drops to less than half. Applying our framework to multiple extrasolar planetary systems that have two planets locked in resonance, we estimate lower limits for the outer planet's migration rate, allowing resonance capture of the inner planet.
Mean motion resonances are comprised of multiple resonant subterms. We find that the corotation subterm can reduce the probability of capture when the planet eccentricity is above a critical value. We present factors that can be used to estimate this critical planet eccentricity. Applying our framework to the migration of Neptune, we find that Neptune's eccentricity is near the critical value that would make its 2 : 1 resonance fail to capture twotinos. The capture probability is affected by the separation between resonant subterms and so is also a function of the precession rates of the longitudes of periapse of both planet and particle near resonance.