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.     

Spin rate distribution of small asteroids

The spin rate distribution of main belt/Mars crossing (MB/MC) asteroids with diameters 3-15 km is uniform in the range from f=1 to 9.5 d⁻¹, and there is an excess of slow rotators with f<1 d⁻¹. The observed distribution appears to be controlled by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. The magnitude of the excess of slow rotators is related to the residence time of slowed down asteroids in the excess and the rate of spin rate change outside the excess. We estimated a median YORP spin rate change of ≈0.022 d⁻¹/Myr for asteroids in our sample (i.e., a median time in which the spin rate changes by 1 d⁻¹ is ≈45 Myr), thus the residence time of slowed down asteroids in the excess is ≈110 Myr. The spin rate distribution of near-Earth asteroids (NEAs) with sizes in the range 0.2-3 km (∼5 times smaller in median diameter than the MB/MC asteroids sample) shows a similar excess of slow rotators, but there is also a concentration of NEAs at fast spin rates with f=9-10 d⁻¹. The concentration at fast spin rates is correlated with a narrower distribution of spin rates of primaries of binary systems among NEAs; the difference may be due to the apparently more evolved population of binaries among MB/MC asteroids.     

Photometric survey of binary near-Earth asteroids

Photometric data on 17 binary near-Earth asteroids (15 of them are certain detections, two are probables) were analysed and characteristic properties of the near-Earth asteroid (NEA) binary population were inferred. We have found that binary systems with a secondary-to-primary mean diameter ratio Ds/Dp?0.18 concentrate among NEAs smaller than 2 km in diameter; the abundance of such binaries decreases significantly among larger NEAs. Secondaries show an upper size limit of Ds=0.5-1 km. Systems with Ds/Dp?0.5 are abundant but larger satellites are significantly less common. Primaries have spheroidal shapes and they rotate rapidly, with periods concentrating between 2.2 to 2.8 h and with a tail of the distribution up to ∼4 h. The fast rotators are close to the critical spin for rubble piles with bulk densities about 2 g/cm³. Orbital periods show an apparent cut-off at Porb∼11 h; closer systems with shorter orbital periods have not been discovered, which is consistent with the Roche limit for strengthless bodies. Secondaries are more elongated on average than primaries. Most, but not all, of their rotations appear to be synchronized with the orbital motion; nonsynchronous secondary rotations may occur especially among wider systems with Porb>20 h. The specific total angular momentum of most of the binary systems is similar to within ±20% and close to the angular momentum of a sphere with the same total mass and density, rotating at the disruption limit; this suggests that the binaries were created by mechanism(s) related to rotation near the critical limit and that they neither gained nor lost significant amounts of angular momentum during or since formation. A comparison with six small asynchronous binaries detected in the main belt of asteroids suggests that the population extends beyond the region of terrestrial planets, but with characteristics shifted to larger sizes and longer periods. The estimated mean proportion of binaries with Ds/Dp?0.18 among NEAs larger than 0.3 km is 15±4%. Among fastest rotating NEAs larger than 0.3 km with periods between 2.2 and 2.8 h, the mean proportion of such binaries is (66⁺¹⁰₋₁₂)%.     

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.     

The Yarkovsky-driven origin of near-Earth asteroids

We investigate the relevance of the Yarkovsky effect for the origin of kilometer and multikilometer near-Earth asteroids (NEAs). The Yarkovsky effect causes a slow migration in semimajor axis of main belt asteroids, some of which are therefore captured into powerful resonances and transported to the NEA space. With an innovative simulation scheme, we determine that in the current steady-state situation 100-160 bodies with H < 18 (roughly larger than 1 km) enter the 3/1 resonance per million years and 40-60 enter the ν₆ resonance. The ranges are due to uncertainties on relevant simulation parameters such as the time scales for collisional disruption and reorientation, their size dependence, and the strength of the Yarkovsky and YORP effects. These flux rates to the resonances are consistent with those independently derived by Bottke et al. (2002, Icarus 156, 399-433) with considerations based only on the NEA orbital distribution and dynamical lifetime. Our results have been obtained assuming that the main belt contains 1,300,000 asteroids with H < 18 and linearly scale with this number. Assuming that the cumulative magnitude distribution of main belt asteroids is N(< H) ∝ 10^γ′H with γ′ = 0.25 in the 15.5 < H < 18 range (consistent with the results of the SDSS survey), we obtain that the bodies captured into the resonances should have a similar magnitude distribution, but with exponent coefficient γ = 0.33-0.40. The lowest value is obtained taking into account the YORP effect, while higher values correspond to a weakened YORP or to YORP-less cases. These values of γ are all compatible with the debiased magnitude distributions of the NEAs according to Rabinowitz et al. (2000, Nature 403, 165-166), Bottke et al. (2000b, Science 288, 2190-2194), and Stuart (2001, Science 294, 1691-1693). Hence the Yarkovsky and YORP effects allow us to understand why the magnitude distribution of NEAs is only moderately steeper than that of the main belt population. The steepest main belt distribution that would still be compatible with the NEA distribution has exponent coefficient γ′ ∼ 0.3.     

A steady-state model of NEA binaries formed by tidal disruption of gravitational aggregates

We present results of a simulation of a steady-state binary near-Earth asteroid (NEA) population. This study combines previous work on tidal disruption of gravitational aggregates [Walsh, K.J., Richardson, D.C., 2006. Icarus 180, 201-216] with a Monte Carlo simulation of NEA planetary encounters. Evolutionary effects include tidal evolution and binary disruption from close planetary encounters. The results show that with the best known progenitor (small Main Belt asteroids) shape and spin distributions, and current estimates of NEA lifetime and encounter probabilities, that tidal disruption should account for approximately 1-2% of NEAs being binaries. Given the best observed estimate of a ∼15% binary NEA fraction, we conclude that there are other formation mechanisms that contribute significantly to this population. We also present the expected distribution of binary orbital and physical properties for the steady-state binary NEAs formed by tidal disruption. We discuss the effects on binary fraction and properties due to changes in the least constrained parameters, and other possible effects on our model that could account for differences between the presented results and the observed binary population. Finally, we model possible effects of a significant population of binaries migrating to the near-Earth population from the Main Belt.     

The cool surfaces of binary near-Earth asteroids

Here we show results from thermal-infrared observations of km-sized binary near-Earth asteroids (NEAs). We combine previously published thermal properties for NEAs with newly derived values for three binary NEAs. The η value derived from the near-Earth asteroid thermal model (NEATM) for each object is then used to estimate an average thermal inertia for the population of binary NEAs and compared against similar estimates for the population of non-binaries. We find that these objects have, in general, surface temperatures cooler than the average values for non-binary NEAs as suggested by elevated η values. We discuss how this may be evidence of higher-than-average surface thermal inertia. This latter physical parameter is a sensitive indicator of the presence or absence of regolith: bodies covered with fine regolith, such as the Earth’s moon, have low thermal inertia, whereas a surface with little or no regolith displays high thermal inertia. Our results are suggestive of a binary formation mechanism capable of altering surface properties, possibly removing regolith: an obvious candidate is the YORP effect.We present also newly determined sizes and geometric visible albedos derived from thermal-infrared observations of three binary NEAs: (5381) Sekhmet, (153591) 2001 SN₂₆₃, and (164121) 2003 YT₁. The diameters of these asteroids are 1.41 ± 0.21 km, 1.56 ± 0.31 km, and 2.63 ± 0.40 km, respectively. Their albedos are 0.23 ± 0.13, 0.24 ± 0.16, and 0.048 ± 0.015, respectively.     

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.     

Do planetary encounters reset surfaces of near Earth asteroids?

Processes such as the solar wind sputtering and micrometeorite impacts can modify optical properties of surfaces of airless bodies. This explains why spectra of the main belt asteroids, exposed to these ‘space weathering’ processes over eons, do not match the laboratory spectra of ordinary chondrite (OC) meteorites. In contrast, an important fraction of Near Earth Asteroids (NEAs), defined as Q-types in the asteroid taxonomy, display spectral attributes that are a good match to OCs. Here we study the possibility that the Q-type NEAs underwent recent encounters with the terrestrial planets and that the tidal gravity (or other effects) during these encounters exposed fresh OC material on the surface (thus giving it the Q-type spectral properties). We used numerical integrations to determine the statistics of encounters of NEAs to planets. The results were used to calculate the fraction and orbital distribution of Q-type asteroids expected in the model as a function of the space weathering timescale, t_sw (see main text for definition), and maximum distance, r^∗, at which planetary encounters can reset the surface. We found that t_sw ∼ 10⁶ yr (at 1 AU) and r^∗ ∼ 5R_pl, where R_pl is the planetary radius, best fit the data. Values t_sw < 10⁵ yr would require that r^∗ > 20R_pl, which is probably implausible because these very distant encounters should be irrelevant. Also, the fraction of Q-type NEAs would be probably much larger than the one observed if t_sw > 10⁷ yr. We found that t_sw ∝ q², where q is the perihelion distance, expected if the solar wind sputtering controls t_sw, provides a better match to the orbital distribution of Q-type NEAs than models with fixed t_sw. We also discuss how the Earth magnetosphere and radiation effects such as YORP can influence the spectral properties of NEAs.     

Axial rotation of near-earth asteroids: The influence of the YORP effect

The distribution of axial rotation velocities of near-Earth asteroids (NEAs) substantially differs from that of the Main-Belt asteroids by an excess of both quickly and slowly rotating objects. Among the possible causes of this difference is the influence of the solar radiation—the so-called YORP effect—that arises from the absorption of solar energy and its reemission in the thermal range by a rotating body of irregular shape. It is known that the magnitude of this effect depends on the asteroid size and the quantity of received solar energy (the insolation). Analysis of the observational data showed that the mean diameter of NEAs decreases from the middle of the distribution to the edges, i.e., the excess of both slowly (ω ≤ 2 rev/day) and quickly (ω = 8–11 rev/day) rotating objects is formed due to the asteroids with sizes smaller than those in the middle of the distribution, which agrees well with the influence of the YORP effect. Moreover, the dependence of the axial rotation velocity of NEAs on the relative insolation shows that, for the NEAs referred to, both excesses are found in orbits where, on average, they receive 8–10% more solar energy than the NEAs in the middle of the distribution. This result also agrees with the character of the influence of the YORP effect and can be considered as an additional argument in its support. Thus, the study showed that one can infer that the currently available observational data suggest the possible influence of the YORP effect on the axial rotation of the near-Earth asteroids having sizes of D ～ 2 km and less. This is the first attempt to find the influence of the YORP effect on the axial rotation of the NEA family as a whole.     

Influence of the YORP effect on rotation rates of near‐Earth asteroids

The distribution of near‐Earth asteroid (NEA) rotation rates differs considerably from the similar distribution of Main Belt asteroids (MBAs) by the presence of excesses of fast and slow rotators, which are not observed or not so prominent in the distribution for MBAs. Among possible reasons for the difference, there can be influence of solar radiation on spin rate of small NEAs, the so‐called “YORP effect,” which appears due to reflection, absorption, and IR re‐emission of the sunlight by an irregularly shaped rotating asteroid. It is known that the YORP‐effect action strongly depends on the amount of solar energy obtained by the body (insolation), its size, and albedo. The analysis of observation data has shown that: (1) the mean diameter of NEAs decreases from the middle of the distribution to its ends, that is, the excesses of slow rotators (ω ≤ 2 rev day^?1) and fast rotators (ω ≥ 8 rev day^?1) are composed of smaller NEAs than in the middle of the distribution; (2) NEAs of both excesses are in the orbits where their insolation is about 8–10% larger than that of NEAs in the middle of the distribution; and (3) the objects in both excesses have a little lower albedo on average than that of objects in the middle of the distribution. All these results qualitatively agree well with the YORP‐effect action and may be considered as independent arguments in favor of it.     

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.     

Thermal inertia of near-Earth asteroids and implications for the magnitude of the Yarkovsky effect

Thermal inertia determines the temperature distribution over the surface of an asteroid and therefore governs the magnitude the Yarkovsky effect. The latter causes gradual drifting of the orbits of km-sized asteroids and plays an important role in the delivery of near-Earth asteroids (NEAs) from the main belt and in the dynamical spreading of asteroid families. At present, very little is known about the thermal inertia of asteroids in the km size range. Here we show that the average thermal inertia of a sample of NEAs in the km-size range is . Furthermore, we identify a trend of increasing thermal inertia with decreasing asteroid diameter, D. This indicates that the dependence of the drift rate of the orbital semimajor axis on the size of asteroids due to the Yarkovsky effect is a more complex function than the generally adopted D⁻¹ dependence, and that the size distribution of objects injected by Yarkovsky-driven orbital mobility into the NEA source regions is less skewed to smaller sizes than generally assumed. We discuss how this fact may help to explain the small difference in the slope of the size distribution of km-sized NEAs and main-belt asteroids.     

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.     

Radar and infrared observations of binary near-Earth Asteroid 2002 CE26

We observed near-Earth Asteroid (NEA) 2002 CE26 in August and September 2004 using the Arecibo S-band (2380-MHz, 12.6-cm) radar and NASA's Infrared Telescope Facility (IRTF). Shape models obtained based on inversion of our delay-Doppler images show the asteroid to be 3.5±0.4 km in diameter and spheroidal; our corresponding nominal estimates of its visual and radar albedos are 0.07 and 0.24, respectively. Our IRTF spectrum shows the asteroid to be C-class with no evidence of hydration. Thermal models from the IRTF data provide a size and visual albedo consistent with the radar-derived estimate. We estimate the spin-pole to be within a few tens of degrees of λ=317°, β=−20°. Our radar observations reveal a secondary approximately 0.3 km in diameter, giving this binary one of the largest size differentials of any known NEA. The secondary is in a near-circular orbit with period 15.6±0.1 h and a semi-major axis of 4.7±0.2 km. Estimates of the binary orbital pole and secondary rotation rate are consistent with the secondary being in a spin-locked equatorial orbit. The orbit corresponds to a primary mass of M=1.95±0.25×¹⁰¹³ kg, leading to a primary bulk density of , one of the lowest values yet measured for a main-belt or near-Earth asteroid.     

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.     

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.     

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.     

Binary near-Earth asteroid formation: Rubble pile model of tidal disruptions

We present numerical simulations of near-Earth asteroid (NEA) tidal disruption resulting in bound, mutually orbiting systems. Using a rubble pile model we have constrained the relative likelihoods for possible physical and dynamical properties of the binaries created. Overall 110,500 simulations were run, with each body consisting of ∼1000 particles. The encounter parameters of close approach distance and velocity were varied, as were the bodies' spin, elongation and spin axis direction. The binary production rate increases for closer encounters, at lower speeds, for more elongated bodies, and for bodies with greater spin. The semimajor axes for resultant binaries are peaked between 5 to 20 primary radii, and there is an overall trend for high eccentricity, with 97% of binaries having e > 0.1. The secondary-to-primary size ratios of the simulated binaries are peaked between 0.1 and 0.2, similar to trends among observed asteroid binaries. The spin rates of the primary bodies are narrowly distributed between 3.5- and 6-h periods, whereas the secondaries' periods are more evenly distributed and can exceed 15-h periods. The spin axes of the primary bodies are very closely aligned with the angular momenta of the binary orbits, whereas the secondary spin axes are nearly random. The shapes of the primaries show a large distribution of axis ratios, where those with low elongation (ratio of long and short axis) are both oblate and prolate, and nearly all with large elongation are prolate. This work presents results that suggest tidal disruption of gravitational aggregates can make binaries physically similar to those currently observed in the NEA population. As well, tidal disruption may create an equal number of binaries with qualities different from those observed, mostly binaries with large separation and with elongated primaries.     

Radar observations of E-class Asteroids 44 Nysa and 434 Hungaria

We observed the E-class main-belt Asteroids (MBAs) 44 Nysa and 434 Hungaria with Arecibo Observatory's S-band (12.6 cm) radar. Both asteroids exhibit polarization ratios higher than those measured for any other MBA: Nysa, μc=0.50±0.02 and Hungaria, μc=0.8±0.1. This is consistent with the high polarization ratios measured for every E-class near-Earth asteroid (NEA) observed by Benner et al. [Benner, L.A.M., and 10 collegues, 2008. Icarus, submitted for publication] and suggests a common cause. Our estimates of radar albedo are 0.19±0.06 for Nysa and 0.22±0.06 for Hungaria. These values are higher than those of most MBAs and, when combined with their high polarization ratios, suggest that the surface bulk density of both asteroids is high. We model Nysa as an ellipsoid of dimension 113×67×65 km (±15%) giving an effective diameter Deff=79±10 km, consistent with previous estimates. The echo waveforms are not consistent with a contact binary as suggested by Kaasalainen et al. [Kaasalainen, M., Torppa, J., Piironen, J., 2002. Astron. Astrophys. 383, L19-L22]. We place a constraint on Hungaria's maximum diameter, Dmax?11 km consistent with previous size estimates.