Rotational fission of contact binary asteroids

The energetics and dynamics of contact binary asteroids as they approach and pass the rotational fission limit is studied. We presume that the asteroids are subject to an external torque, such as from the YORP effect, that increases their angular momentum. Furthermore, we assume the asteroids can be described by a fairly minimal model comprised of a sphere and ellipsoid resting on each other. The minimum energy configurations for contact binary asteroids at different levels of angular momentum are computed and discussed. We find distinct transitions between different configurations as the angular momentum of the system is increased. These indicate that rapidly rotating contact binary asteroids may seek out clearly different relative configurations than slowly rotating systems. We find a single end state of the systems prior to rotational fission, and distinct dynamical outcomes as a function of mass distribution and shape when the rotational fission limit is exceeded. Our theoretical results agree qualitatively with observed properties of near-Earth asteroids, and can be used to help explain the spin-rate barrier, contact binaries, and the observed morphology of most NEO binaries.     

Tidal disruptions: A continuum theory for solid bodies

Although the theory of Roche 1847 for the tidal disruption limits of orbiting satellites assumes a fluid body, a length to diameter of exactly 2.07:1, and a particular body orientation, the theory is commonly applied to the satellites of the Solar System and to small asteroids and comets passing nearby a planet. Clearly these bodies are neither fluid nor generally are that elongated, so a more appropriate theory is needed. Here we present exact analytical results for the distortion and disruption limits of solid spinning ellipsoidal bodies subjected to tidal forces, using the Drucker-Prager strength model with zero cohesion. It is the appropriate model for dry granular materials such as sands and rocks, for rubble-pile asteroids and comets, and for larger satellites, asteroids and comets where the cohesion can be ignored. This study uses the same approach as the studies of spin limits for solid ellipsoidal bodies given in [Holsapple, K.A., 2001. Icarus 154, 432-448; Holsapple, K.A., 2004. Icarus 172, 272-303]. It is a static theory that predicts conditions for breakup and predicts the nature of the deformations at the limit state, but does not track the dynamics of the body as it comes apart. The strength is characterized by a single material parameter associated with an angle of friction, which can range from zero to 90°. The case with zero friction angle has no shear strength whatsoever, so it is then the model of a fluid or gas. The case of 90° represents a material that cannot fail in shear, but still has zero tensile strength. Typical dry soils have angles of friction of 30°-40°. Since the static fluid case is included in the theory as a special case, the classical results of Roche [Roche, E.A., 1847. Acad. Sci. Lett. Montpelier. Mem. Section Sci. 1, 243-262] and Jeans [Jeans, J.H., 1917. Mem. R. Astron. Soc. London 62, 1-48] are included and re-derived in their entirety; but the general solid case has much more variety and applicability. We consider both the spin-locked case, appropriate for most satellites of the Solar System; and the zero spin case, a possible case for a passing stray body. Detailed plots of many special cases are presented, in terms of shape, orientation and mass densities. A very typical result gives a closest approach d=1.5(ρ/ρP)^1/3R in terms of the planet radius R, and the satellite and planet mass densities ρ and ρP. We also use the theory to distinguish between conditions allowing global shape changes leading to new equilibrium states, or those leading to complete disruption. We apply the theory to the potentially hazardous Asteroid 99942 Apophis due to pass very near the Earth in 2029, and conclude it is extremely unlikely to experience any tidal readjustments during its passage. The states of many of the satellites of the Solar System are compared to the theory, and we find that all are well within their tidal disruption limits for expected values of the internal friction.     

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

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

Binary asteroid population: 1. Angular momentum content

We compiled a list of estimated parameters of binary systems among asteroids from near-Earth to trojan orbits. In this paper, we describe the construction of the list, and we present results of our study of angular momentum content in binary asteroids. The most abundant binary population is that of close binary systems among near-Earth, Mars-crossing, and main belt asteroids that have a primary diameter of about 10 km or smaller. They have a total angular momentum very close to, but not generally exceeding, the critical limit for a single body in a gravity regime. This suggests that they formed from parent bodies spinning at the critical rate (at the gravity spin limit for asteroids in the size range) by some sort of fission or mass shedding. The Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect is a candidate to be the dominant source of spin-up to instability. Gravitational interactions during close approaches to the terrestrial planets cannot be a primary mechanism of formation of the binaries, but it may affect properties of the NEA part of the binary population.     

Spin limits of Solar System bodies: From the small fast-rotators to 2003 EL61

Holsapple [Holsapple, K.A., 2001. Icarus 154, 432-448; Holsapple, K.A., 2004. Icarus 172, 272-303] determined the spin limits of bodies using a model for solid bodies without tensile or cohesive strength, but with the pressure-induced shear strengths characteristic of dry sands and gravels. That theory included the classical analyses for fluid bodies given by Maclaurin, Jacobi and others as a special case. For the general solid bodies, it was shown that there exists a very wide range of permissible shapes and spin limits; and explicit algebraic results for those limits were given. This paper gives an extension of those analyses to include geological-like materials that also have tensile and cohesive strength. Those strengths are necessary to explain the smaller, fast-rotating asteroids discovered in the last few years. I find that the spin limits for these more general solids have two limiting regimes: a strength regime for bodies with a diameter <3 km, and a gravity regime for the larger bodies with a diameter >10 km (which is the case covered by the earlier papers). I derive explicit algebraic forms for the dependence of the spin limits on shape, mass density and material strength properties. The comparison of the theory to the database for the spins of asteroids and trans-neptunian objects (TNO's) objects shows excellent agreement. For large bodies (diameter D>10 km), the presence of cohesive and/or tensile strength does not permit higher spin rates than would be allowed for rubble pile bodies. Thus, the fact that the spin rates of all large bodies is limited to periods greater than about 2 h does not imply that they are rubble piles. In contrast, for small bodies (D<10 km) the presence of even a very small amount of strength allows much more rapid spins. Small bodies might then be rubble piles but require a small amount of bonding. Finally, I make some remarks about the application of the theory to the TNO's and large asteroids, and question whether a common assumption by researchers that those bodies must take on relaxed fluid shapes is warranted. If not, then the densities and shapes required by that assumption are not valid. I use 2003 EL61 as a prime example.     

Dynamical passage to approximate equilibrium shapes for spinning, gravitating rubble asteroids

Many asteroids are thought to be particle aggregates held together principally by self-gravity. Here we study — for static and dynamical situations — the equilibrium shapes of spinning asteroids that are permitted for rubble piles. As in the case of spinning fluid masses, not all shapes are compatible with a granular rheology. We take the asteroid to always be an ellipsoid with an interior modeled as a rigid-plastic, cohesion-less material with a Drucker-Prager yield criterion. Using an approximate volume-averaged procedure, based on the classical method of moments, we investigate the dynamical process by which such objects may achieve equilibrium. We first collapse our dynamical approach to its statical limit to derive regions in spin-shape parameter space that allow equilibrium solutions to exist. At present, only a graphical illustration of these solutions for a prolate ellipsoid following the Drucker-Prager failure law is available [Sharma, I., Jenkins, J.T., Burns, J.A., 2005a. Bull. Am. Astron. Soc. 37, 643; Sharma, I., Jenkins, J.T., Burns, J.A., 2005b. Equilibrium shapes of ellipsoidal soil asteroids. In: García-Rojo, R., Hermann, H.J., McNamara, S. (Eds.), Proceedings of the 5th International Conference on Micromechanics of Granular Media, vol. 1. A.A. Balkema, UK; Holsapple, K.A., 2007. Icarus 187, 500-509]. Here, we obtain the equilibrium landscapes for general triaxial ellipsoids, as well as provide the requisite governing formulae. In addition, we demonstrate that it may be possible to better interpret the results of Richardson et al. [Richardson, D.C., Elankumaran, P., Sanderson, R.E., 2005. Icarus 173, 349-361] within the context of a Drucker-Prager material. The graphical result for prolate ellipsoids in the static limit is the same as those of Holsapple [Holsapple, K.A., 2007. Icarus 187, 500-509] because, when worked out, his final equations will match ours. This is because, though the formalisms to reach these expressions differ, in statics, at the lowest level of approximation, volume-averaging and the approach of Holsapple [Holsapple, K.A., 2007. Icarus 187, 500-509] coincide. We note that the approach applied here was obtained independently [Sharma, I., Jenkins, J.T., Burns, J.A., 2003. Bull. Am. Astron. Soc. 35, 1034; Sharma, I., 2004. Rotational Dynamics of Deformable Ellipsoids with Applications to Asteroids. Ph.D. thesis, Cornell University] and it provides a general, though approximate, framework that is amenable to systematic improvements and is flexible enough to incorporate the dynamical effects of a changing shape, different rheologies and complex rotational histories. To demonstrate our technique, we investigate the non-equilibrium dynamics of rigid-plastic, spinning, prolate asteroids to examine the simultaneous histories of shape and spin rate for rubble piles. We have succeeded in recovering most results of Richardson et al. [Richardson, D.C., Elankumaran, P., Sanderson, R.E., 2005. Icarus 173, 349-361], who obtained equilibrium shapes by studying numerically the passage into equilibrium of aggregates containing discrete, interacting, frictionless, spherical particles. Our mainly analytical approach aids in understanding and quantifying previous numerical simulations.     

Generalized YORP evolution: Onset of tumbling and new asymptotic states

Asteroids have a wide range of rotation states. While the majority spin a few times to several times each day in principal axis rotation, a small number spin so slowly that they have somehow managed to enter into a tumbling rotation state. Here we investigate whether the Yarkovsky-Radzievskii-O'Keefe-Paddack (YORP) thermal radiation effect could have produced these unusual spin states. To do this, we developed a Lie-Poisson integrator of the orbital and rotational motion of a model asteroid. Solar torques, YORP, and internal energy dissipation were included in our model. Using this code, we found that YORP can no longer drive the spin rates of bodies toward values infinitely close to zero. Instead, bodies losing too much rotation angular momentum fall into chaotic tumbling rotation states where the spin axis wanders randomly for some interval of time. Eventually, our model asteroids reach rotation states that approach regular motion of the spin axis in the body frame. An analytical model designed to describe this behavior does a good job of predicting how and when the onset of tumbling motion should take place. The question of whether a given asteroid will fall into a tumbling rotation state depends on the efficiency of its internal energy dissipation and on the precise way YORP modifies the spin rates of small bodies.     

Equilibrium figures of spinning bodies with self-gravity

The study of the equilibrium and stability of spinning ellipsoidal fluid bodies with gravity began with Newton in 1687, and continues to the present day. However, no smaller bodies of the Solar System are fluid. Here I model those bodies as elastic-plastic solids using a cohesionless Mohr-Coulomb yield envelope characterized by an angle of friction. This study began in Holsapple 2001. Here new closed-form algebraic formulas for the spin limits of ellipsoidal shapes are derived using an energy method. The fluid results of Maclaurin and Jacobi are again recovered as special cases. I then consider the stability of those equilibrium states. For elastic-plastic solids the common methods cannot be used, because the constitutive equations lack sufficient smoothness at the limiting plastic states. Therefore, I propose and study a new measure of the stability of dynamic processes in general bodies. An energy-based approach is introduced which is shown to include stability approaches used in the statics of nonlinear elastic and elastic-plastic bodies, spectral definitions and the Liapunov methods used for finite-dimensional dynamical systems. The method is applied to spinning, solid, strained bodies. In contrast to the special fluid case, it is found that the strain energy term of solid materials generally induces stability of all equilibrium shapes, except for two possible cases. First, strain softening in the elastic-plastic law can result in instability at the plastic limit spin. Second, a loss of shear stiffness can give unstable states at specific spins less than the limit equilibrium spins. In the latter case, a solid spinning ellipsoidal body without elastic shear stiffness can spin no faster than with a period of about 3.7 hr, else it will fail by shearing deformations. That is distinctly slower than the oft-quoted limit of 2.1 hr at which material would be flung off the equator by tensile forces. However, the final conclusion is that neither cohesion nor tensile strength is required for the shapes and spins of almost all of the larger observed asteroids: we cannot rule out rubble-pile structures.     

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.     

On the shapes and spins of “rubble pile” asteroids

We examine the shape of a “rubble pile” asteroid as it slowly gains angular momentum by YORP torque, to the point where “landsliding” occurs. We find that it evolves to a “top” shape with constant angle of repose from the equator up to mid-latitude, closely resembling the shapes of several nearly critically spinning asteroids imaged by radar, most notably (66391) 1999 KW4 [Ostro, S.J., Margot, J.-L., Benner, L.A.M., Giorgini, J.D., Scheeres, D.J., Fahnestock, E.G., Broschart, S.B., Bellerose, J., Nolan, M.C., Magri, C., Pravec, P., Scheirich, P., Rose, R., Jurgens, R.F., De Jong, E.M., Suzuki, S., 2006. Science 314, 1276-1280]. Similar calculations for non-spinning extremely prolate or oblate “rubble piles” show that even loose rubble can sustain shapes far from fluid equilibrium, thus inferences based on fluid equilibrium are generally useless for inferring bulk properties such as density of small bodies. We also investigate the tidal effects of a binary system with a “top shape” primary spinning at near the critical limit for stability. We find that very close to the stability limit, the tide from the secondary can actually levitate loose debris from the surface and re-deposit it, in a process we call “tidal saltation.” In the process, angular momentum is transferred from the primary spin to the satellite orbit, thus maintaining the equilibrium of near-critical spin as YORP continues to add angular momentum to the system. We note that this process is in fact dynamically related to the process of “shepherding” of narrow rings by neighboring satellites.     

Dynamics of rotationally fissioned asteroids: Source of observed small asteroid systems

We present a model of near-Earth asteroid (NEA) rotational fission and ensuing dynamics that describes the creation of synchronous binaries and all other observed NEA systems including: doubly synchronous binaries, high-e binaries, ternary systems, and contact binaries. Our model only presupposes the Yarkovsky-O’Keefe-Radzievskii-Paddack (YORP) effect, “rubble pile” asteroid geophysics, and gravitational interactions. The YORP effect torques a “rubble pile” asteroid until the asteroid reaches its fission spin limit and the components enter orbit about each other (Scheeres, D.J. [2007]. Icarus 189, 370-385). Non-spherical gravitational potentials couple the spin states to the orbit state and chaotically drive the system towards the observed asteroid classes along two evolutionary tracks primarily distinguished by mass ratio. Related to this is a new binary process termed secondary fission - the secondary asteroid of the binary system is rotationally accelerated via gravitational torques until it fissions, thus creating a chaotic ternary system. The initially chaotic binary can be stabilized to create a synchronous binary by components of the fissioned secondary asteroid impacting the primary asteroid, solar gravitational perturbations, and mutual body tides. These results emphasize the importance of the initial component size distribution and configuration within the parent asteroid. NEAs may go through multiple binary cycles and many YORP-induced rotational fissions during their approximately 10 Myr lifetime in the inner Solar System. Rotational fission and the ensuing dynamics are responsible for all NEA systems including the most commonly observed synchronous binaries.     

The equilibrium of rubble-pile satellites: The Darwin and Roche ellipsoids for gravitationally held granular aggregates

Many new small moons of the giant planets have been discovered recently. In parallel, satellites of several asteroids, e.g., Ida, have been found. Strikingly, a majority of these new-found planetary moons are estimated to have very low densities, which, along with their hypothesized accretionary origins, suggests a rubble internal structure. This, coupled to the fact that many asteroids are also thought to be particle aggregates held together principally by self-gravity, motivates the present investigation into the possible ellipsoidal shapes that a rubble-pile satellite may achieve as it orbits an aspherical primary. Conversely, knowledge of the shape will constrain the granular aggregate's orbit—the closer it gets to a primary, both primary's tidal effect and the satellite's spin are greater. We will assume that the primary body is sufficiently massive so as not to be influenced by the satellite. However, we will incorporate the primary's possible ellipsoidal shape, e.g., flattening at its poles in the case of a planet, and the proloidal shape of asteroids. In this, the present investigation is an extension of the first classical Darwin problem to granular aggregates. General equations defining an ellipsoidal rubble pile's equilibrium about an ellipsoidal primary are developed. They are then utilized to scrutinize the possible granular nature of small inner moons of the giant planets. It is found that most satellites satisfy constraints necessary to exist as equilibrated granular aggregates. Objects like Naiad, Metis and Adrastea appear to violate these limits, but in doing so, provide clues to their internal density and/or structure. We also recover the Roche limit for a granular satellite of a spherical primary, and employ it to study the martian satellites, Phobos and Deimos, as well as to make contact with earlier work of Davidsson [Davidsson, B., 2001. Icarus 149, 375–383]. The satellite's interior will be modeled as a rigid-plastic, cohesion-less material with a Drucker–Prager yield criterion. This rheology is a reasonable first model for rubble piles. We will employ an approximate volume-averaging procedure that is based on the classical method of moments, and is an extension of the virial method [Chandrasekhar, S., 1969. Ellipsoidal Figures of Equilibrium. Yale Univ. Press, New Haven] to granular solid bodies.     

YORP-Induced Long-Term Evolution of the Spin State of Small Asteroids and Meteoroids: Rubincam's Approximation

The rotation states of small asteroids and meteoroids are determined primarily by their collisions, gravitational torques due to the Sun and planets (in the case of close encounters), and internal dissipative effects (that relax the free-precession energy toward the fundamental state of principal-axis rotation). Rubincam has recently pointed out that thermal reemission on irregular-shaped bodies also results in a torque that may secularly change both the rotation rate and the orientation of the spin axis (the so-called YORP effect). Here we pursue investigation of this effect. Keeping the zero thermal-relaxation approximation of Rubincam and the assumption of the principal-axis rotation, we study the YORP effect both for precisely determined shapes of near-Earth asteroids and also for a large statistical sample of automatically generated shapes by the Gaussian-sphere technique of Muinonen. We find that the asymptotic state of the YORP evolution is characterized by an arbitrary value of the obliquity, with higher but nearly equal likelihood of 0°/180° and 90° states. At the adopted approximation, the most typical feature of this end state of the YORP evolution is secular deceleration of the rotation rate, which means that at some instant collisions will randomize the rotation state. In a minority of cases, the final state of the obliquity evolution leads to a permanent acceleration of the body's rotation, eventually resulting in rotational fission. The YORP-induced slow evolution may also play an important role in driving the rotation state of small asteroids toward the resonances between the forced precession due to the solar torque and perturbations of the orbital node and inclination. We find that for small Themis asteroids these resonances are isolated in the relevant range of frequencies, and the YORP evolving rotation may be either temporarily captured or rapidly jump across these resonances. In contrast, the possible values of the forced precession for small Flora asteroids may be resonant with clustered, nonisolated lines of the orbital perturbation. The individual rotation histories of small Flora asteroids may be thus very complicated and basically unpredictable. We comment on possible astronomical consequences of these results.     

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.     

Tidal disruptions: II. A continuum theory for solid bodies with strength, with applications to the Solar System

We present a comprehensive theory for the breakup conditions for ellipsoidal homogeneous secondary bodies subjected to the tidal forces from a nearby larger primary: for materials ranging from purely fluid ones, to granular rubble-pile gravel-like ones, and to those with either cohesive or granular strength including cohesive rocks and metals. The theory includes but greatly extends the classical analyses given by Roche in 1847, which dealt only with fluids, and also our previous analysis [Holsapple, K.A., Michel, P., 2006. Icarus 183, 331-348], which dealt only with solid but non-cohesive bodies. The results here give the distance inside of which breakup must occur, for both a steadily orbiting satellite and for a passing or impacting object. For the fluid bodies there is a single specific shape (a “Roche Ellipsoid”) that can be in equilibrium at any given distance from a primary, and especially only one shape that can exist at the overall minimum distance (d/R)(ρ/ρp)^1/3=2.455, the classical well-known “Roche limit.” In contrast, solid bodies can exist at a given distance from a primary with a range of shapes. Here we give multiple plots of the minimum distances for various important combinations of body shape, spin, mass density, and the strength parameters characterized by an angle of friction and cohesive strength. Such results can be used in different ways. They can be used to estimate limits on strengths and mass densities for orbiting bodies at a known distance and shape. They can be used to determine breakup distances for passing bodies with an assumed strength and shape. They can be used to constraint physical properties such as bulk density of bodies with a known shape that were known to breakup at a given distance. A collection of approximately 40 satellites of the Solar System is used for comparison to the theory. About half of those bodies are closer than the Roche fluid limit and must have some cohesion and/or friction angle to exist at their present orbital distance. The required solid strength for those states is determined. Finally, we apply the theory to the break up of the SL9 comet at close approach with Jupiter. Our results make clear that the literature estimates of its bulk density depend markedly on unknown parameters such as shape, orientation and spin, and most importantly, material strength characterization.     

Secular spin dynamics of inner main-belt asteroids

Understanding the evolution of asteroid spin states is challenging work, in part because asteroids have a variety of orbits, shapes, spin states, and collisional histories but also because they are strongly influenced by gravitational and non-gravitational (YORP) torques. Using efficient numerical models designed to investigate asteroid orbit and spin dynamics, we study here how several individual asteroids have had their spin states modified over time in response to these torques (i.e., 951 Gaspra, 60 Echo, 32 Pomona, 230 Athamantis, 105 Artemis). These test cases which sample semimajor axis and inclination space in the inner main belt, were chosen as probes into the large parameter space described above. The ultimate goal is to use these data to statistically characterize how all asteroids in the main belt population have reached their present-day spin states. We found that the spin dynamics of prograde-rotating asteroids in the inner main belt is generally less regular than that of the retrograde-rotating ones because of numerous overlapping secular spin-orbit resonances. These resonances strongly affect the spin histories of all bodies, while those of small asteroids (?40 km) are additionally influenced by YORP torques. In most cases, gravitational and non-gravitational torques cause asteroid spin axis orientations to vary widely over short (?1 My) timescales. Our results show that (951) Gaspra has a highly chaotic rotation state induced by an overlap of the s and s₆ spin-orbit resonances. This hinders our ability to investigate its past evolution and infer whether thermal torques have acted on Gaspra's spin axis since its origin.     

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.     

Effect of density inhomogeneity on YORP: The case of Itokawa

The effect of density inhomogeneity on the YORP effect for a given shape model is investigated. A density inhomogeneity will cause an offset between the center of figure and the center of mass and a re-orientation of the principal axes away from those associated with the shape alone. Both of these effects can alter the predicted YORP rate of change in angular velocity and obliquity. We apply these corrections to the Itokawa shape model and find that its YORP angular velocity rate is sensitive to offsets between its center of mass and center of figure, with a shift on the order of 15 m being able to change the sign of the YORP effect for that asteroid. Given the non-detection of YORP for Itokawa as of 2008, this can shed light on the density distribution within that body. The theory supports a shift of the asteroid center of mass towards Itokawa's neck region, where there is an accumulation of finer gravels, or towards the asteroid's “Head” region. Detection of the YORP effect for Itokawa should provide some strong constraints on its density distribution. This theory could also be applied to asteroids visited by future spacecraft to constrain density inhomogeneities.     

Angular momentum of binary asteroids: Implications for their possible origin

We describe in this work a thorough study of the physical and orbital characteristics of extensively observed main-belt and trojan binaries, mainly taken from the LAOSA (Large Adaptive Optics Survey of Asteroids [Marchis, F., Baek, M., Berthier, J., Descamps, P., Hestroffer, D., Kaasalainen, M., Vachier, F., 2006c. In: Workshop on Spacecraft Reconnaissance of Asteroid and Comet Interiors. Abstract #3042]) database, along with a selection of bifurcated objects. Dimensionless quantities, such as the specific angular momentum and the scaled primary spin rate, are computed and discussed for each system. They suggest that these asteroidal systems might be the outcome of rotational fission or mass shedding of a parent body presumably subjected to an external torque. One of the most striking features of separated binaries composed of a large primary (Rp>100 km) with a much smaller secondary (Rs<20 km) is that they all have total angular momentum of ∼0.27. This value is quite close to the Maclaurin-Jacobi bifurcation (0.308) of a spinning fluid body. Alternatively, contact binaries and tidally locked double asteroids, made of components of similar size, have an angular momentum larger than 0.48. They compare successfully with the fission equilibrium sequence of a rotating fluid mass. In conclusion, we find that total angular momentum is a useful proxy to assess the internal structure of such systems.