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

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

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

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

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.     

Photometry of Aten asteroids—More than a handful of binaries

Photometric observations were conducted on eight Aten near-Earth asteroids, with the goal of building physical models for the objects (85989) 1999 JD₆, (86450) 2000 CK₃₃, (86667) 2000 FO₁₀, (137170) 1999 HF₁, 1999 MN, 2000 PJ₅, 2002 JC and 2003 NZ₆. The results show rotation periods from 2.3 to almost 26 h. Some objects exhibit amplitudes higher than one magnitude on their lightcurves (1999 JD₆, 2000 CK₃₃ and 2003 NZ₆). Phase curve values (β, H^′, see below) were derived for four Atens and H-G values were found for two. Five Atens were classified by using their B-V and V-R color indices. This taxonomy was compared with the phase coefficient-albedo correlation defined by Belskaya and Shevchenko [Belskaya, I.N., Shevchenko, V.G., 2000. Icarus 147, 94-105]. Color variations during rotation and phase angle change were searched for. Our study demonstrates the high variety among Atens. Five out of the eight Atens are binaries or possible binaries, a significantly higher fraction than the expected 15% [Bottke, W.F., Melosh, H.J., 1996. Nature 281, 51-53]. The lightcurve of 2000 PJ₅ exhibits a binary character with a probable highly eccentric orbital rotation of the secondary component. The different periods of the known binary 1999 HF₁'s are easily detected. Other Atens have lightcurve with features such as high amplitude, V-shaped minima and U-shaped maxima that can be interpreted as a binary asteroid with a synchronous rotation (2003 NZ₆, 2000 CK₃₃ and perhaps also 1999 JD₆). The very red colors of 2000 CK₃₃ suggest a unique surface composition for this near-Earth object.     

Linear stability of collinear equilibrium points around an asteroid as a two-connected-mass: Application to fast rotating Asteroid 2000EB14

This note discusses the stability of collinear equilibrium points around a rotating system composed of two masses rigidly connected by a massless rod in the case, where the centripetal force outweighs the gravitational force. It is found that a stable region appears at L₁ when the ratio of gravitational to centripetal acceleration is less than 0.125, and that there is always no stable area at L₂ and L₃; the result is applied to the fast rotating Asteroid 2000EB₁₄.     

Orbital and thermal evolutions of four potential targets for a sample return space mission to a primitive near-Earth asteroid

In this paper, we present our study of the orbital and thermal evolutions, due to solar radiative heating, of four near-Earth asteroids (NEAs) considered as potential target candidates for sample return space missions to primitive asteroids. We used a dynamical model of the NEA population to estimate the most likely source region and orbital history of these objects. Then, for each asteroid, we integrated numerically over their entire lifetime a set of 14 initially indistinguishable orbit (clones), obtained by small variations of the nominal initial conditions. Using a thermal model, we then computed surface and sub-surface temperatures of these bodies during their dynamical history. Our aim is to determine whether these bodies are likely to have experienced high temperature level, and whether great temperature changes can be expected due to the orbital changes as well as their maximum and minimum values. Such information is important in the framework of sample return space missions whose goal is to bring back pristine materials. The knowledge of the temperature range of materials at different depth over the orbital evolution of potential targets can help defining sampling strategies that ensure the likelihood that unaltered material will be brought back. Our results suggest that for all the considered potential targets, the surface has experienced for some time temperatures greater than 400 K and at most 500 K with 50% probability. This probability drops rapidly with increasing temperature. Sub-surface materials at a depth of only 3 cm are much more protected from high temperature and generally do not reach temperatures exceeding 450 K (with 50% probability). They should thus be unaltered at this depth at least from a Sun-driven heating point of view. On the other hand, surface material for some of the considered objects can have a range of temperature which can make them less reliable as pristine materials. However, it is assumed here that the same material is constantly exposed to solar heat, while regolith turnover may occur. The latter can be caused by different processes such as seismic shaking and/or impact cratering. This would reduce the total time that materials are exposed to a certain temperature. Thus, it is very likely that a sample collected from any of the four considered targets, or any primitive NEA with similar dynamical properties, will have components that will be thermally unaltered as long as some of it comes from only 3 to 5 cm depth. Such a depth is not considered difficult to reach with some of the current designs of sampling devices.     

The orbit and size distribution of small Solar System objects orbiting the Sun interior to the Earth's orbit

We present the first observational measurement of the orbit and size distribution of small Solar System objects whose orbits are wholly interior to the Earth's (Inner Earth Objects, IEOs, with aphelion <0.983 AU). We show that we are able to model the detections of near-Earth objects (NEO) by the Catalina Sky Survey (CSS) using a detailed parameterization of the CSS survey cadence and detection efficiencies as implemented within the Jedicke et al. [Jedicke, R., Morbidelli, A., Spahr, T., Petit, J.M., Bottke, W.F., 2003. Icarus 161, 17-33] survey simulator and utilizing the Bottke et al. [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] model of the NEO population's size and orbit distribution. We then show that the CSS detections of 4 IEOs are consistent with the Bottke et al. [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] IEO model. Observational selection effects for the IEOs discovered by the CSS were then determined using the survey simulator in order to calculate the corrected number and H distribution of the IEOs. The actual number of IEOs with H<18 (21) is 36±26 (530±240) and the slope of the H magnitude distribution (∝¹⁰αH) for the IEOs is . The slope is consistent with previous measurements for the NEO population of αNEO=0.35±0.02 [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.-M., Levison, H.F., Michel, P., Metcalfe, T.S., 2002. Icarus 156, 399-433] and αNEO=0.39±0.013 [Stuart, J.S., Binzel, R.P., 2004. Icarus 170, 295-311]. Based on the agreement between the predicted and observed IEO orbit and absolute magnitude distributions there is no indication of any non-gravitational effects (e.g. Yarkovsky, tidal disruption) affecting the known IEO population.     

On the detection of the Yarkovsky effect on near-Earth asteroids by means of Gaia

The Yarkovsky effect, which causes a slow drifting of the orbital elements (mainly the semimajor axis) of km-sized asteroids and meteors, is the weak non-gravitational force experienced by these bodies due to the emission of thermal photons. This effect is believed to play a role in the delivery of near-Earth asteroids (NEAs) from the main belt, in the spreading of the orbital elements of asteroid families, and in the orbital evolution of potentially hazardous asteroids.Here we present preliminary results of simulationing indicating that the perturbations induced by the Yarkovsky effect on the positions of some tens of NEAs can be observed by means of the high-precision astrometric observations that will be provided by the ESA mission Gaia.     

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.     

Interacting ellipsoids: a minimal model for the dynamics of rubble-pile bodies

A simple mechanical model is formulated to study the dynamics of rubble-pile asteroids, formed by the gravitational re-accumulation of fragments after the collisional breakup of a parent body. In this model, a rubble-pile consists of N interacting fragments represented by rigid ellipsoids, and the equations of motion explicitly incorporate the minimal degrees of freedom necessary to describe the attitude and rotational state of each fragment. In spite of its simplicity, our numerical examples indicate that the overall behavior of our model is in line with several known properties of collisional events, like the energy and angular momentum partition during high velocity impacts. Therefore, it may be considered as a well defined minimal model.     

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.     

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.     

Dynamical evolution of Earth’s quasi-satellites: 2004 GU9 and 2006 FV35

We study the dynamical evolution of Asteroids (164207) 2004 GU₉ and 2006 FV₃₅, which are currently Earth quasi-satellites (QS). Our analysis is based on numerical computation of their orbits, and we also applied the theory of co-orbital motion developed in Wajer (Wajer, P. [2009]. Icarus 200, 147-153) to describe and analyze the objects’ dynamics. 2004 GU₉ stays as an Earth QS for about a 1000 years. In the present epoch it is in the middle of its stay in this regime. After leaving the QS orbit near 2600 this asteroid will move inside the Earth’s co-orbital region on a regular horseshoe (HS) orbit for a few 1000 years. Later, either HS-QS or HS-P transitions are possible, where P means “passing”. Although 2004 GU₉ moves primarily under the influence of the Sun and Earth, Venus plays a significant role in destabilizing the object’s orbit. Our analysis showed that the guiding center of 2006 FV₃₅ moves deep inside the averaged potential well, and since the asteroid’s argument of perihelion precesses at a rate of approximately , it prevents the QS state begin left for a long period of time; consequently the asteroid has occupied this state for about 10⁴ years and will stay in this orbit for about 800 more years. Near 2800 the asteroid’s close approach with Venus will cause it to exit the QS state, but probably it will still be moving inside the Earth’s co-orbital region and will experience transitions between HS, TP (tadpole) and P types of motion.     

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.     

Physical modeling of near-Earth Asteroid (29075) 1950 DA

Near-Earth Asteroid (29075) 1950 DA may closely encounter Earth in 2880. The probability of Earth impact may be as high as 1/300, but the outcome of the encounter depends critically on the physical properties of the asteroid [Giorgini et al., 2002. Science 196, 132-136]. We have used Arecibo and Goldstone radar data and optical lightcurves to estimate the shape, spin state, and surface structure of 1950 DA. The data allow two distinct models. One rotates prograde and is roughly spheroidal with mean diameter 1.16±0.12 km. The other rotates retrograde and is oblate and about 30% larger. Both models suggest a nickel-iron or enstatite chondritic composition. Ground-based observations should be able to determine which model is correct within the next several decades.     

Rotation and cometary activity of KBO (29981) 1999 TD10

We describe observations of the scattered Kuiper Belt object (29981) 1999 TD₁₀ performed during five observing runs at two observatories, over 370 days from 2000 September to 2001 September. They show significant brightness variations that fit a double-peaked lightcurve with period 15.448±0.012 h in V and R bands. The phase effect in V band, 0.09±0.01 mag deg⁻¹, is smaller than that of Pluto but larger than that of several KBOs, while in R band it is 0.030±0.005 mag deg⁻¹. We find color variation between the two bands, which implies a non-homogeneous albedo distribution on the surface. Evidence of surface activity near perihelion in the form of a coma/tail is presented using radial image profiles and a 2D contour map.     

Albedo and size determination of potentially hazardous asteroids: (99942) Apophis

The near-Earth object (99942) Apophis will make an extremely close approach to the Earth in 2029, and currently has approximately a one-in-45,000 chance of impacting our planet in 2036 (JPL Sentry, November 2006). Computation of the orbital evolution of this object is limited by insufficient knowledge of physical properties required to determine the role played by non-gravitational effects. Using polarimetric observations, we have obtained the first reliable determination of the albedo of Apophis, obtaining 0.33±0.08. We also derive an updated estimate of the asteroid's absolute magnitude: H=19.7±0.4. Using this albedo and H, we find that Apophis has a diameter of 270±60 m, slightly smaller than preliminary estimates based upon an assumed albedo. Our observations demonstrate the feasibility of polarimetric observations aimed at obtaining albedos and sizes of small, potentially hazardous asteroids.     

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.     

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.     

The topographic limits of gravitationally bound, rotating sand piles

Rubble pile asteroids can attain shapes that are dramatically different from those of rotating, self-gravitating equilibrium fluids. A new numerical technique, called “seed growth,” is demonstrated for calculating three-dimensional bodies that are self-gravitating and rotating, and whose every surface is approximately at a constant angle, ?, with respect to the local horizontal. By altering the configuration of cusps, which are points along a constant longitude path where the surface angle changes sign but not magnitude, multiple solution shapes that satisfy the condition that all surface slopes are at a constant angle are possible. Five different cusp configurations are explored here, three of which yield solutions for 20°???30°. Rotational effects are explored, and it is found that for some solution shapes, the ratios of their shortest to longest dimensions, c/a, can fall outside the limits published in the literature for rotating, cohesionless, spheroidal bodies. Solution shapes show some similarities to observed small bodies, such as the saturnian satellite Atlas, the near-Earth Asteroid 1999 KW4, and some contact binary asteroids.