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.     

Binary asteroid systems: Tidal end states and estimates of material properties

The locations of the fully despun, double synchronous end states of tidal evolution, where the rotation rates of both the primary and secondary components in a binary system synchronize with the mean motion about the center of mass, are derived for spherical components. For a given amount of scaled angular momentum J/J′, the tidal end states are over-plotted on a tidal evolution diagram in terms of mass ratio of the system and the component separation (semimajor axis in units of primary radii). Fully synchronous orbits may not exist for every combination of mass ratio and angular momentum; for example, equal-mass binary systems require J/J′ > 0.44. When fully synchronous orbits exist for prograde systems, tidal evolution naturally expands the orbit to the stable outer synchronous solution. The location of the unstable inner synchronous orbit is typically within two primary radii and often within the radius of the primary itself. With the exception of nearly equal-mass binaries, binary asteroid systems are in the midst of lengthy tidal evolutions, far from their fully synchronous tidal end states. Of those systems with unequal-mass components, few have even reached the stability limit that splits the fully synchronous orbit curves into unstable inner and stable outer solutions.Calculations of material strength based on limiting the tidal evolution time to the age of the Solar System indicate that binary asteroids in the main belt with 100-km-scale primary components are consistent with being made of monolithic or fractured rock as expected for binaries likely formed from sub-catastrophic impacts in the early Solar System. To tidally evolve in their dynamical lifetime, near-Earth binaries with km-scale primaries or smaller created via a spin-up mechanism must be much weaker mechanically than their main-belt counterparts even if formed in the main belt prior to injection into the near-Earth region. Small main-belt binaries, those having primary components less than 10 km in diameter, could bridge the gap between the large main-belt binaries and the near-Earth binaries, as, depending on the age of the systems, small main-belt binaries could either be as strong as the large main-belt binaries or as weak as the near-Earth binaries. The inherent uncertainty in the age of a binary system is the leading source of error in calculation of material properties, capable of affecting the product of rigidity μ and tidal dissipation function Q by orders of magnitude. Several other issues affecting the calculation of μQ are considered, though these typically affect the calculation by no more than a factor of two. We also find indirect evidence within all three groups of binary asteroids that the semimajor axis of the mutual orbit in a binary system may evolve via another mechanism (or mechanisms) in addition to tides with the binary YORP effect being a likely candidate.     

Simulation and analysis of the dynamics of binary near-Earth Asteroid (66391) 1999 KW4

Near-Earth Asteroid (66391) 1999 KW4 was the subject of the recently published first extensive radar imaging, shape and mutual orbit modeling, and physical and dynamical characterization of a binary asteroid. In this paper we present in detail our numerical simulation of KW4 behind that work. Our propagations of the system with some variation in estimated parameters cover the set of KW4's possible current dynamical states consistent with the body models and other information obtained directly from the observations. We also apply our implementation of this simulation capability to address some of the dynamical mechanisms by which KW4 may be moved into the more energetically excited of those possible current states, particularly solar gravity interaction. Through comparison of the results with certain features of the observation data, we conclude that the actual KW4 system is not in the most energetically relaxed configuration but must be moderately excited. The system occupies a generalized Cassini state 2 which is different from that considered in most previously published treatments of Cassini states in that it involves co-precession of the primary's spin frame and the mutual orbit rather than co-precession of a satellite's spin frame and that satellite's orbit about the primary. We present a simple analytical theory describing the system's dynamics, which should be applicable to any other binary systems, of which KW4 is representative, in which a massive, roughly oblate primary is spinning rapidly relative to the rate of its mutual orbit with an on-average synchronous, elongated secondary. We examine separately both the effect of the larger binary component's oblateness, and the effect of the smaller component's roughly triaxial ellipsoid shape. The simple analytical formulae obtained agree with full-detail numerical simulation results, and can be used for remote estimation of binary mass properties from observed system motion.     

Orbital evolution of small binary asteroids

About 15% of both near-Earth and main-belt asteroids with diameters below 10 km are now known to be binary. These small asteroid binaries are relatively uniform and typically contain a fast-spinning, flattened primary and a synchronously rotating, elongated secondary that is 20-40% as large (in diameter) as the primary. The principal formation mechanism for these binaries is now thought to be YORP (Yarkovsky-O’Keefe-Radzievskii-Paddack) effect induced spin-up of the primary followed by mass loss and accretion of the secondary from the released material. It has previously been suggested (?uk, M. [2007]. Astrophys. J. 659, L57-L60) that the present population of small binary asteroids is in a steady state between production through YORP and destruction through binary YORP (BYORP), which should increase or decrease secondary’s orbit, depending on the satellite’s shape. However, BYORP-driven evolution has not been directly modeled until now. Here we construct a simple numerical model of the binary’s orbital as well the secondary’s rotational dynamics which includes BYORP and selected terms representing main solar perturbations. We find that many secondaries should be vulnerable to chaotic rotation even for relatively low-eccentricity mutual orbits. We also find that the precession of the mutual orbit for typical small binary asteroids might be dominated by the perturbations from the prolate and librating secondary, rather than the oblate primary. When we evolve the mutual orbit by BYORP we find that the indirect effects on the binary’s eccentricity (through the coupling between the orbit and the secondary’s spin) dominate over direct ones caused by the BYORP acceleration. In particular, outward evolution causes eccentricity to increase and eventually triggers chaotic rotation of the secondary. We conclude that the most likely outcome will be reestablishing of the synchronous lock with a “flipped” secondary which would then evolve back in. For inward evolution we find an initial decrease of eccentricity and secondary’s librations, to be followed by later increase. We think that it is likely that various forms of dissipation we did not model may damp the secondary’s librations close to the primary, allowing for further inward evolution and a possible merger. We conclude that a merger or a tidal disruption of the secondary are the most likely outcomes of the BYORP evolution. Dissociation into heliocentric pairs by BYORP alone should be very difficult, and satellite loss might be restricted to the minority of systems containing more than one satellite at the time.     

Detailed prediction for the BYORP effect on binary near-Earth Asteroid (66391) 1999 KW4 and implications for the binary population

A previous theory by the authors for detailed modeling of the binary YORP effect is reviewed and expanded to accommodate doubly-synchronous binary systems, as well as a method for non-dimensionalizing the coefficients for application to binary systems where a shape model to compute its own coefficients is not available. The theory is also expanded to account for the effects of primary J₂ and the Sun’s 3rd body perturbation on the secular orbit evolution. The newly expanded theory is applied to the binary near-Earth Asteroid 1999 KW4, for which a detailed shape model is available. The result of simulation of the secular evolutionary equations shows that the KW4 orbit will be double in size in approximately 22,000 years, and will reach the Hill radius in approximately 54,000 years. The simulation also shows that the eccentricity will alternate growing and shrinking in magnitude, depending on the location of the solar node in the body-fixed frame. Therefore the eccentricity is not fixed to evolve in the opposite sign as the semi-major axis unless the circulation of the node (with a period of 500 years) is averaged out as well. The current orbit expansion rate for KW4 of 7 cm per year is shown to be detectable with observations of the mean anomaly which grows quadratically in time with an expanding orbit. Finally, the KW4 results are scaled for application to a number of other binary systems for which detailed shape models are not available. This application shows that the orbits considered can expand to their Hill radius in the range of 10⁴-10⁶ years. This implies rapid formation of binary systems is necessary to support the large percentage of binaries observed in the NEA population.     

On YORP-induced spin deformations of asteroids

The alteration of the spin states of small Solar-System bodies by the YORP thermal effect has recently become a plausible and, for some, the favorite candidate for the formation of binary asteroids. The idea is that if an asteroid is slowly spun up to a state where some strength measure is exceeded; it can no longer remain rigid and adjusts to a new configuration. Such a process might involve global fission, global shape changes without fission, or gradual surface mass loss with subsequent mass re-accumulations forming a secondary body.Here I analyze the changes in the shape, spin, and state during slowly increasing angular momentum of rubble-pile, self-gravitating, homogeneous ellipsoidal bodies undergoing homogeneous motions. I use, as appropriate for rubble-pile asteroids, the strength models of granular materials with zero tensile strength (cohesionless but arbitrary dilatancy); those are characterized by the “angle of friction” material constant. There are distinct limit spins depending on that angle of friction and the shape, which were previously presented [Holsapple, K.A., 2001. Icarus 154, 432-448; Holsapple, K.A., 2004. Icarus 172, 272-303]. Here the deformations and state changes when the angular momentum is slowly increased from that of a limit spin state are determined, to study the YORP processes. When a body is at its limit spin and the angular momentum increases further, the body deforms in a unique way along definite paths in the ellipsoidal shape space: it evolves as an elongating shape with an increasing rotational inertia, which in most cases produces a decreasing spin. I give exact analytical solutions for those shape and spin histories, as well as the histories of the mass density, angular momentum and energy. Comparison to other approaches is made.     

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.     

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.     

Effects of thermal radiation on the dynamics of binary NEAs

The Yarkovsky force, produced when thermal radiation is re-emitted asymmetrically, causes significant orbital evolution of asteroids in the 10 m-10 km size range. When acting on a non-spherical body, the momentum carried by this radiation generally produces a torque, called the YORP effect, which may be important in re-orienting asteroidal spins. Here we explore a related effect, the “binary YORP” (BYORP), that can modify the orbit of a synchronously rotating secondary in a binary system. It arises because a locked secondary is effectively an asymmetric appendage of the primary. It should be particularly important for Near-Earth Asteroids (NEAs) owing to their small sizes, proximity to the Sun, and benign collisional environment. To estimate BYORP's strength, we subjected 100 random Gaussian spheroids to the thermal radiation model of Rubincam (2000, Radiative spin-up and spin-down of small asteroids, Icarus, 148, 2-11). For most shapes, a significant torque arose on the secondary's orbit, typically modifying the orbit's size, eccentricity and inclination in less than 10⁵ years, for components of 1 and 0.3 km radii, separated by 2 km, at 1 AU, each of density 1750 kg m⁻³. Together YORP and BYORP are capable of synchronizing secondaries and circularizing orbits, making tidal dissipation unnecessary to explain the evolved state of observed NEA pairs. However, BYORP's rapid timescale poses a problem for the abundance of observed NEA binaries, since their formation rate is thought to be much slower. We consider and reject the following resolutions of this quandary: (i) the approximation using Gaussian spheroids inadequately models YORP; (ii) most secondaries are not synchronous, but inhabit other spin-orbit resonances (very unlikely); (iii) tidal dissipation is much more efficient than previously estimated, and thus capable of stabilizing observed systems; and (iv) moderately close encounters with planets can re-orient secondaries, turning BYORP into a slower, random-walk process. Finally, we speculate that most observed binary NEAs are in a stable state in which the obliquity-changing torques of YORP (acting on the primary) and BYORP cancel out, and that those systems must be close to 55° inclination, where the momentum-changing torques of both YORP and BYORP tend to be very small. Some retrograde systems might develop such that the nodes precess at a Sun-synchronous rate, while some prograde ones might move into the “evection” resonance. All three of these hypotheses can be tested directly by comparison with the i, Ω and ? observed for NEA binaries.     

Secular orbit variation due to solar radiation effects: a detailed model for BYORP

A detailed derivation of the effect of solar radiation pressure on the orbit of a body about a primary orbiting the Sun is given. The result is a set of secular equations that can be used for long-term predictions of changes in the orbit. Solar radiation pressure is modeled as a Fourier series in the body’s rotation state, where the coefficients are based on the shape and radiation properties of the body as parameters. In this work, the assumption is made that the body is in a synchronous orbit about the primary and rotates at a constant rate. This model is used to write explicit variational equations of the energy, eccentricity vector, and angular momentum vector for an orbiting body. Given that the effect of the solar radiation pressure and the orbit are periodic functions, they are readily averaged over an orbit. Furthermore, the equations can be averaged again over the orbit of the primary about the Sun to give secular equations for long-term prediction. This methodology is applied to both circular and elliptical orbits, and the full equations for secular changes to the orbit in both cases are presented. These results can be applied to natural systems, such as the binary asteroid system 1999 KW4, to predict their evolution due to the Binary YORP effect, or to artificial Earth orbiting, nadir-pointing satellites to enable more precise models for their orbital evolution.     

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.     

Orbital YORP and asteroid orbit evolution, with application to Apophis

Photon thrust from shape alone can produce quasi-secular changes in an asteroid's orbital elements. An asteroid in an elliptical orbit with a north–south shape asymmetry can steadily alter its elements over timescales longer than one orbital trip about the Sun. This thrust, called here orbital YORP (YORP = Yarkovsky–O'Keefe–Radzievskii–Paddack), operates even in the absence of thermal inertia, which the Yarkovsky effects require. However, unlike the Yarkovsky effects, which produce secular orbital changes over millions or billions of years, the change in an asteroid's orbital elements from orbital YORP operates only over the precession timescale of the orbit or of the asteroid's spin axis; this is generally only thousands or tens of thousands of years. Thus while the orbital YORP timescale is too short for an asteroid to secularly journey very far, it is long enough to warrant investigation with respect to 99942 Apophis, which might conceivably impact the Earth in 2036. A near-maximal orbital YORP effect is found by assuming Apophis is without thermal inertia and is shaped like a hemisphere, with its spin axis lying in the orbital plane. With these assumptions orbital YORP can change its along-track position by up to ±245 km, which is comparable to Yarkovsky effects. Though Apophis' shape, thermal properties, and spin axis orientation are currently unknown, the practical upper and lower limits are liable to be much less than the ±245 km extremes. Even so, the uncertainty in position is still likely to be much larger than the 0.5 km “keyhole” Apophis must pass through during its close approach in 2029 in order to strike the Earth in 2036.     

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.     

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.     

Extreme sensitivity of the YORP effect to small-scale topography

Radiation recoil (YORP) torques are shown to be extremely sensitive to small-scale surface topography, using numerical simulations. Starting from a set of “base objects” representative of the near-Earth object population, random realizations of three types of small-scale topography are added: Gaussian surface fluctuations, craters, and boulders. For each, the expected relative errors in the spin and obliquity components of the YORP torque caused by the observationally unresolved small-scale topography are computed. Gaussian power, at angular scales below an observational limit, produces expected errors of order 100% if observations constrain the surface to a spherical harmonic order l?10. For errors under 10%, the surface must be constrained to at least l=20. A single crater with diameter roughly half the object's mean radius, placed at random locations, results in expected errors of several tens of percent. The errors scale with crater diameter D as D² for D>0.3 and as D³ for D<0.3 mean radii. Objects that are identical except for the location of a single large crater can differ by factors of several in YORP torque, while being photometrically indistinguishable at the level of hundredths of a magnitude. Boulders placed randomly on identical base objects create torque errors roughly 3 times larger than do craters of the same diameter, with errors scaling as the square of the boulder diameter. A single boulder comparable to Yoshinodai on 25143 Itokawa, moved by as little as twice its own diameter, can alter the magnitude of the torque by factors of several, and change the sign of its spin component at all obliquities. Most of the total torque error produced by multiple unresolved craters is contributed by the handful of largest craters; but both large and small boulders contribute comparably to the total boulder-induced error. A YORP torque prediction derived from groundbased data can be expected to be in error by of order 100% due to unresolved topography. Small surface changes caused by slow spin-up or spin-down may have significant stochastic effects on the spin evolution of small bodies. For rotation periods between roughly 2 and 10 h, these unpredictable changes may reverse the sign of the YORP torque. Objects in this spin regime may random-walk up and down in spin rate before the rubble-pile limit is exceeded and fissioning or loss of surface objects occurs. Similar behavior may be expected at rotation rates approaching the limiting values for tensile-strength dominated objects.     

Radar observations and the shape of near-Earth asteroid 2008 EV5

We observed the near-Earth ASTEROID 2008 EV5 with the Arecibo and Goldstone planetary radars and the Very Long Baseline Array during December 2008. EV5 rotates retrograde and its overall shape is a 400 ± 50 m oblate spheroid. The most prominent surface feature is a ridge parallel to the asteroid’s equator that is broken by a concavity about 150 m in diameter. Otherwise the asteroid’s surface is notably smooth on decameter scales. EV5’s radar and optical albedos are consistent with either rocky or stony-iron composition. The equatorial ridge is similar to structure seen on the rubble-pile near-Earth asteroid (66391) 1999 KW4 and is consistent with YORP spin-up reconfiguring the asteroid in the past. We interpret the concavity as an impact crater. Shaking during the impact and later regolith redistribution may have erased smaller features, explaining the general lack of decameter-scale surface structure.     

Computing the effects of YORP on the spin rate distribution of the NEO population

The overall change of NEO spin rate due to planetary encounters and YORP is evaluated by using a Monte Carlo model. A large sample of test objects mimicking a source population is evolved over a timescale comparable with the Solar System age until they reach a steady state spin distribution that should reproduce the current NEO distribution. The spin change due to YORP is computed for each body according to a simplified model based on Scheeres [Scheeres, D.J., 2007a. Icarus 188, 430-450].The steady state cumulative distribution of NEO spin rates obtained from our simulation nicely reproduces the observed one, once our results are biased to match the diameter distribution of the sample of objects included in the observational database. The excellent agreement strongly suggests that YORP is responsible for the concentration of spin at low rotation rates. In fact, in the absence of YORP the steady state population significantly deviates from the observed one. The spin evolution due to YORP is also so rapid for NEOs that the initial rotation rate distribution of any source population is quickly relaxed to that of the observed population. This has profound consequences for the study of NEO origin since we cannot trace the sources of NEOs from their rotation rate only.     

Zero secular torque on asteroids from impinging solar photons in the YORP effect: A simple proof

YORP torques, where “YORP” stands for “Yarokovsky-O’Keefe-Radzievskii-Paddack,” arise mainly from sunlight reflected off a Solar System object and the infrared radiation emitted by it. We show here, through the most elementary demonstration that we can devise, that secular torques from impinging solar photons are generally negligible and thus cause little secular evolution of an asteroid’s obliquity or spin rate.     

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.     

The spin state of 433 Eros and its possible implications

In this paper, we show that Asteroid (433) Eros is currently residing in a spin-orbit resonance, with its spin axis undergoing a small-amplitude libration about the Cassini state 2 of the proper mode in the nonsingular orbital element sinI/2exp(?Ω), where I the orbital inclination and Ω the longitude of the node. The period of this libration is ?53.4 kyr. By excluding these libration wiggles, we find that Eros' pole precesses with the proper orbital plane in inertial space with a period of ?61.4 kyr. Eros' resonant state forces its obliquity to oscillate with a period of ?53.4 kyr between ?76° and ?89.5°. The observed value of ?89° places it near the latter extreme of this cycle. We have used these results to probe Eros' past orbit and spin evolution. Our computations suggest that Eros is unlikely to have achieved its current spin state by solar and planetary gravitational perturbations alone. We hypothesize that some dissipative process such as thermal torques (e.g., the so-called YORP effect) may be needed in our model to obtain a more satisfactory match with data. A detailed study of this problem is left for future work.