Planet X and the Extended Scattered Disk

The effects that a hypothetical trans-Plutonian planet would produce on theorbital distribution of the Classical Edgeworth-Kuiper-Belt, has beensurveyed for different physical and orbital parameters of the hypotheticalbody in Melita et al. (2003a). The best fits were obtained by a moderatelyeccentric and inclined Earth-sized object with a semimajor axis of ～ 70AU. However the history of some objects in the `Extended Scattered disk’still represent a puzzle. One possibility is that they can be `extracted’from the Scattered disk by the planetoid. In this work we confirm that such anhypothesis would not explain the present orbit of 2000 CR105, given theconditions for a gap as observed to be formed in the Classical EKB.     

The plane of the Edgeworth-Kuiper belt

We examine possible locations for the primordial disk of the Edgeworth-Kuiper Belt (EKB), using several subsets of the known objects as markers of the total mass distribution. Using a secular perturbation theory, we find that the primordial plane of the EKB could have remained thin enough to escape detection only if it is clustered very closely about the invariable plane of the Solar System.     

Consequences of tidal interaction between disks and orbiting protoplanets for the evolution of multi-planet systems with architecture resembling that of Kepler 444

We study orbital evolution of multi-planet systems with masses in the terrestrial planet regime induced through tidal interaction with a protoplanetary disk assuming that this is the dominant mechanism for producing orbital migration and circularization. We develop a simple analytic model for a system that maintains consecutive pairs in resonance while undergoing orbital circularization and migration. This model enables migration times for each planet to be estimated once planet masses, circularization times and the migration time for the innermost planet are specified. We applied it to a system with the current architecture of Kepler 444 adopting a simple protoplanetary disk model and planet masses that yield migration times inversely proportional to the planet mass, as expected if they result from torques due to tidal interaction with the protoplanetary disk. Furthermore the evolution time for the system as a whole is comparable to current protoplanetary disk lifetimes. In addition we have performed a number of numerical simulations with input data obtained from this model. These indicate that although the analytic model is inexact, relatively small corrections to the estimated migration rates yield systems for which period ratios vary by a minimal extent. Because of relatively large deviations from exact resonance in the observed system of up to 2 %, the migration times obtained in this way indicate only weak convergent migration such that a system for which the planets did not interact would contract by only \({\sim }1\,\%\) although undergoing significant inward migration as a whole. We have also performed additional simulations to investigate conditions under which the system could undergo significant convergent migration before reaching its final state. These indicate that migration times have to be significantly shorter and resonances between planet pairs significantly closer during such an evolutionary phase. Relative migration rates would then have to decrease allowing period ratios to increase to become more distant from resonances as the system approached its final state in the inner regions of the protoplanetary disk.     

Kuiper Belt structure as a reflection of the migration process of a planet

One of the main particular features of the structure of the Kuiper Belt is that it contains clusters of objects of small orbital eccentricity and inclination (“cold population”). In order to solve the problem of the origin of the objects, we considered the process of the gravitational interaction of a comparatively small-mass planet with a planetesimal disk. We found that one particular property of the process is that the planet changes its direction of migration. The interaction with the planet results in the transportation of a considerable portion of planetesimals from the inner zone out to the Kuiper Belt. After such a transition of the objects, the planet returns to the inner regions of the planetesimal disk. Numerical simulations show that the reversible migration of a planet of a mass similar to that of the Earth can explain the main properties of the Kuiper Belt’s cold population orbit distribution.     

The effect on the Edgeworth–Kuiper Belt of a large distant tenth planet

We investigate the orbital evolution of both real and hypothetical Edgeworth–Kuiper Objects in order to determine whether any conclusions can be drawn regarding the existence, or otherwise, of the tenth planet postulated by Murray. We find no qualitative difference in the orbital evolution, and so conclude that the hypothetical planet has been placed on an orbit at such a large heliocentric distance that no evidence for the existence, or non-existence, can be found from a study of the known Edgeworth–Kuiper Objects.     

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.     

Exploring the 7:4 mean motion resonance—II: Scattering evolutionary paths and resonance sticking

In our preliminary study, we have investigated basic properties and dynamical evolution of classical TNOs around the 7:4 mean motion resonance with Neptune (a∼43.7 AU), motivated by observational evidences that apparently present irregular features near this resonance (see [Lykawka and Mukai, 2005a. Exploring the 7:4 mean motion resonance—I. Dynamical evolution of classical trans-Neptunian objects (TNOs). Space Planet. Sci. 53, 1175-1187]; hereafter “Paper I”). In this paper, we aim to explore the dynamical long-term evolution in the scattered disk (but not its early formation) based on the computer simulations performed in Paper I together with extra computations. Specifically, we integrated the orbital motion of test particles (totalizing a bit more than 10,000) placed around the 7:4 mean motion resonance under the effect of the four giant planets for the age of the Solar System. In order to investigate chaotic diffusion, we also conducted a special simulation with on-line computation of proper elements following tracks in phase space over 4-5 Gyr. We found that: (1) A few percent (1-2%) of the test particles survived in the scattered disk with direct influence of other Neptunian mean motion resonances, indicating that resonance sticking is an extremely common phenomenon and that it helps to enhance scattered objects longevity. (2) In the same region, the so-called extended scattered TNOs are able to form via very long resonance trapping under certain conditions. Namely, if the body spends more than about 80% of its dynamical lifetime trapped in mean motion resonance(s) and there is the action of a k+1 or (k+2)/2 mean motion resonance (e.g., external mean motion resonances with Neptune described as (j+k)/j with j=1 and 2, respectively). According to this hypothetical mechanism, 5-15% of current scattered TNOs would possess thus probably constituting a significant part of the extended scattered disk. (3) Moreover, considering hot orbital initial conditions, it is likely that the trans-Neptunian belt (or Edgeworth-Kuiper belt) has been providing members to the scattered disk, so that scattered TNOs observed today would consist of primordial scattered bodies mixed with TNOs that came from unstable regions of the trans-Neptunian belt in the past.Considering the three points together, our results demonstrated that the scattered disk has been evolving continuously since early times until present.     

Planet X and the origins of the shower and steady state flux of short-period comets

In the planet X model periodic comet showers are associated with passages of the planet's perihelion and aphelion points through a primordial disk of comets believed to lie beyond the orbit of Neptune. A strong feature of this model is that the required orbital elements and mass of planet X are consistent with independently predicted values based on the residuals in the motions of Uranus and Neptune. Here we present a more extensive analysis of the model taking into account the fact that only those comets scattered directly into the zones of influence of Saturn and Jupiter can contribute to a shower whose duration is consistent with observation (≲ 15 myr). These requirements impose a minimum planetary inclination of ≈25°, which in turn restricts the semimajor axis to be ≲100 AU. A fraction of the comets scattered directly into the zones of influence of Uranus and Neptune will evolve on time scales of ∼10⁸ years into the steady state flux of short-period comets. We find that the absolute numbers of shower and steady state are comparable and compatible with the known terrestrial cratering rate, assuming the existence of long-lived extinct comet cores. Canonical planet X model parameters, deduced in part from the scattering dynamics analysis, are: semimajor axis ≈80 AU, eccentricity ≈0.3, inclination ≈45°, and mass ≈5m_⊕. An analysis is given which suggests that planet X, in its present orbit, can create the requisite density gradient of comets near perihelion and aphelion during the lifetime of the Solar System. The required inclination of planet X's orbit (≳25°) may explain the failure of previous surveys to discover the planet as its present latitude is not likely to be near the ecliptic. It it exists, the best immediate hope of finding planet X is the ongoing IRAS search in the 100-μm band and the full sky optical survey by Shoemaker and Shoemaker. Independent of the question of periodic comet showers, the existence of planet X and the comet disk can readily explain the origin of the steady state flux of short-period comets over a wide range of parameters.     

Evolution of a Keplerian disk of colliding and fragmenting particles: a kinetic model with application to the Edgeworth-Kuiper belt

We present a kinetic model of a disk of solid particles, orbiting a primary and experiencing inelastic collisions. In distinction to other collisional models that use a 2D (mass-semimajor axis) binning and perform a separate analysis of the velocity (eccentricity, inclination) evolution, we choose mass and orbital elements as independent variables of a phase space. The distribution function in this space contains full information on the combined mass, spatial, and velocity distributions of particles. A general kinetic equation for the distribution function is derived, valid for any set of orbital elements and for any collisional outcome, specified by a single kernel function. The first implementation of the model utilizes a 3D phase space (mass-semimajor axis-eccentricity) and involves averages over the inclination and all angular elements. We assume collisions to be destructive, simulate them with available material- and size-dependent scaling laws, and include collisional damping. A closed set of kinetic equations for a mass-semimajor axis-eccentricity distribution is written and transformation rules to usual mass and spatial distributions of the disk material are obtained. The kinetic “core” of our approach is generic. It is possible to add inclination as an additional phase space variable, to include cratering collisions and agglomeration, dynamical friction and viscous stirring, gravity of large perturbers, drag forces, and other effects into the model. As a specific application, we address the collisional evolution of the classical population in the Edgeworth-Kuiper belt (EKB). We run the model for different initial disk's masses and radial profiles and different impact strengths of objects. Our results for the size distribution, collisional timescales, and mass loss are in agreement with previous studies. In particular, collisional evolution is found to be most substantial in the inner part of the EKB, where the separation size between the survivors over EKB's age and fragments of earlier collisions lies between a few and several tens of km. The size distribution in the EKB is not a single Dohnanyi-type power law, reflecting the size dependence of the critical specific energy in both strength and gravity regimes. The net mass loss rate of an evolved disk is nearly constant and is dominated by disruption of larger objects. Finally, assuming an initially uniform distribution of orbital eccentricities, we show that an evolved disk contains more objects in orbits with intermediate eccentricities than in nearly circular or more eccentric orbits. This property holds for objects of any size and is explained in terms of collisional probabilities. The effect should modulate the eccentricity distribution shaped by dynamical mechanisms, such as resonances and truncation of perihelia by Neptune.     

Tidal interactions in multi-planet systems

We study systems of close orbiting planets evolving under the influence of tidal circularization. It is supposed that a commensurability forms through the action of disk induced migration and orbital circularization. After the system enters an inner cavity or the disk disperses the evolution continues under the influence of tides due to the central star which induce orbital circularization. We derive approximate analytic models that describe the evolution away from a general first order resonance that results from tidal circularization in a two planet system and which can be shown to be a direct consequence of the conservation of energy and angular momentum. We consider the situation when the system is initially very close to resonance and also when the system is between resonances. We also perform numerical simulations which confirm these models and then apply them to two and four planet systems chosen to have parameters related to the GJ 581 and HD 10180 systems. We also estimate the tidal dissipation rates through effective quality factors that could result in evolution to observed period ratios within the lifetimes of the systems. Thus the survival of, or degree of departure from, close commensurabilities in observed systems may be indicative of the effectiveness of tidal disipation, a feature which in turn may be related to the internal structure of the planets involved.     

Neptune migration model with one extra planet

We explore conventional Neptune migration model with one additional planet of mass at 0.1-2.0M_⊕. This planet inhabited in the 3:2 mean motion resonance with Neptune during planet migration epoch, and then escaped from the Kuiper belt when jovian planets parked near the present orbits. Adding this extra planet and assuming the primordial disk truncated at about 45 AU in the conventional Neptune migration model, it is able to explain the complex structure of the observed Kuiper belt better than the usual Neptune migration model did in several respects, which are the following. (1) High-inclination Plutinos with i?15-35° are produced. (2) Generating the excitation of the classical Kuiper belt objects, which have moderate eccentricities and inclinations. (3) Producing the larger ratio of Neptune’s 3:2 to 2:1 resonant particles, and the lower ratio of particles in the 3:2 resonance to those in the classical belt, which may be more consistent with observations. (4) Finally, several Neptune’s 5:2 resonant particles are obtained. However, numerical experiments imply that this model is a low-probability event. In addition to the low probability, two features produced by this model may be inconsistent with the observations. They are small number of low-inclination particles in the classical belt, and the production of a remnant population with near-circular and low-inclination orbit within . According to our present study, including one extra planet in the conventional Neptune migration model as the scenario we explored here may be unsuitable because of the low probability, and the two drawbacks mentioned above, although this model can explain better several features which is hard to produce by the conventional Neptune migration model. The issues of low-probability event and the lack of low-inclination KBOs in the classical belt are interesting and may be studied further under a more realistic consideration.     

Searches for transplutonian planets using long-period comets

We discuss the dynamical connection of long-period and nearly parabolic comets with hypothetical transplutonian planets. The statistics includes 792 comets with periods P > 200 years. The orbital plane of the parent planet can be determined from the observed distribution of the perihelia and poles of cometary orbits. The radius of a planetary orbit can be calculated using the Radzievsky-Tisseran criterion. We calculated the minimum distance of each of the 792 orbits to 11 hypothetical planetary orbits. Testing for the kinematic connection of comets with transplutonian planets yielded a negative result. The presence of the nodes of cometary orbits in the transplutonian region is shown to be the result of a geometric effect. We found a high concentration of the nodes and perihelia of cometary orbits in the zone of the terrestrial planets.     

Eccentricities of Planets in Binary Systems

The most puzzling property of the extrasolar planets discovered by recent radial velocity surveys is their high orbital eccentricities, which are very difficult to explain within our current theoretical paradigm for planet formation. Current data reveal that at least 25% of these planets, including some with particularly high eccentricities, are orbiting a component of a binary star system. The presence of a distant companion can cause significant secular perturbations in the orbit of a planet. At high relative inclinations, large-amplitude, periodic eccentricity perturbations can occur. These are known as “Kozai cycles” and their amplitude is purely dependent on the relative orbital inclination. Assuming that every planet host star also has a (possibly unseen, e.g., substellar) distant companion, with reasonable distributions of orbital parameters and masses, we determine the resulting eccentricity distribution of planets and compare it to observations? We find that perturbations from a binary companion always appear to produce an excess of planets with both very high (?0.6) and very low (e ? 0.1) eccentricities. The paucity of near-circular orbits in the observed sample implies that at least one additional mechanism must be increasing eccentricities. On the other hand, the overproduction of very high eccentricities observed in our models could be combined with plausible circularization mechanisms (e.g., friction from residual gas) to create more planets with intermediate eccentricities (e? 0.1–0.6).     

The evidence of an early stellar encounter in Edgeworth-Kuiper belt

We investigate the influence of a stellar fly-by encounter on the Edgeworth-Kuiper belt objects through numerical orbital calculations, in order to explain both mass depletion and high orbital inclinations of the classical Edgeworth-Kuiper belt (CEKB) objects, which have semimajor axis of 42-48 AU and perihelia beyond 35 AU. The observationally inferred total mass of the CEKB is ∼1/10 Earth masses, which is only ∼0.02 of that extrapolated from the minimum-mass solar nebula model. The CEKB consists of bimodal population: “hot population” with inclinations i?0.2-0.6 radians and “cold population” with i?0.1. The observationally suggested difference in size and color of objects between the two populations may imply different origins of the two populations. We find that both the depletion of solid materials in the CEKB and the formation of the hot population are accounted for by a single close stellar encounter with pericenter distance of 80-100 AU and inclination relative to the initial protoplanetary disk ?50°-70°. Such a stellar encounter highly pumps up eccentricities of most objects in the CEKB and then their perihelia migrate within 35 AU. These objects would be removed by Neptune's perturbations after Neptune is formed at or migrates to the current position (30 AU). Less than 10% of the original objects remain in stable orbits with small eccentricities and perihelion distances larger than 35 AU, in the CEKB, which is consistent with the observation. We find that i of the remaining objects are as large as that of the observed hot population. The only problem is how to stop Neptune's migration at ∼30 AU, which is addressed in a separate paper. The depletion by the stellar encounter extends deeply into ∼30-35 AU, which provides the basis of the formation model for the cold population through Neptune's outward migration by Levison and Morbidelli (2003, Nature, 426, 419-421). The combination of our model with Levison and Morbidelli's model could consistently explain the mass depletion, truncation at 50 AU, bimodal distribution in i, and differences in size and color between the hot and the cold populations in the CEKB.     

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

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

The Meudon Multicolor Survey (2MS) of Centaurs and trans-neptunian objects: extended dataset and status on the correlations reported

We present here the latest B−V, V−R, and R−I color measurements obtained with the CFH12K mosaic camera of the 3.6-m Canada-France-Hawaii Telescope (CFHT). This work is the latest extension of the Meudon Multicolor Survey (2MS) and extends the total number of Centaurs and trans-neptunian objects (TNOs) in the dataset to 71. With this large and homogeneous dataset, we performed relevant statistical analyses to search for correlations with physical and orbital parameters and interrelations with related populations (cometary nuclei and irregular satellites). With a larger dataset, we confirm the correlations found for the Classical TNOs in our previous survey: some colors are significantly correlated with perihelion distance and inclination. The only exception is with the eccentricity. However, results strongly depend on which objects are considered Classicals, and with a dynamically more restricted definition these correlations are no longer present. We also find that strongly significant trends with orbital parameters are not detected for Centaurs, Plutinos or scattered disk objects (SDOs). We also make for the first time reliable statistical comparison between TNOs and related populations (e.g., Centaurs, irregular satellites, short period comets—i.e., SPCs). We find that (1) the colors of SPCs do not match either their TNO or Centaur precursors, and this suggests that some process modifies the surface of SPCs at entry into the inner Solar System. The only exception concerns colors of SDOs from which we could statistically assess that SPCs and SDOs could be drawn from a same single parent distribution. (2) Not surprisingly, Centaurs are compatible with each of the Edgeworth-Kuiper belt dynamical groups at a highly significant level except with the SDOs. (3) Centaurs' colors still present a strong dichotomy between a neutral/slightly red group (e.g., Chiron) and a very red group (e.g., Pholus). (4) The irregular satellite population is not compatible with any of the Centaur, Plutino or Classical populations; however, the similarity of their color properties with SDOs suggests that both groups can be extracted from the same parent distribution. However, due to the small number of Centaurs and SDOs these conclusions cannot be taken as definitive.     

Migration and evolution of extrasolar planets

Giant planets in circumstellar disks can migrate inward from their initial (formation) positions at several AUs. Inward radial migration of the planet is caused by torques between the planet and the disk; outward radial migration of the planet is caused by torques between the planet and the spinning star, and by torques due to Roche lobe overflow and consequent mass loss from the planet. We present self-consistent numerical considerations of the problem of migrating giant planets by summing torques on planets for various physical parameters of the disk and of planets. We find that Jupiter-mass planets can stably arrive and survive at small heliocentric distances, thus reproducing observed properties of some of the recently discovered extra-solar planets. The range of fates of massive planets is broad, and some perish by losing all their mass onto the central star during Roche lobe overflow, while others survive for the lifetime of the central star. Surviving planets cluster into two groups when examined in terms of final mass and final heliocentric distance: those which have lost mass and those which have not. Some of the observed extrasolar planets fall into each of these two exclusive classes. We also find that there is an inner boundary for planets' final heliocentric distances, caused by tidal torques with the central star. Planets in small orbits are shown to be stable against atmospheric loss.     

The Kuiper Belt and its Primordial Sculpting

We discuss the structure of the Kuiper belt as it can be inferred from the first decade of observations. In particular, we focus on its most intriguing properties – the mass deficit, inclination distribution, the apparent existence of an outer edge and of a correlation among inclinations, colours and sizes – which clearly show that the belt has lost its pristine structureof a dynamically cold proto-planetary disk. Understanding how the Kuiperbelt acquired its present structure will provide insight into the formationof the outer planetary system and on its early evolution. We outline ascenario of primordial sculpting – issued from a combination of mechanismsproposed by various authors – that seems to explain most of the observedproperties of the Kuiper belt. Several aspects are not yet totallyclear. But, for the first time, we have a view – if not of the detailedsculpture – at least of its rough cast.     

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.