Mass dispersal and angular momentum transfer during collisions between rubble-pile asteroids

We perform N-body simulations of impacts between initially non-rotating rubble-pile asteroids, and investigate mass dispersal and angular momentum transfer during such collisions. We find that the fraction of the dispersed mass (Mdisp) is approximately proportional to , where Qimp is the impact kinetic energy; the power index α is about unity when the impactor is much smaller than the target, and 0.5?α<1 for impacts with a larger impactor. Mdisp is found to be smaller for more dissipative impacts with small values of the restitution coefficient of the constituent particles. We also find that the efficiency of transfer of orbital angular momentum to the rotation of the largest remnant depends on the degree of disruption. In the case of disruptive oblique impacts where the mass of the largest remnant is about half of the target mass, most of the orbital angular momentum is carried away by the escaping fragments and the efficiency becomes very low (<0.05), while the largest remnant acquires a significant amount of spin angular momentum in moderately disruptive impacts. These results suggest that collisions likely played an important role in rotational evolution of small asteroids, in addition to the recoil force of thermal re-radiation.     

The asteroids as outcomes of catastrophic collisions

The role of catastrophic collisions in the evolution of the asteroids is discussed in detail, employing extrapolations of experimental results on the outcomrs of high-velocity impacts. We determine the range of the probable largest collision for target asteroids of different sizes during the solar system's lifetime, and we conclude that all the asteroids have undergone collisional events capable of overcoming the material's solid-state cohesion. Such events do not lead inescapably to complete disruption of the targets, because (i) for a previously unfractured target, experiments show that fragments of significant size can survive breakup, depending on the energy and geometry of the collision; (ii) self-gravitation can easily cause a reaccumulation of fragments for targets exceeding a critical size, which seems to be of the order of 100 km. In the intermediate diameter range 100?D ?300 km, where formation of gravitationally bound “rubble piles” is frequent, the transfer of angular momentum can be large enough to produce objects with triaxial equilibrium shapes (Jacobi ellipsoids) or to cause fission into binary systems. In the same size range, low-velocity escape of collisional fragments can also occur, leading to the formation of dynamical families. Asteroids smaller than ～100 km are mostly multigeneration fragments, while for $\text{D?300}\text{km}$ the collisional process produces nearly spheroidal objects covered by megaregoliths; whether their rotation is “primordial” or collisionally generated depends critically on the past flux of colliders. The complex and size-dependent phenomenology predicted by the theory compares satisfactorily with the observational evidence, as derived both by a classification of asteroids in terms of their size, spin rate, and lightcurve amplitude, and by a comparison between the rotational properties of family and nonfamily asteroids. The fundamental result of this investigation is that almost all asteroids are outcomes of catastrophic collisions, and that these events cause either complete fragmentation of the target bodies or, at least, drastic readjustments of their internal structure, shape, and spin rate.     

Lunar-forming collisions with pre-impact rotation

Prior models of lunar-forming impacts assume that both the impactor and the target protoearth were not rotating prior to the Moon-forming event. However, planet formation models suggest that such objects would have been rotating rapidly during the late stages of terrestrial accretion. In this paper I explore the effects of pre-impact rotation on impact outcomes through more than 100 hydrodynamical simulations that consider a range of impactor masses, impact angles and impact speeds. Pre-impact rotation, particularly in the target protoearth, can substantially alter collisional outcomes and leads to a more diverse set of final planet-disk systems than seen previously. However, the subset of these impacts that are also lunar-forming candidates—i.e. that produce a sufficiently massive and iron-depleted protolunar disk—have properties similar to those determined for collisions of non-rotating objects [Canup, R.M., Asphaug, E., 2001. Nature 412, 708-712; Canup, R.M., 2004a. Icarus 168, 433-456]. With or without pre-impact rotation, a lunar-forming impact requires an impact angle near 45 degrees, together with a low impact velocity that is not more than 10% larger than the Earth's escape velocity, and produces a disk containing up to about two lunar masses that is composed predominantly of material originating from the impactor. The most significant differences in the successful cases involving pre-impact spin occur for impacts into a retrograde rotating protoearth, which allow for larger impactors (containing up to 20% of Earth's mass) and provide an improved match with the current Earth-Moon system angular momentum compared to prior results. The most difficult state to reconcile with the Moon is that of a rapidly spinning, low-obliquity protoearth before the giant impact, as these cases produce disks that are not massive enough to yield the Moon.     

Collapse and explosion in a rotating star

The collapse, bounce, shock wave and expansion of the envelope of a rotating star have been analysed in the adiabatic approximation using the particle-in-cell method. The bounce takes place first in the equatorial plane and a shock wave arises there which shortly afterwards crosses the surface of the star. In the envelope, and to a less extent in the remainder of the star, there is a fast and lasting meridional motion the direction of which changes. As a consequence of the fast meridional motion in the envelope, mass and angular momentum are transported towards the axis of rotation. If the initial star rotates fast enough this will cause a secondary radial expansion in the polar region and a mass ejection. These motions reduce the strong anisotropy caused originally by the equatorial expansion. Strong whirls may arise along the axis of rotation. In the remainder of the star the meridional motion becomes supersonic. The temperature in the envelope depends to a high degree on the choice of the equation of state. Massloss is proportional to the energy initially added. The final loss of angular momentum and of energy is quite large, both losses being about 25%.     

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.     

On YORP-induced spin deformations of asteroids

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

Collapse and fragmentation of rotating magnetized clouds – II. Binary formation and fragmentation of first cores

Subsequent to Paper I, the evolution and fragmentation of a rotating magnetized cloud are studied with use of three-dimensional magnetohydrodynamic nested grid simulations. After the isothermal runaway collapse, an adiabatic gas forms a protostellar first core at the centre of the cloud. When the isothermal gas is stable for fragmentation in a contracting disc, the adiabatic core often breaks into several fragments. Conditions for fragmentation and binary formation are studied. All the cores which show fragmentation are geometrically thin, as the diameter-to-thickness ratio is larger than 3. Two patterns of fragmentation are found. (1) When a thin disc is supported by centrifugal force, the disc fragments into a ring configuration (ring fragmentation). This is realized in a rapidly rotating adiabatic core as  Ω > 0.2τ⁻¹_ff  , where Ω and  τ_ff  represent the angular rotation speed and the free-fall time of the core, respectively. (2) On the other hand, the disc is deformed to an elongated bar in the isothermal stage for a strongly magnetized or rapidly rotating cloud. The bar breaks into 2–4 fragments (bar fragmentation). Even if a disc is thin, the disc dominated by the magnetic force or thermal pressure is stable and forms a single compact body. In either ring or bar fragmentation mode, the fragments contract and a pair of outflows is ejected from the vicinities of the compact cores. The orbital angular momentum is larger than the spin angular momentum in the ring fragmentation. On the other hand, fragments often quickly merge in the bar fragmentation, since the orbital angular momentum is smaller than the spin angular momentum in this case. Comparison with observations is also shown.     

Impact evolution of asteroid shapes: 1. random mass redistribution

We explore whether the cumulative effect of small-scale meteoroid bombardment can drive asteroids into nonaxisymmetric shapes comparable to those of known objects (elongated prolate forms, twin-lobed binaries, etc). We simulate impact cratering as an excavation followed by the launch, orbit, and reimpact of ejecta. Orbits are determined by the gravity and rotation of the evolving asteroid, whose shape and spin change as cratering occurs repeatedly. For simplicity we consider an end-member evolution where impactors are all much smaller than the asteroid and where all ejecta remain bound. Given those assumptions, we find that cumulative small impacts on rotating asteroids lead to oblate shapes, irrespective of the chosen value for angle of repose or for initial angular momentum. The more rapidly a body is spinning, the more flattened the outcome, but oblateness prevails. Most actual asteroids, by contrast, appear spherical to prolate. We also evaluate the timescale for reshaping by small impacts and compare it to the timescale for catastrophic disruption. For all but the steepest size distributions of impactors, reshaping from small impacts takes more than an order of magnitude longer than catastrophic disruption. We conclude that small-scale cratering is probably not dominant in shaping asteroids, unless our assumptions are naive. We believe we have ruled out the end-member scenario; future modeling shall include angular momentum evolution from impacts, mass loss in the strength regime, and craters with diameters up to the disruption threshold. The ultimate goal is to find out how asteroids get their shapes and spins and whether tidal encounters in fact play a dominant role.     

A comparison between rubble-pile and monolithic targets in impact simulations: Application to asteroid satellites and family size distributions

Collisions are a fundamental process in the creation of asteroid families and in satellite formation. For this reason, understanding the outcome of impacts is fundamental to the accurate modeling of the formation and evolution of such systems. Smoothed-Particle Hydrodynamics/N-body codes have become the techniques of choice to study large-scale impact outcomes, including both the fragmentation of the parent body and the gravitational interactions between fragments. It is now possible to apply this technique to targets with either monolithic or rubble-pile internal structures. In this paper we apply these numerical techniques to rubble-pile targets, extending previous investigations by Durda et al. (Durda, D.D., Bottke, W.F., Enke, B.L., Merline, W.J., Asphaug, E., Richardson, D.C., Leinhardt, Z.M. [2004]. Icarus 170, 243–257; Durda, D.D., Bottke, W.F., Nesvorný, D., Enke, B.L., Merline, W.J., Asphaug, E., Richardson, D.C. [2007]. Icarus 186, 498–516). The goals are to study asteroid–satellite formation and the morphology of the size–frequency distributions (SFDs) from 175 impact simulations covering a range of collision speeds, impact angles, and impactor sizes. Our results show that low-energy impacts into rubble-pile and monolithic targets produce different features in the resulting SFDs and that these are potentially diagnostic of the initial conditions for the impact and the internal structure of the parent bodies of asteroid families. In contrast, super-catastrophic events (i.e., high-energy impacts with large specific impact energy) result in SFDs that are similar to each other. We also find that rubble-pile targets are less efficient in producing satellites than their monolithic counterparts. However, some features, such as the secondary-to-primary diameter ratio and the relative separation of components in binary systems, are similar for these two different internal structures of parent bodies.     

Secular acceleration of solar rotation from 1943 to 1986

The rotation of the surface layer of the Sun is found to have been accelerated secularly from the sunspot data of 1943 to 1986. To represent the overall state of rotation of the differentially rotating Sun, we define an indexM, by integrating the angular momentum density over the whole surface of the Sun, and call it the angular momentum layer density. The indexM increased monotonically or secularly from 1943 to 1986. This period corresponds to solar cycles 18, 19, 20, and 21. The monotonic increase ofM indicates that a net angular momentum must have steadily been coming from the layer down below the surface. The differential rotation latitudinal dependence profile did not change much from cycle 18 to cycle 20, but at cycle 21 the degree of equatorial acceleration dropped. This aspect is discussed in the context of the 55-year grand cycle. Cycle 21 is the start of grand cycle VI. The latitudinal dependence is less steep at cycle 21. The time scale of secular change of the indexM reflects the time scale of change of linkage of the surface and the deep layer in form of the angular momentum transfer, and that the time scale of the profile change of the differential rotation reflects the time scale of the angular momentum transfer within the surface layer.     

A 2D study of the effects of the size of a solid core on the equatorial flow in giant planets

Zonal winds simulated in two-dimensional computer models of turbulent convection in the equatorial plane of giant planets have greater surface amplitudes for cases with smaller solid cores, and therefore larger buoyancy driving, all other properties being equal. This differential rotation in radius is maintained by the convergence of angular momentum flux, which occurs because of the convective flow that develops due to the effects of planetary rotation and density stratification. The superposition of the convective flow and the stronger zonal flow produces wave-like, instead of cellular convection.     

Mass Loss from Rotating Hot-stars: Inhibition of Wind Compressed Disks by Nonradial Line-forces

We review the dynamics of radiatively driven mass loss from rapidly rotating hot-stars. We first summarize the angular momentum conservation process that leads to formation of a Wind Compressed Disk(WCD) when material from a rapidly rotating star is driven gradually outward in the radial direction. We next describe how stellar oblateness and asymmetries in the Sobolev line-resonance generally leads to nonradialcomponents of the driving force is a line-driven wind, including an azimuthal spin-down force acting against the sense of the wind rotation, and a latitudinal force away from the equator. We summarize results from radiation-hydrodynamical simulations showing that these nonradial forces can lead to an effective suppressionof the equatorward flow needed to form a WCD, as well as a modest (∼ 25%) spin-downof the wind rotation. Furthermore, contrary to previous expectations that the wind mass flux should be enhanced by the reduced effective gravity near the equator, we show here that gravity darkening effects can actually lead to a reducedmass loss, and thus lower density, in the wind from the equatorial region. Finally, we examine the equatorial bistability model, and show that a sufficiently strong jump in wind driving parameters can, in principle, overcome the effect of reduced radiative driving flux, thus still allowing moderate enhancements in density in an equatorial, bistability zone wind. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the angular velocity distribution in a nonhomogeneously rotating star

Some results of the calculations of the angular velocity distribution in the equatorial plane of a nonhomogeneously rotating star, in particular of the Sun, are presented. Hence the hypothesis is advanced of the turbulent viscosity as one of the main factors of the angular momentum transport from a rapidly rotating core through the radiative equilibrium region.     

Creating an isotopically similar Earth–Moon system with correct angular momentum from a giant impact

The giant impact hypothesis is the dominant theory explaining the formation of our Moon. However, the inability to produce an isotopically similar Earth–Moon system with correct angular momentum has cast a shadow on its validity. Computer-generated impacts have been successful in producing virtual systems that possess many of the observed physical properties. However, addressing the isotopic similarities between the Earth and Moon coupled with correct angular momentum has proven to be challenging. Equilibration and evection resonance have been proposed as means of reconciling the models. In the summer of 2013, the Royal Society called a meeting solely to discuss the formation of the Moon. In this meeting, evection resonance and equilibration were both questioned as viable means of removing the deficiencies from giant impact models. The main concerns were that models were multi-staged and too complex. We present here initial impact conditions that produce an isotopically similar Earth–Moon system with correct angular momentum. This is done in a single-staged simulation. The initial parameters are straightforward and the results evolve solely from the impact. This was accomplished by colliding two roughly half-Earth-sized impactors, rotating in approximately the same plane in a high-energy, off-centered impact, where both impactors spin into the collision.     

On the angular momentum in star formation

We discuss the rotation of interstellar clouds which are in a stage immediately before star formation. Cloud collisions seem to be the principal cause of the observed rotation of interstellar clouds. The rotational motion of the clouds is strongly influenced by turbulence.Theories dealing with the resolution of the angular momentum problem in star formation are classified into five major groups. We develop the old idea that the angular momentum of an interstellar cloud passes during star formation into the angular momentum of double star systems and/or circumstellar clouds.It is suggested that a rotating gas cloud contracts into a ring-like structure which fragments into self-gravitating subcondensations. By collisions and gas accretion these subcondensations accrete into binary systems surrounded by circumstellar clouds. Using some rough approximations we find analytical expressions for the semi-major axis of the binary system and for the density of the circumstellar clouds as a function of the initial density and of the initial angular velocity of an interstellar cloud. The obtained values are well within the observational limits.     

On equilibrium rotation of the central object of an accretion disc

The hydrodynamic interaction of an accretion disc with its central object is reanalysed within the framework of the slim-disc approximation. Arguments are presented against an interpretation of the total angular momentum flux as an eigenvalue of the system. A simple intuitive consideration is provided, which shows that the central object may be in a state of stationary rotation even if the disc imposes the constraint of a finite angular momentum flux into it. It is argued that equilibrium rotation is characterized by vanishing viscous torque rather than by zero total angular momentum flux. As a consequence, the central object can be in a state of stationary rotation below the break-up limit, although its angular momentum increases. Despite accretion, even for positive total angular momentum flux and subcritical rotation, central objects are spun down within a considerable range of their parameters. The results are illustrated by application to FU Orionis systems.     

The possible role of r-modes in post-glitch relaxation of the Crab pulsar

The loss of angular momentum through gravitational radiation, driven by the excitation of r-modes, is considered for neutron stars that have rotation frequencies lower than the associated critical frequency. We find that for reasonable values of the initial amplitudes of such pulsation modes of the star, being excited at the event of a glitch in a pulsar, the total post-glitch losses correspond to a negligible fraction of the initial rise of the spin frequency in the case of Vela and older pulsars. However, for the Crab pulsar the same effect would result, within a few months, in a decrease in its spin frequency by an amount larger than its glitch-induced frequency increase. This could provide an explanation for the peculiar behaviour observed in the post-glitch relaxations of the Crab pulsar.     

Cumulative damage in strength-dominated collisions of rocky asteroids: Rubble piles and brick piles

Laboratory impact experiments were performed to investigate the conditions that produce large-scale damage in rock targets. Aluminum cylinders (6.3 mm diameter) impacted basalt cylinders (69 mm diameter) at speeds ranging from 0.7 to 2.0 km/s. Diagnostics included measurements of the largest fragment mass, velocities of the largest remnant and large fragments ejected from the periphery of the target, and X-ray computed tomography imaging to inspect some of the impacted targets for internal damage. Significant damage to the target occurred when the kinetic energy per unit target mass exceeded roughly of the energy required for catastrophic shattering (where the target is reduced to one-half its original mass). Scaling laws based on a rate-dependent strength were developed that provide a basis for extrapolating the results to larger strength-dominated collisions. The threshold specific energy for widespread damage was found to scale with event size in the same manner as that for catastrophic shattering. Therefore, the factor of four difference between the two thresholds observed in the lab also applies to larger collisions. The scaling laws showed that for a sequence of collisions that are similar in that they produce the same ratio of largest fragment mass to original target mass, the fragment velocities decrease with increasing event size. As a result, rocky asteroids a couple hundred meters in diameter should retain their large ejecta fragments in a jumbled rubble-pile state. For somewhat larger bodies, the ejection velocities are sufficiently low that large fragments are essentially retained in place, possibly forming ordered “brick-pile” structures.     

Westward rotation of the atmospheric angular momentum vector of Titan by thermal tides

The orientation of the atmospheric angular momentum vector of Titan and its temporal variation predicted by a general circulation model are analysed and interpreted. The atmospheric angular momentum vector is tilted by a few degrees from the polar axis and the vector rotates (precesses) westward with a constant period of 1 Titan day. The fast westward rotation is likely to be caused by migrating diurnal thermal tides. The tilt is almost cancelled out in the troposphere by the wavenumber 2 pattern of Saturn's gravitational tide, but is more pronounced in the stratosphere, where thermal tides are significant. The predicted tilt angle and the equatorial angular momentum vary with season and maximize when the hemispheric asymmetry of the axial angular momentum or superrotation attains its peak.     

Fragmentation of a Molecular Cloud Core versus Fragmentation of the Massive Protoplanetary Disk in the Main Accretion Phase

We present three-dimensional numerical simulations on binary formation through fragmentation. The simulations follow gravitational collapse of a molecular cloud core up to growth of the first core by accretion. At the initial stage, the gravity is only slightly dominant over the gas pressure. We made various models by changing initial velocity distribution (rotation speed, rotation law, and bar-mode perturbation). The cloud fragments whenever the cloud rotates sufficiently slowly to allow collapse but faster enough to form a disk before first-core formation. The latter condition is equivalent to Ω₀ t _ff ? 0.05, where Ω₀ and t _ff f denote the initial central angular velocity and the freefall time measured from the central density, and the condition is independent of the initial rotation law and bar-mode perturbation. Fragmentation is classified into six types. When the initial cloud rotates rigidly the cloud collapses to form a adiabatic disk supported by rotation. When the bar-mode perturbation is very minor, the disk deforms to a rotating bar, and the bar fragments. Otherwise, the adiabatic disk evolves into a central core surrounded by a circumstellar disk, and the the circumstellar disk fragments. When the initial cloud rotates differentially, the cloud deforms to a ring or bar in the isothermal collapse phase. The ring fragments into free or more cores, while the bar fragments into only two cores. In the latter case, the core merges due to low orbital angular momentum and new satellite cores form in the later stages.