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.     

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.     

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.     

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.     

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.     

Binary asteroid population. 2. Anisotropic distribution of orbit poles of small,inner main-belt binaries

Our photometric observations of 18 main-belt binary systems in more than one apparition revealed a strikingly high number of 15 having positively re-observed mutual events in the return apparitions. Our simulations of the survey showed that it cannot be due to an observational selection effect and that the data strongly suggest that poles of mutual orbits between components of binary asteroids in the primary size range 3–8 km are not distributed randomly: The null hypothesis of an isotropic distribution of the orbit poles is rejected at a confidence level greater than 99.99%. Binary orbit poles concentrate at high ecliptic latitudes, within 30° of the poles of the ecliptic. We propose that the binary orbit poles oriented preferentially up/down-right are due to either of the two processes: (i) the YORP tilt of spin axes of their parent bodies toward the asymptotic states near obliquities 0° and 180° (pre-formation mechanism) or (ii) the YORP tilt of spin axes of the primary components of already formed binary systems toward the asymptotic states near obliquities 0° and 180° (post-formation mechanism). The alternative process of elimination of binaries with poles closer to the ecliptic by dynamical instability, such as the Kozai effect due to gravitational perturbations from the Sun, does not explain the observed orbit pole concentration. This is because for close binary asteroid systems, the gravitational effects of primary’s irregular shape dominate the solar-tide effect.     

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

Orbit Determination of Binary Asteroids

In addition to the detection of an asteroid moon or a binary asteroid, the knowledge of the satellite’s true orbit is of high importance to derive fundamental physical parameters of the binary system such as its mass and to shed light on its possible formation history and dynamical evolution (prograde/retrograde orbit, large/small eccentricity or inclination, etc.). A new methodology for preliminary orbit determination of binary asteroids – and visual binaries in general – is proposed. It is based on Thiele–Innes method combined with a ‘trial and error’ Monte-Carlo technique. This method provides the full set of solutions (bundle of orbits, with the 7 orbital elements) even for a reduced number of observations. The mass is a direct by-product of this orbit determination, from which one can next infer the bulk-density and porosity. In addition to the bundle of orbits, the method provides the marginal probability densities of the foreseen parameters. Such error analysis – since it avoids linear approximation – can be of importance for the prediction of the satellite’s position in the plane-of-sky during future stellar occultations or subsequent observations, but also for the analysis of the orbit’s secular evolution. After briefly describing the method, we present the algorithm and its application to some practical cases, with particular emphasis on asteroids binaries and applications on orbital evolution.     

The dynamical evolution of massive black hole binaries - II. Self-consistent N-body integrations

We use a hybrid N-body program to study the evolution of massive black hole binaries in the centers of galaxies, mainly to understand the factors affecting the binary eccentricity, the response of the galaxy to the binary merger, and the effect of loss-cone depletion on the merger time. The scattering experiments from paper I (Quinlan, 1996)[NewA, 1, 35] showed that the merger time is not sensitive to the eccentricity growth unless a binary forms with at least a moderate eccentricity. We find here that the eccentricity can become large under some conditions if a binary forms in a galaxy with a flat core or with a radial bias in its velocity distribution, especially if the dynamical friction is enhanced by resonances as suggested by Rauch & Tremaine (1996)[NewA, 1, 149]. But the necessary conditions seem unlikely, and our prediction from paper I remains unchanged: in most cases the eccentricity will start and remain small. The ejection of stars caused by the hardening of a binary may explain why large elliptical galaxies have weaker density cusps than smaller galaxies. If so, the central velocity distributions in those galaxies should have strong tangential anisotropies. The wandering of a binary from the center of a galaxy counteracts the effects of loss-cone depletion and helps the binary merge.     

Aegaeon (Saturn LIII), a G-ring object

Aegaeon (Saturn LIII, S/2008 S1) is a small satellite of Saturn that orbits within a bright arc of material near the inner edge of Saturn’s G-ring. This object was observed in 21 images with Cassini’s Narrow-Angle Camera between June 15 (DOY 166), 2007 and February 20 (DOY 051), 2009. If Aegaeon has similar surface scattering properties as other nearby small saturnian satellites (Pallene, Methone and Anthe), then its diameter is approximately 500 m. Orbit models based on numerical integrations of the full equations of motion show that Aegaeon’s orbital motion is strongly influenced by multiple resonances with Mimas. In particular, like the G-ring arc it inhabits, Aegaeon is trapped in the 7:6 corotation eccentricity resonance with Mimas. Aegaeon, Anthe and Methone therefore form a distinctive class of objects in the Saturn system: small moons in corotation eccentricity resonances with Mimas associated with arcs of debris. Comparisons among these different ring-arc systems reveal that Aegaeon’s orbit is closer to the exact resonance than Anthe’s and Methone’s orbits are. This could indicate that Aegaeon has undergone significant orbital evolution via its interactions with the other objects in its arc, which would be consistent with the evidence that Aegaeon’s mass is much smaller relative to the total mass in its arc than Anthe’s and Methone’s masses are.     

Equilibrium shapes of rubble-pile binaries: The Darwin ellipsoids for gravitationally held granular aggregates

Binaries are in vogue; many minor-planets like asteroids are being found to be binary or contact-binary systems. Even ternaries like 87 Sylvia have been discovered. The densities of these binaries are often estimated to be very low, and this, along with suspected accretionary origins, hints at a rubble interior. As in the case of fluid objects, a rubble-pile is unable to sustain all manners of spin, self-gravitation, and tidal interactions. This motivates the present study of the possible ellipsoidal shapes and mutual separations that members of a rubble-pile binary system may achieve. Conversely, knowledge of a granular binary’s shape and separation will constrain its internal structure - the ability of the binary’s members to sustain elongated shapes and/or maintain contact will hint at appreciable internal frictional strength. Because the binary’s members are allowed to be of comparable mass, the present investigation constitutes an extension of the second classical Darwin problem to granular aggregates.General equations defining the ellipsoidal rubble-pile binary system’s equilibrium are developed. These are then specialized to a pair of spin-locked, possibly unequal, prolate ellipsoidal granular aggregates aligned along their long axes. We observe that contact rubble-pile binaries can indeed exist. Further, depending on the binary’s geometry, an equilibrium contact binary’s members may, in fact, disrupt if separated. These results are applied to four suspected or known binaries: 216 Kleopatra, 25143 Itokawa, 624 Hektor and 90 Antiope. This exercise helps to bound the shapes and/or provide information about the interiors of these binaries.The binary’s interior will be modeled as a rigid-plastic, cohesionless material with a Drucker-Prager yield criterion. This rheology is a reasonable first model for rubble piles. We 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 University Press, New Haven, CT) to granular solid bodies. The present approach also helps us present an incrementally consistent approach to investigate the equilibrium shapes of fluid binaries, while highlighting the inconsistencies and errors inherent in the popular “Roche binary approximation”.     

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.     

Grooves on small saturnian satellites and other objects: Characteristics and significance

High-resolution images from the Cassini Imaging Science Subsystem (ISS) show parallel sets of grooves on Epimetheus and Pandora. Grooves have previously been observed on other satellites and asteroids, including Phobos, Gaspra, Ida, Eros, and minor occurrences on Phoebe. Sets of parallel grooves are so far observed only on satellites known or likely to be subject to significant tidal stresses, such as forced librations. Grooves on asteroids and on satellites not subject to significant forced librations occur in more globally disorganized patterns that may reflect impacts, varying internal structures, or even thermal stresses. The patterns and individual morphologies of grooves on the tidally-affected satellites suggest fracturing in weak materials due to tidal stresses and forced librations.     

Reassessing the origin of Triton

The origin of Neptune’s large, circular but retrograde satellite Triton has remained largely unexplained. There is an apparent consensus that its origin lies in it being captured, but until recently no successful capture mechanism has been found. Agnor and Hamilton (Agnor, C.B., Hamilton, D.P. [2006]. Nature 441, 192-194) demonstrated that the disruption of a trans-neptunian binary object which had Triton as a member, and which underwent a very close encounter with Neptune, was an effective mechanism to capture Triton while its former partner continued on a hyperbolic orbit. The subsequent evolution of Triton’s post-capture orbit to its current one could have proceeded through gravitational tides (Correia, A.C.M. [2009]. Astrophys. J. 704, L1-L4), during which time Triton was most likely semi-molten (McKinnon, W.B. [1984]. Nature 311, 355-358). However, to date, no study has been performed that considered both the capture and the subsequent tidal evolution. Thus it is attempted here with the use of numerical simulations. The study by Agnor and Hamilton (Agnor, C.B., Hamilton, D.P. [2006]. Nature 441, 192-194) is repeated in the framework of the Nice model (Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F. [2005]. Nature 435, 459-461) to determine the post-capture orbit of Triton. After capture Triton is then subjected to tidal evolution using the model of Mignard (Mignard, F. [1979]. Moon Planets 20, 301-315; Mignard, F. [1980]. Moon Planets 23, 185-201). The perturbations from the Sun and the figure of Neptune are included. The perturbations from the Sun acting on Triton just after its capture cause it to spend a long time in its high-eccentricity phase, usually of the order of 10 Myr, while the typical time to circularise to its current orbit is some 200 Myr, consistent with earlier studies. The current orbit of Triton is consistent with an origin through binary capture and tidal evolution, even though the model prefers Triton to be closer to Neptune than it is today. The probability of capturing Triton in this manner is approximately 0.7%. Since the capture of Triton was at most a 50% event - since only Neptune has one, but Uranus does not - we deduce that in the primordial trans-neptunian disc there were some 100 binaries with at least one Triton-sized member. Morbidelli et al. (Morbidelli, A., Levison, H.F., Bottke, W.F., Dones, L., Nesvorný, D. [2009]. Icarus 202, 310-315) concludes there were some 1000 Triton-sized bodies in the trans-neptunian proto-planetary disc, so the primordial binary fraction with at least one Triton-sized member is 10%. This value is consistent with theoretical predictions, but at the low end. If Triton was captured at the same time as Neptune’s irregular satellites, the far majority of these, including Nereid, would be lost. This suggests either that Triton was captured on an orbit with a small semi-major axisa ? 50R_N (a rare event), or that it was captured before the dynamical instability of the Nice model, or that some other mechanism was at play. The issue of keeping the irregular satellites remains unresolved.     

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.     

Discovery of a binary Centaur

We have identified a binary companion to (42355) 2002 CR₄₆ in our ongoing deep survey using the Hubble Space Telescope's High Resolution Camera. It is the first companion to be found around an object in a non-resonant orbit that crosses the orbits of giant planets. Objects in orbits of this kind, the Centaurs, have experienced repeated strong scattering with one or more giant planets and therefore the survival of binaries in this transient population has been in question. Monte Carlo simulations suggest, however, that binaries in (42355) 2002 CR₄₆-like heliocentric orbits have a high probability of survival for reasonable estimates of the binary's still-unknown system mass and separation. Because Centaurs are thought to be precursors to short period comets, the question of the existence of binary comets naturally arises; none has yet been definitively identified. The discovery of one binary in a sample of eight observed by HST suggests that binaries in this population may not be uncommon.     

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.     

The correlated colors of transneptunian binaries

We report resolved photometry of the primary and secondary components of 23 transneptunian binaries obtained with the Hubble Space Telescope. V-I colors of the components range from 0.7 to 1.5 with a median uncertainty of 0.06 magnitudes. The colors of the primaries and secondaries are correlated with a Spearman rank correlation probability of 99.99991%, 5 sigma for a normal distribution. Fits to the primary vs. secondary colors are identical to within measurement uncertainties. The color range of binaries as a group is indistinguishable from that of the larger population of apparently single transneptunian objects. Whatever mechanism produced the colors of apparently single TNOs acted equally on binary systems. The most likely explanation is that the colors of transneptunian objects and binaries alike are primordial and indicative of their origin in a locally homogeneous, globally heterogeneous protoplanetary disk.