Study of dynamical processes at the final stages of planetary system formation: Resonance motion of giant planets

According to current observational data, planets of many exoplanetary systems have resonant motion. The formation of resonance configurations is studied within a unified model of planetary migration. Planets in the observed systems 24 Sex, HD 37124, HD 73526, HD 82943, HD 128311, HD 160691, Kepler 9, NN Ser, which are moving in the 2: 1 resonance, could have been captured into this resonance due to both the Type I and II migration with a wide range of parameters. The migration conditions are defined for the formation of HD 45364 and HD 200964 that are in the 3: 2 and 4: 3 first-order resonances, correspondingly. The results obtained for HD 200964 show that planets can be captured in the first-order resonances, when the outer-to-inner orbital period ratios for the planets are less than 3: 2, only if Type I migration rates are large, and the mass of at least one planet is substantially less than the modern masses of the observed giant planets. The formation of the HD 102272, HD 108874, HD 181433 and HD 202206 systems with planets in high-order resonances is considered. The capture into these resonances can be realized with very slow Type II migration. Possible bounds for migration parameters are considered. In particular, it has been found that the capture of HD 108874 into the 4: 1 resonance is possible only if the angle between the plane of planetary orbits and the plane of sky is appreciably less than 90°, i.e., the planetary masses are a few times larger than the minimum values. The capture of HD 202206 into the 5: 1 resonance is possible at low migration rates; however, another mechanism is required to explain the high observed eccentricity of the inner planet (for example, strong gravitational interaction between the planets). Resonant configurations can be disrupted due to the interaction between planets and remaining fragments of the planetesimal disk as, for example, may occur in the three-planet system 47 UMa. The specific orbital features observed for this system are explained.     

Orbital Migration Of Giant Planets: Using Numerical Integration To Investigate Consequences For Other Bodies

A number of extrasolar planets have been detected in close orbits around nearby stars. It is probable that these planets did not form in these orbits but migrated from their formation locations beyond the ice line. Orbital migration mechanisms involving angular momentum transfer through tidal interactions between the planets and circumstellar gas-dust disks or by gravitational interaction with a residual planetesimal disk together with several means of halting inward migration have been identified. These offer plausible schemes to explain the orbits of observed extrasolar giant planets and giant planets within the Solar System. Recent advances in numerical integration methods and in the power of computer workstations have allowed these techniques to be applied to modelling directly the mechanisms and consequences of orbital migration in the Solar System. There is now potential for these techniques also to be applied to modelling the consequences of the orbital migration of planets in the observed exoplanetary systems. In particular the detailed investigation of the stability of terrestrial planets in the habitable zone of these systems and the formation of terrestrial planets after the dissipation of the gas disk is now possible. The stability of terrestrial planets in the habitable zone of selected exoplanetary systems has been established and the possibility of the accretion of terrestrial planets in these systems is being investigated by the author in collaboration with Barrie W. Jones (Open University), and with John Chambers (NASA-Ames) and Mark Bailey of Armagh Observatory, using numerical integration. The direct simulation of orbital migration by planetesimal scattering must probably await faster hardware and/or more efficient algorithms. This revised version was published online in July 2006 with corrections to the Cover Date.     

Multiple planetary systems: Properties of the current sample

We carry out analyses on stellar and planetary properties of multiple exoplanetary systems in the currently available sample. With regards to the stars, we study their temperature, distance from the Sun, and metallicity distributions, finding that the stars that harbour multiple exoplanets tend to have subsolar metallicities, in contrast to metal-rich Hot Jupiter hosts; while non-Hot Jupiter single planet hosts form an intermediate group between these two, with approximately solar metallicities. With regards to the planetary systems, we select those with four or more planets and analyse their configurations in terms of stability (via Hill radii), compactness, and size variations. We find that most planetary pairs are stable, and that the compactness correlates to the size variation: More compact systems have more similarly sized planets and vice versa. We also investigate the spectral energy distributions of the stars hosting multiple exoplanetary systems, seeking infra-red excesses that could indicate the presence of debris disks. These disks would be leftovers from the planetary formation process, and could be considered as analogues of the Solar System’s Asteroid or Kuiper belts. We identify potential candidates for disks that are good targets for far infra-red follow-up observations to confirm their existence.     

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.     

Jupiter and Super-Earth embedded in a gaseous disc

In this paper we investigate the evolution of a pair of interacting planets – a Jupiter-mass planet and a Super-Earth with a mass of  5.5 M _⊕  – orbiting a Solar-type star and embedded in a gaseous protoplanetary disc. We focus on the effects of type I and II orbital migrations, caused by the planet–disc interaction, leading to the capture of the Super-Earth in first-order mean-motion resonances by the Jupiter. The stability of the resulting resonant system in which the Super-Earth is on the internal orbit relative to the Jupiter is studied numerically by means of full 2D hydrodynamical simulations. Our main aim is to determine the Super-Earth behaviour in the presence of the gas giant in the system. It is found that the Jupiter captures the Super-Earth into the interior 3:2 or 4:3 mean-motion resonance, and that the stability of such configurations depends on the initial positions of the planets and on the evolution of the eccentricity. If the initial separation of the orbits of the planets is larger than or close to that required for the exact resonance, the final outcome is the migration of the pair of planets at a rate similar to that of the gas giant, at least for the time of our simulations. Otherwise, we observe a scattering of the Super-Earth from the disc. The evolution of planets immersed in a gaseous disc is compared with their behaviour in the case of the classical three-body problem when the disc is absent.     

Capture into first-order resonances and long-term stability of pairs of equal-mass planets

Massive planets form within the lifetime of protoplanetary disks, and therefore, they are subject to orbital migration due to planet–disk interactions. When the first planet reaches the inner edge of the disk, its migration stops and consequently the second planet ends up locked in resonance with the first one. We detail how the resonant trapping works comparing semi-analytical formulae and numerical simulations. We restrict to the case of two equal-mass coplanar planets trapped in first-order resonances, but the method can be easily generalized. We first describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using series expansions up to different orders in eccentricity as well as a non-expanded Hamiltonian. Then we show that during convergent migration the planets evolve along these families of equilibrium points. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict precisely analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that in most cases the resonant configurations achieved by migration become unstable after the removal of the protoplanetary disk. Here we probe the stability of the resonances as a function of planetary mass. For this purpose, we fictitiously increase the masses of resonant planets, adiabatically maintaining the low-amplitude libration regime until instability occurs. We discuss two hypotheses for the instability, that of a low-order secondary resonance of the libration frequency with a fast synodic frequency of the system, and that of minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter nevertheless play the decisive role in the destabilization of resonant pairs. We show evidence that as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.     

Close encounters in young stellar clusters: implications for planetary systems in the solar neighbourhood

The stars that populate the solar neighbourhood were formed in stellar clusters. Through N -body simulations of these clusters, we measure the rate of close encounters between stars. By monitoring the interaction histories of each star, we investigate the singleton fraction in the solar neighbourhood. A singleton is a star which formed as a single star, has never experienced any close encounters with other stars or binaries, or undergone an exchange encounter with a binary. We find that, of the stars which formed as single stars, a significant fraction is not singletons once the clusters have dispersed. If some of these stars had planetary systems, with properties similar to those of the Solar System, the planets' orbits may have been perturbed by the effects of close encounters with other stars or the effects of a companion star within a binary. Such perturbations can lead to strong planet–planet interactions which eject several planets, leaving the remaining planets on eccentric orbits. Some of the single stars exchange into binaries. Most of these binaries are broken up via subsequent interactions within the cluster, but some remain intact beyond the lifetime of the cluster. The properties of these binaries are similar to those of the observed binary systems containing extrasolar planets. Thus, dynamical processes in young stellar clusters will alter significantly any population of Solar System-like planetary systems. In addition, beginning with a population of planetary systems exactly resembling the Solar System around single stars, dynamical encounters in young stellar clusters may produce at least some of the extrasolar planetary systems observed in the solar neighbourhood.     

Inner Solar System dynamical analogs of plutinos

By studying orbits of asteroids potentially in 3:2 exterior mean motion resonance with Earth, Venus, and Mars, we have found plutino analogs. We identify at least 27 objects in the inner Solar System dynamically protected from encounter through this resonance. These are four objects associated with Venus, six with Earth, and seventeen with Mars. Bodies in the 3:2 exterior resonance (including those in the plutino resonance associated with Neptune) orbit the Sun twice for every three orbits of the associated planet, in such a way that with sufficiently low libration amplitude close approaches to the planet are impossible. As many as 15% of Kuiper Belt objects share the 3:2 resonance, but are poorly observed. One of several resonance sweeping mechanisms during planetary migration is likely needed to explain the origin and properties of 3:2 resonant Kuiper Belt objects. Such a mechanism likely did not operate in the inner Solar System. We suggest that scattering by the next planet out allows entry to, and exit from, 3:2 resonance for objects associated with Venus or Earth. 3:2 resonators of Mars, on the other hand, do not cross the paths of other planets, and have a long lifetime. There may exist some objects trapped in the 3:2 Mars resonance which are primordial, with our tests on the most promising objects known to date indicating lifetimes of at least tens of millions of years. Identifying 3:2 resonant systems in the inner Solar System permits this resonance to be studied on shorter timescales and with better determined orbits than has been possible to date, and introduces new mechanisms for entry into the resonant configuration.     

The contribution of the ARIEL space mission to the study of planetary formation

The study of extrasolar planets and of the Solar System provides complementary pieces of the mosaic represented by the process of planetary formation. Exoplanets are essential to fully grasp the huge diversity of outcomes that planetary formation and the subsequent evolution of the planetary systems can produce. The orbital and basic physical data we currently possess for the bulk of the exoplanetary population, however, do not provide enough information to break the intrinsic degeneracy of their histories, as different evolutionary tracks can result in the same final configurations. The lessons learned from the Solar System indicate us that the solution to this problem lies in the information contained in the composition of planets. The goal of the Atmospheric Remote-Sensing Infrared Exoplanet Large-survey (ARIEL), one of the three candidates as ESA M4 space mission, is to observe a large and diversified population of transiting planets around a range of host star types to collect information on their atmospheric composition. ARIEL will focus on warm and hot planets to take advantage of their well-mixed atmospheres, which should show minimal condensation and sequestration of high-Z materials and thus reveal their bulk composition across all main cosmochemical elements. In this work we will review the most outstanding open questions concerning the way planets form and the mechanisms that contribute to create habitable environments that the compositional information gathered by ARIEL will allow to tackle.     

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.     

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.     

On the migration of a system of protoplanets

The evolution of a system consisting of a protoplanetary disc with two embedded Jupiter-sized planets is studied numerically. The disc is assumed to be flat and non-self-gravitating; this is modelled by the planar (two-dimensional) Navier–Stokes equations. The mutual gravitational interaction of the planets and the star, and the gravitational torques of the disc acting on the planets and the central star are included. The planets have an initial mass of one Jupiter mass M _Jup each, and the radial distances from the star are one and two semimajor axes of Jupiter, respectively.
During the evolution a joint wide annular gap is created by the planets. Both planets increase their mass owing to accretion of gas from the disc: after about 2500 orbital periods of the inner planet it has reached a mass of 2.3  M _Jup, while the outer planet has reached a mass of 3.2  M _Jup. The net gravitational torques exerted by the disc on the planets result in an inward migration of the outer planet on time-scales comparable to the viscous evolution time of the disc. The semimajor axis of the inner planet remains constant as there is very little gas left in its vicinity to induce any migration. When the distance of close approach eventually becomes smaller than the mutual Hill radius, the eccentricities increase strongly and the system may become unstable.
If disc depletion occurs rapidly enough before the planets come too close to each other, a stable system similar to our own Solar system may remain. Otherwise the orbits may become unstable and produce systems like υ And.     

Disc–Planet Interactions: Migration and Resonances in Extrasolar Planetary Systems

We consider orbital resonances in multiplanet systems. These are expected to arise during or just after formation in a gaseous disc. Disc–planet interaction naturally produces orbital migration and circularization through the action of tidal torques which in turn may lead to an orbital resonance. The mass and angular momentum content of the disc is likely to be comparable to that in the planets so that it is essential to fully incorporate the disc in the analysis.We study the orbital evolution of two planets locked in 2:1 commensurability through migration tidally induced by the disc using both analytic methods and numerical hydrodynamic simulations. The planets are assumed to orbit in an inner cavity containing at most only a small amount of disc material. Results are found to be sensitive to initial surface density profile, planet masses and disc parameters. The evolution may range between attaining and subsequently maintaining a resonance lock with two angles librating to divergent migration with no commensurability formed. In the former case eccentricities increase monotonically with time while the system undergoes inward migration. If the migration is halted by loss of the disc leaving the planets in a final configuration, there is likely to be a low probability of seeing resonant planets at small radii as well as a sensitive dependence on past history.We have also considered a multiplanet system in secular apsidal resonance. We consider the system as being in just one secular normal mode and include the effects of a gaseous disc. It is suggested that a normal mode may be selected by adding in some weak dissipative process in the disc and that it may remain, involving only the planets, when the disc is slowly removed.     

Physical significance and the role of Lagrangian points and the Oort clouds of planets in the Solar System

Here, the role played by Lagrangian points and the Oort clouds of planets in the evolution and structure of the Solar System has been discussed. It is revealed that the Lagrangian points are not mere isolated points in space associated with the orbit of a planet at which the resultant gravitational force of the system of three bodies is zero as thought previously to be, but their existence has much deeper physical significance as regards the origin of the Solar System and those of satellite systems of planets.     

Planet formation around M-dwarfs: the moving snow line and super-Earths

Planets result from a series of processes within a circumstellar disk. Evidence comes from the near planar orbits in the Solar System and other planetary systems, observations of newly formed disks around young stars, and debris disks around main-sequence stars. As planet-hunting techniques improve, we approach the ability to detect systems like the Solar System, and place ourselves in context with planetary systems in general. Along the way, new classes of planets with unexpected characteristics are discovered. One of the most recent classes contains super Earth-mass planets orbiting a few AU from low-mass stars. In this contribution, we outline a semi-analytic model for planet formation during the pre-main sequence contraction phase of a low-mass star. As the star contracts, the “snow line”, which separates regions of rocky planet formation from regions of icy planet formation, moves inward. This process enables rapid formation of icy protoplanets that collide and merge into super-Earths before the star reaches the main sequence. The masses and orbits of these super-Earths are consistent with super-Earths detected in recent microlensing experiments.     

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)     

Predictions for the frequency and orbital radii of massive extrasolar planets

We investigate the migration of massive extrasolar planets caused by gravitational interaction with a viscous protoplanetary disc. We show that a model in which planets form at 5 au at a constant rate, before migrating, leads to a predicted distribution of planets that is a steeply rising function of log( a ), where a is the orbital radius. Between 1 and 3 au, the expected number of planets per logarithmic interval in a roughly doubles. We demonstrate that, once selection effects are accounted for, this is consistent with current data, and then extrapolate the observed planet fraction to masses and radii that are inaccessible to current observations. In total, approximately 15 per cent of stars targeted by existing radial velocity searches are predicted to possess planets with masses  0.3< M _p sin( i )<10 M _J  and radii  0.1< a <5 au  . A third of these planets (around 5 per cent of the target stars) lie at the radii most amenable to detection via microlensing. A further  5–10  per cent of stars could have planets at radii of  5< a <8 au  that have migrated outwards. We discuss the probability of forming a system (akin to the Solar system) in which significant radial migration of the most massive planet does not occur. Approximately  10–15  per cent of systems with a surviving massive planet are estimated to fall into this class. Finally, we note that a smaller fraction of low-mass planets than high-mass planets is expected to survive without being consumed by the star. The initial mass function for planets is thus predicted to rise more steeply towards small masses than the observed mass function.     

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

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

The total number of giant planets in debris discs with central clearings

Infrared spectra from the Spitzer Space Telescope ( SSC ) of many debris discs are well fit with a single blackbody temperature which suggest clearings within the disc. We assume that clearings are caused by orbital instability in multiple planet systems with similar configurations to our own. These planets remove dust-generating planetesimal belts as well as dust generated by the outer disc that is scattered or drifts into the clearing. From numerical integrations, we estimate a minimum planet spacing required for orbital instability (and so planetesimal and dust removal) as a function of system age and planet mass. We estimate that a 10⁸ yr old debris disc with a dust disc edge at a radius of 50 au hosted by an A star must contain approximately five Neptune mass planets between the clearing radius and the iceline in order to remove all primordial objects within it. We infer that known debris disc systems contain at least a fifth of a Jupiter mass in massive planets. The number of planets and spacing required is insensitive to the assumed planet mass. However, an order of magnitude higher total mass in planets could reside in these systems if the planets are more massive.