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.     

Close encounters of nearly parabolic comets and planets

An overview is given of close encounters of nearly parabolic comets (NPCs; with periods of P > 200 years and perihelion distances of q > 0.1 AU; the number of the comets is N = 1041) with planets. The minimum distances Δ_min between the cometary and planetary orbits are calculated to select comets whose Δ_min are less than the radius of the planet’s sphere of influence. Close encounters of these comets with planets are identified by numerical integration of the comets’ equations of motion over an interval of ±50 years from the time of passing the perihelion. Close encounters of NPCs with Jupiter in 1663–2011 are reported for seven comets. An encounter with Saturn is reported for comet 2004 F2 (in 2001).     

Observable consequences of planet formation models in systems with close-in terrestrial planets

To date, two planetary systems have been discovered with close-in, terrestrial-mass planets     . Many more such discoveries are anticipated in the coming years with radial velocity and transit searches. Here we investigate the different mechanisms that could form 'hot Earths' and their observable predictions. Models include: (1) in situ accretion; (2) formation at larger orbital distance followed by inward 'type 1' migration; (3) formation from material being 'shepherded' inward by a migrating gas giant planet; (4) formation from material being shepherded by moving secular resonances during dispersal of the protoplanetary disc; (5) tidal circularization of eccentric terrestrial planets with close-in perihelion distances and (6) photoevaporative mass-loss of a close-in giant planet. Models 1–4 have been validated in previous work. We show that tidal circularization can form hot Earths, but only for relatively massive planets     with very close-in perihelion distances (≲0.025 au), and even then the net inward movement in orbital distance is at most only 0.1–0.15 au. For planets of less than     , photoevaporation can remove the planet's envelope and leave behind the solid core on a Gyr time-scale, but only for planets inside 0.025–0.05 au. Using two quantities that are observable by current and upcoming missions, we show that these models each produce unique signatures, and can be observationally distinguished. These observables are the planetary system architecture (detectable with radial velocities, transits and transit timing) and the bulk composition of transiting close-in terrestrial planets (measured by transits via the planet's radius).     

Orbits of short period comets captured from the Oort cloud

Oort cloud comets occasionally obtain orbits which take them through the planetary region. The perturbations by the planets are likely to change the orbit of the comet. We model this process by using a Monte Carlo method and cross sections for orbital changes, i.e. changes in energy, inclination and perihelion distance, in a single planet-comet encounter. The influence of all major planets is considered. We study the distributions of orbital parameters of observable comets, i.e. those which have perihelion distance smaller than a given value. We find that enough comets are captured from the Oort cloud in order to explain the present populations of short period comets. The median value of cos i for the Jupiter family is 0.985 while it is 0.27 for the Halley types. The results may explain the orbital features of short period comets, assuming that the active lifetime of a comet is not much greater than 400 orbital revolutions.     

On the process of accretion in the formation of the planets and comets

The physically reasonable assumption that the seed bodies which initiated the accretion of the individual asteroids, planets, and comets (subsequently these objects are collectively called planetoids) formed by stochastic processes requires a radius distribution function which is unique except for two scaling parameters: the total number of planetoids and their most probable radius. The former depends on the ease of formation of the seed bodies while the second is uniquely determined by the average pre-encounter velocity, V, of the accretable material relative to an individual planetoid. This theoretical radius function can be fit to the initial asteroid radius distribution which Anders (1965) derived from the present-day distribution by allowing for fragmentation collisions among the asteroids since their formation. Normalizing the theoretical function to this empirical distribution reveals that there were about 10² precollision asteroids and that V = (2?4) × 10^?2 km/sec which was presumably the turbulent velocity in the Solar Nebula. Knowing V we can determine the scale height of the dust in the Solar Nebula and consequently its space density. The density of accretable material determines the rate of accretion of the planetoids. From this we find, for example, that the Earth formed in about 8 × 10⁶ yr and it attained a maximum temperature through accretion of about 3 × 10³°K. From the total mass of the terrestrial planets and the theoretical radius function we find that about 2 × 10³ planetoids formed in the vicinity of the terrestrial planets. Except for the asteroids the smaller planetoids have since been accreted by the terrestrial planets. About 15% of the present mass of the terrestrial planets was accumulated by the secondary accretion of these smaller primary planetoids. There are far fewer primary planetoids than craters on the Moon or Mars. The craters were likely produced by the collisional breakup of a few primary planetoids with masses between one-tenth and one lunar mass. This deduction comes from comparing the collision cross sections of the planetoids in this mass range to that of the terrestrial planets. This comparison shows that two to three collisions leading to the breakup of four to six objects likely occurred among these objects before their accretion by the terrestrial planets. The number of these fragments is quite adequate to explain the lunar and Martin craters. Furthermore the mass spectrum of such fragments is a power-law distribution which results in a power-law distribution of crater radii of just the type observed on the Moon and Mars. Applying the same analysis to the planetoids which formed in the vicinity of the giant planets reveals that it is unlikely that any fragmentation collisions took place among them before they were accreted by these planets due to the integrated collision cross section of the giant planets being about three orders of magnitude greater than that of the terrestrial planets. We can thus anticipate a marked scarcity of impact craters on the satellites of these outer planets. This prediction can be tested by future space probes. Our knowledge of the radius function of the comets is consistent with their being primary planetoids. The primary difference between the radius function of the planetoids which formed in the inner part of the solar system and that of the comets results from the fact that the seed bodies which grew into the comets formed far more easily than those which grew into the asteroids and the terrestrial planets. Thus in the outer part of the Solar Nebula the principal solid material (water and ammonia snow) accreted into a huge (～10¹²⁺) number of relatively small objects (comets) while in the inner part of the nebula the solid material (hard-to-stick refractory substances) accumulated into only a few (～10³) large objects (asteroids and terrestrial planets). Uranus and Neptune presumably formed by the secondary accretion of the comets.     

Orbits of short period comets captured from the Oort cloud

Oort cloud comets occasionally obtain orbits which take them through the planetary region. The perturbations by the planets are likely to change the orbit of the comet. We model this process by using a Monte Carlo method and cross sections for orbital changes, i.e. changes in energy, inclination and perihelion distance, in a single planet-comet encounter. The influence of all major planets is considered. We study the distributions of orbital parameters of observable comets, i.e. those which have perihelion distance smaller than a given value. We find that enough comets are captured from the Oort cloud in order to explain the present populations of short period comets. The median value of cos i for the Jupiter family is 0.985 while it is 0.27 for the Halley types. The results may explain the orbital features of short period comets, assuming that the active lifetime of a comet is not much greater than 400 orbital revolutions.     

Orbital evolution of long-period comets I. Trivariate Monte Carlo simulation

The mean orbital evolution of long-period comets for 16 representative initial orbits to short-period comets is calculated by a Monte Carlo method. First, trivariate perturbation distributions of barycentric Kepler energy, total angular momentum, and its z component in single encounters of comets with Jupiter are obtained numerically. Their characteristics are examined in detail and the distributions are found to be simple, symmetric, and easy to handle. Second, utilizing these distributions, we have done trivariate Monte Carlo simulations of the orbital evolution of long-period comets, with special emphasis on high-inclination orbits. About half of the 16 initial orbits are traced up to 5000 returns. For each of these orbits, the mean values of semimajor axis, perihelion distance, and inclination; their standard deviations, survival, and capture rates; as well as time scales of orbital evolution are calculated as functions of return number. Survival rates of the initial orbits with high inclination (～90°) and small perihelion distance (～1–2 AU) have been found to be only two or three times smaller than those of the main-source orbits of short-period comets established quantitatively by Everhart. The time scales of orbitsl evolution of the former, however, are nearly 10 times longer than the latter. There is a general trend that, for smaller perihelion distance, the survival efficiency becomes higher. The results of this paper should be considered a basis for a succeeding paper (Paper II) in which the physical lifetime of comets will be determined, and a comparison with the orbital data will be done.     

Evolution of comet orbits under the perturbing influence of the giant planets and nearby stars

The orbital evolution of 500 hypothetical comets during 10⁹ years is studied numerically. It is assumed that the birthplace of such comets was the region of Uranus and Neptune from where they were deflected into very elongated orbits by perturbations of these planets. Then, we adopted the following initial orbital elements: perihelion distances between 20 and 30 AU, inclinations to the ecliptic plane smaller than 20°, and semimajor axes from 5 × 10³ to 5 × 10⁴ AU. Gravitational perturbations by the four giant planets and by hypothetical stars passing at distances from the Sun smaller than 5 × 10⁵ AU are considered. During the simulation, somewhat more than 50% of the comets were lost from the solar system due to planetary or stellar perturbations. The survivors were removed from the planetary region and left as members of what is generally known as the cometary cloud. At the end of the studied period, the semimajor axes of the surviving comets tend to be concentrated in the interval 2 × 10⁴ < a < 3 × 10⁴ AU. The orbital planes of the comets with initial $\text{a ≧ 3 × 10}{}^4\text{AU}$ acquired a complete randomization while the others still maintain a slight predominance of direct orbits. In addition, comet orbits with final a < 6 × 10⁴AU preserve high eccentricities with an average value greater than 0.8 Most “new” comets from the sample entering the region interior to Jupiter's orbit had already registered earlier passages through the planetary region. By scaling up the rate of paritions of hypothetical new comets with the observed one, the number of members of the cometary cloud is estimated to be about 7 × 10¹⁰ and the conclusion is drawn that Uranus and Neptune had to remove a number of comets ten times greater.     

Relationship of comets and planets

We analyze our earlier data on the numerical integration of the equations of motion for 274 short-period comets (with the period P<200 yr) on a time interval of 6000 yr. As many as 54 comets had no close approaches to planets, 13 comets passed through the Saturnian sphere of action, and one comet passed through the Uranian sphere of action. The orbital elements of these 68 comets changed by no more than ±3 percent in a space of 6000 yr. As many as 206 comets passed close to Jupiter. We confirm Everhart’s conclusion that Jupiter can capture long-period comets with q = 4–6 AU and i < 9° into short-period orbits. We show that nearly parabolic comets cross the solar system mainly in the zone of terrestrial planets. No relationship of nearly parabolic comets and terrestrial planets was found for the epoch of the latest apparition of comets. Guliev’s conjecture about two trans-Plutonian planets is based on the illusory excess of cometary nodes at large heliocentric distances. The existence of cometary nodes at the solar system periphery turns out to be a solely geometrical effect.     

A co-accretional model of satellite formation

Numerical calculation of a simple accretion model including the effects of tidal friction indicate that coformation is tenable only if the planet's Q is less than about 10³. The parameter which most strongly affects the final mass ratio of the pair is the time at which the secondary embryo is introduced. Our model yields the proper Moon-Earth mass ratio if the Moon embryo is introduced when the Earth is only about $\text{1}\text{10}$ of its final mass. The lunar orbit remains at about 10 Earth radii throughout most of the growth.This model of satellite formation overcomes two difficulties of the “circumterrestrial cloud” model of Ruskol (1960, 1963, 1972): (1) The difficulty of accumulating a mass as great as the entire Moon before gravitational instability reduces the cloud to a small number of moonlets is removed. (2) The differences between terrestrial and outer planet satellite systems is easily understood in terms of the differences in Q between these planets. The high Q of the outer planets does not allow a satellite embryo to survive a significant portion of the accretion process, thus only small bodies which formed very late in the accumulation of the planet remain as satellites. The low Q of the terrestrial planets allows satellite embryos of these planets to survive during accretion, thus massive satellites such as the Earth's Moon are expected. The present lack of such satellites of the other terrestrial planets may be the result of tidal evolution, either infall following primary despinning (Burns, 1973) or escape due to increase in orbit eccentricity.     

Long-term evolution of Oort Cloud comets: methods and comparisons

Two long-term simulation methods for cometary orbits, a Monte Carlo method and a direct integration method, are compared with each other. The comparison is done in seven inclination and perihelion distance intervals, and shows differences in dynamical lifetime and capture probabilities for the following main reasons. We use a finite energy step approximation in the Monte Carlo method and the method considers only close approaches with the planets. The differences can be taken into account statistically and it is possible to calculate the correction factors for the capture probability and dynamical lifetime in the Monte Carlo method. Both corrections depend on the inclination and on the value of the minimum energy step. The capture probabilities of the short-period comets originating in the Oort Cloud are calculated by the corrected Monte Carlo method and compared with published results.     

On the relationship of long-period comets to known and unknown planets

Statistical aspects of the question of the dynamical relationship of long-period comets to giant planets (Saturn, Uranus, Neptune), dwarf planets (Pluto, 2003 EL₆₁, 2003 UB₃₁₃), and five hypothetical planets are investigated. Data for 859 and 888 comets (2005, 2006) with periods P > 200 years are used in the statistics. No comets of the Kreutz, Marsden, Kracht, and Meyer groups are considered. The minimum interorbital distances of comets from the listed planetary bodies are mainly analyzed. Effects testifying to the kinematical relationship of some of the comets to planets have been established by testing the cometary data relative to 67 planes. The distribution of the perihelia of long-period comets is also discussed.     

Evolution of cometary perihelion distances in oort cloud: Another statistical approach

The evolution of the perihelion distance distribution in the Oort cloud was studied over the age of the solar system, under the gravitational perturbations of random passing stars, using a statistical approach. These perturbations are accounted for through an empirical relation relating the change in cometary perihelion distance to the closest-approach comet-star distance; this relation is deduced from a previous study [H. Scholl, A. Cazenave, and A. Brahic, Astron. Astrophys.112, 157–166 (1982)]. Two kinds of initial perihelion distances are considered: (a) perihelion distances <2500 AU, associated with an origin of comets as icy planetesimals in the region of the giant planets, and (b) larger perihelion distances (up to 5 × 10⁴ AU), possibly representative of comet formation as satellite fragments in the accretion disk of the primitive solar nebula. Distant star-comet encounters, as well as rare close encounters, are considered. Several quantities are estimated: (i) number of “new” comets entering into the planetary region, (ii) number of comets escaping the Sun sphere of influence or lost by hyperbolic ejection and (iii) percentage of total comet loss over the age of the solar system. From these quantities, the current and original cloud populations are deduced, as well as the corresponding cloud mass, for the two types of formation scenarios.     

Asteroids with large secular orbital variations

As the classical linear theory of secular perturbations for asteroids is known not to be adequate for computing the perturbations of asteroids with high eccentricities and/or inclinations, a seminumerical method to calculate the secular perturbations by including higher-degree terms in the disturbing function has been developed. It is here applied to asteroids with small values of $\text{(1 ? e}{}^2\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{cos}\text{i}$, since the secular variations as well as their deviations from the results derived by the classical linear theory are generally large for such asteroids. It is found that the arguments of perihelion for five of the numbered asteroids are librating around 90 or 270°. For asteroids with $\text{(1 ? e}{}^2\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext><mtext>cos</mtext><mtext>\; i</mtext>}}$ less than 0.85 the results of the secular variations are tabulated. Also the stability of such orbits is discussed by comparing the orbital properties of short-periodic comets with them. Generally speaking, orbits of the asteroids are more stable than those of the short-periodic comets, and asteroids with librating arguments of perihelion are more stable than those with circular coplanar orbits although their orbital elements are changed more by secular perturbations.     

The magnitude distribution, perihelion distribution and flux of long-period comets

The mass distribution and perihelion distribution of long-period comets are re-assessed. The mass distribution index is found to be 1.598±0.016 , indicating that the distribution is somewhat steeper than was obtained by previous analyses of an amalgam of all the available historical data. The number of long-period comets that have orbital perihelion distances, q , that fall in a specific q to q +d q range is found to be independent of q . It is also noted that the flux of long-period comets to the inner Solar system has remained constant throughout recorded history.
The number of long-period comets, , per 1-au interval of perihelion distance, per year, brighter than H , entering the inner Solar system is found to be given by log₁₀ =−2.607+0.359 H . It is therefore estimated that, for example, about 0.5, 30 and 2000 long-period comets with absolute magnitudes brighter than 0, 5 and 10 respectively pass the Sun on orbits with perihelion distances less than 2.0 au, every century.     

Dynamical capture and physical decay of short-period comets

The brightness evolution of short-period comets is discussed in connection with their physical lifetimes. It is shown that changes in the fraction of the free-subliming area of the nuclear surface may be more important than mass decrease in determining brightness variations. The decrease in the activity of short-period comets caused by the buildup of a dust mantle may be interrupted—and partially reversed—by dust blowoffs that leave exposed areas of fresh ices. Short-period comets may thus be subject to random brightness fluctuations that make quite uncertain any derivation of their physical lifetime based on comparisons of their absolute brightness at different apparitions. As an alternate procedure, the numerical integration of the whole sample of short-period comet orbits carried out by A. Carusi, L.Kresák, E. Perozzi and G. B. Valsecchi (1984, Long-Term Evolution of Short-Period Comets. Istituto Astrofisica Spaziale Internal Report 12, Rome) is used to draw conclusions about the transfer rate of their perihelia from Jupiter's region to the region of the terrestrial planets (heliocentric distances<1.5 AU). It is found that about one short-period comet per century reaches the region of the terrestrial planets. From this result and under the assumption of a steady-state comet population, an average lifetime of the order of 6 × 10³ years (～10³ revolutions) is derived for a typical kilometer-sized short-period comet of perihelion distance q ～ 1 AU. Such a rather long comet lifetime, as compared to some previous derivations, is consistent with the survival of some periodic comets on small-q orbits of long dynamical time scales.     

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.     

Bidirectional reflectance spectroscopy: 3. Correction for macroscopic roughness

A mathematically rigorous formalism is derived by which an arbitrary photometric function for the bidirectional reflectance of a smooth surface may be corrected to include effects of general macroscopic roughness. The correction involves only one arbitrary parameter, the mean slope angle $\text{θ}$, and is applicable to surfaces of any albedo. Using physically reasonable assumptions and mathematical approximations the correction expressions are evaluated analytically to second order in $\text{θ}$. The correction is applied to the bidirectional reflectance function of B. Hapke (1981, J. Geophys. Res.86, 3039–3054). Expressions for both the differential and integral brightnesses are obtained. Photometric profiles on hypothetical smooth and rough planets of low and high albedo are shown to illustrate the effects of macroscopic roughness. The theory is applied to observations of Mercury and predicts the integral phase function, the apparent polar darkening, and the lack of limb brightness surge on the planet. The roughness-corrected bidirectional reflectance function is sufficiently simple that it can be conveniently evaluated on a programmable hand-held calculator.     

Dynamical evolution of comets and the problem of cometary fading

Possibilities to explain the observed 1/a-distribution are discussed in the light of improved understanding of the dynamical evolution of long-period comets. It appears that the ‘fading problem’ applies both to single-injection and continuous-injection models. Although uncertainties due to nongravitational effects do not allow detailed results to be drawn from the observed 1/a-distribution at small perihelion distance q, that for q ? 1.5 AU shows that a constant fading probability cannot explain all the features of the observed distribution. Assuming that comets can reappear following a period of fading, values for the assumed constant fading and renewal probabilities, and the total cometary flux have been estimated for q > 1.5 AU.     

Origin and evolution of the long-period comets

We consider the changes of cometary perihelion distances as a process of diffusion in the value of q, due to perturbations by stars. We find more comets at large q values than is observed. This suggests that a large number of long-period comets is not observed.