Micrometeorite impact annealing of ice in the outer Solar System

The spectra of water ice on the surfaces of icy satellites and Kuiper Belt Objects (KBOs) indicate that the surface ice on these bodies is in a crystalline state. This conflicts with theoretical models, which predict that radiation (galactic cosmic rays and solar ultraviolet) should damage the crystalline structure of ice on geologically short timescales. Temperatures are too low in the outer Solar System for the ice to anneal, and reflectance spectra of these bodies should match those of amorphous solid water (ASW). We assess whether the kinetic energy deposited as heat by micrometeorite impacts on outer Solar System bodies is sufficient to anneal their surface ice down to a near-infrared optical depth . We calculate the kinetic energy flux from interplanetary micrometeorite impacts, including gravitational focusing. We also calculate the thermal diffusion of impact heat in various surfaces and the rate of annealing of ice. We conclude that the rate of annealing from micrometeorite impacts is sufficient to explain the crystallinity of ice on nearly all the surfaces of the saturnian, uranian and neptunian satellites. We discuss how the model can be used in conjunction with spectra of KBOs to probe dust fluxes in the Kuiper Belt.     

Colors and collision rates within the Kuiper belt: problems with the collisional resurfacing scenario

We present a numerical check of the collisional resurfacing (CR) hypothesis proposed to explain the observed color diversity within the Kuiper Belt (where surface reddening due to space weathering is counteracted by regular resurfacing of neutral material after mutual collisions). Deterministic simulations are performed in order to estimate the relative spatial distribution of kinetic energy received by collisions, , for a population of target Kuiper Belt objects (KBOs) embedded in a swarm of impactors distributed within the belt. Four different impactor disks have been considered, depending on the excitation and the external limit of the belt and the density of the scattered KBOs (SKBOs) population. The obtained results are compared to the relative color index distribution within the observed Kuiper Belt, in order to derive possible similarities between the high vs low objects spatial distribution in our simulations and the bluer vs redder KBOs distribution in the “real” Kuiper Belt. Such similarities are found for several important features, in particular the general correlations between highly impacted objects and high rms excitation and low perihelion q values that are in good agreement with equivalent correlations found for the bluest objects of the observed belt. Nevertheless, simulations disagree with observations on two crucial points. (1) The plutinos are significantly more collisionally affected than the rest of our test KBO population, whereas there is no observed tendency toward bluer plutinos. (2) There is always a much stronger correlation between and eccentricities than inclinations, whereas observations show just the opposite feature. The presence of numerous SKBO impactors could significantly damp these problematic features, but cannot erase them. Whether these contradictions invalidate the whole CR scenario or not remains yet uncertain, since the physical processes at play are still far from being fully understood and the sample of available observational data is still relatively limited. But it seems nevertheless that the scenario might not hold in its simple present form.     

Fragmentation model dependence of collision cascades

Mass depletion of bodies through successive collisional disruptions (i.e., collision cascade) is one of the most important processes in the studies of the asteroids belt, the Edgeworth-Kuiper belt, debris disks, and planetary formation. The collisional disruption is divided into two types, i.e., catastrophic disruption and cratering. Although some studies of the collision cascades neglected the effect of cratering, it is unclear which type of disruption makes a dominant contribution to the collision cascades. In the present study, we construct a simple outcome model describing both catastrophic disruption and cratering, which has some parameters characterizing the total ejecta mass, the mass of the largest fragment, and the power-law exponent of the size distribution of fragments. Using this simple outcome model with parameters, we examine the model dependence of the mass depletion time in collision cascades for neglect of coalescence of colliding bodies due to high collisional velocities. We find the cratering collisions are much more effective in collision cascades than collisions with catastrophic disruption in a wide region of the model parameters. It is also found that the mass depletion time in collision cascades is mainly governed by the total ejecta mass and almost insensitive to the mass of the largest fragment and the power-law exponent of fragments for a realistic parameter region. The total ejecta mass is usually determined by the ratio of the impact energy divided by the target mass (i.e. Q-value) to its threshold value for catastrophic disruption, as well as in our simple model. We derive a mass depletion time in collision cascades, which is determined by of the high-mass end of collision cascades. The mass depletion time derived with our model would be applicable to debris disks and planetary formation.     

The Kozai resonance in the outer solar system and the dynamics of long-period comets

We consider the secular evolution of the orbits of bodies in the Outer Solar System under the perturbations of the jovian planets assumed on coplanar and circular orbits. Through the approach used for asteroidal belt by Yoshihide Kozai in 1962, we obtain that the Kozai resonance do not affect the behavior of bodies belonging to the Kuiper belt but concerns the long-timescale evolution of long-period comets. In particular this resonance appears as a process contributing to produce Sun-grazer comets.     

Mass of the Kuiper belt

The Kuiper belt includes tens of thousand of large bodies and millions of smaller objects. The main part of the belt objects is located in the annular zone between 39.4 and 47.8 au from the Sun; the boundaries correspond to the average distances for orbital resonances 3:2 and 2:1 with the motion of Neptune. One-dimensional, two-dimensional, and discrete rings to model the total gravitational attraction of numerous belt objects are considered. The discrete rotating model most correctly reflects the real interaction of bodies in the Solar system. The masses of the model rings were determined within EPM2017—the new version of ephemerides of planets and the Moon at IAA RAS—by fitting spacecraft ranging observations. The total mass of the Kuiper belt was calculated as the sum of the masses of the 31 largest trans-Neptunian objects directly included in the simultaneous integration and the estimated mass of the model of the discrete ring of TNO. The total mass is \((1.97 \pm 0.35)\times 10^{-2} \ m_{\oplus }\). The gravitational influence of the Kuiper belt on Jupiter, Saturn, Uranus, and Neptune exceeds at times the attraction of the hypothetical 9th planet with a mass of \(\sim 10 \ m_{\oplus }\) at the distances assumed for it. It is necessary to take into account the gravitational influence of the Kuiper belt when processing observations and only then to investigate residual discrepancies to discover a possible influence of a distant large planet.     

Numerical simulations of impacts involving porous bodies: II. Comparison with laboratory experiments

In this paper, we compare the outcome of high-velocity impact experiments on porous targets, composed of pumice, with the results of simulations by a 3D SPH hydrocode in which a porosity model has been implemented. The different populations of small bodies of our Solar System are believed to be composed, at least partially, of objects with a high degree of porosity. To describe the fragmentation of such porous objects, a different model is needed than that used for non-porous bodies. In the case of porous bodies, the impact process is not only driven by the presence of cracks which propagate when a stress threshold is reached, it is also influenced by the crushing of pores and compaction. Such processes can greatly affect the whole body's response to an impact. Therefore, another physical model is necessary to improve our understanding of the collisional process involving porous bodies. Such a model has been developed recently and introduced successfully in a 3D SPH hydrocode [Jutzi, M., Benz, W., Michel, P., 2008. Icarus 198, 242-255]. Basic tests have been performed which already showed that it is implemented in a consistent way and that theoretical solutions are well reproduced. However, its full validation requires that it is also capable of reproducing the results of real laboratory impact experiments. Here we present simulations of laboratory experiments on pumice targets for which several of the main material properties have been measured. We show that using the measured material properties and keeping the remaining free parameters fixed, our numerical model is able to reproduce the outcome of these experiments carried out under different impact conditions. This first complete validation of our model, which will be tested for other porous materials in the future, allows us to start addressing problems at larger scale related to small bodies of our Solar System, such as collisions in the Kuiper Belt or the formation of a family by the disruption of a porous parent body in the main asteroid belt.     

Coupling dynamical and collisional evolution of small bodies: II. Forming the Kuiper belt, the Scattered Disk and the Oort Cloud

The Oort Cloud, the Kuiper belt and the Scattered Disk are dynamically distinct populations of small bodies evolving in the outer regions of the Solar System. Whereas their collisional activity is now quiet, gravitational interactions with giant planets may have shaped these populations both dynamically and collisionally during their formation. Using a hybrid approach [Charnoz, S., Morbidelli, A., 2003. Icarus 166, 141-166], the present paper tries to couple the primordial collisional and dynamical evolution of these three populations in a self-consistent way. A critical parameter is the primordial size-distribution. We show that the initial planetesimal size distribution that allows an effective mass depletion of the Kuiper belt by collisional grinding, would decimate also the population of comet-size bodies that end in the Oort Cloud and, in particular, in the Scattered Disk. As a consequence, the Oort Cloud and the Scattered Disk would be too anemic, by a factor 20 to 100, relative to the estimates achieved from the observation of the fluxes of long period and Jupiter family comets, respectively. For these two reservoirs to have a sufficient number of comets, the initial size distribution in the planetesimal disk had to be such that the mass depletion by collisional erosion of the Kuiper belt was negligible. Consequently the current mass deficit of the Kuiper belt needs to be explained by dynamical mechanisms.     

Low-speed impacts between rubble piles modeled as collections of polyhedra

We present results of modeling rubble piles as collections of polyhedra. The use of polyhedra allows more realistic (irregular) shapes and interactions (e.g. collisions), particularly for objects of different sizes. Rotational degrees of freedom are included in the modeling, which may be important components of the motion. We solved the equations of rigid-body dynamics, including frictional/inelastic collisions, for collections of up to several hundred elements. As a demonstration of the methods and to compare with previous work by other researchers, we simulated low-speed collisions between km-scale bodies with the same general parameters as those simulated by Leinhardt et al. [Leinhardt, Z.M., Richardson, D.C., Quinn, T., 2000. Icarus 146, 133-151]. High-speed collisions appropriate to present-day asteroid encounters require additional treatment of shock effects and fragmentation and are the subject of future work; here we study regimes appropriate to planetesimal accretion and re-accretion in the aftermath of catastrophic events. Collisions between equal-mass objects at low speeds () were simulated for both head-on and off-center collisions between rubble piles made of a power-law mass spectrum of sub-elements. Very low-speed head-on collisions produce single objects from the coalescence of the impactors. For slightly higher speeds, extensive disruption occurs, but re-accretion produces a single object with most of the total mass. For increasingly higher speeds, the re-accreted object has smaller mass, finally resulting in complete catastrophic disruption with all sub-elements on escape trajectories and only small amounts of mass in re-accreted bodies. Off-center collisions at moderately low speeds produce two re-accreted objects of approximately equal mass, separating at greater than escape speed. At high speed, complete disruption occurs as with the high-speed head-on collisions. Head-on collisions at low to moderate speeds result in objects of mostly oblate shape, while higher speed collisions produce mostly prolate objects, as do off-center collisions at moderate and high speeds. Collisions carried out with the same dissipative coefficients (coefficient of restitution ?n=0.8, zero friction) as used by Leinhardt et al. [Leinhardt, Z.M., Richardson, D.C., Quinn, T., 2000. Icarus 146, 133-151] result in a value for specific energy for disruption , somewhat lower than the value of 2 J/kg found by them, while collisions with a lower coefficient of restitution and friction [?n=0.5, ?t=0, μ=0.5, similar to those used by Michel, et al. [Michel, P., Benz, W., Richardson, D.C., 2004. Planet. Space Sci. 52, 1109-1117] for SPH + N-body calculations] yield .     

A derivation of the luminosity function of the Kuiper belt from a broken power-law size distribution

We have derived a model of the Kuiper belt luminosity function exhibited by a broken power-law size distribution. This model allows direct comparison of the observed luminosity function to the underlying size distribution. We discuss the importance of the radial distribution model in determining the break diameter. We determine a best-fit break-diameter of the Kuiper belt size-distribution of 30<Db<90 km via a maximum-likelihood fit of our model to the observed luminosity function. We also confirm that the observed luminosity function for m(R)∼21-28 is consistent with a broken power-law size distribution, and exhibits a break at .     

Masses of the Main Asteroid Belt and the Kuiper Belt from the Motions of Planets and Spacecraft

Dynamicalmass estimates for the main asteroid belt and the trans-Neptunian Kuiper belt have been found from their gravitational influence on the motion of planets. Discrete rotating models consisting ofmovingmaterial points have been used tomodel the total attraction fromsmall or as yet undetected bodies of the belts. The masses of the model belts have been included in the set of parameters being refined and determined and have been obtained by processing more than 800 thousand modern positional observations of planets and spacecraft. We have processed the observations and determined the parameters based on the new EPM2017 version of the IAA RAS planetary ephemerides. The large observed radial extent of the belts (more than 1.2 AU for the main belt and more than 8 AU for the Kuiper belt) and the concentration of bodies in the Kuiper belt at a distance of about 44 AU found from observations have been taken into account in the discrete models. We have also used individual mass estimates for large bodies of the belts as well as for objects that spacecraft have approached and for bodies with satellites. Our mass estimate for the main asteroid belt is (4.008 ± 0.029) × 10^?4/m_⊕ (3σ). The bulk of the Kuiper belt objects are in the ring zone from 39.4 to 47.8 AU. The estimate of its total mass together with the mass of the 31 largest trans-Neptunian Kuiper belt objects is (1.97 ± 0.30) × 10^?2m_⊕ (3σ), which exceeds the mass of the main asteroid belt almost by a factor of 50. The mass of the 31 largest trans-Neptunian objects (TNOs) is only about 40% of the total one.     

Color diversity among Kuiper belt objects: The collisional resurfacing model revisited

A re-evaluation of the collisional resurfacing model based on up-to-date Kuiper belt objects size distribution and a more precise treatment of the cosmic-ray environment at the outer Solar System is presented. The result of the irradiation due to cosmic-rays with different energies altered in a different way the material of the objects, producing, under certain conditions, a thick irradiation mantle. Since the collisional resurfacing model is based in the competition between darkening by cosmic-rays and resurfacing due to impacts, the color of objects in different regions of the belt could vary if the projectile populations in those regions are truncated at a different radius.     

Asteroids were born big

How big were the first planetesimals? We attempt to answer this question by conducting coagulation simulations in which the planetesimals grow by mutual collisions and form larger bodies and planetary embryos. The size frequency distribution (SFD) of the initial planetesimals is considered a free parameter in these simulations, and we search for the one that produces at the end objects with a SFD that is consistent with Asteroid belt constraints. We find that, if the initial planetesimals were small (e.g. km-sized), the final SFD fails to fulfill these constraints. In particular, reproducing the bump observed at diameter in the current SFD of the asteroids requires that the minimal size of the initial planetesimals was also ∼100 km. This supports the idea that planetesimals formed big, namely that the size of solids in the proto-planetary disk “jumped” from sub-meter scale to multi-kilometer scale, without passing through intermediate values. Moreover, we find evidence that the initial planetesimals had to have sizes ranging from 100 to several 100 km, probably even 1000 km, and that their SFD had to have a slope over this interval that was similar to the one characterizing the current asteroids in the same size range. This result sets a new constraint on planetesimal formation models and opens new perspectives for the investigation of the collisional evolution in the Asteroid and Kuiper belts as well as of the accretion of the cores of the giant planets.     

Numerical modelling of heating in porous planetesimal collisions

Collisions between planetesimals at speeds of several kilometres per second were common during the early evolution of our Solar System. However, the collateral effects of these collisions are not well understood. In this paper, we quantify the efficiency of heating during high-velocity collisions between planetesimals using hydrocode modelling. We conducted a series of simulations to test the effect on shock heating of the initial porosity and temperature of the planetesimals, the relative velocity of the collision and the relative size of the two colliding bodies. Our results show that while heating is minor in collisions between non-porous planetesimals at impact velocities below 10 km s⁻¹, in agreement with previous work, much higher temperatures are reached in collisions between porous planetesimals. For example, collisions between nearly equal-sized, porous planetesimals can melt all, or nearly all, of the mass of the bodies at collision velocities below 7 km s⁻¹. For collisions of small bodies into larger ones, such as those with an impactor-to-target mass ratio below 0.1, significant localised heating occurs in the target body. At impact velocities as low as 5 km s⁻¹, the mass of melt will be nearly double the mass of the impactor, and the mass of material shock heated by 100 K will be nearly 10 times the mass of the impactor. We present a first-order estimate of the cumulative effects of impact heating on a porous planetesimal parent body by simulating the impact of a population of small bodies until a disruptive event occurs. Before disruption, impact heating is volumetrically minor and highly localised; in no case was more than about 3% of the parent body heated by more than 100 K. However, heating during the final disruptive collision can be significant; in about 10% of cases, almost all of the parent body is heated to 700 K (from an initial temperature of ∼300 K) and more than a tenth of the parent body mass is melted. Hence, energetic collisions between planetesimals could have had important effects on the thermal evolution of primitive materials in the early Solar System.     

Colours and Composition of the Centaurs

Centaurs are widely believed to come from the Edgeworth-Kuiper belt, located beyond the orbit of Neptune. From here they can be injected into the inner part of the Solar System through planetary perturbations or mutual collisions. Due to their origin and dynamical evolution, Centaurs are supposed to constitute a transition population of objects from the large reservoir of Trans-Neptunian Objects (TNOs) to the active bodies of the inner Solar System. On the basis of the present knowledge of the physical properties of Centaurs and TNOs a similarity between the two populations appears evident. This is the strongest observational constraint supporting the theory of common origin.     

Fugitives from the Vesta family

Spectroscopic observations of Asteroid (4) Vesta and numerous members of the Vesta family located in the inner asteroid belt have determined that these objects have reflectance properties of basaltic material. A plausible hypothesis is that the surface of Vesta was punctured by large impacts in the past which dispersed fragments of its basaltic crust into space and produced one of the most prominent asteroid families ever created in the belt. Until recently, Vesta was the only known object in the asteroid belt which underwent differentiation and survived to the present epoch. Since 2000, many new small basaltic asteroids have been discovered in the inner and outer parts of the asteroid belt, possibly representing fragments from distinct differentiated bodies. These discoveries may help us to better understand the number and nature of objects in the inner Solar System that underwent geological differentiation. To investigate these issues we performed extensive numerical simulations whose aim was to reproduce, as precisely as possible, the dynamical evolution of Vesta's ejected fragments over timescales comparable to the family's age. Specifically, we numerically integrated the orbital evolution of 6600 test bodies with orbits that started within the Vesta family and dynamically evolved over 2 Gy. Our model included gravitational perturbation of all planets (except Mercury) and the Yarkovsky effect. The results show that a relatively large fraction of the original Vesta family members may have evolved out of the family borders defined by clustering algorithms and are now dispersed over the inner asteroid belt. We compared the orbital distribution of our model fragments with the orbital locations of known basaltic asteroids in various parts of the inner main belt to find that: (i) Most basaltic asteroids with semimajor axis located outside the Vesta family's borders in the inner main belt, including (809) Lundia and (956) Elisa, are most likely fugitives from the Vesta family that have evolved to their current orbits via various identified dynamical pathways. Our results also suggest that the Vesta family is at least ∼1 Gy old. (ii) Interestingly, orbits of many basaltic asteroids with , like those of (4796) Lewis and (5379) Abehiroshi, are displaced from the Vesta family to low inclinations and are not obtained in our simulations with sufficient efficiency. We propose that: (i) these small basaltic asteroids may be fragments of differentiated bodies other than (4) Vesta; or (ii) they were liberated from the Vesta's surface before (or during) the Late Heavy Bombardment epoch ∼3.8 Gy ago and their orbital inclinations separated from that of Vesta when secular resonances swept through the region.     

Asteroids with Satellites: Inventory, Properties, and Prospects for Future Discoveries

The population of binary asteroids numbers over 160 systems, and they can be found amongst near-Earth asteroids (NEAs), Main-Belt asteroids (MBAs), Jupiter Trojans, Centaurs and trans-Neptunian objects (TNOs). The discoveries have been made with space missions, radar observations, photometric lightcurves, and high resolution imaging from the ground and space. The properties of each population are widely different due to varying formation mechanisms and discovery techniques for each group. Future large-aperture telescopes will be capable of imaging both components for nearly all known systems and will drastically improve prospects for discovery of smaller and more tightly bound systems throughout the Solar System. The study of binary asteroids has provided valuable estimates on asteroid density and structure, a better understanding of the radiative YORP-effect, insights on catastrophic collisions, and may prove to be a key diagnostic for understanding the formation and evolution of the Kuiper Belt population.     

The orbit, mass, and albedo of transneptunian binary (66652) 1999 RZ253

We have observed (66652) 1999 RZ₂₅₃ with the Hubble Space Telescope at seven separate epochs and have fit an orbit to the observed relative positions of this binary. Two orbital solutions have been identified that differ primarily in the inclination of the orbit plane. The best fit corresponds to an orbital period, days, semimajor axis a=4660±170 km and orbital eccentricity e=0.460±0.013 corresponding to a system mass m=3.7±0.4×¹⁰¹⁸ kg. For a density of the albedo at 477 nm is p₄₇₇=0.12±0.01, significantly higher than has been commonly assumed for objects in the Kuiper belt. Multicolor, multiepoch photometry shows this pair to have colors typical for the Kuiper belt with a spectral gradient of 0.35 per 100 nm in the range between 475 and 775 nm. Photometric variations at the four epochs we observed were as large as 12±3% but the sampling is insufficient to confirm the existence of a lightcurve.     

Resolution of the kuiper belt object color controversy: two distinct color populations

Four years ago, we reported that the surface colors of ancient, icy bodies at and beyond the orbit of Neptune—Kuiper belt objects— divide into two distinct color populations. Our report has proven quite controversial. Specifically, every other research group looking with large telescopes at Kuiper belt objects finds a continuous range of colors rather than two distinct populations. Here we report new color measurements of 18 objects, primarily from the Keck I 10-m telescope, that confirm the existence of two populations. We have combined the color measurements of the other groups to create a data set comparable in size to our data set. We have carried out a Monte Carlo statistical analysis and found that both data sets are consistent with two color populations and our data set, which has smaller uncertainties, rules out a continuum of colors. In addition, our new observations and those in the literature confirm our earlier report that classical KBOs with perihelion distances beyond 40 AU exhibit extremely red surface colors. Our results rule out a continuous color distribution for both our complete sample and subsamples with perihelion distances greater than or less than 40 AU. We suspect the color patterns will result in a better understanding of the formation and evolution of the outer Solar System.     

Embedded star clusters and the formation of the Oort Cloud

Observations suggest most stars originate in clusters embedded in giant molecular clouds [Lada, C.J., Lada, E.A., 2003. Annu. Rev. Astron. Astrophys. 41, 57-115]. Our Solar System likely spent 1-5 Myrs in such regions just after it formed. Thus the Oort Cloud (OC) possibly retains evidence of the Sun's early dynamical history and of the stellar and tidal influence of the cluster. Indeed, the newly found objects (90377) Sedna and 2000 CR₁₀₅ may have been put on their present orbits by such processes [Morbidelli, A., Levison, H.F., 2004. Astron. J. 128, 2564-2576]. Results are presented here of numerical simulations of the orbital evolution of comets subject to the influence of the Sun, Jupiter and Saturn (with their current masses on orbits appropriate to the period before the Late Heavy Bombardment (LHB) [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461]), passing stars and tidal force associated with the gas and stars of an embedded star cluster. The cluster was taken to be a Plummer model with 200-400 stars, with a range of initial central densities. The Sun's orbit was integrated in the cluster potential together with Jupiter and Saturn and the test particles. Stellar encounters were incorporated by directly integrating the effects of stars passing within a sphere centred on the Sun of radius equal to the Plummer radius for low-density clusters and half a Plummer radius for high-density clusters. The gravitational influence of the gas was modeled using the tidal force of the cluster potential. For a given solar orbit, the mean density, 〈ρ〉, was computed by orbit-averaging the density of material encountered. This parameter proved to be a good measure for predicting the properties of the OC. On average 2-18% of our initial sample of comets end up in the OC after 1-3 Myr. A comet is defined to be part of the OC if it is bound and has q>35 AU. Our models show that the median distance of an object in the OC scales approximately as 〈ρ〉^−1/2 when . Our models easily produce objects on orbits like that of (90377) Sedna [Brown, M.E., Trujillo, C., Rabinowitz, D., 2004. Astrophys. J. 617, 645-649] within ∼1 Myr in cases where the mean density is or higher; one needs mean densities of order to create objects like 2000 CR₁₀₅ by this mechanism, which are reasonable (see, e.g., Guthermuth, R.A., Megeath, S.T., Pipher, J.L., Williams, J.P., Allen, L.E., Myers, P.C., Raines, S.N., 2005. Astrophys. J. 632, 397-420). Thus the latter object may also be part of the OC. Close stellar passages can stir the primordial Kuiper Belt to sufficiently high eccentricities (e?0.05; Kenyon, S.J., Bromley, B.C., 2002. Astron. J. 123, 1757-1775) that collisions become destructive. From the simulations performed it is determined that there is a 50% or better chance to stir the primordial Kuiper Belt to eccentricities e?0.05 at 50 AU when . The orbit of the new object 2003 UB₃₁₃ [Brown, M.E., Trujillo, C.A., Rabinowitz, D.L., 2005. Astrophys. J. 635, L97-L100] is only reproduced for mean cluster densities of the order of , but in the simulations it could not come to be on its current orbit by this mechanism without disrupting the formation of bodies in the primordial Kuiper Belt down to 20 AU. It is therefore improbable that the latter object is created by this mechanism.