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.     

Interacting ellipsoids: a minimal model for the dynamics of rubble-pile bodies

A simple mechanical model is formulated to study the dynamics of rubble-pile asteroids, formed by the gravitational re-accumulation of fragments after the collisional breakup of a parent body. In this model, a rubble-pile consists of N interacting fragments represented by rigid ellipsoids, and the equations of motion explicitly incorporate the minimal degrees of freedom necessary to describe the attitude and rotational state of each fragment. In spite of its simplicity, our numerical examples indicate that the overall behavior of our model is in line with several known properties of collisional events, like the energy and angular momentum partition during high velocity impacts. Therefore, it may be considered as a well defined minimal model.     

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.     

Impact cratering and break up of the small bodies of the Solar System

We consider the largest impact craters observed on small satellites and asteroids and the impact disruption of such bodies. Observational data are considered from 21 impact-like structures on 13 satellites and 8 asteroids (target body radii in the range 0.7-265 km). If the radius of the target body is R and the diameter of the largest crater observed on this body D, the ratio D/R is then the main observational parameter of interest. This is found on the observed bodies and compared to data obtained in the laboratory. Taking the largest observed value for D/R as a proxy for the ratio Dc/R (where Dc is the diameter of the largest crater that can be formed on a body without shattering it) it was found that for the observed icy satellites Dc,icy≈1.2R and for the asteroids and the rocky satellites Dc,rocky≈1.6R. In laboratory experiments with ice targets at impactor speeds of 1 to 3 km s⁻¹ we obtained Dc,icy≈1.64R.     

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.     

Experimental study on the impact fragmentation of core-mantle bodies: Implications for collisional disruption of rocky planetesimals with sintered core covered with porous mantle

Impact experiments of inhomogeneous targets such as layered bodies consisting of a dense core and porous mantle were conducted to clarify the effect of the layered structure on impact strength. The layered structure of small bodies could be the result of the thermal evolution of planetesimals in the solar nebula. So, the impact disruption of thermally evolved bodies with core-mantle structure is important for the origin of small bodies such as asteroids. We investigated the impact strength of rocky-layered bodies with porous mantle-sintered cores, which could be formed at an initial stage of thermal evolution. Spherical targets composed of soda-lime glass or quartz core and porous gypsum mantle were prepared as an analog of small bodies with a core-mantle structure, and the internal structure was changed. A nylon projectile was impacted at the impact velocity from 1 to 5 km/s. The impact strength of the core-mantle targets decreases with the increase of the core/target mass ratio (RCM) in the specific energy range from 1×¹⁰³ to 4×¹⁰⁴ J/kg. We observed two distinct destruction modes characterized by the damage to the core: one shows a damaged core and fractured mantle, and the other shows an intact core and broken mantle. The former mode was usually observed with increasing RCM, and the boundary condition of the core destruction () was experimentally found to be , where is the specific energy required to disrupt a glass core. From this empirical equation, it might be possible to discuss the destruction conditions of a thermally evolved body with a porous mantle-sintered core structure. We speculate that the impact strength of the body could be significantly reduced with the progress of internal evolution at the initial stage of thermal evolution.     

Mass dispersal and angular momentum transfer during collisions between rubble-pile asteroids. II. Effects of initial rotation and spin-down through disruptive collisions

Expanding on our previous N-body simulation of impacts between initially non-rotating rubble-pile objects [Takeda, T., Ohtsuki, K., 2007. Icarus 189, 256-273], we examine effects of initial rotation of targets on mass dispersal and change of spin rates. Numerical results show that the collisional energy needed to disrupt a rubble-pile object is not sensitive to initial rotation of the target, in most of the parameter range studied in our simulations. We find that initial rotation of targets is slowed down through disruptive impacts for a wide range of parameters. The spin-down is caused by escape of high-velocity ejecta and asymmetric re-accumulation of fragments. When these effects are significant, rotation is slowed down even when the angular momentum added by an impactor is in the same direction as the initial rotation of the target. Spin-down is most efficient when the impact occurs in the equatorial plane of the target, because in this case most of the ejected fragments originate from the equatorial region of the target and a significant amount of angular momentum can be easily removed. In the case of impacts from directions inclined relative to the target's equatorial plane, spin-down still occurs with reduced degree, unless impacts occur onto the pole region from the vertical direction. Our results suggest that such spin-down through disruptive impacts may have played an important role in spin evolution of asteroids through collisions in the gravity-dominated regime.     

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].     

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.”     

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.     

The Asteroid Veritas: An intruder in a family named after it?

The Veritas family is located in the outer main belt and is named after its apparent largest constituent, Asteroid (490) Veritas. The family age has been estimated by two independent studies to be quite young, around 8 Myr. Therefore, current properties of the family may retain signatures of the catastrophic disruption event that formed the family. In this paper, we report on our investigation of the formation of the Veritas family via numerical simulations of catastrophic disruption of a 140-km-diameter parent body, which was considered to be made of either porous or non-porous material, and a projectile impacting at 3 or 5 km/s with an impact angle of 0° or 45°. Not one of these simulations was able to produce satisfactorily the estimated size distribution of real family members. Based on previous studies devoted to either the dynamics or the spectral properties of the Veritas family, which already treated (490) Veritas as a special object that may be disconnected from the family, we simulated the formation of a family consisting of all members except that asteroid. For that case, the parent body was smaller (112 km in diameter), and we found a remarkable match between the simulation outcome, using a porous parent body, and the real family. Both the size distribution and the velocity dispersion of the real reduced family are very well reproduced. On the other hand, the disruption of a non-porous parent body does not reproduce the observed properties very well. This is consistent with the spectral C-type of family members, which suggests that the parent body was porous and shows the importance of modeling the effect of this porosity in the fragmentation process, even if the largest members are produced by gravitational reaccumulation during the subsequent gravitational phase. As a result of our investigations, we conclude that it is very likely that the Asteroid (490) Veritas and probably several other small members do not belong to the family as originally defined, and that the definition of this family should be revised. Further investigations will be performed to better constrain the definitions and properties of other asteroid families of different types, using the appropriate model of fragmentation. The identification of very young families in turn will continue to serve as a tool to check the validity of numerical models.     

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.     

Mass dispersal and angular momentum transfer during collisions between rubble-pile asteroids

We perform N-body simulations of impacts between initially non-rotating rubble-pile asteroids, and investigate mass dispersal and angular momentum transfer during such collisions. We find that the fraction of the dispersed mass (Mdisp) is approximately proportional to , where Qimp is the impact kinetic energy; the power index α is about unity when the impactor is much smaller than the target, and 0.5?α<1 for impacts with a larger impactor. Mdisp is found to be smaller for more dissipative impacts with small values of the restitution coefficient of the constituent particles. We also find that the efficiency of transfer of orbital angular momentum to the rotation of the largest remnant depends on the degree of disruption. In the case of disruptive oblique impacts where the mass of the largest remnant is about half of the target mass, most of the orbital angular momentum is carried away by the escaping fragments and the efficiency becomes very low (<0.05), while the largest remnant acquires a significant amount of spin angular momentum in moderately disruptive impacts. These results suggest that collisions likely played an important role in rotational evolution of small asteroids, in addition to the recoil force of thermal re-radiation.     

Can planetesimals left over from terrestrial planet formation produce the lunar Late Heavy Bombardment?

The lunar Late Heavy Bombardment (LHB) defines a time between ∼3.8 to possibly 4.1 Gy ago when the Nectarian and early-Imbrium basins on the Moon with reasonably well-constrained ages were formed. Some have argued that these basins were produced by a terminal cataclysm that caused a spike in the inner Solar System impactor flux during this interval. Others have suggested the basins were formed by the tail end of a monotonically decreasing impactor population originally produced by planet formation processes in the inner Solar System. Here we investigate whether this so-called declining bombardment scenario of the LHB is consistent with constraints provided by planet formation models as well as the inferred ages of Nectaris, Serenitatis, Imbrium, and Orientale. We did this by modeling the collisional and dynamical evolution of the post-planet formation population (PPP) for a range of starting PPP masses. Using a Monte Carlo code, we computed the probability that the aforementioned basins were created at various times after the Moon-forming event approximately 4.54 Ga. Our results indicate that the likelihood that the declining bombardment scenario produced Nectaris, Serenitatis, Imbrium, and Orientale (or even just Imbrium and Orientale) at any of their predicted ages is extremely low and can be ruled out at the 3σ confidence level, regardless of the PPP's starting mass. The reason is that collisional and dynamical evolution quickly depletes the PPP, leaving behind a paucity of large projectiles capable of producing the Moon's youngest basins between 3.8-4.1 Gy ago. If collisions are excluded from our model, we find that the PPP produces numerous South Pole-Aitken-like basins during the pre-Nectarian period. This is inconsistent with our understanding of lunar topography. Accordingly, our results lead us to conclude that the terminal cataclysm scenario is the only existing LHB paradigm at present that is both viable from a dynamical modeling perspective and consistent with existing constraints.     

Impact experiments of porous gypsum-glass bead mixtures simulating parent bodies of ordinary chondrites: Implications for re-accumulation processes related to rubble-pile formation

Laboratory impact experiments were conducted for gypsum-glass bead targets simulating the parent bodies of ordinary chondrites. The effects of the chondrules included in the parent bodies on impact disruption were experimentally investigated in order to determine the impact conditions for the formation of rubble-pile bodies after catastrophic disruption. The targets included glass beads with a diameter ranging from 100 μm to 3 mm and the volume fraction was 0.6, similar to that of ordinary chondrites, which is about 0.65-0.75. Nylon projectiles with diameters of 10 mm and 2 mm were impacted at 60-180 m s⁻¹ by a single-stage gas gun and at 4 km s⁻¹ by a two-stage light gas gun, respectively. The impact strength of the gypsum-glass bead target was found to range from 56 to 116 J kg⁻¹ depending on the glass bead size, and was several times smaller than that of the porous gypsum target, 446 J kg⁻¹ in low-velocity collisions. The impact strengths of the 100 μm bead target and the porous gypsum target strongly depended on the impact velocity: those obtained in high-velocity collisions were several times greater than those obtained in low-velocity collisions. The velocities of fragments ejected from two corners on the impact surface of the target, measured in the center of the mass system, were slightly dependent on the target materials, irrespective of impact velocity. These results suggest that chondrule-including planetesimals (CiPs) can reconstruct rubble-pile bodies in catastrophic disruptions at the size of the planetesimal smaller than that of planetesimals without chondrules.     

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.     

Compaction and fragmentation of porous gypsum targets from low-velocity impacts

We performed low-velocity impact experiments of gypsum spheres with porosity ranging from 0 to 61% and diameter ranging from 25 to 83 mm. The impact velocity was from 0.2 to 22 m/s. The target was an iron plate. The outcome of gypsum spheres with porosity 31-61% was different from those of non-porous ice [Higa M., Arakawa, M., Maeno, N., 1996. Planet. Space Sci. 44, 917-925; Higa M., Arakawa, M., Maeno, N., 1998. Icarus 133, 310-320] and non-porous gypsum. In between the intact and fragmentation modes, the outcome of the non-porous ice and gypsum was crack growth at the impact point. However, the outcome of the porous gypsum was compaction. We found that the restitution coefficients of the porous gypsum spheres were all in a similar range, in spite of the difference of the porosity and size at impact velocities up to about 10 m/s where they begin to be fragmented in pieces. Moreover, there is not a large difference between the restitution coefficient of porous and non-porous gypsum. These results collectively indicate that restitution coefficient of gypsum spheres of cm-size is not strongly dependent upon the porosity and compaction process.     

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.