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.     

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.     

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

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.     

On the asteroid belt's orbital and size distribution

For absolute magnitudes greater than the current completeness limit of H-magnitude ∼15 the main asteroid belt's size distribution is imperfectly known. We have acquired good-quality orbital and absolute H-magnitude determinations for a sample of small main-belt asteroids in order to study the orbital and size distribution beyond H=15, down to sub-kilometer sizes (H>18). Based on six observing nights over a 11-night baseline we have detected, measured photometry for, and linked observations of 1087 asteroids which have one-week time baselines or more. The linkages allow the computation of full heliocentric orbits (as opposed to statistical distances determined by some past surveys). Judged by known asteroids in the field the typical uncertainty in the (a/e/i) orbital elements is less than 0.03 AU/0.03/0.5°. The distances to the objects are sufficiently well known that photometric uncertainties (of 0.3 magnitudes or better) dominate the error budget of their derived H-magnitudes. The detected asteroids range from HR=12-22 and provide a set of objects down to sizes below 1 km in diameter. We find an on-sky surface density of 210 asteroids per square degree in the ecliptic with opposition magnitudes brighter than mR=23, with the cumulative number of asteroids increasing by a factor of 10^0.27/mag from mR=18 down to the mR?23.5 limit of our survey. In terms of absolute H magnitudes, we find that beyond H=15 the belt exhibits a constant power-law slope with the number increasing proportional to ¹⁰0.30H from H?15 to 18, after which incompleteness begins in the survey. Examining only the subset of detections inside 2.5 AU, we find weak evidence for a mildly shallower slope for H=15-19.5. We provide the information necessary such that anyone wishing to model the main asteroid belt can compare a detailed model to our detected sample.     

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 primordial excitation and clearing of the asteroid belt—Revisited

We have performed new simulations of two different scenarios for the excitation and depletion of the primordial asteroid belt, assuming Jupiter and Saturn on initially circular orbits as predicted by the Nice Model of the evolution of the outer Solar System [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465]. First, we study the effects of sweeping secular resonances driven by the depletion of the solar nebula. We find that these sweeping secular resonances are incapable of giving sufficient dynamical excitation to the asteroids for nebula depletion timescales consistent with estimates for solar-type stars, and in addition cannot cause significant mass depletion in the asteroid belt or produce the observed radial mixing of different asteroid taxonomic types. Second, we study the effects of planetary embryos embedded in the primordial asteroid belt. These embedded planetary embryos, combined with the action of jovian and saturnian resonances, can lead to dynamical excitation and radial mixing comparable to the current asteroid belt. The mass depletion driven by embedded planetary embryos alone, even in the case of an eccentric Jupiter and Saturn, is roughly 10-20× less than necessary to explain the current mass of the main belt, and thus a secondary depletion event, such as that which occurs naturally in the Nice Model, is required. We discuss the implications of our new simulations for the dynamical and collisional evolution of the main belt.     

Catastrophic disruption of pre-shattered parent bodies

In this paper, we analyze the effect of the internal structure of a parent body on its fragment properties following its disruption in different impact energy regimes. To simulate an asteroid breakup, we use the same numerical procedure as in our previous studies, i.e., a 3D SPH hydrocode to compute the fragmentation phase and the parallel N-body code pkdgrav to compute the subsequent gravitational re-accumulation phase. To explore the importance of the internal structure in determining the collisional outcome, we consider two different parent body models: (1) a purely monolithic one and (2) a pre-shattered one which consists of several fragments separated by damaged zones and small voids. We present here simulations spanning two different impact energy regimes—barely disruptive and highly catastrophic—corresponding to the formation of the Eunomia and Koronis families, respectively. As we already found for the intermediate energy regime represented by the Karin family, pre-shattered parent bodies always lead to outcome properties in better agreement with those of real families. In particular, the fragment size distribution obtained by disrupting a monolithic body always contains a large gap between the largest fragment and the next largest ones, whereas it is much more continuous in the case of a pre-shattered parent body. In the latter case, the ejection speeds of large fragments are also higher and a smaller impact energy is generally required to achieve a similar degree of disruption. Hence, unless the internal structure of bodies involved in a collision is known, predicting accurately the outcome is impossible. Interestingly, disrupting a pre-shattered parent body to reproduce the Koronis family yields a fragment size distribution characterized by four almost identical largest objects, as observed in the real family. This peculiar outcome has been found before in laboratory experiments but is obtained for the first time following gravitational re-accumulation. Finally, we show that material belonging to the largest fragments of a family originates from well-defined regions inside the parent body (the extent and location of which are dependent upon internal structure), despite the many gravitational interactions that occur during the re-accumulation process. Hence fragment formation does not proceed stochastically but results directly from the velocity field imparted during the impact.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

V-type asteroids in the middle main belt

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