Asteroid families: evidence of ageing of the proper elements

In recent times it has been emphasized that the present kinematical structures of asteroid families should be evolved with respect to the original post-impact situations, according to numerical simulations performed taking into account also the previously neglected Yarkovsky effect. In this paper we show that also a “classical” approach based on an analysis of the current kinematical properties of families leads to conclude that the distributions of proper eccentricities and semimajor axes of family members exhibit evidence of an evolution. The importance of this approach is that it yields a fully independent and quantitative estimate of an evolutionary spreading of the proper elements. In particular, we find that the original post-impact families had to be on the average about twice more compact than the families we observe now, when considering family members down to about 5 km in size. This result can be used in future analyses to derive estimates of the ages of different families, and to better constrain the typical values of the ejection velocities of the fragments in family-forming events.     

Orbital evolution of the Gefion and Adeona asteroid families: close encounters with massive asteroids and the Yarkovsky effect

Asteroid families are groupings of minor planets identified by clustering in their proper orbital elements; these objects have spectral signatures consistent with an origin in the break-up of a common parent body. From the current values of proper semimajor axes a of family members one might hope to estimate the ejection velocities with which the fragments left the putative break-up event (assuming that the pieces were ejected isotropically). However, the ejection velocities so inferred are consistently higher than N-body and hydro-code simulations, as well as laboratory experiments, suggest. To explain this discrepancy between today’s orbital distribution of asteroid family members and their supposed launch velocities, we study whether asteroid family members might have been ejected from the collision at low speeds and then slowly drifted to their current positions, via one or more dynamical processes. Studies show that the proper a of asteroid family members can be altered by two mechanisms: (i) close encounters with massive asteroids, and (ii) the Yarkovsky non-gravitational effect. Because the Yarkovsky effect for kilometer-sized bodies decreases with asteroid diameter D, it is unlikely to have appreciably moved large asteroids (say those with D > 15 km) over the typical family age (1-2 Gyr).For this reason, we numerically studied the mobility of family members produced by close encounters with main-belt, non-family asteroids that were thought massive enough to significantly change their orbits over long timescales. Our goal was to learn the degree to which perturbations might modify the proper a values of all family members, including those too large to be influenced by the Yarkovsky effect. Our initial simulations demonstrated immediately that very few asteroids were massive enough to significantly alter relative orbits among family members. Thus, to maximize gravitational perturbations in our 500-Myr integrations, we investigated the effect of close encounters on two families, Gefion and Adeona, that have high encounter probabilities with 1 Ceres, by far the largest asteroid in the main belt. Our results show that members of these families spreads in a of less than 5% since their formation. Thus gravitational interactions cannot account for the large inferred escape velocities.The effect of close encounters with massive asteroids is, however, not entirely negligible. For about 10% of the simulated bodies, close encounters increased the “inferred” ejection velocities from sub-100 m/s to values greater than 100 m/s, beyond what hydro-code and N-body simulations suggest are the maximum possible initial ejection velocity for members of Adeona and Gefion with D > 15 km. Thus this mechanism of mobility may be responsible for the unusually high inferred ejection speeds of a few of the largest members of these two families.To understand the orbital evolution of the entire family, including smaller members, we also performed simulations to account for the drift of smaller asteroids caused by the Yarkovsky effect. Our two sets of simulations suggest that the two families we investigated are relatively young compared to larger families like Koronis and Themis, which have estimated ages of about 2 Byr. The Adeona and Gefion families seems to be no more than 600 and 850 Myr old, respectively.     

Asteroid Mean Elements: Higher Order and Iterative Theories

Mean orbital elements are obtained from osculating ones by removing the short periodic perturbations. Large catalogues of asteroid mean elements need to be computed, as a first step in the computation of proper elements, used to study asteroid families. The algorithms for this purpose available so far are only accurate to first order in the masses of the perturbing planet; the mean elements have satisfactory accuracy for most of the asteroid belt, but degraded accuracy in the neighbourhoods of the main mean motion resonances, especially the 2:1. We investigate a number of algorithms capable of improving this approximation; they belong to the two classes of Breiter-type methods and iterative methods. The former are obtained by applying some higher order numerical integration scheme, such as Runge–Kutta, to the differential equation whose solution is a transformation removing the fast angular variables from the equations; they can be used to compute a full second order theory, however, only if the full second order determining function is explicitly computed, and this is computationally too cumbersome for a complicated problem such as the N-body. The latter are fixed point iterative schemes, with the first order theory as an iteration step, used to compute the inverse map from mean to osculating elements; formally the method is first order, but because they implement a fixed frequency perturbation theory, they are more accurate than conventional single iteration methods; a similar method is already in use in our computation of proper from mean elements. Many of these methods are tested on a sample of asteroid orbits taken from the Themis family, up to the edge of the 2:1 resonance, and the dispersion of the values of the computed mean semimajor axis over 100 000 years is used as quality control. The results of these tests indicate that the iterative methods are superior, in this specific application, to the Breiter methods, in accuracy and reliability. This is understood as the result of the cancellations occurring between second order perturbation terms: the incomplete second order theory, resulting from the use of a Breiter method with the first order determining function only, can be less accurate than complete, fixed frequency theories of the first order. We have therefore computed new catalogues of asteroid mean and proper elements, incorporating an iterative algorithm in both steps (osculating to mean and mean to proper elements). This new data set, significantly more reliable even in the previously degraded regions of Themis and Cybele, is in the public domain. This revised version was published online in July 2006 with corrections to the Cover Date.     

Asteroid taxonomic signatures from photometric phase curves

We explore the correlation between an asteroid’s taxonomy and photometric phase curve using the H, G₁₂ photometric phase function, with the shape of the phase function described by the single parameter G₁₂. We explore the usability of G₁₂ in taxonomic classification for individual objects, asteroid families, and dynamical groups. We conclude that the mean values of G₁₂ for the considered taxonomic complexes are statistically different, and also discuss the overall shape of the G₁₂ distribution for each taxonomic complex. Based on the values of G₁₂ for about half a million asteroids, we compute the probabilities of C, S, and X complex membership for each asteroid. For an individual asteroid, these probabilities are rather evenly distributed over all of the complexes, thus preventing meaningful classification. We then present and discuss the G₁₂ distributions for asteroid families, and predict the taxonomic complex preponderance for asteroid families given the distribution of G₁₂ in each family. For certain asteroid families, the probabilistic prediction of taxonomic complex preponderance can clearly be made. In particular, the C complex preponderant families are the easiest to detect, the Dora and Themis families being prime examples of such families. We continue by presenting the G₁₂-based distribution of taxonomic complexes throughout the main asteroid belt in the proper element phase space. The Nysa–Polana family shows two distinct regions in the proper element space with different G₁₂ values dominating in each region. We conclude that the G₁₂-based probabilistic distribution of taxonomic complexes through the main belt agrees with the general view of C complex asteroid proportion increasing towards the outer belt. We conclude that the G₁₂ photometric parameter cannot be used in determining taxonomic complex for individual asteroids, but it can be utilized in the statistical treatment of asteroid families and different regions of the main asteroid belt.     

The Flora Family: A Case of the Dynamically Dispersed Collisional Swarm?

Asteroid families are believed to originate by catastrophic disruptions of large asteroids. They are nowadays identified as clusters in the proper orbital elements space. The proper elements are analytically defined as constants of motion of a suitably simplified dynamical system. Indeed, they are generally nearly constant on a 10⁷-10⁸-year time scale. Over longer time intervals, however, they may significantly change, reflecting the accumulation of the tiny nonperiodic evolutions provided by chaos and nonconservative forces. The most important effects leading to a change of the proper orbital elements are (i) the chaotic diffusion in narrow mean motion resonances, (ii) the Yarkovsky nongravitational force, and (iii) the gravitational impulses received at close approaches with large asteroids. A natural question then arises: How are the size and shape of an asteroid family modified due to evolution of the proper orbital elements of its members over the family age? In this paper, we concentrate on the dynamical dispersion of the proper eccentricity and inclination, which occurs due to (i), but with the help of (ii) and (iii). We choose the Flora family as a model case because it is unusually dispersed in eccentricity and inclination and, being located in the inner main belt, is intersected by a large number of effective mean motion resonances with Mars and Jupiter. Our results suggest that the Flora family dynamically disperses on a few 10⁸-year time scale and that its age may be significantly less than 10⁹ years. We discuss the possibility that the parent bodies of the Flora family and of the ordinary L chondrite meteorites are the same object. In a broader sense, this work suggests that the common belief that the present asteroid families are simple images of their primordial dynamical structure should be revised.     

The size distributions of asteroid families in the SDSS Moving Object Catalog 4

Asteroid families, traditionally defined as clusters of objects in orbital parameter space, often have distinctive optical colors. We show that the separation of family members from background interlopers can be improved with the aid of SDSS colors as a qualifier for family membership. Based on an ∼88,000 object subset of the Sloan Digital Sky Survey Moving Object Catalog 4 with available proper orbital elements, we define 37 statistically robust asteroid families with at least 100 members (12 families have over 1000 members) using a simple Gaussian distribution model in both orbital and color space. The interloper rejection rate based on colors is typically ∼10% for a given orbital family definition, with four families that can be reliably isolated only with the aid of colors. About 50% of all objects in this data set belong to families, and this fraction varies from about 35% for objects brighter than an H magnitude of 13 and rises to 60% for objects fainter than this. The fraction of C-type objects in families decreases with increasing H magnitude for H>13, while the fraction of S-type objects above this limit remains effectively constant. This suggests that S-type objects require a shorter timescale for equilibrating the background and family size distributions via collisional processing. The size distribution varies significantly among families, and is typically different from size distributions for background populations. The size distributions for 15 families display a well-defined change of slope and can be modeled as a “broken” double power-law. Such “broken” size distributions are twice as likely for S-type familes than for C-type families (73% vs. 36%), and are dominated by dynamically old families. The remaining families with size distributions that can be modeled as a single power law are dominated by young families (<1 Gyr). When size distribution requires a double power-law model, the two slopes are correlated and are steeper for S-type families. No such slope-color correlation is discernible for families whose size distribution follows a single power law. For several very populous families, we find that the size distribution varies with the distance from the core in orbital-color space, such that small objects are more prevalent in the family outskirts. This “size sorting” is consistent with predictions based on the Yarkovsky effect.     

Reanalysis of asteroid families structure through visible spectroscopy

The taxonomic properties of the main asteroid families are analyzed and discussed in the light of an updated definition of the families using a large proper elements database and the asteroids taxonomy derived from reflectance spectra recently obtained by two large visible spectroscopic surveys: the SMASS II and the S3OS2. Our analysis indicates that most families are quite homogeneous taxonomically and mineralogically—whenever there exists a mineralogical constraint—, being probably originated from homogeneous parent bodies. The exceptions are the Nysa family, that should likely be considered a clan, and the Eos family that encompasses a broad range of taxonomies, whose mineralogical relations cannot be completely ruled out. Only in a few cases the families may be taxonomically distinguished from the background population. That is the case of the Minerva/Gefion, Adeona, Dora, Merxia, Hoffmeister, Koronis, Eos, and Veritas families. Some of the families presented in this work show a larger spectral diversity than previously reported, as it is the case for the Maria and Koronis families. On the other hand, the Veritas family is found to be homogeneous, in sharp contrast with previous works. Mineralogical relations are reported whenever they could be found in the literature and we examine the possible constraints posed by the presence of different taxonomies in certain families.     

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.     

Reconstructing the orbital history of the Veritas family

The family of (490) Veritas is a young, dynamically heterogeneous asteroid family, located in the outer main belt. As such, it represents a valuable example for studying the effects of chaotic diffusion on the shape of asteroid families. The Veritas family can be decomposed into several groups, in terms of the principal mechanisms that govern the local dynamics, which are analyzed here. A relatively large spread in proper eccentricity is observed, for the members of two chaotic groups. We show that different types of chaos govern the motion of bodies within each group, depending on the extent of overlap among the components of the corresponding resonant multiplets. In particular, one group appears to be strongly diffusive, while the other is not. Studying the evolution of the diffusive group and applying statistical methods, we estimate the age of the family to be τ=(8.7±1.7) Myr. This value is statistically compatible with that of 8.3 Myr previously derived by Nesvorný et al. [Nesvorný, D., Bottke, W.F., Levison, H.F., Dones, L., 2003. Astrophys. J. 591, 486-497], who analyzed the secular evolution of family members on regular orbits. Our methodology, applied here in the case of the Veritas family, can be used to reconstruct the orbital history of other, dynamically complex, asteroid families and derive approximate age estimates for young asteroid families, located in diffusive regions of the main belt. Possible refinements of the method are also discussed.     

Catastrophic fragmentation and formation of families: Preliminary results from a new numerical model

Preliminary results of an improved version of the semiempirical model for catastrophic break up processes developed by Paolicchi et al., (1989) are presented. Among the several changes with respect to the old version, the most important seem to be related to the new treatment of gravitational effects, including self-compression and reaccumulation of fragments. In particular, the new model is able to analyze processes involving both cm-sized objects, like those studied by means of laboratory experiments, as well as much larger bodies, for which self-gravitational effects are dominant; moreover, in this latter case the model seems in principle adequate to describe with the same physics very different phenomena, like the formation of plausible asteroid families and the creation of single, rapidly spinning, objects. This fact, if confirmed by refined analyses, may be of high importance for our general understanding of asteroid collisional evolution.     

Synthetic Proper Elements for Outer Main Belt Asteroids

For the orbits with low to moderate inclination and eccentricity, in the asteroid main belt, the analytically computed proper elements are accurate to a level very close to the best result achievable by any analytical theory. This fundamental limitation results from the infinite web of resonances and because of the occurrence of chaotic motions. Still, there are some regions of the belt in which these proper elements are of degraded accuracy, thus preventing a reliable definition of asteroid families and detailed studies of the dynamical structure. We have used a different method to compute asteroid proper elements, following the approach introduced in the LONGSTOP project to describe the secular dynamics of the major outer planets. By applying purely numerical techniques, we produced so-called synthetic proper elements for a catalog of 10,256 asteroids with osculating semimajor axes between 2.5 and 4.0AU.The procedure consisted of simultaneous integration of asteroid and planetary orbits for 2Myr, with online filtering of the short-periodic perturbations. The output of the integration was spectrally resolved, and the principal harmonics (proper values) extracted from the time series. For each asteroid we have also tested the accuracy and stability in time of the proper elements, and estimated the maximum Lyapunov Characteristic Exponent to monitor the chaotic behaviors. This provided information on the reliability of the data for each orbit, in particular allowing to select 1,852 cases for an extended integration (10Myr) of the orbits showing instability. The results indicate that for more than half of the cases the proper elements have a time stability improved by more than a factor 3 with respect to the elements computed by the previous analytical theory. But of course there are also unstable cases for which the proper elements are less accurate and reliable, the extreme examples being 23 orbits exhibiting hyperbolic escape from the solar system. This form of escape from the asteroid belt could be responsible for a significant mass loss over the age of the solar system.     

OpenOrb: Open‐source asteroid orbit computation software including statistical ranging

Abstract— We are making an open‐source asteroid orbit computation software package called OpenOrb publicly available. OpenOrb is built on a well‐established Bayesian inversion theory, which means that it is to a large part complementary to orbit‐computation packages currently available. In particular, OpenOrb is the first package that contains tools for rigorously estimating the uncertainties resulting from the inverse problem of computing orbital elements using scarce astrometry. In addition to the well‐known least‐squares method, OpenOrb also contains both Monte‐Carlo (MC) and Markov‐Chain MC (MCMC; Oszkiewicz et al. [2009]) versions of the statistical ranging method. Ranging allows the user to obtain sampled, non‐Gaussian orbital‐element probability‐density functions and is therefore optimized for cases where the amount of astrometry is scarce or spans a relatively short time interval. Ranging‐based methods have successfully been applied to a variety of different problems such as rigorous ephemeris prediction, orbital element distribution studies for transneptunian objects, the computation of invariant collision probabilities between near‐Earth objects and the Earth, detection of linkages between astrometric asteroid observations within an apparition as well as between apparitions, and in the rigorous analysis of the impact of orbital arc length and/or astrometric uncertainty on the uncertainty of the resulting orbits. Tools for making ephemeris predictions and for classifying objects based on their orbits are also available in OpenOrb. As an example, we use OpenOrb in the search for candidate retrograde and/or high‐inclination objects similar to 2008 KV₄₂ in the known population of transneptunian objects that have an observational time span shorter than 30 days.     

The accuracy of proper orbital elements and the properties of asteroid families: Comparison with the linear theory

The accuracy and reliability of the proper orbital elements used to define asteroid families are investigated by simulating numerically the dynamical evolution of families assumed to arise from the “explosion” of a parent object. The orbits of the simulated family asteroids have then been integrated in the frame of the elliptic restricted three-body problem Sun-Jupiter-asteroid, for times of the order of the circulation periods of perihelia and nodes. By filtering out short-periodic perturbations, we have monitored the behavior of the proper eccentricities and inclinations, computed according to the linear secular perturbation theory. Significant long-period variations have been found especially for families having nonnegligible eccentricities and/or inclinations (like the Eos family), and strong disturbances due to the proximity of mean motion commensurabilities with Jupiter have been evidenced (for instance, in the case of the Themis family). These phenomena can cause a significant “noise” on the proper eccentricities and inclinations, probably affecting in some cases the derived family memberships. They can also give rise to a spurious anisotropy in the fragment ejection velocity fields computed from the dispersion in proper elements observed in each family, and this could explain the puzzling anisotropies of this kind actually found in real families by D. Brouwer (1951, Astron. J. 56, 9–32) and by V. Zappalà, P. Farinella, Z. Knežević, and P. Paolicchi (1984), Icarus 59, 261–285).     

The peculiar case of the Agnia asteroid family

The Agnia asteroid family, a cluster of asteroids located near semimajor axis a=2.79 AU, has experienced significant dynamical evolution over its lifetime. The family, which was likely created by the breakup of a diameter D∼50 km parent body, is almost entirely contained within the high-order secular resonance z₁. This means that unlike other families, Agnia's full extent in proper eccentricity and inclination is a byproduct of the large-amplitude resonant oscillations produced by this resonance. Using numerical integration methods, we found that the spread in orbital angles observed among Agnia family members would have taken at least 40 Myr to create; this sets a lower limit on the family's age. To determine the upper bound on Agnia's age, we used a Monte Carlo model to track how the small members in the family evolve in semimajor axis by Yarkovsky thermal forces. Our results indicate the family is no more than 140 Myr old, with a best-fit age of 100⁺³⁰₋₂₀ Myr. Using two independent methods, we also determined that the D∼5 km fragments were ejected from the family-forming event at a velocity near 15 m/s. This velocity is consistent with results from numerical hydrocode simulations of asteroid impacts and observations of other similarly sized asteroid families. Finally, we found that 57% of known Agnia fragments were initially prograde rotators. The reason for this limited asymmetry is unknown, though we suspect it is a fluke produced by the stochastic nature of asteroid disruption events.     

The origin of the Flora family

Evidence is presented indicating that the Flora family is of common origin. The distribution of proper elements and physical properties of Flora-family asteroids are compared with those of families believed to have formed from the catastrophic disruption of parent bodies. Differences in these orbital and physical properties suggest that the creation of the Flora family was more complex. Available evidence concerning the Flora family, together with recent models for the collisional evolution of the asteroids, suggests that this family may have originated from a binary or multiple asteroid. A mechanism in which the Flora family may have been produced by the disruption of a former major satellite of 8 Flora is presented and compared with other possible modes of formation.     

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.     

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.     

Yarkovsky footprints in the Eos family

The Eos asteroid family is the third most populous, after Themis and Koronis, and one of the largest non-random groups of asteroids in the main belt. It has been known and studied for decades, but its structure and history still presented difficulties to understand. We first revise the Eos family identification as a statistical cluster in the space of proper elements. Using the most to-date catalogue of proper elements we determine a nominal Eos family, defined by us using the hierarchical-clustering method with the cut-off velocity of 55 m/s, contains some 4400 members. This unforeseen increase in known Eos asteroids allows us to perform a much more detailed study than was possible so far. We show, in particular, that most of the previously thought peculiar features are explained within the following model: (i) collisional disruption of the parent body leads to formation of a compact family in the proper element space (with characteristic escape velocities of the observed asteroids of tens of meters per second, compatible with hydrocode simulations), and (ii) as time goes, the family dynamically evolves due to a combination of the thermal effects and planetary perturbations. This model allows us to explain sharp termination of the family at the J7/3 mean motion resonance with Jupiter, uneven distribution of family members about the J9/4 mean motion resonance with Jupiter, semimajor axis distribution of large vs small members in the family and anomalous residence of Eos members inside the high-order secular resonance z₁. Our dynamical method also allows us to estimate Eos family age to . Several formal members of the Eos family are in conflict with our model and these are suspected interlopers. We use spectroscopic observations, whose results are also reported here, and results of 5-color wide-band Sloan Digital Sky Survey photometry to prove some of them are indeed spectrally incompatible with the family.     

Yarkovsky/YORP chronology of asteroid families

Asteroid families are the byproducts of catastrophic collisions whose fragments form clusters in proper semimajor axis, eccentricity, and inclination space. Although many families have been observed in the main asteroid belt, only two very young families, Karin and Veritas, have well-determined ages. The ages of other families are needed, however, if we hope to infer information about their ejection velocity fields, space weathering processes, etc. In this paper, we developed a method that allows us to estimate the ages of moderately young asteroid families (approximately in between 0.1 and 1 Gyr). We apply it to four suitable cases—Erigone, Massalia, Merxia, and Astrid—and derive their likely ages and approximate ejection velocity fields. We find that Erigone and Merxia were produced by large catastrophic disruption events (i.e., parent body ?100 km) that occurred approximately 280 and 330 Myr ago, respectively. The Massalia family was likely produced by a cratering event on Asteroid (20) Massalia less than 200 Myr ago. Finally, the Astrid family, which was produced by the disruption of a 60-70 km asteroid, is 100-200 Myr old, though there is considerable uncertainty in this result. We estimate that the initial ejection velocities for these families were only a few tens of meters per second, consistent with numerical hydrocode models of asteroid impacts. Our results help to verify that asteroid families are constantly undergoing dynamical orbital evolution from thermal (Yarkovsky) forces and spin vector evolution from thermal (YORP) torques.     

A stability study of asteroid families near the 3:1 and 5:2 resonance with Jupiter

By using theD-criterion Lindblad (1992) has identified 14 asteroid families from a sample of 4100 numbered asteroids with proper elements from Milani and Kneevi (1990). Taxonomic types and other physical properties for a significant number of objects in five of the families show strong homogeneity within each family, further strengthening their internal relationship.To test the hypothesis of a common origin in, e.g., a catastrophic collision event, we have set out to integrate the orbits of the members of the Maria, Dora and Oppavia-Gefion families over some 10⁶ years. The mean distance for the Maria family is close to the 3:1 mean-motion resonance with Jupiter, while the other two families lie close to the 5:2 resonance.We used a simplified solar system model which included the perturbations by Jupiter and Saturn only and implemented Everhart's variable stepsize integrator RA15. All close encounters between the family members (within 0.1 AU) were recorded as well. Preliminary results from integrations over 4×10⁵ years are presented here.The statistics of close encounters show pronounced peaks for several members within each family, while for others no significant levels above the background of random encounters or even very low frequencies were found. This indicates a subclustering within the families. Quite a lot of very close (<0.005 AU) mutual encounters are found, which suggest that, at least for the larger members in a family, the mutual gravitational interactions could be of some importance for the real orbital evolutions.The encounter statistics between the Dora and Oppavia family members suggest a possible interrelationship between this two groups.