The shallow magnitude distribution of asteroid families

It is well known that asteroid families have steeper absolute magnitude (H) distributions for H < 12-13 values than the background population. Beyond this threshold, the shapes of the absolute magnitude distributions in the family/background populations are difficult to determine, primarily because both populations are not yet observationally complete. Using a recently generated catalog containing the proper elements of 106,284 main belt asteroids and an innovative approach, we debiased the absolute magnitude distribution of the major asteroid families relative to the local background populations. Our results indicate that the magnitude distributions of asteroid families are generally not steeper than those of the local background populations for H > 13 (i.e., roughly for diameters smaller than 10 km). In particular, most families have shallower magnitude distributions than the background in the range 15-17 mag. Thus, we conclude that, contrary to previous speculations, the population of kilometer-size asteroids in the main belt is dominated by background bodies rather than by members of the most prominent asteroid families. We believe this result explains why the Spacewatch, Sloan Digital Sky Survey, and Subaru asteroid surveys all derived a shallow magnitude distribution for the dimmer members of the main belt population.We speculate on a few dynamical and collisional scenarios that can explain this shallow distribution. One possibility is that the original magnitude distributions of the families (i.e., at the moment of the formation event) were very shallow for H larger than ∼ 13, and that most families have not yet had the time to collisionally evolve to the equilibrium magnitude distribution that presumably characterizes the background population. A second possibility is that family members smaller than about 10 km, eroded over time by collisional and dynamical processes, have not yet been repopulated by the break-up of larger family members. For this same reason, the older (and possibly characterized by a weaker impact strength) background population shows a shallow distribution in the range 15-60 km.     

The fossilized size distribution of the main asteroid belt

Planet formation models suggest the primordial main belt experienced a short but intense period of collisional evolution shortly after the formation of planetary embryos. This period is believed to have lasted until Jupiter reached its full size, when dynamical processes (e.g., sweeping resonances, excitation via planetary embryos) ejected most planetesimals from the main belt zone. The few planetesimals left behind continued to undergo comminution at a reduced rate until the present day. We investigated how this scenario affects the main belt size distribution over Solar System history using a collisional evolution model (CoEM) that accounts for these events. CoEM does not explicitly include results from dynamical models, but instead treats the unknown size of the primordial main belt and the nature/timing of its dynamical depletion using innovative but approximate methods. Model constraints were provided by the observed size frequency distribution of the asteroid belt, the observed population of asteroid families, the cratered surface of differentiated Asteroid (4) Vesta, and the relatively constant crater production rate of the Earth and Moon over the last 3 Gyr. Using CoEM, we solved for both the shape of the initial main belt size distribution after accretion and the asteroid disruption scaling law . In contrast to previous efforts, we find our derived function is very similar to results produced by numerical hydrocode simulations of asteroid impacts. Our best fit results suggest the asteroid belt experienced as much comminution over its early history as it has since it reached its low-mass state approximately 3.9-4.5 Ga. These results suggest the main belt's wavy-shaped size-frequency distribution is a “fossil” from this violent early epoch. We find that most diameter D?120 km asteroids are primordial, with their physical properties likely determined during the accretion epoch. Conversely, most smaller asteroids are byproducts of fragmentation events. The observed changes in the asteroid spin rate and lightcurve distributions near D∼100-120 km are likely to be a byproduct of this difference. Estimates based on our results imply the primordial main belt population (in the form of D<1000 km bodies) was 150-250 times larger than it is today, in agreement with recent dynamical simulations.     

Karin cluster formation by asteroid impact

Insights into collisional physics may be obtained by studying the asteroid belt, where large-scale collisions produced groups of asteroid fragments with similar orbits and spectra known as the asteroid families. Here we describe our initial study of the Karin cluster, a small asteroid family that formed 5.8±0.2 Myr ago in the outer main belt. The Karin cluster is an ideal ‘natural laboratory’ for testing the codes used to simulate large-scale collisions because the observed fragments produced by the 5.8-Ma collision suffered apparently only limited dynamical and collisional erosion. To date, we have performed more than 100 hydrocode simulations of impacts with non-rotating monolithic parent bodies. We found good fits to the size-frequency distribution of the observed fragments in the Karin cluster and to the ejection speeds inferred from their orbits. These results suggest that the Karin cluster was formed by a disruption of an ≈33-km-diameter asteroid, which represents a much larger parent body mass than previously estimated. The mass ratio between the parent body and the largest surviving fragment, (832) Karin, is ≈0.15-0.2, corresponding to a highly catastrophic event. Most of the parent body material was ejected as fragments ranging in size from yet-to-be-discovered sub-km members of the Karin cluster to dust grains. The impactor was ≈5.8 km across. We found that the ejections speeds of smaller fragments produced by the collision were larger than those of the larger fragments. The mean ejection speeds of >3-km-diameter fragments were . The model and observed ejection velocity fields have different morphologies perhaps pointing to a problem with our modeling and/or assumptions. We estimate that ∼5% of the large asteroid fragments created by the collision should have satellites detectable by direct imaging (separations larger than 0.1 arcsec). We also predict a large number of ejecta binary systems with tight orbits. These binaries, located in the outer main belt, could potentially be detected by lightcurve observations. Hydrocode modeling provides important constraints on the interior structure of asteroids. Our current work suggests that the parent asteroid of the Karin cluster may have been an unfractured (or perhaps only lightly fractured) monolithic object. Simulations of impacts into fractured/rubble pile targets were so far unable to produce the observed large gap between the first and second largest fragment in the Karin cluster, and the steep slope at small sizes (≈6.3 differential index). On the other hand, the parent asteroid of the Karin cluster was produced by an earlier disruptive collision that created the much larger, Koronis family some 2-3 Gyr ago. Standard interpretation of hydrocode modeling then suggests that the parent asteroid of the Karin cluster should have been formed as a rubble pile from Koronis family debris. We discuss several solutions to this apparent paradox.     

The shape distribution of asteroid families: Evidence for evolution driven by small impacts

A statistical analysis of brightness variability of asteroids reveals how their shapes evolve from elongated to rough spheroidal forms, presumably driven by impact-related phenomena. Based on the Sloan Digital Sky Survey Moving Object Catalog, we determined the shape distribution of 11,735 asteroids, with special emphasis on eight prominent asteroid families. In young families, asteroids have a wide range of shape elongations, implying fragmentation-formation. In older families we see an increasing number of rough spheroids, in agreement with the predictions of an impact-driven evolution. Old families also contain a group of moderately elongated members, which we suggest correspond to higher-density, more impact-resistant cores of former fragmented asteroids that have undergone slow shape erosion. A few percent of asteroids have very elongated shapes, and can either be young fragments or tidally reshaped bodies. Our results confirm that the majority of asteroids are gravitationally bound “rubble piles.”     

Size-frequency distributions of fragments from SPH/N-body simulations of asteroid impacts: Comparison with observed asteroid families

We investigate the morphology of size-frequency distributions (SFDs) resulting from impacts into 100-km-diameter parent asteroids, represented by a suite of 161 SPH/N-body simulations conducted to study asteroid satellite formation [Durda, D.D., Bottke, W.F., Enke, B.L., Merline, W.J., Asphaug, E., Richardson, D.C., Leinhardt, Z.M., 2004. Icarus 170, 243-257]. The spherical basalt projectiles range in diameter from 10 to 46 km (in equally spaced mass increments in logarithmic space, covering six discrete sizes), impact speeds range from 2.5 to 7 km/s (generally in 1 km/s increments), and impact angles range from 15° to 75° (nearly head-on to very oblique) in 15° increments. These modeled SFD morphologies match very well the observed SFDs of many known asteroid families. We use these modeled SFDs to scale to targets both larger and smaller than 100 km in order to gain insights into the circumstances of the impacts that formed these families. Some discrepancies occur for families with parent bodies smaller than a few tens of kilometers in diameter (e.g., 832 Karin), however, so due caution should be used in applying our results to such small families. We find that ∼20 observed main-belt asteroid families are produced by the catastrophic disruption of D>100 km parent bodies. Using these data as constraints, collisional modeling work [Bottke Jr., W.F., Durda, D.D., Nesvorný, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H.F., 2005b. Icarus 179, 63-94] suggests that the threshold specific energy, , needed to eject 50% of the target body's mass is very close to that predicted by Benz and Asphaug [Benz, W., Asphaug, E., 1999. Icarus 142, 5-20].     

Physical properties of asteroid dust bands and their sources

Disruptive collisions in the main belt can liberate fragments from parent bodies ranging in size from several micrometers to tens of kilometers in diameter. These debris bodies group at initially similar orbital locations. Most asteroid-sized fragments remain at these locations and are presently observed as asteroid families. Small debris particles are quickly removed by Poynting-Robertson drag or comminution but their populations are replenished in the source locations by collisional cascade. Observations from the Infrared Astronomical Satellite (IRAS) showed that particles from particular families have thermal radiation signatures that appear as band pairs of infrared emission at roughly constant latitudes both above and below the Solar System plane. Here we apply a new physical model capable of linking the IRAS dust bands to families with characteristic inclinations. We use our results to constrain the physical properties of IRAS dust bands and their source families. Our results indicate that two prominent IRAS bands at inclinations ≈2.1° and ≈9.3° are byproducts of recent asteroid disruption events. The former is associated with a disruption of a ≈30-km asteroid occurring 5.8 Myr ago; this event gave birth to the Karin family. The latter came from the breakup of a large >100-km-diameter asteroid 8.3 Myr ago that produced the Veritas family. Using an N-body code, we tracked the dynamical evolution of ≈10⁶ particles, 1 μm to 1 cm in diameter, from both families. We then used these results in a Monte Carlo code to determine how small particles from each population undergo collisional evolution. By computing the thermal emission of particles, we were able to compare our results with IRAS observations. Our best-fit model results suggest the Karin and Veritas family particles contribute by 5-9% in 10-60-μm wavelengths to the zodiacal cloud's brightness within 50° latitudes around the ecliptic, and by 9-15% within 10° latitudes. The high brightness of the zodiacal cloud at large latitudes suggests that it is mainly produced by particles with higher inclinations than what would be expected for asteroidal particles produced by sources in the main belt. From these results, we infer that asteroidal dust represents a smaller fraction of the zodiacal cloud than previously thought. We estimate that the total mass accreted by the Earth in Karin and Veritas particles with diameters 20-400 μm is ≈15,000-20,000 tons per year (assuming 2 g cm⁻³ particles density). This is ≈30-50% of the terrestrial accretion rate of cosmic material measured by the Long Duration Exposure Facility. We hypothesize that up to ≈50% of our collected interplanetary dust particles and micrometeorites may be made up of particle species from the Veritas and Karin families. The Karin family IDPs should be about as abundant as Veritas family IDPs though this ratio may change if the contribution of third, near-ecliptic source is significant. Other sources of dust and/or large impact speeds must be invoked to explain the remaining ≈50-70%. The disproportional contribution of Karin/Veritas particles to the zodiacal cloud (only 5-9%) and to the terrestrial accretion rate (30-50%) suggests that the effects of gravitational focusing by the Earth enhance the accretion rate of Karin/Veritas particles relative to those in the background zodiacal cloud. From this result and from the latitudinal brightness of the zodiacal cloud, we infer that the zodiacal cloud emission may be dominated by high-speed cometary particles, while the terrestrial impactor flux contains a major contribution from asteroidal sources. Collisions and Poynting-Robertson drift produce the size-frequency distribution (SFD) of Karin and Veritas particles that becomes increasingly steeper closer to the Sun. At 1 AU, the SFD is relatively shallow for small particle diameters D (differential slope exponent of particles with D?100 μm is ≈2.2-2.5) and steep for D?100 μm. Most of the mass at 1 AU, as well as most of the cross-sectional area, is contributed by particles with D≈100-200 μm. Similar result has been found previously for the SFD of the zodiacal cloud particles at 1 AU. The fact that the SFD of Karin/Veritas particles is similar to that of the zodiacal cloud suggests that similar processes shaped these particle populations. We estimate that there are ≈5×¹⁰²⁴ Karin and ≈10²⁵ Veritas family particles with D>30 μm in the Solar System today. The IRAS observation of the dust bands may be satisfactorily modeled using ‘averaged’ SFDs that are constant with semimajor axis. These SFDs are best described by a broken power-law function with differential power index α≈2.1-2.4 for D?100 μm and by α?3.5 for 100 μm?D?1 cm. The total cross-sectional surface area of Veritas particles is a factor of ≈2 larger than the surface area of the particles producing the inner dust bands. The total volumes in Karin and Veritas family particles with 1 μm<D<1 cm correspond to D=11 km and D=14 km asteroids with equivalent masses ≈1.5×¹⁰¹⁸ g and ≈3.0×¹⁰¹⁸ g, respectively (assuming 2 g cm⁻³ bulk density). If the size-frequency and radial distribution of particles in the zodiacal cloud were similar to those in the asteroid dust bands, we estimate that the zodiacal cloud represents ∼3×¹⁰¹⁹ g of material (in particles with 1 μm<D<1 cm) at ±10° around the ecliptic and perhaps as much as ∼¹⁰²⁰ g in total. The later number corresponds to about a 23-km-radius sphere with 2 g cm⁻³ density.     

Linking the collisional history of the main asteroid belt to its dynamical excitation and depletion

The main belt is believed to have originally contained an Earth mass or more of material, enough to allow the asteroids to accrete on relatively short timescales. The present-day main belt, however, only contains ∼5×¹⁰⁻⁴ Earth masses. Numerical simulations suggest that this mass loss can be explained by the dynamical depletion of main belt material via gravitational perturbations from planetary embryos and a newly-formed Jupiter. To explore this scenario, we combined dynamical results from Petit et al. [Petit, J. Morbidelli, A., Chambers, J., 2001. The primordial excitation and clearing of the asteroid belt. Icarus 153, 338-347] with a collisional evolution code capable of tracking how the main belt undergoes comminution and dynamical depletion over 4.6 Gyr [Bottke, W.F., Durda, D., Nesvorny, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H., 2005. The fossilized size distribution of the main asteroid belt. Icarus 175, 111-140]. Our results were constrained by the main belt's size-frequency distribution, the number of asteroid families produced by disruption events from diameter D>100 km parent bodies over the last 3-4 Gyr, the presence of a single large impact crater on Vesta's intact basaltic crust, and the relatively constant lunar and terrestrial impactor flux over the last 3 Gyr. We used our model to set limits on the initial size of the main belt as well as Jupiter's formation time. We find the most likely formation time for Jupiter was 3.3±2.6 Myr after the onset of fragmentation in the main belt. These results are consistent with the estimated mean disk lifetime of 3 Myr predicted by Haisch et al. [Haisch, K.E., Lada, E.A., Lada, C.J., 2001. Disk frequencies and lifetimes in young clusters. Astrophys. J. 553, L153-L156]. The post-accretion main belt population, in the form of diameter D?1000 km planetesimals, was likely to have been 160±40 times the current main belt's mass. This corresponds to 0.06-0.1 Earth masses, only a small fraction of the total mass thought to have existed in the main belt zone during planet formation. The remaining mass was most likely taken up by planetary embryos formed in the same region. Our results suggest that numerous D>200 km planetesimals disrupted early in Solar System history, but only a small fraction of their fragments survived the dynamical depletion event described above. We believe this may explain the limited presence of iron-rich M-type, olivine-rich A-type, and non-Vesta V-type asteroids in the main belt today. The collisional lifetimes determined for main belt asteroids agree with the cosmic ray exposure ages of stony meteorites and are consistent with the limited collisional evolution detected among large Koronis family members. Using the same model, we investigated the near-Earth object (NEO) population. We show the shape of the NEO size distribution is a reflection of the main belt population, with main belt asteroids driven to resonances by Yarkovsky thermal forces. We used our model of the NEO population over the last 3 Gyr, which is consistent with the current population determined by telescopic and satellite data, to explore whether the majority of small craters (D<0.1-1 km) formed on Mercury, the Moon, and Mars were produced by primary impacts or by secondary impacts generated by ejecta from large craters. Our results suggest that most small craters formed on these worlds were a by-product of secondary rather than primary impacts.     

The formation of asteroid satellites in large impacts: results from numerical simulations

We present results of 161 numerical simulations of impacts into 100-km diameter asteroids, examining debris trajectories to search for the formation of bound satellite systems. Our simulations utilize a 3-dimensional smooth-particle hydrodynamics (SPH) code to model the impact between the colliding asteroids. The outcomes of the SPH models are handed off as the initial conditions for N-body simulations, which follow the trajectories of the ejecta fragments to search for the formation of satellite systems. Our results show that catastrophic and large-scale cratering collisions create numerous fragments whose trajectories can be changed by particle-particle interactions and by the reaccretion of material onto the remaining target body. Some impact debris can enter into orbit around the remaining target body, which is a gravitationally reaccreted rubble pile, to form a SMAshed Target Satellite (SMATS). Numerous smaller fragments escaping the largest remnant may have similar trajectories such that many become bound to one another, forming Escaping Ejecta Binaries (EEBs). Our simulations so far seem to be able to produce satellite systems qualitatively similar to observed systems in the main asteroid belt. We find that impacts of 34-km diameter projectiles striking at 3 km s⁻¹ at impact angles of ∼30° appear to be particularly efficient at producing relatively large satellites around the largest remnant as well as large numbers of modest-size binaries among their escaping ejecta.     

The formation of asteroid satellites in large impacts: results from numerical simulations

We present results of 161 numerical simulations of impacts into 100-km diameter asteroids, examining debris trajectories to search for the formation of bound satellite systems. Our simulations utilize a 3-dimensional smooth-particle hydrodynamics (SPH) code to model the impact between the colliding asteroids. The outcomes of the SPH models are handed off as the initial conditions for N-body simulations, which follow the trajectories of the ejecta fragments to search for the formation of satellite systems. Our results show that catastrophic and large-scale cratering collisions create numerous fragments whose trajectories can be changed by particle-particle interactions and by the reaccretion of material onto the remaining target body. Some impact debris can enter into orbit around the remaining target body, which is a gravitationally reaccreted rubble pile, to form a SMAshed Target Satellite (SMATS). Numerous smaller fragments escaping the largest remnant may have similar trajectories such that many become bound to one another, forming Escaping Ejecta Binaries (EEBs). Our simulations so far seem to be able to produce satellite systems qualitatively similar to observed systems in the main asteroid belt. We find that impacts of 34-km diameter projectiles striking at 3 km s⁻¹ at impact angles of ∼30° appear to be particularly efficient at producing relatively large satellites around the largest remnant as well as large numbers of modest-size binaries among their escaping ejecta.     

Fragment properties at the catastrophic disruption threshold: The effect of the parent body’s internal structure

Numerical simulations of asteroid breakups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by , which results in the escape of half of the target’s mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of , and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of on the impact speed of the projectile.In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. Indeed, decreasing the strength of non-porous targets make them easier to disrupt in this regime, while increasing the strength of porous targets has much less influence on the value of . Therefore, one cannot say that non-porous targets are systematically easier or more difficult to disrupt than porous ones, as the outcome highly depends on the assumed strength values. In the gravity regime, we also confirm that the process of gravitational reaccumulation is at the origin of the largest remnant’s mass in both cases. We then propose some power-law relationships between and both target’s size and impact speed that can be used in collisional evolution models. The resulting fragment size distributions can also be reasonably fitted by a power-law whose exponent ranges between −2.2 and −2.7 for all target diameters in both cases and independently on the impact velocity (at least in the small range investigated between 3 and 5 km/s). Then, although ejection velocities in the gravity regime tend to be higher from porous targets, they remain on the same order as the ones from non-porous targets.     

Chemical and mineralogical size segregation in the impact disruption of inhomogeneous, anhydrous meteorites

We performed impact disruption experiments on pieces from eight different anhydrous chondritic meteorites—four weathered ordinary chondrite finds from North Africa (NWA791, NWA620, NWA869 and MOR001), three almost unweathered ordinary chondrite falls (Mbale, Gao, and Saratov), and an almost unweathered carbonaceous chondrite fall (Allende). In each case the impactor was a small (1/8 or 1/4 in) aluminum sphere fired at the meteorite target at , comparable to the mean collision speed in the main-belt. Some of the ∼5 to debris from each disruption was collected in aerogel capture cells, and the captured particles were analyzed by in situ synchrotron-based X-ray fluorescence. For each meteorite, many of the smallest particles ( up to in size, depending on the meteorite) exhibit very high Ni/Fe ratios compared to the Ni/Fe ratios measured in the larger particles , a composition consistent with the smallest debris being dominated by matrix material while the larger debris is dominated by fragments from olivine chondrules. These results may explain why the interplanetary dust particles (IDPs) collected from the Earth's stratosphere are C-rich and volatile-rich compared to the presumed solar nebula composition. The IDPs may simply sample the matrix of an inhomogeneous parent body, structurally and mineralogically similar to the chondritic meteorites, which are inhomogeneous assemblages of compact, strong, C- and volatile-poor chondrules that are distributed in a more porous, C- and volatile-rich matrix. In addition, these results may explain why the micrometeorites, which are to millimeters in size, recovered from the polar ices are Ni- and S-poor compared to chondritic meteorites, since these polar micrometeorites may preferentially sample fragments from the Ni- and S-poor olivine chondrules. These results indicate that the average composition of the IDPs may be biased towards the composition of the matrix of the parent body while the average composition of the polar micrometeorites may be more heavily weighted towards the composition of the chondrules and clasts. Thus, neither the IDPs nor the polar micrometeorites may sample the bulk composition of their respective parent bodies.We determined the threshold collisional specific energy for these chondritic meteorites to be 1419 J/kg, about twice the value for terrestrial basalt. Comparison of the mass of the largest fragment produced in the disruption of an sample of the porous ordinary chondrite Saratov with the largest fragment produced in the disruption of an sample of the compact ordinary chondrite MOR001 when each was struck by an impactor having approximately the same kinetic energy confirms that it requires significantly more energy to disrupt a porous target than a non-porous target.These results may also have important implications for the design of spacecraft missions intended to sample the composition and mineralogy of the chondritic asteroids and other inhomogeneous bodies. A Stardust-like spacecraft intended to sample asteroids by collecting only the small debris from a man-made impact onto the asteroid may collect particles that over-sample the matrix of the target and do not provide a representative sample of the bulk composition. The impact collection technique to be employed by the Japanese HAYABUSA (formerly MUSES-C) spacecraft to sample the asteroid Itokawa may result in similar mineral segregation.     

Yarkovsky depletion and asteroid collisional evolution

The orbital parameters of small asteroids change with time, as a consequence of the so-called Yarkovsky effect. This leads to a steady removal of objects from the Main Belt, which takes place when the objects reach one of the major resonant regions in the orbital elements space. The process may influence the evolution of the inventory and size distribution of Main Belt asteroids, but it has not been taken into account by classical models of the collisional evolution of the asteroid population. In this paper we discuss the role of the Yarkovsky effect in producing the current observed size distribution. We show that adding Yarkovsky effect to purely collisional mechanisms may increase the removal of objects at sizes around 1 km by a factor of about 2 with respect to a purely collisional scenario. Moreover, waves in the size distribution may also be produced. However, taking also into account current uncertainties in the efficiency of purely collisional mechanisms, the role of the Yarkovsky effect seems not dominant, and cannot be unambiguously determined.     

Families among high-inclination asteroids

We present a new classification of families identified among the population of high-inclination asteroids. We computed synthetic proper elements for a sample of 18,560 numbered and multi-opposition objects having sine of proper inclination greater than 0.295. We considered three zones at different heliocentric distances (inner, intermediate and outer region) and used the standard approach based on the Hierarchical Clustering Method (HCM) to identify families in each zone. In doing so, we used slightly different approach with respect to previously published methodologies, to achieve a more reliable and robust classification. We also used available SDSS color data to improve membership and identify likely family interlopers. We found a total of 38 families, as well as a significant number of clumps and clusters deserving further investigation.     

Dynamics of the Hungaria asteroids

To try to understand the dynamical and collisional evolution of the Hungaria asteroids we have built a large catalog of accurate synthetic proper elements. Using the distribution of the Hungaria, in the spaces of proper elements and of proper frequencies, we can study the dynamical boundaries and the internal structure of the Hungaria region, both within a purely gravitational model and also showing the signature of the non-gravitational effects. We find a complex interaction between secular resonances, mean motion resonances, chaotic behavior and Yarkovsky-driven drift in semimajor axis. We also find a rare occurrence of large scale instabilities, leading to escape from the region. This allows to explain the complex shape of a grouping which we suggest is a collisional family, including most Hungaria but by no means all; we provide an explicit list of non-members of the family. There are finer structures, of which the most significant is a set of very close asteroid couples, with extremely similar proper elements. Some of these could have had, in a comparatively recent past, very close approaches with low relative velocity. We argue that the Hungaria, because of the favorable observing conditions, may soon become the best known sub-group of the asteroid population.     

Craters on asteroids: Reconciling diverse impact records with a common impacting population

O'Brien and Greenberg [O'Brien, D.P., Greenberg, R., 2005. Icarus 178, 179-212] developed a self-consistent numerical model of the collisional and dynamical evolution of the main-belt and NEA populations that was tested against a diverse range of observational and theoretical constraints. In this paper, we use those results to update the asteroid cratering model of Greenberg et al. [Greenberg, R., Nolan, M.C., Bottke, W.F., Kolvoord, R.A., Veverka, J., 1994. Icarus 107, 84-97; Greenberg, R., Bottke, W.F., Nolan, M., Geissler, P., Petit, J., Durda, D.D., Asphaug, E., Head, J., 1996. Icarus 120, 106-118], and show that the main-belt asteroid population from the O'Brien and Greenberg collisional/dynamical evolution modeling is consistent with the crater records on Gaspra, Ida, Mathilde, and Eros, the four asteroids that have been observed by spacecraft.     

Combined effect of YORP and collisions on the rotation rate of small Main Belt asteroids

The rotation rate distribution of small Main Belt asteroids is dominated by YORP and collisions. These mechanism act differently depending on the size of the bodies and give rise to non-linear effects when they both operate. Using a Monte Carlo method we model the formation of a steady state population of small asteroids under the influence of both mechanisms and the rotation rate distribution is compared to the observed one as derived from Pravec et al. (Pravec, P. et al. [2008]. Icarus 197, 497-504). A better match to observations is obtained with respect to the case in which only YORP is considered. In particular, an excess of slow rotators is produced in the model with both collisions and YORP because bodies driven to slow rotation by YORP have a random walk-like evolution of the spin induced by repeated collisions with small projectiles. This is a dynamical evolution different from tumbling and it lasts until a large impact takes the body to a faster rotation rate. According to our model, the rotational fission of small asteroids is a very frequent event and might explain objects like P/2010 A2 and its associated tail of millimeter-sized dust particles. The mass loss during fission of small asteroids might significantly influence the overall collisional evolution of the belt. Fission can in fact be considered as an additional erosion mechanism, besides cratering and fragmentation, acting only at small diameters.     

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 .     

The collisional and dynamical evolution of the main-belt and NEA size distributions

The size distribution of main belt of asteroids is determined primarily by collisional processes. Large asteroids break up and form smaller asteroids in a collisional cascade, with the outcome controlled by the strength-size relationship for asteroids. In addition to collisional processes, the non-collisional removal of asteroids from the main belt (and their insertion into the near-Earth asteroid (NEA) population) is critical, and involves several effects: strong resonances increase the orbital eccentricity of asteroids and cause them to enter the inner planet region; chaotic diffusion by numerous weak resonances causes a slow leak of asteroids into the Mars- and Earth-crossing populations; and the Yarkovsky effect, a radiation force on asteroids, is the primary process that drives asteroids into these resonant escape routes. Yarkovsky drift is size-dependent and can modify the main-belt size distribution. The NEA size distribution is primarily determined by its source, the main-belt population, and by the size-dependent processes that deliver bodies from the main belt. All of these effects are simulated in a numerical collisional evolution model that incorporates removal by non-collisional processes. We test our model against a wide range of observational constraints, such as the observed main-belt and NEA size distributions, the number of asteroid families, the preserved basaltic crust of Vesta and its large south-pole impact basin, the cosmic ray exposure ages of meteorites, and the cratering records on asteroids. We find a strength-size relationship for main-belt asteroids and non-collisional removal rates from the main belt such that our model fits these constraints as best as possible within the parameter space we explore. Our results are consistent with other independent estimates of strength and removal rates.     

Catastrophic disruption of asteroids and family formation: a review of numerical simulations including both fragmentation and gravitational reaccumulations

In the last few years, thanks to the development of sophisticated numerical codes, a major breakthrough has been achieved in our understanding of the processes involved in small body collisions. In this review, we summarize the most recent results provided by numerical simulations, accounting for both the fragmentation of an asteroid and the gravitational interactions of the generated fragments. These studies have greatly improved our knowledge of the mechanisms that are at the origin of some observed features in the asteroid belt. In particular, the simulations have demonstrated that, for bodies larger than several kilometers, the collisional process not only involves the fragmentation of the asteroid but also the gravitational interactions between the ejected fragments. This latter mechanism can lead to the formation of large aggregates by gravitational reaccumulation of smaller fragments, and helps explain the presence of large members within asteroid families. Numerical simulations of the complete process have thus reproduced successfully for the first time the main properties of asteroid families, each formed by the disruption of a large parent body, and provided information on the possible internal structure of the parent bodies. A large amount of work remains necessary, however, to understand in deeper detail the physical process as a function of material properties and internal structures that are relevant to asteroids, and to determine in a more quantitative way the outcome properties such as fragment shapes and rotational states.     

Non-power law behavior in fragmentation cascades

Collisions resulting in fragmentation are important in shaping the mass spectrum of minor bodies in the asteroid belt, the Kuiper Belt, and debris disks. Models of fragmentation cascades typically find that in steady-state, the solution for the particle mass distribution is a power law in the mass. However, previous studies have typically assumed that the mass of the largest fragment produced in a collision with just enough energy to shatter the target and disperse half its mass to infinity is directly proportional to the target mass. We show that if this assumption is not satisfied, then the power law solution for the steady-state particle mass distribution is modified by a multiplicative factor, which is a slowly varying function of the mass. We derive analytic solutions for this correction factor and confirm our results numerically. We find that this correction factor proves important when extrapolating over many orders of magnitude in mass, such as when inferring the number of large objects in a system based on infrared observations. In the course of our work, we have also discovered an unrelated type of non-power law behavior: waves can persist in the mass distribution of objects even in the absence of upper or lower cutoffs to the mass distribution or breaks in the strength law.