Long term dynamical evolution and classification of classical TNOs

Classical trans-Neptunian objects (TNOs) are believed to represent the most dynamically pristine population in the trans-Neptunian belt (TNB) offering unprecedented clues about the formation of our Solar System. The long term dynamical evolution of classical TNOs was investigated using extensive simulations. We followed the evolution of more than 17000 particles with a wide range of initial conditions taking into account the perturbations from the four giant planets for 4 Gyr. The evolution of objects in the classical region is dependent on both their inclination and semimajor axes, with the inner (a<45 AU) and outer regions (a>45 AU) evolving differently. The reason is the influence of overlapping secular resonances with Uranus and Neptune (40–42 AU) and the 5:3 (a∼ ∼42.3 AU), 7:4 (a∼ ∼43.7 AU), 9:5 (a∼ ∼44.5 AU) and 11:6 (a∼ ∼ 45.0 AU) mean motion resonances strongly sculpting the inner region, while in the outer region only the 2:1 mean motion resonance (a∼ ∼47.7 AU) causes important perturbations. In particular, we found: (a) A substantial erosion of low-i bodies (i<10°) in the inner region caused by the secular resonances, except those objects that remained protected inside mean motion resonances which survived for billion of years; (b) An optimal stable region located at 45 AU<a<47 AU, q>40 AU and i>5° free of major perturbations; (c) Better defined boundaries for the classical region: 42–47.5 AU (q>38 AU) for cold classical TNOs and 40–47.5 AU (q>35 AU) for hot ones, with i=4.5° as the best threshold to distinguish between both populations; (d) The high inclination TNOs seen in the 40–42 AU region reflect their initial conditions. Therefore they should be classified as hot classical TNOs. Lastly, we report a good match between our results and observations, indicating that the former can provide explanations and predictions for the orbital structure in the classical region.     

Origin of scattered disk resonant TNOs: Evidence for an ancient excited Kuiper belt of 50 AU radius

Resonance occupation of trans-neptunian objects (TNOs) in the scattered disk (>48 AU) was investigated by integrating the orbits of 85 observed members for 4 Gyr. Twenty seven TNOs were locked in the 9:4, 16:7, 7:3, 12:5, 5:2, 8:3, 3:1, 4:1, 11:2, and 27:4 resonances. We then explored mechanisms for the origin of the resonant structure in the scattered disk, in particular the long-term 9:4, 5:2, and 8:3 resonant TNOs (median 4 Gyr), by performing large scale simulations involving Neptune scattering and planetary migration over an initially excited planetesimals disk (wide range of eccentricities and inclinations). To explain the formation of Gyr-resident populations in such distant resonances, our results suggest the existence of a primordial planetesimal disk of at least 45-50 AU radius that suffered a dynamical perturbation leading to 0.1-0.3 or greater eccentricities and a range of inclinations up to ∼20° during early stages of the Solar System history, before planetary migration.     

The scattered disk population as a source of Oort cloud comets: evaluation of its current and past role in populating the Oort cloud

We have integrated the orbits of the 76 scattered disk objects (SDOs), discovered through the end of 2002, plus 399 clones for 5 Gyr to study their dynamical evolution and the probability of falling in one of the following end states: reaching Jupiter's influence zone, hyperbolic ejection, or transfer to the Oort cloud. We find that nearly 50% of the SDOs are transferred to the Oort cloud (i.e., they reach heliocentric distances greater than 20,000 AU in a barycentric elliptical orbit), from which about 60% have their perihelia beyond Neptune's orbit (31 AU<q<36 AU) at the moment of reaching the Oort cloud. This shows that Neptune acts as a dynamical barrier, scattering most of the bodies to near-parabolic orbits before they can approach or cross Neptune's orbit in non-resonant orbits (that may allow their transfer to the planetary region as Centaurs via close encounters with Neptune). Consequently, Neptune's dynamical barrier greatly favors insertion in the Oort cloud at the expense of the other end states mentioned above. We found that the current rate of SDOs with radii R>1 km incorporated into the Oort cloud is about 5 yr⁻¹, which might be a non-negligible fraction of comet losses from the Oort cloud (probably around or even above 10%). Therefore, we conclude that the Oort cloud may have experienced and may be even experiencing a significant renovation of its population, and that the trans-neptunian belt—via the scattered disk—may be the main feeding source.     

A distant planetary-mass solar companion may have produced distant detached objects

Most known trans-neptunian objects (TNO's) are either on low eccentricity orbits or could have been perturbed to their current trajectories via gravitational interactions with known bodies. However, one or two recently-discovered TNO's are distant detached objects (DDO's) (perihelion, q>40 AU and semimajor axis, a>50 AU) whose origins are not as easily understood. We investigate the parameter space of a hypothetical distant planetary-mass solar companion which could detach the perihelion of a Neptune-dominated TNO into a DDO orbit. Perturbations of the giant planets are also included. The problem is analyzed using two models. In the first model, we start with a distribution of undetached, low-inclination TNO's having a wide range of semimajor axes. The planetary perturbations and the companion perturbation are treated in the adiabatic, secularly averaged tidal approximation. This provides a starting point for a more detailed analysis by providing insights as to the companion parameter space likely to create DDO's. The second model includes the companion and the planets and numerically integrates perturbations on a sampling that is based on the real population of scattered disk objects (SDO's). A single calculation is performed including the mutual interactions and migration of the planets. By comparing these models, we distinguish the distant detached population that can be attributable to the secular interaction from those that require additional planetary perturbations. We find that a DDO can be produced by a hypothetical Neptune-mass companion having semiminor axis, bc?2000 AU or a Jupiter-mass companion with bc?5000 AU. DDO's produced by such a companion are likely to have small inclinations to the ecliptic only if the companion's orbit is significantly inclined. We also discuss the possibility that the tilt of the planets' invariable plane relative to the solar equatorial plane has been produced by such a hypothetical distant planetary-mass companion. Perturbations of a companion on Oort cloud comets are also considered.     

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.     

Hilda asteroids among Jupiter family comets

Hilda asteroids and comets are similar from the compositional point of view. The D-taxonomic class prevailing among Hildas has all the characteristics found in cometary spectra. Jupiter Family Comets (JFCs) coming from the trans-neptunian region are under the gravitational control of Jupiter, making them a dynamically unstable population with a mean dynamical lifetime of 10⁴ to ¹⁰⁵ years. In contrast, Hilda asteroids residing in the 3:2 mean motion resonance with Jupiter are a very stable population. But once they escape from the resonance, they are dynamically controlled by Jupiter, and in this sense their behavior resembles that of JFC. We performed a numerical simulation to analyze the dynamical evolution that Hildas follow after escaping from the resonance, and their contribution to the JFC population. We found that 8% of the particles leaving the resonance end up impacting Jupiter. 98.7% of the escaped Hildas live at least 1000 years as a JFC, with a mean lifetime of 1.4×¹⁰⁶ years. In particular, escaped Hildas stay mainly in the region of perihelion distances greater than 2.5 AU. On the other hand, the number of escaped Hildas reaching the inner Solar System (q<2.5 AU) is negligible. So, there are almost no Hilda asteroids among the NEO population. We also analyzed the possibility that the Shoemaker-Levy 9 were an escaped Hilda asteroid. In this case, it would be possible to give stronger constraints to its pre-capture orbital elements.     

Trans-neptunian objects

Summary. The trans-neptunian objects (TNOs) constitute a new class of solar system object that was discovered only recently to exist beyond the orbit of Neptune. About 400 trans-neptunian objects have been detected over the past nine years and more than ten new objects are being discovered every month. All of the TNOs known to date fit into three dynamical classes: the classical, the resonant and the scattered objects. The total mass of the TNOs currently orbiting the Sun is estimated from the observed luminosity distribution to be of the order of 10–20% of the Earth's mass. However, theoretical investigations of the formation and evolution of the trans-neptunian belt into its currently observed shape suggest that it was much more massive in the past. The physical characterisation of TNOs starts to reveal some of the basic properties of these objects, such as size, shape and rotation and provides a first glance into the diversity of their surfaces. TNOs cover a very diverse range of colours, possibly reflecting different surface compositions. First evidence for the presence of water ice was found in a spectrum of one TNO while others do not show the characteristic absorption bands. The TNOs are now regarded as the likely source of some short-period comets. Owing to giant-planet and collisional perturbations, some TNOs may evolve into Centaurs, i.e. objects orbiting the Sun in the region between Jupiter and Neptune, which are further perturbed to become Jupiter-family short-period comets. Together with smaller debris generated by collisional shattering, the TNOs might represent a belt that has evolved from a more massive circumstellar disc into its present structure. Received 15 May 2001 / Published online 5 October 2001     

Are There Many Inactive Jupiter-Family Comets among the Near-Earth Asteroid Population?

We analyze the dynamical evolution of Jupiter-family (JF) comets and near-Earth asteroids (NEAs) with aphelion distances Q>3.5 AU, paying special attention to the problem of mixing of both populations, such that inactive comets may be disguised as NEAs. From numerical integrations for 2×10⁶ years we find that the half lifetime (where the lifetime is defined against hyperbolic ejection or collision with the Sun or the planets) of near-Earth JF comets (perihelion distances q<1.3 AU) is about 1.5×10⁵ years but that they spend only a small fraction of this time (∼ a few 10³ years) with q<1.3 AU. From numerical integrations for 5×10⁶ years we find that the half lifetime of NEAs in “cometary” orbits (defined as those with aphelion distances Q>4.5 AU, i.e., that approach or cross Jupiter's orbit) is 4.2×10⁵ years, i.e., about three times longer than that for near-Earth JF comets. We also analyze the problem of decoupling JF comets from Jupiter to produce Encke-type comets. To this end we simulate the dynamical evolution of the sample of observed JF comets with the inclusion of nongravitational forces. While decoupling occurs very seldom when a purely gravitational motion is considered, the action of nongravitational forces (as strong as or greater than those acting on Encke) can produce a few Enckes. Furthermore, a few JF comets are transferred to low-eccentricity orbits entirely within the main asteroid belt (Q<4 AU and q>2 AU). The population of NEAs in cometary orbits is found to be adequately replenished with NEAs of smaller Q's diffusing outward, from which we can set an upper limit of ∼20% for the putative component of deactivated JF comets needed to maintain such a population in steady state. From this analysis, the upper limit for the average time that a JF comet in near-Earth orbit can spend as a dormant, asteroid-looking body can be estimated to be about 40% of the time spent as an active comet. More likely, JF comets in near-Earth orbits will disintegrate once (or shortly after) they end their active phases.     

Higher albedos and size distribution of large transneptunian objects

Transneptunian objects (TNOs) orbit beyond Neptune and do offer important clues about the formation of our solar system. Although observations have been increasing the number of discovered TNOs and improving their orbital elements, very little is known about elementary physical properties such as sizes, albedos and compositions. Due to TNOs large distances (>40 AU) and observational limitations, reliable physical information can be obtained only from brighter objects (supposedly larger bodies). According to size and albedo measurements available, it is evident the traditionally assumed albedo p=0.04 cannot hold for all TNOs, especially those with approximately absolute magnitudes H?5.5. That is, the largest TNOs possess higher albedos (generally >0.04) that strongly appear to increase as a function of size. Using a compilation of published data, we derived empirical relations which can provide estimations of diameters and albedos as a function of absolute magnitude. Calculations result in more accurate size/albedo estimations for TNOs with H?5.5 than just assuming p=0.04. Nevertheless, considering low statistics, the value p=0.04 sounds still convenient for H>5.5 non-binary TNOs as a group. We also discuss about physical processes (e.g., collisions, intrinsic activity and the presence of tenuous atmospheres) responsible for the increase of albedo among large bodies. Currently, all big TNOs (>700 km) would be capable to sustain thin atmospheres or icy frosts composed of CH₄, CO or N₂ even for body bulk densities as low as 0.5 g cm⁻³. A size-dependent albedo has important consequences for the TNOs size distribution, cumulative luminosity function and total mass estimations. According to our analysis, the latter can be reduced up to 50% if higher albedos are common among large bodies.Lastly, by analyzing orbital properties of classical TNOs (), we confirm that cold and hot classical TNOs have different concentration of large bodies. For both populations, distinct absolute magnitude distributions are maximized for an inclination threshold equal to 4.5° at >99.63% confidence level. Furthermore, more massive classical bodies are anomalously present at , a result statistically significant and apparently not caused by observational biases. This feature would provide a new constraint for transneptunian belt formation models.     

Resonance sticking in the scattered disk

We investigate the dynamical evolution of trans-neptunian objects (TNOs) in typical scattered disk orbits (scattered TNOs) by performing simulations using several thousand particles lying initially on Neptune-encountering orbits. We explore the role of resonance sticking in the scattered disk, a phenomenon characterized by multiple temporary resonance captures (‘resonances’ refers to external mean motion resonances with Neptune, which can be described in the form r:s, where the arguments r and s are integers). First, all scattered TNOs evolve through intermittent temporary resonance capture events and gravitational scattering by Neptune. Each scattered TNO experiences tens to hundreds of resonance captures over a period of 4 Gyr, which represents about 38% of the object's lifetime (mean value). Second, resonance sticking plays an important role at semimajor axes , where the great majority of such captures occurred. It is noteworthy that the stickiest (i.e., dominant) resonances in the scattered disk are located within this distance range and are those possessing the lowest argument s. This was evinced by r:1, r:2 and r:3 resonances, which played the greatest role during resonance sticking evolution, often leading to captures in several of their neighboring resonances. Finally, the timescales and likelihood of temporary resonance captures are roughly proportional to resonance strength. The dominance of low s resonances is also related to the latter. In sum, resonance sticking has an important impact on the evolution of scattered TNOs, contributing significantly to the longevity of these objects.     

On the dynamic connection of comets with Uranus

We consider the connection with Uranus for: (1) 945 near-parabolic comets (the period P > 200 years, the perihelion distance q > 0.1 AU), (2) 1277 Kreutz comets (P > 200 years, q < 0.01 AU), and (3) 414 short-period comets (P < 200 years). It turns out that none of near-parabolic comets passed through Uranus’s activity sphere, none of the Kreutz comets approach Uranus closer than 11 AU, and only two short-period comets, C/2006 U7 and C/2006 F2, could have a close approach to Uranus during 5000 years.     

Evidence for an Extended Scattered Disk

By telescopic tracking, we have established that the transneptunian object (TNO) 2000 CR₁₀₅ has a semimajor axis of 220±1 AU and perihelion distance of 44.14±0.02 AU, beyond the domain which has heretofore been associated with the “scattered disk” of Kuiper Belt objects interacting via gravitational encounters with Neptune. We have also firmly established that the TNO 1995 TL₈ has a high perihelion (of 40.08±0.02 AU). These objects, and two other recent discoveries which appear to have perihelia outside 40 AU, have probably been placed on these orbits by a gravitational interaction which is not strong gravitational scattering off of any of the giant planets on their current orbits. Their existence may thus have profound cosmogonic implications for our understanding of the formation of the outer Solar System. We discuss some viable scenarios which could have produced these objects, including long-term diffusive chaos and scattering off of other massive bodies in the outer Solar System. This discovery implies that there must be a large population of TNOs in an “extended scattered disk” with perihelia above the previously suggested 38 AU boundary. The total population is difficult to estimate due to the ease with which such objects would have been lost. This illustrates the great value of frequent and well time-sampled recovery observations of trans-neptunian objects within their discovery opposition.     

V-type asteroids in the middle main belt

V-type asteroids are bodies whose surfaces are constituted of basalt. In the Main Asteroid Belt, most of these asteroids are assumed to come from the basaltic crust of Asteroid (4) Vesta. This idea is mainly supported by (i) the fact that almost all the known V-type asteroids are in the same region of the belt as (4) Vesta, i.e., the inner belt (semi-major axis 2.1<a<2.5 AU), (ii) the existence of a dynamical asteroid family associated to (4) Vesta, and (iii) the observational evidence of at least one large craterization event on Vesta's surface. One V-type asteroid that is difficult to fit in this scenario is (1459) Magnya, located in the outer asteroid belt, i.e., too far away from (4) Vesta as to have a real possibility of coming from it. The recent discovery of the first V-type asteroid in the middle belt (2.5<a<2.8 AU), (21238) 1995WV7 [Binzel, R.P., Masi, G., Foglia, S., 2006. Bull. Am. Astron. Soc. 38, 627; Hammergren, M., Gyuk, G., Puckett, A., 2006. ArXiv e-print, astro-ph/0609420], located at ∼2.54 AU, raises the question of whether it came from (4) Vesta or not. In this paper, we present spectroscopic observations indicating the existence of another V-type asteroid at ∼2.53 AU, (40521) 1999RL95, and we investigate the possibility that these two asteroids evolved from the Vesta family to their present orbits by a semi-major axis drift due to the Yarkovsky effect. The main problem with this scenario is that the asteroids need to cross the 3/1 mean motion resonance with Jupiter, which is highly unstable. Combining N-body numerical simulations of the orbital evolution, that include the Yarkovsky effect, with Monte Carlo models, we compute the probability that an asteroid of a given diameter D evolves from the Vesta family and crosses over the 3/1 resonance, reaching a stable orbit in the middle belt. Our results indicate that an asteroid like (21238) 1995WV7 has a low probability (∼1%) of having evolved through this mechanism due to its large size (D∼5 km), because the Yarkovsky effect is not sufficiently efficient for such large asteroids. However, the mechanism might explain the orbits of smaller bodies like (40521) 1999RL95 (D∼3 km) with ∼70-100% probability, provided that we assume that the Vesta family formed ?3.5 Gy ago. We estimate the debiased population of V-type asteroids that might exist in the same region as (21238) and (40521) (2.5<a?2.62 AU) and conclude that about 10 to 30% of the V-type bodies with D>1 km may come from the Vesta family by crossing over the 3/1 resonance. The remaining 70-90% must have a different origin.     

Injection of Oort Cloud comets: the fundamental role of stellar perturbations

We present Monte Carlo simulations of the dynamical evolution of the Oort cloud over the age of the Solar System, using an initial sample of one million test comets without any cloning. Our model includes perturbations due to the Galactic tide (radial and vertical) and passing stars. We present the first detailed analysis of the injection mechanism into observable orbits by comparing the complete model with separate models for tidal and stellar perturbations alone. We find that a fundamental role for injecting comets from the region outside the loss cone (perihelion distance q > 15 AU) into observable orbits (q < 5 AU) is played by stellar perturbations. These act in synergy with the tide such that the total injection rate is significantly larger than the sum of the two separate rates. This synergy is as important during comet showers as during quiescent periods and concerns comets with both small and large semi-major axes. We propose different dynamical mechanisms to explain the synergies in the inner and outer parts of the Oort Cloud. We find that the filling of the observable part of the loss cone under normal conditions in the present-day Solar System rises from <1% for a < 20 000 AU to about 100% for a ? 100 000 AU.     

HST Photometry of trans-Neptunian Objects

From July 2001 to June 2002, an HST snapshot program obtained V, R and I photometry for 72 TNOs. The TNOs were sorted by dynamical class, and Spearman rank correlation statistics were calculated for each combination of color and orbital parameter. No strong correlations were found for the combined sample of TNOs, the resonant TNOs, or the non-resonant TNOs (classical). The results presented here suggest that if correlations reported by other authors are real, they are evident only at shorter wavelengths than observed in our survey.     

Resonant motion of trans-Neptunian objects in high-eccentricity orbits

A substantial fraction of the Edgeworth-Kuiper belt objects are presently known to move in resonance with Neptune (the principal commensurabilities are 1/2, 3/5, 2/3, and 3/4). We have found that many of the distant (with orbital semimajor axes a > 50 AU) trans-Neptunian objects (TNOs) also execute resonant motions. Our investigation is based on symplectic integrations of the equations of motion for all multiple-opposition TNOs with a > 50 AU with allowance made for the uncertainties in their initial orbits. Librations near such commensurabilities with Neptune as 4/9, 3/7, 5/12, 2/5, 3/8, 4/27, and others have been found. The largest number of distant TNOs move near the 2/5 resonance with Neptune: 12 objects librate with a probability higher than 0.75. The multiplicity of objects moving in 2/5 resonance and the longterm stability of their librations suggest that this group of resonant objects was formed at early formation stages of the Solar system. For most of the other resonant objects, the librations are temporary. We also show the importance of asymmetric resonances in the large changes in TNO perihelion distances.     

Secular light curve of Comet 9P/Tempel 1

In support of the Deep Impact Mission, we have updated the secular light curve of 9P/Tempel 1 presented in Paper I [Ferrín, I., 2005. Icarus 178, 493-516], with new data sets. The secular light curves (SLC) of the comet are presented in the log and time plots (Figs. 1 and 2) and provide a clear profile of the overall shape of the envelope. We arrive at the following conclusions: (1) Improved values of 18 photometric parameters are derived including the turn on and turn off points, RON=−3.47±0.05 AU, ROFF=+4.20±0.05 AU, and TON=−410±25 d, TOFF=+555±25 d. (2) The improved SLC shows a most interesting and peculiar shape, with a linear power law of slope n=7.7±0.1 from RON=−3.47 AU to RBP=−2.08±0.05 AU, and then converts to a law with curvature. The break point of the power law at RBP=−2.08 AU, mV(1,R)=14.0±0.1 mag, is interpreted as a change in sublimating something more volatile than water ice (most probably CO₂), to water ice sublimation. In other words, the comet's sublimation is controlled by two different substances. (3) The photometric-age (defined in Paper I) and the time-age of the comet [Ferrín, I., 2006. Icarus. In press] are recomputed, and results in a value P-AGE=21±2 and T-AGE=11±2 comet years. Thus 9P is a young comet. (4) The comet is active almost up to aphelion since the turn off point has been determined at ROFF=+4.20±0.05 AU while aphelion takes place at Q=+4.74 AU. (5) The comet exhibits activity post-aphelion which is not understood. Two hypothesis are advanced to explain this behavior.     

Considerations on the magnitude distributions of the Kuiper belt and of the Jupiter Trojans

By examining the absolute magnitude (H) distributions (hereafter HD) of the cold and hot populations in the Kuiper belt and of the Trojans of Jupiter, we find evidence that the Trojans have been captured from the outer part of the primordial trans-neptunian planetesimal disk. We develop a sketch model of the HDs in the inner and outer parts of the disk that is consistent with the observed distributions and with the dynamical evolution scenario known as the ‘Nice model’. This leads us to predict that the HD of the hot population should have the same slope of the HD of the cold population for 6.5<H<9, both as steep as the slope of the Trojans' HD. Current data partially support this prediction, but future observations are needed to clarify this issue. Because the HD of the Trojans rolls over at H∼9 to a collisional equilibrium slope that should have been acquired when the Trojans were still embedded in the primordial trans-neptunian disk, our model implies that the same roll-over should characterize the HDs of the Kuiper belt populations, in agreement with the results of Bernstein et al. [Bernstein, G.M., and 5 colleagues, 2004. Astron. J. 128, 1364-1390] and Fuentes and Holman [Fuentes, C.I., Holman, M.J., 2008. Astron. J. 136, 83-97]. Finally, we show that the constraint on the total mass of the primordial trans-neptunian disk imposed by the Nice model implies that it is unlikely that the cold population formed beyond 35 AU.     

The population of faint Jupiter family comets near the Earth

We study the population of faint Jupiter family comets (JFCs) that approach the Earth (perihelion distances q<1.3 AU) by applying a debiasing technique to the observed sample. We found for the debiased cumulative luminosity function (CLF) of absolute total magnitudes H₁₀ a bimodal distribution in which brighter comets (H₁₀?9) follow a linear relation with a steep slope α=0.65±0.14, while fainter comets follow a much shallower slope α=0.25±0.06 down to H₁₀∼18. The slope can be pushed up to α=0.35±0.09 if a second break in the H₁₀ distribution to a much shallower slope is introduced at H₁₀∼16. We estimate a population of about 10³ faint JFCs with q<1.3 AU and 10<H₁₀<15 (radii ∼0.1-0.5 km). The shallowness of the CLF for faint near-Earth JFCs may be explained either as: (i) the source population (the scattered disk) has an equally very shallow distribution in the considered size range, or (ii) the distribution is flattened by the disintegration of small objects before that they have a chance of being observed. The fact that the slope of the magnitude distribution of the faint active JFCs is very similar to that found for a sample of dormant JFCs candidates suggests that for a surviving (i.e., not disintegrated) object, the probability of becoming dormant versus keeping some activity is roughly size independent.     

Recent orbital evolution and the internal structures of Enceladus and Dione

We study the orbital behavior of Saturn’s satellites Enceladus and Dione during their passage through the 2:1 mean-motion resonances to constrain their interior structures, parameterized by the quantity k₂/Q (assumed constant). Enceladus’ evolution after escape from the second-order e-Enceladus e-Dione resonance requires that (k₂/Q)_Enceladus<8×10^-4, for that Q_Saturn>18,000. This result is in agreement with [Meyer, J., Wisdom, J., 2008b. Icarus 193, 213-223]. The present-day libration amplitude of Enceladus requires that (k₂/Q)_Enceladus>1.2×10^-4, assuming that Q_Saturn<10⁵. Dione’s present-day eccentricity indicates that (k₂/Q)_Dione?3×10^-4 for Q_Saturn>18,000. Assuming Maxwellian viscoelastic behavior, we find that for Enceladus a convective ice shell overlying an ocean is too dissipative to match the orbital constraints. We conclude that a conductive shell overlying an ocean is more likely, and discuss the implications of this result. Dione’s ice shell is also likely to be conductive, but our results are less constraining.