Over one quarter of white dwarfs contain observable metallic debris from the breakup of exo-asteroids. Understanding the physical and orbital history of this debris would enable us to self-consistently link planetary system formation and fate. One major debris reservoir is generated by YORP-induced rotational fission during the giant branch phases of stellar evolution, where the stellar luminosity can exceed the Sun's by four orders of magnitude. Here, we determine the efficacy of the giant branch YORP effect for asteroids with nonzero internal strength, and model post-fission evolution by imposing simple analytic fragmentation prescriptions. We find that even the highest realistic internal strengths cannot prevent the widespread fragmentation of asteroids and the production of a debris field over 100 au in size. We compute the number of successive fission events as they occur in progressively smaller time intervals as the star ascends the giant branches, providing a way to generate size distributions of asteroid fragments. The results are highly insensitive to progenitor stellar mass. We also conclude that the ease with which giant branch YORP breakup can generate binary asteroid subsystems is strongly dependent on internal strength. Formed binary subsystems in turn could be short-lived due to the resulting luminosity-enhanced BYORP effect.

This paper considers the energy required for collections of finite density bodies to undergo escape under internal gravitational interactions alone. As the level of the system energy is increased there are different combinations of components that can escape, until the total energy becomes positive, when the entire system can undergo mutual disruption. The results are also defined for bodies modeled as a continuum. These results provide rigorous constraints for the disruption of rubble pile asteroids when only considering gravitational interaction effects, with the energy provided by rotation of an initial collection of the system. These issues are considered for discrete particles in the N body problem and for size distributions of discrete particles in the continuum limit.

We obtained thorough photometric observations of two binary near-Earth asteroids (66391) Moshup = 1999 KW4 and (88710) 2001 SL9 from June 2000 to June 2019 and derived physical and dynamical properties of the binary systems. We found that the data for 1999 KW4 are inconsistent with a constant orbital period and we obtained unique solution with a quadratic drift of the mean anomaly of the satellite of -0.65 +/- 0.16 deg/yr2 (all quoted uncertainties are 3sigma). This means that the semimajor axis of the mutual orbit of the components of this binary system increases in time with a mean rate of 1.2 +/- 0.3 cm/yr. The data for 2001 SL9 are also inconsistent with a constant orbital period and we obtained two solutions for the quadratic drift of the mean anomaly: 2.8 +/- 0.2 and 5.2 +/- 0.2 deg/yr2, implying that the semimajor axis of the mutual orbit of the components decreases in time with a mean rate of -2.8 +/- 0.2 or -5.1 +/- 0.2 cm/yr for the two solutions, respectively. The expanding orbit of 1999 KW4 may be explained by mutual tides interplaying with binary YORP (BYORP) effect (McMahon and Scheeres, 2010). However, a modeling of the BYORP drift using radar-derived shapes of the binary components predicted a much higher value of the orbital drift than the observed one. It suggests that either the radar-derived shape model of the secondary is inadequate for computing the BYORP effect, or the present theory of BYORP overestimates it. It is possible that the BYORP coefficient has instead an opposite sign than predicted; in that case, the system may be moving into an equilibrium between the BYORP and the tides. In the case of 2001 SL9, the BYORP effect is the only known physical mechanism that can cause the inward drift of its mutual orbit.

On 15 February 2013, the asteroid 367943 Duende (2012 DA14) experienced a near-Earth encounter at an altitude of 27,700 km or 4.2 Earth radii. We present here the results of an extensive, multi-observatory campaign designed to probe for spectral and/or rotational changes to Duende due to gravitational interactions with the Earth during the flyby. Our spectral data reveal no changes within the systematic uncertainties of the data. Post-flyby lightcurve photometry places strong constraints on the rotation state of Duende, showing that it is in non-principal axis rotation with fundamental periods of P_1 = 8.71 +/- 0.03 and P_2 = 23.7 +/- 0.2 hours. Multiple lightcurve analysis techniques, coupled with theoretical considerations and delay-doppler radar imaging, allows us to assign these periods to specific rotational axes of the body. In particular we suggest that Duende is now in a non-principal, short axis mode rotation state with a precessional period equal to P_1 and oscillation about the symmetry axis at a rate equal to P_2. Temporal and signal-to-noise limitations inherent to the pre-flyby photometric dataset make it difficult to definitively diagnose whether these periods represent a change imparted due to gravitational torques during the flyby. However, based on multiple analysis techniques and a number of plausibility arguments, we suggest that Duende experienced a rotational change during the planetary encounter with an increase in its precessional rotation period. Our preferred interpretation of the available data is that the precession rate increased from 8.4 hours prior to the flyby to 8.7 hours afterwards. A companion paper by Benson et al. (2019) provides a more detailed dynamical analysis of this event and compares the data to synthetic lightcurves computed from a simple shape model of Duende. (abbreviated abstract)

The full two-body problem (F2BP) is often used to model binary asteroid systems, representing the bodies as two finite mass distributions whose dynamics are influenced by their mutual gravity potential. The emergent behavior of the F2BP is highly coupled translational and rotational mutual motion of the mass distributions. A large fraction of characterized binary asteroids appear to be at, or near, the doubly synchronous equilibrium, which occurs when both bodies are tidally-locked and in a circular co-orbit. Stable oscillations about this equilibrium can be shown, for the nonplanar system, to be combinations of seven fundamental frequencies of the system and the mutual orbit rate. The fundamental frequencies arise as the linear periods of center manifolds identified about the equilibrium which are heavily influenced by each body's mass parameters. We leverage these eight dynamical constraints to investigate the observability of binary asteroid mass parameters via dynamical observations. This is accomplished by deriving a relationship between the fundamental frequencies and mass parameters for doubly synchronous systems. This relationship allows us to show the sensitivity of the dynamics to changes in the mass parameters, first for the planar dynamics, and then for the nonplanar dynamics. In so doing we are able to predict the idealized estimation covariance of the mass parameters based on observation quality and define idealized observation accuracies for desired mass parameter certainties. We apply these tools to 617 Patroclus, a doubly synchronous Trojan binary and flyby target of the LUCY mission, as well as the Pluto and Charon system in order to predict mutual behavior

The migration of cohesive regolith on the surface of an otherwise monolithic or strong asteroid is studied using theoretical and simulation models. The theory and simulations show that under an increasing spin rate (such as due to the YORP effect), the regolith covering is preferentially lost across certain regions of the body. For regolith with little or no cohesive strength, failure occurs by landsliding from the mid latitudes of the body at high enough spin rates. As the cohesive strength of the regolith increases, failure occurs by fission of grains (or coherent chunks of grains) across a greater extent of latitudes and eventually will first occur at the equator. As the spin rate is further increased, failure regions migrate from the first failure point to higher and lower latitudes. Eventually failure will encompass the equatorial region, however there always remains a region of high latitudes (around the poles) that will not undergo failure for arbitrarily high spin rates (unless disturbed by some other phenomenon). With these results a scaling law is derived that can be used to determine whether observed asteroids could retain surface regolith grains of a given size. The implications of this for the interpretation of spectral observations of small asteroids and boulder migration on large asteroids are discussed.

Asteroids and other Small Solar System Bodies (SSSBs) are of high general and scientific interest in many aspects. The origin, formation, and evolution of our Solar System (and other planetary systems) can be better understood by analysing the constitution and physical properties of small bodies in the Solar System. Currently, two space missions (Hayabusa2, OSIRIS-REx) have recently arrived at their respective targets and will bring a sample of the asteroids back to Earth. Other small body missions have also been selected by, or proposed to, space agencies. The threat posed to our planet by near-Earth objects (NEOs) is also considered at the international level, and this has prompted dedicated research on possible mitigation techniques. The DART mission, for example, will test the kinetic impact technique. Even ideas for industrial exploitation have risen during the last years. Lastly, the origin of water and life on Earth appears to be connected to asteroids. Hence, future space mission projects will undoubtedly target some asteroids or other SSSBs. In all these cases and research topics, specific knowledge of the structure and mechanical behaviour of the surface as well as the bulk of those celestial bodies is crucial. In contrast to large telluric planets and dwarf planets, a large proportion of such small bodies is believed to consist of gravitational aggregates ('rubble piles') with no -- or low -- internal cohesion, with varying macro-porosity and surface properties (from smooth regolith covered terrain, to very rough collection of boulders), and varying topography (craters, depressions, ridges) [...].

162173 Ryugu, the target of Hayabusa2, has a round shape with an equatorial ridge, which is known as a spinning top-shape. A strong centrifugal force is a likely contributor to Ryugu's top-shaped features. Observations by Optical Navigation Camera onboard Hayabusa2 show a unique longitudinal variation in geomorphology; the western side of this asteroid, later called the western bulge, has a smooth surface and a sharp equatorial ridge, compared to the other side. Here, we propose a structural deformation process that generated the western bulge. Applying the mission-derived shape model, we employ a finite element model technique to analyze the locations that experience structural failure within the present shape. Assuming that materials are uniformly distributed, our model shows the longitudinal variation in structurally failed regions when the spin period is shorter than ~3.75 h. Ryugu is structurally intact in the subsurface region of the western bulge while other regions are sensitive to structural failure. We infer that this variation is indicative of the deformation process that occurred in the past, and the western bulge is more relaxed structurally than the other region. Our analysis also shows that this deformation process might occur at a spin period between ~3.5 h and ~3.0 h, providing the cohesive strength ranging between ~4 Pa and ~10 Pa.

We observed the near-Earth asteroid (101955) Bennu from the ground in 1999 and 2005, and with the Hubble Space Telescope in 2012, to constrain its rotation rate. The data reveal an acceleration of $2.64 \pm 1.05 \times 10^{-6} \mathrm{deg\ day}^{-2}$, which could be due to a change in the moment of inertia of Bennu or to spin up from the YORP effect or other source of angular momentum. The best solution is within 1 sigma of the period determined by Nolan et al. (2013). The OSIRIS-REx mission will determine the rotation state independently in 2019. Those measurements should show whether the change in rotation rate is a steady increase (due, for example, to the YORP effect) or some other phenomenon. The precise shape and surface properties measured by the OSIRIS-REx science team will allow for a better understanding of variations in rotation rate of small asteroids.

We studied 93 asteroid pairs. We estimated times elapsed since separation of pair members that are between 7*10^3 and a few 10^6 yr. We derived the rotation periods for all the primaries and a sample of secondaries. We derived the absolute magnitude differences of the asteroid pairs that provide their mass ratios. We refined their WISE geometric albedos and estimated their taxonomic classifications. For 17 pairs, we determined their pole positions. In 2 pairs where we obtained the spin poles for both components, we saw the same sense of rotation for both components and constrained the angles between their original spin vectors at the time of their separation. We found that the primaries of 13 pairs are actually binary or triple systems, i.e., they have one or two bound secondaries (satellites). As by-product, we found 3 new young asteroid clusters (each of them consisting of three known asteroids on highly similar orbits). We compared the obtained asteroid pair data with theoretical predictions and discussed their implications. We found that 86 of the 93 studied pairs follow the trend of primary rotation period vs mass ratio that was found by Pravec et al. (2010). Of the 7 outliers, 3 appear insignificant (may be due to our uncertain or incomplete knowledge), but 4 are high mass ratio pairs that were unpredicted by the theory of asteroid pair formation by rotational fission. We discuss a (remotely) possible way that they could be created by rotational fission of flattened parent bodies followed by re-shaping of the formed components. The 13 pairs with binary primaries are particularly interesting systems that place important constraints on formation and evolution of asteroid pairs. We present two hypotheses for their formation: The pairs having both bound and unbound secondaries could be `failed asteroid clusters', or they could be formed by a cascade primary spin fission process.

There is a correlation between the components of the YORP effect of most asteroids, which drives the obliquity and spin rate of the affected bodies in a consistent pattern. This allows for a clear and unambiguous picture for how the spin rates and poles of asteroids affected by YORP will evolve and simplifies the overall picture for how populations will migrate on average. The YORP effect can also lead to a previously unexplored equilibrium state for affected bodies. This equilibrium state is a function of the usual "normal YORP" effect (which arises due to the global shape asymmetry of the asteroid) and the "tangential YORP" (which arises due the transport of thermal energy through rocky surface features). Estimates from current shape models show that 10-20\% of asteroids have the proper condition to be captured in this equilibrium state, indicating that the occurrence of this state may be significant. The existence of this attractor for the asteroid population means that objects affected by YORP may leave their usual YORP cycles and maintain a constant spin rate over long time periods -- this has significant implications for our interpretation of asteroid spin rate evolution and related theories for their physical evolution.

Many asteroids are rubble piles with irregular shapes. While the irregular shapes of large asteroids may be attributed to collisional events, those of small asteroids may result from not only impact events but also rotationally induced failure, a long-term consequence of small torques caused by, for example, solar radiation pressure. A better understanding of shape deformation induced by such small torques will allow us to give constraints on the evolution process of an asteroid and its structure. However, no quantitative study has been reported to provide the relationship between an asteroid's shape and its failure mode due to its fast rotation. Here, we use a finite element model (FEM) technique to analyze the failure modes and conditions of 24 asteroids with diameters less than 30 - 40 km, which were observed at high resolution by ground radar or asteroid exploration missions. Assuming that the material distribution is uniform, we investigate how these asteroids fail structurally at different spin rates. Our FEM simulations describe the detailed deformation mode of each irregularly shaped asteroid at fast spin. The failed regions depend on the original shape. Spheroidal objects structurally fail from the interior, while elongated objects experience structural failure on planes perpendicular to the minimum moment of inertia axes in the middle of their structure. Contact binary objects have structural failure across their most sensitive cross sections. We further investigate if our FEM analysis is consistent with earlier works that theoretically explored a uniformly rotating triaxial ellipsoid. The results show that global shape variations may significantly change the failure condition of an asteroid. Our work suggests that it is critical to take into account the actual shapes of asteroids to explore their failure modes in detail.

The evolution of rotation states of small asteroids is governed by the YORP effect, nonetheless some asteroids can stop their YORP evolution by attaining a stable equilibrium. The same is true for binary asteroids subjected to the BYORP effect. Here we discuss a new type of equilibrium that combines these two, which is possible in a singly synchronous binary system. This equilibrium occurs when the normal YORP, the tangential YORP and the binary YORP compensate each other, and tidal torques distribute the angular momentum between the components of the system and dissipate energy. Such a system if unperturbed would remain singly synchronous in perpetuity with constant spin and orbit rates, as the tidal torques dissipate the incoming energy from impinging sunlight at the same rate. The probability of the existence of this kind of equilibrium in a binary system is found to be on the order of a few percent.

The stability of a system of $N$ equal sized mutually gravitating spheres resting on each other in a straight line and rotating in inertial space is considered. This is a generalization of the "Euler Resting" configurations previously analyzed in the finite density 3 and 4 body problems. Specific questions for the general case are how rapidly the system must spin for the configuration to stabilize, how rapidly it can spin before the components separate from each other, and how these results change as a function of $N$. This paper shows that the Euler Resting configuration can only be stable for up to 5 bodies, and that for 6 or more bodies the configuration can never be stable. This places an ideal limit of 5:1 on the aspect ratio of a rubble pile body's shape.

We propose an analytic model for the normal YORP and diurnal Yarkovsky effects experienced by a convex asteroid. Both the YORP torque and the Yarkovsky force are expressed as integrals of a universal function over the surface of an asteroid. Although in general this function can only be calculated numerically from the solution of the heat conductivity equation, approximate solutions can be obtained in quadratures for important limiting cases. We consider 3 such simplified models: Rubincam's approximation (zero heat conductivity), low thermal inertia limit (including the next order correction and thus valid for small heat conductivity), and high thermal inertia limit (valid for large heat conductivity). All three simplified models are compared with the exact solution.

A doubly synchronous binary asteroid simultaneously experiences YORP and BYORP, the former being independent of the radius of the orbit, while the latter linearly dependent on the radius. In many systems YORP and BYORP can compensate each other at some radius, causing an equilibrium state, to which other orbits converge in the course of their evolution. We derive mathematical formalism of such equilibria in terms of dimensionless YORP and BYORP coefficients. We compute YORP and BYORP coefficients for a set of photometric and radar shape models, and find equilibria to be relatively common.

Binary near-Earth asteroid (65803) Didymos is the target of the proposed NASA Double Asteroid Redirection Test (DART), part of the Asteroid Impact & Deflection Assessment (AIDA) mission concept. In this mission, the DART spacecraft is planned to impact the secondary body of Didymos, perturbing mutual dynamics of the system. The primary body is currently rotating at a spin period close to the spin barrier of asteroids, and materials ejected from the secondary due to the DART impact are likely to reach the primary. These conditions may cause the primary to reshape, due to landslides, or internal deformation, changing the permanent gravity field. Here, we propose that if shape deformation of the primary occurs, the mutual orbit of the system would be perturbed due to a change in the gravity field. We use a numerical simulation technique based on the full two-body problem to investigate the shape effect on the mutual dynamics in Didymos after the DART impact. The results show that under constant volume, shape deformation induces strong perturbation in the mutual motion. We find that the deformation process always causes the orbital period of the system to become shorter. If surface layers with a thickness greater than ~0.4 m on the poles of the primary move down to the equatorial region due to the DART impact, a change in the orbital period of the system and in the spin period of the primary will be detected by ground-based measurement.

The implication of small amounts of cohesion within relatively small rubble pile asteroids is investigated with regard to their evolution under the persistent presence of the YORP effect. We find that below a characteristic size, which is a function of cohesive strength, density and other properties, rubble pile asteroids can enter a "disaggregation phase" in which they are subject to repeated fissions after which the formation of a stabilizing binary system is not possible. Once this threshold is passed rubble pile asteroids may be disaggregated into their constituent components within a finite time span. These constituent components will have their own spin limits -- albeit potentially at a much higher spin rate due to the greater strength of a monolithic body. The implications of this prediction are discussed and include modification of size distributions, prevalence of monolithic bodies among meteoroids and the lifetime of small rubble pile bodies in the solar system. The theory is then used to place constraints on the strength of binary asteroids characterized as a function of their type.

In May of 2011, NASA selected the Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) asteroid sample return mission as the third mission in the New Frontiers program. The other two New Frontiers missions are New Horizons, which explored Pluto during a flyby in July 2015 and is on its way for a flyby of Kuiper Belt object 2014 MU69 on Jan. 1, 2019, and Juno, an orbiting mission that is studying the origin, evolution, and internal structure of Jupiter. The spacecraft departed for near-Earth asteroid (101955) Bennu aboard an United Launch Alliance Atlas V 411 evolved expendable launch vehicle at 7:05 p.m. EDT on September 8, 2016, on a seven-year journey to return samples from Bennu. The spacecraft is on an outbound-cruise trajectory that will result in a rendezvous with Bennu in August 2018. The science instruments on the spacecraft will survey Bennu to measure its physical, geological, and chemical properties, and the team will use these data to select a site on the surface to collect at least 60 g of asteroid regolith. The team will also analyze the remote-sensing data to perform a detailed study of the sample site for context, assess Bennus resource potential, refine estimates of its impact probability with Earth, and provide ground-truth data for the extensive astronomical data set collected on this asteroid. The spacecraft will leave Bennu in 2021 and return the sample to the Utah Test and Training Range (UTTR) on September 24, 2023.

When gravitational aggregates are spun to fission they can undergo complex dynamical evolution, including escape and reconfiguration. Previous work has shown that a simple analysis of the full 2-body problem provides physically relevant insights for whether a fissioned system can lead to escape of the components and the creation of asteroid pairs. In this paper we extend the analysis to the full 3-body problem, utilizing recent advances in the understanding of fission mechanics of these systems. Specifically, we find that the full 3-body problem can eject a body with as much as 0.31 of the total system mass, significantly larger than the 0.17 mass limit previously calculated for the full 2-body problem. This paper derives rigorous limits on a fissioned 3-body system with regards to whether fissioned system components can physically escape from each other and what other stable relative equilibria they could settle in. We explore this question with a narrow focus on the Spherical Full Three Body Problem studied in detail earlier.

Context. Tangential YORP is a thermophysical effect that can alter the rotation rate of asteroids and is distinct from the "normal" YORP effect, but to date has only been studied for asteroids with zero obliquity. Aims. The tangential YORP force produced by spherical boulders on the surface of an asteroid with an arbitrary obliquity is studied. Methods. A finite element method is used to simulate heat conductivity inside a boulder, to find the recoil force experienced by it. Then an ellipsoidal asteroid uniformly covered by such boulders is considered and the torque is numerically integrated over its surface. Results. Tangential YORP is found to operate on non-zero obliquities and decreases by a factor of 2 for increasing obliquity.

From the results of a comprehensive asteroid population evolution model, we conclude that the YORP-induced rotational fission hypothesis can be consistent with the observed population statistics of small asteroids in the main belt including binaries and contact binaries. The foundation of this model is the asteroid rotation model of Marzari et al. (2011), which incorporates both the YORP effect and collisional evolution. This work adds to that model the rotational fission hypothesis and the binary evolution model of Jacobson & Scheeres (2011). The asteroid population evolution model is highly constrained by these and other previous works, and therefore it has only two significant free parameters: the ratio of low to high mass ratio binaries formed after rotational fission events and the mean strength of the binary YORP (BYORP) effect. We successfully reproduce characteristic statistics of the small asteroid population: the binary fraction, the fast binary fraction, steady-state mass ratio fraction and the contact binary fraction. We find that in order for the model to best match observations, rotational fission produces high mass ratio ($> 0.2$) binary components with four to eight times the frequency as low mass ratio ($< 0.2$) components, where the mass ratio is the mass of the secondary component divided by the mass of the primary component. This is consistent with post-rotational fission binary system mass ratio being drawn from either a flat or a positive and shallow distribution, since the high mass ratio bin is four times the size of the low mass ratio bin; this is in contrast to the observed steady-state binary mass ratio, which has a negative and steep distribution. This can be understood in the context of the BYORP-tidal equilibrium hypothesis, which predicts that low mass ratio binaries survive for a significantly longer period of time than high mass ratio systems.

The relative equilibria for the spherical, finite density 3 body problem are identified. Specifically, there are 28 distinct relative equilibria in this problem which include the classical 5 relative equilibria for the point-mass 3-body problem. None of the identified relative equilibria exist or are stable over all values of angular momentum. The stability and bifurcation pathways of these relative equilibria are mapped out as the angular momentum of the system is increased. This is done under the assumption that they have equal and constant densities and that the entire system rotates about its maximum moment of inertia. The transition to finite density greatly increases the number of relative equilibria in the 3-body problem and ensures that minimum energy configurations exist for all values of angular momentum.

Surface shedding of an asteroid is a failure mode where surface materials fly off due to strong centrifugal forces beyond the critical spin period, while the internal structure does not deform significantly. This paper proposes a possible structure of an asteroid interior that leads to such surface shedding due to rapid rotation rates. A rubble pile asteroid is modeled as a spheroid composed of a surface shell and a concentric internal core, the entire assembly called the test body. The test body is assumed to be uniformly rotating around a constant rotation axis. We also assume that while the bulk density and the friction angle are constant, the cohesion of the surface shell is different from that of the internal core. First, developing an analytical model based on limit analysis, we provide the upper and lower bounds for the actual surface shedding condition. Second, we use a Soft-Sphere Discrete Element Method (SSDEM) to study dynamical deformation of the test body due to a quasi-static spin-up. In this paper we show the consistency of both approaches. Additionally, the SSDEM simulations show that the initial failure always occurs locally and not globally. In addition, as the core becomes larger, the size of lofted components becomes smaller. These results imply that if there is a strong enough core in a progenitor body, surface shedding is the most likely failure mode.

The Yarkovsky effect describes a small but significant force that affects the orbital motion of meteoroids and asteroids smaller than $30-40$ kilometers in diameter. It is caused by sunlight; when these bodies heat up in the Sun, they eventually re-radiate the energy away in the thermal waveband, which in turn creates a tiny thrust. This recoil acceleration is much weaker than solar and planetary gravitational forces, but it can produce measurable orbital changes over decades and substantial orbital effects over millions to billions of years. The same physical phenomenon also creates a thermal torque that, complemented by a torque produced by scattered sunlight, can modify the rotation rates and obliquities of small bodies as well. This rotational variant has been coined the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. During the past decade or so, the Yarkovsky and YORP effects have been used to explore and potentially resolve a number of unsolved mysteries in planetary science dealing with small bodies. Here we review the main results to date, and preview the goals for future work.

Rozitis et al. recently reported that near-Earth asteroid (29075) 1950 DA, whose bulk density ranges from 1.0 g/cm3 to 2.4 g/cm3, is a rubble pile and requires a cohesive strength of at least 44 Pa to 74 Pa to keep from failing due to its fast spin period. Since their technique for giving failure conditions required the averaged stress over the whole volume, it discarded information about the asteroid's failure mode and internal stress condition. This paper develops a finite element model and revisits the stress and failure analysis of 1950 DA. For the modeling, we do not consider material-hardening and softening. Under the assumption of an associated flow rule and uniform material distribution, we identify the deformation process of 1950 DA when its constant cohesion reaches the lowest value that keeps its current shape. The results show that to avoid structural failure the internal core requires a cohesive strength of at least 75 Pa - 85 Pa. It suggests that for the failure mode of this body, the internal core first fails structurally, followed by the surface region. This implies that if cohesion is constant over the whole volume, the equatorial ridge of 1950 DA results from a material flow going outward along the equatorial plane in the internal core, but not from a landslide as has been hypothesized. This has additional implications for the likely density of the interior of the body.

The Design Reference Asteroid (DRA) is a compilation of all that is known about the OSIRIS-REx mission target, asteroid (101955) Bennu. It contains our best knowledge of the properties of Bennu based on an extensive observational campaign that began shortly after its discovery, and has been used to inform mission plan development and flight system design. The DRA will also be compared with post-encounter science results to determine the accuracy of our Earth-based characterization efforts. The extensive observations of Bennu in 1999 has made it one of the best-characterized near-Earth asteroids. Many physical parameters are well determined, and span a number of categories: Orbital, Bulk, Rotational, Radar, Photometric, Spectroscopic, Thermal, Surface Analog, and Environment Properties. Some results described in the DRA have been published in peer-reviewed journals while others have been reviewed by OSIRIS-REx Science Team members and/or external reviewers. Some data, such as Surface Analog Properties, are based on our best knowledge of asteroid surfaces, in particular those of asteroids Eros and Itokawa. This public release of the OSIRIS-REx Design Reference Asteroid is a annotated version of the internal OSIRIS-REx document OREX-DOCS-04.00-00002, Rev 9 (accepted by the OSIRIS-REx project on 2014-April-14). The supplemental data products that accompany the official OSIRIS-REx version of the DRA are not included in this release. We are making this document available as a service to future mission planners in the hope that it will inform their efforts.

Conditions for regolith landslides to occur on spinning, gravitating spheroidal asteroids and their aftermath are studied. These conditions are developed by application of classical granular mechanics stability analysis to the asteroid environment. As part of our study we determine how slopes evolve across the surface of these bodies as a function of spin rate, the dynamical fate of material that exceeds the angle of repose, and an analysis of how the shape of the body may be modified based on these results. We find specific characteristics for body surfaces and shapes when spun near the surface disruption limit and develop what their observable implications are. The small, oblate and rapidly spinning asteroids such as 1999 KW4 Alpha and 2008 EV5 exhibit some of these observable traits. The detailed mechanisms outlined here can also provide insight and constraints on the recently observed active asteroids such as P/2013 P5, and the creation of asteroidal meteor streams.

Tangential YORP, or TYORP, has recently been demonstrated to be an important factor in the evolution of an asteroid's rotation state. It is complementary to normal YORP, or NYORP, which used to be considered previously. While NYORP is produced by non-symmetry in the large-scale geometry of an asteroid, TYORP is due to heat conductivity in stones on the surface of the asteroid. Yet to date TYORP has been studied only in a simplified 1-dimensional model, substituting stones by high long walls. This article for the first time considers TYORP in a realistic 3-dimensional model, also including shadowing and self-illumination effects via ray tracing. TYORP is simulated for spherical stones lying on regolith. The model includes only 5 free parameters, and the dependence of the TYORP on each of them is studied. The TYORP torque appears to be smaller than previous estimates from 1-dimensional model, but still comparable to the NYORP torques. These results can be used to estimate TYORP of different asteroids, and also as a basis for more sophisticated models of TYORP.

Jewitt et al. (2014) recently reported that main belt comet P/2013 R3 experienced a breakup, probably due to rotational disruption, with its components separating on mutually hyperbolic orbits. We propose a technique for constraining physical properties of the proto-body, especially the initial spin period and cohesive strength, as a function of the body's estimated size and density. The breakup conditions are developed by combining mutual orbit dynamics of the smaller components and the failure condition of the proto-body. Given a proto-body with a bulk density ranging from 1000 kg/m3 to 1500 kg/m3 (a typical range of the bulk density of C-type asteroids), we obtain possible values of the cohesive strength (40 - 210 Pa) and the initial spin state (0.48 - 1.9 hr). From this result, we conclude that although the proto-body could have been a rubble pile, it was likely spinning beyond its gravitational binding limit and would have needed cohesive strength to hold itself together. Additional observations of P/2013 R3 will enable stronger constraints on this event, and the present technique will be able to give more precise estimates of its internal structure.

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.

We explore the hypothesis that, due to small van der Waals forces between constituent grains, small rubble pile asteroids have a small but non-zero cohesive strength. The nature of this model predicts that the cohesive strength should be constant independent of asteroid size, which creates a scale dependence with relative strength increasing as size decreases. This model counters classical theory that rubble pile asteroids should behave as scale-independent cohesionless collections of rocks. We explore a simple model for asteroid strength that is based on these weak forces, validate it through granular mechanics simulations and comparisons with properties of lunar regolith, and then explore its implications and ability to explain and predict observed properties of small asteroids in the NEA and Main Belt populations, and in particular of asteroid 2008 TC3. One conclusion is that the population of rapidly rotating asteroids could consist of both distributions of smaller grains (i.e., rubble piles) and of monolithic boulders.

The size distribution of small asteroids in the Main Belt is assumed to be determined by an equilibrium between the creation of new bodies out of the impact debris of larger asteroids and the destruction of small asteroids by collisions with smaller projectiles. However, for a diameter less than 6 km we find that YORP-induced rotational disruption significantly contributes to the erosion even exceeding the effects of collisional fragmentation. Including this additional grinding mechanism in a collision evolution model for the asteroid belt, we generate size-frequency distributions from either an accretional (Weidenschilling, 2011) or an "Asteroids were born big" (Morbidelli, 2009) initial size-frequency distribution that are consistent with observations reported in Gladman et al. (2009). Rotational disruption is a new mechanism that must be included in all future collisional evolution models of asteroids.

Given the spin state by Magnusson (1990), the shape model by Ostro et al. (2000), and the mass by Descamps et al. (2011), this paper evaluates a dynamically and structurally stable size of Asteroid (216) Kleopatra. In particular, we investigate two different failure modes: material shedding from the surface and structural failure of the internal body. We construct zero-velocity curves in the vicinity of this asteroid to determine surface shedding, while we utilize a limit analysis to calculate the lower and upper bounds of structural failure under the zero-cohesion assumption. Surface shedding does not occur at the current spin period (5.385 hr) and cannot directly initiate the formation of the satellites. On the other hand, this body may be close to structural failure; in particular, the neck may be situated near a plastic state. In addition, the neck's sensitivity to structural failure changes as the body size varies. We conclude that plastic deformation has probably occurred around the neck part in the past. If the true size of this body is established through additional measurements, this method will provide strong constraints on the current friction angle for the body.

We propose and analyze a new mechanism for the formation of the wide asynchronous binary population. These binary asteroids have wide semi-major axes relative to most near-Earth and Main Belt asteroid systems. Confirmed members have rapidly rotating primaries and satellites that are not tidally locked. Previously suggested formation mechanisms from impact ejecta, planetary flybys and directly from rotational fission events cannot satisfy all of the observations. The newly hypothesized mechanism works as follows: (i) these systems are formed from rotational fission, (ii) their satellites are tidally locked, (iii) their orbits are expanded by the BYORP effect, (iv) their satellites de-synchronize due to the adiabatic invariance between the libration of the secondary and the mutual orbit, and (v) the secondary avoids resynchronization due to the the YORP effect. This seemingly complex chain of events is a natural pathway for binaries with satellites that have particular shapes, which define the BYORP effect torque that acts on the system. After detailing the theory, we analyze each of the wide asynchronous binary members and candidates to assess their most likely formation mechanism. Finally, we suggest possible future observations to check and constrain our hypothesis.

This paper explores the dynamic properties of the planar system of an ellipsoidal satellite in an equatorial orbit about an oblate primary. In particular, we investigate the conditions for which the satellite is bound in librational motion or when the satellite will circulate with respect to the primary. We find the existence of stable equilibrium points about which the satellite can librate, and explore both the linearized and non-linear dynamics around these points. Absolute bounds are placed on the phase space of the libration-orbit coupling through the use of zero-velocity curves that exist in the system. These zero-velocity curves are used to derive a sufficient condition for when the satellite's libration is bound to less than 90 degrees. When this condition is not satisfied so that circulation of the satellite is possible, the initial conditions at zero libration angle are determined which lead to circulation of the satellite. Exact analytical conditions for circulation and the maximum libration angle are derived for the case of a small satellite in orbits of any eccentricity.

Minimum energy configurations in celestial mechanics are investigated. It is shown that this is not a well defined problem for point-mass celestial mechanics but well-posed for finite density distributions. This naturally leads to a granular mechanics extension of usual celestial mechanics questions such as relative equilibria and stability. This paper specifically studies and finds all relative equilibria and minimum energy configurations for $N=1,2,3$ and develops hypotheses on the relative equilibria and minimum energy configurations for $N\gg 1$ bodies.

Synchronous binary asteroids may exist in a long-term stable equilibrium, where the opposing torques from mutual body tides and the binary YORP (BYORP) effect cancel. Interior of this equilibrium, mutual body tides are stronger than the BYORP effect and the mutual orbit semi-major axis expands to the equilibrium; outside of the equilibrium, the BYORP effect dominates the evolution and the system semi-major axis will contract to the equilibrium. If the observed population of small (0.1 - 10 km diameter) synchronous binaries are in static configurations that are no longer evolving, then this would be confirmed by a null result in the observational tests for the BYORP effect. The confirmed existence of this equilibrium combined with a shape model of the secondary of the system enables the direct study of asteroid geophysics through the tidal theory. The observed synchronous asteroid population cannot exist in this equilibrium if described by the canonical "monolithic" geophysical model. The "rubble pile" geophysical model proposed by \citet{Goldreich2009} is sufficient, however it predicts a tidal Love number directly proportional to the radius of the asteroid, while the best fit to the data predicts a tidal Love number inversely proportional to the radius. This deviation from the canonical and \citet{Goldreich2009} models motivates future study of asteroid geophysics. Ongoing BYORP detection campaigns will determine whether these systems are in an equilibrium, and future determination of secondary shapes will allow direct determination of asteroid geophysical parameters.

This document covers five broad topics from theory of small bodies: planetesimal formation, cosmochemistry, thermal evolution, collisions and dynamics. Each of these topics is described in a separate section, where we prioritize selected main issues over completeness. The text points out the principal unresolved problems in each area, and suggests ways how progress can be made. This includes support for new code development, observations and experimental work that can be used to constrain theory, new research directions, and studies of cross-over regimes where the system's behavior is determined by several competing processes. The suggested development areas are placed in the context of NASA space exploration.

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

Asteroid pairs sharing similar heliocentric orbits were found recently. Backward integrations of their orbits indicated that they separated gently with low relative velocities, but did not provide additional insight into their formation mechanism. A previously hypothesized rotational fission process4 may explain their formation - critical predictions are that the mass ratios are less than about 0.2 and, as the mass ratio approaches this upper limit, the spin period of the larger body becomes long. Here we report photometric observations of a sample of asteroid pairs revealing that primaries of pairs with mass ratios much less than 0.2 rotate rapidly, near their critical fission frequency. As the mass ratio approaches 0.2, the primary period grows long. This occurs as the total energy of the system approaches zero requiring the asteroid pair to extract an increasing fraction of energy from the primary's spin in order to escape. We do not find asteroid pairs with mass ratios larger than 0.2. Rotationally fissioned systems beyond this limit have insufficient energy to disrupt. We conclude that asteroid pairs are formed by the rotational fission of a parent asteroid into a proto-binary system which subsequently disrupts under its own internal system dynamics soon after formation.

The scaling of physical forces to the extremely low ambient gravitational acceleration regimes found on the surfaces of small asteroids is performed. Resulting from this, it is found that van der Waals cohesive forces between regolith grains on asteroid surfaces should be a dominant force and compete with particle weights and be greater, in general, than electrostatic and solar radiation pressure forces. Based on this scaling, we interpret previous experiments performed on cohesive powders in the terrestrial environment as being relevant for the understanding of processes on asteroid surfaces. The implications of these terrestrial experiments for interpreting observations of asteroid surfaces and macro-porosity are considered, and yield interpretations that differ from previously assumed processes for these environments. Based on this understanding, we propose a new model for the end state of small, rapidly rotating asteroids which allows them to be comprised of relatively fine regolith grains held together by van der Waals cohesive forces.

We have developed a theoretical description of how of an asteroid's polarization-phase curve will be affected by the removal of the dust from the surface due to a size-dependent phenomenon such as radiation pressure-driven escape of levitated particles. We test our calculations against new observations of four small (D ~ 1 km) near-Earth asteroids [(85236), (142348), (162900) and 2006 SZ_217] obtained with the Dual Beam Imaging Polarimeter on the University of Hawaii's 2.2 m telescope, as well as previous observations of (25143) Itokawa and (433) Eros. We find that the polarization of the light reflected from an asteroid is controlled by the mineralogical and chemical composition of the surface and is independent of dust particle. The relation between the slope of the polarization-phase curve beyond the inversion angle and the albedo of an asteroid is thus independent of the surface regolith size distribution and is valid for both Main Belt and Near-Earth asteroids.

Given the impressive investment by the nation in observational Astronomy and Astrophysics facilities coming on line now and in the near future, we advocate for an increased investment in applied and fundamental research on Astrophysical and Celestial Dynamics (ACD). Specifically we call for a) continued and expanded support for applied research in ACD, b) creation of support for fundamental research in ACD and its subfields, and c) the creation of a unified program to help scientists coordinate and collaborate in their research in these fields. The benefits of this proposal are threefold. First, it will enable researchers to interpret and understand the implications of newly observed phenomena that will invariably arise from new facilities and surveys. Second, research on fundamentals will foster connections between specialists, leveraging advances found in one sub-field and making them available to others. Third, a coordinated approach for applied and fundamental research in ACD will help academic institutions in the United States to produce future researchers trained and knowledgeable in essential subfields such as Mathematical Celestial Mechanics and able to continue its advancement in conjunction with the increase in observations.

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

The rosette-shaped motion of a particle in a central force field is known to be classically solvable by quadratures. We present a new approach of describing and characterizing such motion based on the eccentricity vector of the two body problem. In general, this vector is not an integral of motion. However, the orbital motion, when viewed from the nonuniformly rotating frame defined by the orientation of the eccentricity vector, can be solved analytically and will either be a closed periodic circulation or libration. The motion with respect to inertial space is then given by integrating the argument of periapsis with respect to time. Finally we will apply the decomposition to a modern central potential, the spherical Hernquist-Newton potential, which models dark matter halos of galaxies with central black holes.