The New Horizons mission has returned stunning images of the bilobate Kuiper belt object (486958) Arrokoth. It is a contact binary, formed from two intact and relatively undisturbed predecessor objects joined by a narrow contact region. We use a version of pkdgrav, an N-body code that allows for soft-sphere collisions between particles, to model a variety of possible merger scenarios with the aim of constraining how Arrokoth may have evolved from two Kuiper belt objects into its current contact binary configuration. We find that the impact must have been quite slow (less than 5 m/s) and grazing (impact angles greater than 75 degrees) in order to leave intact lobes after the merger, in the case that both progenitor objects were rubble piles. A gentle contact between two bodies in a close synchronous orbit seems most plausible.

For problems in astrophysics, planetary science and beyond, numerical simulations are often limited to simulating fewer particles than in the real system. To model collisions, the simulated particles (aka superparticles) need to be inflated to represent a collectively large collisional cross section of real particles. Here we develop a superparticle-based method that replicates the kinetic energy loss during real-world collisions, implement it in an $N$-body code and test it. The tests provide interesting insights into dynamics of self gravitating collisional systems. They show how particle systems evolve over several free fall timescales to form central concentrations and equilibrated outer shells. The superparticle method can be extended to account for the accretional growth of objects during inelastic mergers.

The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

The New Horizons spacecraft's encounter with the cold classical Kuiper belt object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary planetesimal. We investigate how it formed, finding it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collapsing, solid particle cloud. The geometric alignment of the lobes indicates the lobes were a co-orbiting binary that experienced angular momentum loss and subsequent merger, possibly due to dynamical friction and collisions within the cloud or later gas drag. Arrokoth's contact-binary shape was preserved by the benign dynamical and collisional environment of the cold classical Kuiper belt, and so informs the accretion processes that operated in the early Solar System.

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) [...].

Understanding the collisional fragmentation and subsequent reaccumulation of fragments is crucial for studies of the formation and evolution of the small-body populations. Using an SPH / N-body approach, we investigate the size-frequency distributions (SFDs) resulting from the disruption of 100 km-diameter targets consisting of porous material, including the effects of pore-crushing as well as friction. Overall, the porous targets have a significantly higher impact strength (Q*D) than the rubble-pile parent bodies investigated previously (Benavidez et al., 2012) and show a behavior more similar to non-porous monolithic targets (Durda et al., 2007). Our results also confirm that for a given specific impact energy, the SFDs resulting from a parent body disruption are strongly dependent on the size scale.

The optical images of near-Earth asteroid 4179 Toutatis acquired by Chang'e-2 spacecraft show that Toutatis has an elongated contact binary configuration, with the contact point located along the long axis. We speculate that such configuration may have resulted from a low-speed impact between two components. In this work, we performed a series of numerical simulations and compared the results with the optical images, to examine the mechanism and better understand the formation of Toutatis. Herein we propose an scenario that an assumed separated binary precursor could undergo a close encounter with Earth, leading to an impact between the primary and secondary, and the elongation is caused by Earth's tide. The precursor is assumed to be a doubly synchronous binary with a semi-major axis of 4 $R_p$ (radius of primary) and the two components are represented as spherical cohesionless self-gravitating granular aggregates. The mutual orbits are simulated in a Monte Carlo routine to provide appropriate parameters for our \emph{N}-body simulations of impact and tidal distortion. We employ the \emph{pkdgrav} package with a soft-sphere discrete element method (SSDEM) to explore the entire scenarios. The results show that contact binary configurations are natural outcomes under this scenario, whereas the shape of the primary is almost not affected by the impact of the secondary. However, our simulations further provide an elongated contact binary configuration best-matching to the shape of Toutatis at an approaching distance $r_p$ = 1.4 $\sim$ 1.5 $R_e$ (Earth radius), indicative of a likely formation scenario for configurations of Toutatis-like elongated contact binaries.

We report on our study of asteroidal breakups, i.e. fragmentations of targets, subsequent gravitational reaccumulation and formation of small asteroid families. We focused on parent bodies with diameters $D_{\rm pb} = 10$ km. Simulations were performed with a smoothed-particle hydrodynamics (SPH) code combined with an efficient N-body integrator. We assumed various projectile sizes, impact velocities and impact angles (125 runs in total). Resulting size-frequency distributions are significantly different from scaled-down simulations with $D_{\rm pb} = 100$ km targets (Durda et al., 2007). We derive new parametric relations describing fragment distributions, suitable for Monte-Carlo collisional models. We also characterize velocity fields and angular distributions of fragments, which can be used as initial conditions for N-body simulations of small asteroid families. Finally, we discuss a number of uncertainties related to SPH simulations.

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.

Solar system small bodies come in a wide variety of shapes and sizes, which are achieved following very individual evolutional paths through billions of years. This paper focuses on the reshaping process of rubble-pile asteroids driven by meteorite impacts. In our study, numerous possible equilibrium configurations are obtained via Monte Carlo simulation, and the structural stability of these configurations is determined via eigen analysis of the geometric constructions. The eigen decomposition reveals a connection between the cluster's reactions and the types of external disturbance. Numerical simulations are performed to verify the analytical results. The gravitational N-body code pkdgrav is used to mimic the responses of the cluster under intermittent non-dispersive impacts. We statistically confirm that the stability index, the total gravitational potential and the volume of inertia ellipsoid show consistent tendency of variation. A common regime is found in which the clusters tend towards crystallization under intermittent impacts, i.e., only the configurations with high structural stability survive under the external disturbances. The results suggest the trivial non-disruptive impacts might play an important role in the rearrangement of the constituent blocks, which may strengthen these rubble piles and help to build a robust structure under impacts of similar magnitude. The final part of this study consists of systematic simulations over two parameters, the projectile momentum and the rotational speed of the cluster. The results show a critical value exists for the projectile momentum, as predicted by theory, below which all clusters become responseless to external disturbances; and the rotation proves to be significant for it exhibits an "enhancing" effect on loose-packed clusters, which coincides with the observation that several fast-spinning asteroids have low bulk densities.

As the target of the proposed Asteroid Impact & Deflection Assessment (AIDA) mission, the near-Earth binary asteroid 65803 Didymos represents a special class of binary asteroids, those whose primaries are at risk of rotational disruption. To gain a better understanding of these binary systems and to support the AIDA mission, this paper investigates the creep stability of the Didymos primary by representing it as a cohesionless self-gravitating granular aggregate subject to rotational acceleration. To achieve this goal, a soft-sphere discrete element model (SSDEM) capable of simulating granular systems in quasi-static states is implemented and a quasi-static spin-up procedure is carried out. We devise three critical spin limits for the simulated aggregates to indicate their critical states triggered by reshaping and surface shedding, internal structural deformation, and shear failure, respectively. The failure condition and mode, and shear strength of an aggregate can all be inferred from the three critical spin limits. The effects of arrangement and size distribution of constituent particles, bulk density, spin-up path, and interparticle friction are numerically explored. The results show that the shear strength of a spinning self-gravitating aggregate depends strongly on both its internal configuration and material parameters, while its failure mode and mechanism are mainly affected by its internal configuration. Additionally, this study provides some constraints on the possible physical properties of the Didymos primary based on observational data and proposes a plausible formation mechanism for this binary system. With a bulk density consistent with observational uncertainty and close to the maximum density allowed for the asteroid, the Didymos primary in certain configurations can remain geo-statically stable without including cohesion.

We study the B ring's complex optical depth structure. The source of this structure may be the complex dynamics of the Keplerian shear and the self-gravity of the ring particles. The outcome of these dynamic effects depends sensitively on the collisional and physical properties of the particles. Two mechanisms can emerge that dominate the macroscopic physical structure of the ring: self-gravity wakes and viscous overstability. Here we study the interplay between these two mechanisms by using our recently developed particle collision method that allows us to better model the inter-particle contact physics. We find that for a constant ring surface density and particle internal density, particles with rough surfaces tend to produce axisymmetric ring features associated with the viscous overstability, while particles with smoother surfaces produce self-gravity wakes.

In this work, we analyze Jovian Trojans in the space of suitable resonant elements and we identify clusters of possible collisional origin by two independent methods: the hierarchical clustering and a so-called "randombox". Compared to our previous work (Bro\v{z} and Rozehnal 2011), we study a twice larger sample. Apart from Eurybates, Ennomos and $1996\,\rm RJ$ families, we have found three more clusters --- namely families around asteroids (20961)~Arkesilaos, (624)~Hektor in the $L_4$ libration zone and (247341)~$2001\,\rm UV_{209}$ in $L_5$. The families fulfill our stringent criteria, i.e. a high statistical significance, an albedo homogeneity and a steeper size-frequency distribution than that of background. In order to understand their nature, we simulate their long term collisional evolution with the Boulder code (Morbidelli et al. 2009) and dynamical evolution using a modified SWIFT integrator (Levison and Duncan, 1994). Within the framework of our evolutionary model, we were able to constrain the the age of the Hektor family to be either 1 to 4 Gyr or, less likely, 0.1 to 2.5 Gyr, depending on initial impact geometry. Since (624) Hektor itself seems to be a bilobed--shape body with a satellite (Marchis et al. 2014), i.e. an exceptional object, we address its association with the D--type family and we demonstrate that the moon and family could be created during a single impact event. We simulated the cratering event using a Smoothed Particle Hydrodynamics (SPH, Benz and Asphaug, 1994). This is also the first case of a family associated with a D--type parent body.

In order to understand the rate of merger of stellar-mass black hole binaries (BHBs) by gravitational wave (GW) emission it is important to determine the major pathways to merger. We use numerical simulations to explore the evolution of BHBs inside the radius of influence of supermassive black holes (SMBHs) in galactic centers. In this region the evolution of binaries is dominated by perturbations from the central SMBH. In particular, as first pointed out by Antonini and Perets, the Kozai-Lidov (KL) mechanism trades relative inclination of the BHB to the SMBH for eccentricity of the BHB, and for some orientations can bring the BHB to an eccentricity near unity. At very high eccentricities, GW emission from the BHB can become efficient, causing the members of the BHB to coalesce. We use a novel combination of two N-body codes to follow this evolution. We are forced to simulate small systems to follow the behavior accurately. We have completed 400 simulations that range from $\sim$ 300 stars around a $10^{3}$ M$_{\odot}$ black hole to $\sim$ 4500 stars around a $10^{4}$ M$_{\odot}$ black hole. These simulations are the first to follow the internal orbit of a binary near a SMBH while also following the changes to its external orbit self-consistently. We find that this mechanism could produce mergers at a maximum rate per volume of $\sim 100$ Gpc$^{-3}$ yr$^{-1}$ or considerably less if the inclination oscillations of the binary remain constant as the BHB inclination to the SMBH changes, or if the binary black hole fraction is small.

In the last decade, thanks to the development of sophisticated numerical codes, major breakthroughs have been achieved in our understanding of the formation of asteroid families by catastrophic disruption of large parent bodies. In this review, we describe numerical simulations of asteroid collisions that reproduced the main properties of families, accounting for both the fragmentation of an asteroid at the time of impact and the subsequent gravitational interactions of the generated fragments. The simulations demonstrate that the catastrophic disruption of bodies larger than a few hundred meters in diameter leads to the formation of large aggregates due to gravitational reaccumulation of smaller fragments, which helps explain the presence of large members within asteroid families. Thus, for the first time, numerical simulations successfully reproduced the sizes and ejection velocities of members of representative families. Moreover, the simulations provide constraints on the family dynamical histories and on the possible internal structure of family members and their parent bodies.

We present the first self-consistent simulations of the coupled spin-shape evolution of small gravitational aggregates under the influence of the YORP effect. Because of YORP's sensitivity to surface topography, even small centrifugally driven reconfigurations of aggregates can alter the YORP torque dramatically, resulting in spin evolution that can differ qualitatively from the rigid-body prediction. One third of our simulations follow a simple evolution described as a modified YORP cycle. Two-thirds exhibit one or more of three distinct behaviors---stochastic YORP, self-governed YORP, and stagnating YORP---which together result in YORP self-limitation. Self-limitation confines rotation rates of evolving aggregates to far narrower ranges than those expected in the classical YORP cycle, greatly prolonging the times over which objects can preserve their sense of rotation. Simulated objects are initially randomly packed, disordered aggregates of identical spheres in rotating equilibrium, with low internal angles of friction. Their shape evolution is characterized by rearrangement of the entire body, including the deep interior. They do not evolve to axisymmetric top shapes with equatorial ridges. Mass loss occurs in one-third of the simulations, typically in small amounts from the ends of a prolate-triaxial body. We conjecture that YORP self-limitation may inhibit formation of top-shapes, binaries, or both, by restricting the amount of angular momentum that can be imparted to a deformable body. Stochastic YORP, in particular, will affect the evolution of collisional families whose orbits drift apart under the influence of Yarkovsky forces, in observable ways.

Our knowledge of the strengths of small bodies in the Solar System is limited by our poor understanding of their internal structures, and this, in turn, clouds our understanding of the formation and evolution of these bodies. Observations of the rotational states of asteroids whose diameters are larger than a few hundreds of meters have revealed that they are dominated by gravity and that most are unlikely to be monoliths; however, there is a wide range of plausible internal structures. Numerical and analytical studies of shape and spin limits of gravitational aggregates and their collisional evolution show a strong dependence on shear strength. In order to study this effect, we carry out a systematic exploration of the dependence of collision outcomes on dissipation and friction parameters of the material components making up the bodies. We simulate the catastrophic disruption (leading to the largest remnant retaining 50% of the original mass) of km-size asteroids modeled as gravitational aggregates using pkdgrav, a cosmology N-body code adapted to collisional problems and recently enhanced with a new soft-sphere collision algorithm that includes more realistic contact forces. We find that for a range of three different materials, higher friction and dissipation values increase the catastrophic disruption threshold by about half a magnitude. Furthermore, we find that pre-impact rotation systematically increases mass loss on average, regardless of the target's internal configuration. Our results have important implications for the efficiency of planet formation via planetesimal growth, and also more generally to estimate the impact energy threshold for catastrophic disruption, as this generally has only been evaluated for non-spinning bodies without detailed consideration of material properties.

This study is carried out in the framework of sample-return missions to asteroids that use a low-speed projectile as the primary component of its sampling mechanism (e.g., JAXA's Hayabusa and Hayabusa2 missions). We perform numerical simulations of such impacts into granular materials using different projectile shapes under Earth's gravity. We then compare the amounts of ejected mass obtained in our simulations against what was found in experiments that used similar setups, which allows us to validate our numerical approach. For the targets, we consider 2 different monodisperse grain-diameter sizes: 5 mm and 3 mm. The impact speed of the projectile is 11 m s$^{-1}$ directed downward, perpendicular to the surface of the targets. Using an implementation of the soft-sphere discrete element method (SSDEM) in the $N$-Body gravity tree code PKDGRAV, previously validated in the context of low-speed impacts into sintered glass bead agglomerates, we find a noticeable dependence of the amount of ejected mass on the projectile shape. As found in experiments, in the case of the larger target grain size (5 mm), a conically shaped projectile ejects a greater amount of mass than do projectiles of other shapes, including disks and spheres. We then find that numerically the results are sensitive to the normal coefficient of restitution of the grains, especially for impacts into targets comprised of smaller grains (3 mm). We also find that static friction plays a more important role for impacts into targets comprised of the larger grains. As a preliminary demonstration, one of these considered setups is simulated in a microgravity environment. As expected, a reduction in gravity increases both the amount of ejected mass and the timescale of the impact process. Our methodology is also adaptable to the conditions of sampling mechanisms included in specific mission designs.

Asteroid (99942) Apophis' close approach in 2029 will be one of the most significant small-body encounter events in the near future and offers a good opportunity for in situ exploration to determine the asteroid's surface properties and measure any tidal effects that might alter its regolith configuration. Resurfacing mechanics has become a new focus for asteroid researchers due to its important implications for interpreting surface observations, including space weathering effects. This paper provides a prediction for the tidal effects during the 2029 encounter, with an emphasis on whether surface refreshing due to regolith movement will occur. The potential shape modification of the object due to the tidal encounter is first confirmed to be negligibly small with systematic simulations, thus only the external perturbations are taken into account for this work (despite this, seismic shaking induced by shifting blocks might still play a weak role and we will look into this mechanism in future work). A two-stage approach is developed to model the responses of asteroid surface particles (the regolith) based on the soft-sphere implementation of the parallel N-body gravity tree code pkdgrav. A full-body model of Apophis is sent past the Earth on the predicted trajectory to generate the data of all forces acting at a target point on the surface. A sandpile constructed in the local frame is then used to approximate the regolith materials; all the forces the sandpile feels during the encounter are imposed as external perturbations to mimic the regolith's behavior in the full scenario. The local mechanical environment on the asteroid surface is represented in detail, leading to an estimation of the change in global surface environment due to the encounter. Typical patterns of perturbation are presented that depend on the asteroid orientation and sense of rotation at perigee.

Out of the handful of asteroids that have been imaged, some have distributions of blocks that are not easily explained. In this paper, we investigate the possibility that seismic shaking leads to the size sorting of particles in asteroids. In particular, we focus on the so-called Brazil Nut Effect (BNE) that separates large particles from small ones under vibrations. We study the BNE over a wide range of parameters by using the N-body code PKDGRAV, and find that the effect is largely insensitive to the coefficients of restitution, but sensitive to friction constants and oscillation speeds. Agreeing with the previous results, we find that convection drives the BNE, where the intruder rises to the top of the particle bed. For the wide-cylinder case, we also observe a "whale" effect, where the intruder follows the convective current and does not stay at the surface. We show that the non-dimensional critical conditions for the BNE agree well with previous studies. We also show that the BNE is scalable for low-gravity environments and that the rise speed of an intruder is proportional to the square root of the gravitational acceleration. Finally, we apply the critical conditions to observed asteroids, and find that the critical oscillation speeds are comparable to the seismic oscillation speeds that are expected from non-destructive impacts.

We carry out a systematic exploration of the effect of pre-impact rotation on the outcomes of low-speed collisions between planetesimals modeled as gravitational aggregates. We use pkdgrav, a cosmology code adapted to collisional problems and recently enhanced with a new soft-sphere collision algorithm that includes more realistic contact forces. A rotating body has lower effective surface gravity than a non-rotating one and therefore might suffer more mass loss as the result of a collision. What is less well understood, however, is whether rotation systematically increases mass loss on average regardless of the impact trajectory. This has important implications for the efficiency of planet formation via planetesimal growth, and also more generally for the determination of the impact energy threshold for catastrophic disruption (leading to the largest remnant retaining 50% of the original mass), as this has generally only been evaluated for non-spinning bodies. We find that for most collision scenarios, rotation lowers the threshold energy for catastrophic dispersal. For head-on collisions, we develop a semi-analytic description of the change in the threshold description as a function of the target's pre-impact rotation rate, and find that these results are consistent with the "universal law" of catastrophic disruption developed by Leinhardt & Stewart. Using this approach, we introduce re-scaled catastrophic disruption variables that take into account the interacting mass fraction of the target and the projectile in order to translate oblique impacts into equivalent head-on collisions.

The absolute-magnitude distributions of seven supernova types are presented. The data used here were primarily taken from the Asiago Supernova Catalogue, but were supplemented with additional data. We accounted for both foreground and host-galaxy extinction. A bootstrap method is used to correct the samples for Malmquist bias. Separately, we generate volume-limited samples, restricted to events within 100 Mpc. We find that the superluminous events (M_B < -21) make up only about 0.1% of all supernovae in the bias-corrected sample. The subluminous events (M_B > -15) make up about 3%. The normal Ia distribution was the brightest with a mean absolute blue magnitude of -19.25. The IIP distribution was the dimmest at -16.75.

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular "atmosphere" that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate.

Surfaces of planets and small bodies of our Solar System are often covered by a layer of granular material that can range from a fine regolith to a gravel-like structure of varying depths. Therefore, the dynamics of granular materials are involved in many events occurring during planetary and small-body evolution thus contributing to their geological properties. We demonstrate that the new adaptation of the parallel N-body hard-sphere code pkdgrav has the capability to model accurately the key features of the collective motion of bidisperse granular materials in a dense regime as a result of shaking. As a stringent test of the numerical code we investigate the complex collective ordering and motion of granular material by direct comparison with laboratory experiments. We demonstrate that, as experimentally observed, the scale of the collective motion increases with increasing small-particle additive concentration. We then extend our investigations to assess how self-gravity and external gravity affect collective motion. In our reduced-gravity simulations both the gravitational conditions and the frequency of the vibrations roughly match the conditions on asteroids subjected to seismic shaking, though real regolith is likely to be much more heterogeneous and less ordered than in our idealised simulations. We also show that collective motion can occur in a granular material under a wide range of inter-particle gravity conditions and in the absence of an external gravitational field. These investigations demonstrate the great interest of being able to simulate conditions that are to relevant planetary science yet unreachable by Earth-based laboratory experiments.

We present several lines of evidence based on different kinds of observations to conclude that rotational fission has likely occurred for a fraction of the known Trans-Neptunian Objects (TNOs). It is also likely that a number of binary systems have formed from that process in the trans-neptunian belt. We show that Haumea is a potential example of an object that has suffered a rotational fission. Its current fast spin would be a slight evolution of a primordial fast spin, rather than the result of a catastrophic collision, because the percentage of objects rotating faster than 4 hours would not be small in a maxwellian distribution of spin rates that fits the current TNO rotation database. On the other hand, the specific total angular momentum of Haumea and its satellites falls close to that of the high size ratio asteroid binaries, which are thought to be the result of rotational fissions or mass shedding. We also present N-body simulations of rotational fissions applied to the case of Haumea, which show that this process is feasible, might have generated satellites, and might have even created a "family" of bodies orbitally associated to Haumea. The orbitally associated bodies may come from the direct ejection of fragments according to our simulations, or through the evolution of a proto-satellite formed during the fission event. Also, the disruption of an escaped fragment after the fission might create the orbitally related bodies. If any of those mechanisms are correct, other rotational fission families may be detectable in the trans-neptunian belt in the future, and perhaps even TNO pairs might be found (pairs of bodies sharing very similar orbital elements, but not bound together).

Space weathering is the generic term used for processes that modify the optical properties of surfaces of atmosphereless rocky bodies under exposure to the space environment. The general agreement about the relevance of the effects of space weathering on the spectral properties of S-complex asteroids fails when some basic quantitative estimates are attempted. In particular, there is severe disagreement regarding the typical timescales for significant spectral reddening to occur, ranging from 1 Myr to 1 Gyr. Generally speaking, the spectral reddening of an individual object can be considered as the sum of three terms, one (which is relevant for statistical analyses) depending on the exposure of the object to space weathering during its lifetime, a second one due to the original surface composition, and a third one (a "noise" term) due to the combination of poorly constrained effects (e.g., structure and texture of the surface). The surface of an asteroid is usually covered by regolith, and its presence and properties presumably play a critical role in the weathering processes. In this paper we discuss the role played by collisional evolution in affecting the spectral properties of asteroids and refreshing the surfaces due to the formation of ejecta, and the necessity of a simultaneous modeling of collisions and weathering processes. We introduce a new idea, based on the possibility of a sort of saturation of the refreshing process whenever a massive reaccumulation of the impact ejecta takes place. In this case, a dependence of the overall reddening on the asteroid size should naturally come out. We show that this conclusion is indeed supported by available main belt asteroid spectroscopic data.

A large fraction of 100-km-class low-inclination objects in the classical Kuiper Belt (KB) are binaries with comparable mass and wide separation of components. A favored model for their formation was capture during the coagulation growth of bodies in the early KB. Instead, recent studies suggested that large objects can rapidly form in the protoplanetary disks when swarms of locally concentrated solids collapse under their own gravity. Here we examine the possibility that KB binaries formed during gravitational collapse when the excess of angular momentum prevented the agglomeration of available mass into a solitary object. We find that this new mechanism provides a robust path toward the formation of KB binaries with observed properties, and can explain wide systems such as 2001 QW322 and multiples such as (47171) 1999 TC36. Notably, the gravitational collapse is capable of producing 100% binary fraction for a wide range of the swarm's initial angular momentum values. The binary components have similar masses (80% have the secondary-over-primary radius ratio >0.7) and their separation ranges from ~1,000 to ~100,000 km. The binary orbits have eccentricities from e=0 to ~1, with the majority having e<0.6. The binary orbit inclinations with respect to the initial angular momentum of the swarm range from i=0 to ~90 deg, with most cases having i<50 deg. Our binary formation mechanism implies that the primary and secondary components in each binary pair should have identical bulk composition, which is consistent with the current photometric data. We discuss the applicability of our results to the Pluto-Charon, Orcus-Vanth, (617) Patroclus-Menoetius and (90) Antiope binary systems.

Hill's equations are an approximation that is useful in a number of areas of astrophysics including planetary rings and planetesimal disks. We derive a symplectic method for integrating Hill's equations based on a generalized leapfrog. This method is implemented in the parallel N-body code, PKDGRAV and tested on some simple orbits. The method demonstrates a lack of secular changes in orbital elements, making it a very useful technique for integrating Hill's equations over many dynamical times. Furthermore, the method allows for efficient collision searching using linear extrapolation of particle positions.

Numerical simulations of asteroid break-ups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by Q*D, which results in the escape of half of the target's mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize Q*D for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of Q*D, and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of Q*D on the impact speed of the projectile. In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. We then propose some power-law relationships between Q*D and both target's size and impact speed that can be used in collisional evolution models.

We have developed a semi-analytic method of parameterizing N-body simulations of self-gravity wakes in Saturn's rings, describing their properties by means of only 6 numbers: 3 optical depths and 3 weighting factors. These numbers are obtained using a density-estimation procedure that finds the frequencies of various values of local density within a simulated ring patch. Application of our parameterization to a suite of N-body simulations implies that the distribution of local optical depths is trimodal, rather than bimodal as previous authors have assumed. Rings dominated by self-gravity wakes appear to be mostly empty space. The implications of this result for the analysis of occultation data are more conceptual than practical. The only adjustment needed is that the model parameter $\tau_{gap}$ should be interpreted as representing the area-weighted average optical depth within the gaps (or inter-wake regions). The most significant consequence of our results applies to the question of why "propeller" structures observed in the mid-A ring are seen as relative-bright features, even though the most prominent features of simulated propellers are regions of relatively low density. We find preliminary quantitative support for the hypothesis that propellers would be bright if they involve a local and temporary disruption of self-gravity wakes, flooding the region with more "photometrically active" material (i.e., material that can contribute to the rings' local optical depth) even though the overall density is lower, and thus raising their apparent brightnesses in agreement with observations. We suggest that this mechanism be tested by future detailed numerical models.

Asteroids are leftover pieces from the era of planet formation that help us understand conditions in the early Solar System. Unlike larger planetary bodies that were subject to global thermal modification during and subsequent to their formation, these small bodies have kept at least some unmodified primordial material from the solar nebula. However, the structural properties of asteroids have been modified considerably since their formation. Thus, we can find among them a great variety of physical configurations and dynamical histories. In fact, with only a few possible exceptions, all asteroids have been modified or completely disrupted many times during the age of the Solar System. This picture is supported by data from space mission encounters with asteroids that show much diversity of shape, bulk density, surface morphology, and other features. Moreover, the gravitational attraction of these bodies is so small that some physical processes occur in a manner far removed from our common experience on Earth. Thus, each visit to a small body has generated as many questions as it has answered. In this review we discuss the current state of research into asteroid disruption processes, focusing on collisional and rotational mechanisms. We find that recent advances in modeling catastrophic disruption by collisions have provided important insights into asteroid internal structures and a deeper understanding of asteroid families. Rotational disruption, by tidal encounters or thermal effects, is responsible for altering many smaller asteroids, and is at the origin of many binary asteroids and oddly shaped bodies.

We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 solar mass star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 10^5 particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation.

In this paper we extend our numerical method for simulating terrestrial planet formation from Leinhardt and Richardson (2005) to include dynamical friction from the unresolved debris component. In the previous work we implemented a rubble pile planetesimal collision model into direct N-body simulations of terrestrial planet formation. The new collision model treated both accretion and erosion of planetesimals but did not include dynamical friction from debris particles smaller than the resolution limit for the simulation. By extending our numerical model to include dynamical friction from the unresolved debris, we can simulate the dynamical effect of debris produced during collisions and can also investigate the effect of initial debris mass on terrestrial planet formation. We find that significant initial debris mass, 10% or more of the total disk mass, changes the mode of planetesimal growth. Specifically, planetesimals in this situation do not go through a runaway growth phase. Instead they grow concurrently, similar to oligarchic growth. In addition to including the dynamical friction from the unresolved debris, we have implemented particle tracking as a proxy for monitoring compositional mixing. Although there is much less mixing due to collisions and gravitational scattering when dynamical friction of the background debris is included, there is significant inward migration of the largest protoplanets in the most extreme initial conditions.

Photometry data were collected from the literature and analyzed for supernovae that are thought to have a gamma-ray burst association. There are several gamma-ray burst afterglow light curves that appear to have a supernova component. For these light curves, the supernova component was extracted and analyzed. A supernova light curve model was used to help determine the peak absolute magnitudes as well as estimates for the kinetic energy, ejected mass and nickel mass in the explosion. The peak absolute magnitudes are, on average, brighter than those of similar supernovae (stripped-envelope supernovae) that do not have a gamma-ray burst association, but this can easily be due to a selection effect. However, the kinetic energies and ejected masses were found to be considerably higher, on average, than those of similar supernovae without a gamma-ray burst association.

Asteroid shapes and satellites: role of gravitational reaccumulation. Following current evidences, it is widely accepted that many asteroids would be "gravitational aggregates", i.e. bodies lacking internal cohesion. They could mainly be originated during the catastrophic disruption of some parent bodies, through the gravitational reaccumulation of the resulting fragments. The same events produced the dynamical families that we observe. In this work we address the problem of the origin of shapes of gravitational aggregates, that could contain signatures of their origin. We use a N-body code to simulate the collapse of a cloud of fragments, with a variety of initial velocity distributions and total angular momentum. The fragments are treated as inhelastic spheres, that rapidly accumulate to form rotating aggregates. The resulting shapes and rotational properties are compared with theoretical predictions. The results show that only a precise category of shapes (flattened spheroids) are created via this mechanism. This may provide interesting constraints on the evolution of asteroid shapes, in particular for those with one or more satellites.

A description of the light curves of microlensed Type Ia supernovae (SNe Ia) as extended and expanding sources is presented. We give examples of what microlensing by stellar-mass Chang-Refsdal lenses can do to a small percentage of supernova light curves. We find that in addition to overall brightening, significant changes in light--curve shapes can also occur. Peaks can be distorted, plateaus can appear, and even secondary peaks can be formed. The effects of both the relative motion of the lens and the supernova and the expansion of the supernova are given and compared. The effects of relative motion are more pronounced when a distant supernova ($z_s \sim 1$) impacts well within the Einstein ring of a nearby microlens ($z_d \sim 0.05$) and are less important for more distant deflectors. We also find that the increase in shear that comes with increased deflector distance tends to reduce the time variability of microlensing. We briefly discuss the probability of observing these effects.

The absolute visual magnitudes of three Type IIb, 11 Type Ib and 13 Type Ic supernovae (collectively known as stripped-envelope supernovae) are studied by collecting data on the apparent magnitude, distance, and interstellar extinction of each event. Weighted and unweighted mean absolute magnitudes of the combined sample as well as various subsets of the sample are reported. The limited sample size and the considerable uncertainties, especially those associated with extinction in the host galaxies, prevent firm conclusions regarding differences between the absolute magnitudes of supernovae of Type Ib and Ic, and regarding the existence of separate groups of overluminous and normal-luminosity stripped-envelope supernovae. The spectroscopic characteristics of the events of the sample are considered. Three of the four overluminous events are known to have had unusual spectra. Most but not all of the normal luminosity events had typical spectra. Light curves of stripped-envelope supernovae are collected and compared. Because SN 1994I in M51 was very well observed it often is regarded as the prototypical Type Ic supernova, but it has the fastest light curve in the sample. Light curves are modeled by means of a simple analytical technique that, combined with a constraint on E/M from spectroscopy, yields internally consistent values of ejected mass, kinetic energy, and nickel mass.

A detailed description of the apparent light curves of microlensed SNe Ia as extended and expanding sources is presented. We show that microlensing amplification can have significant effects on supernova observation. A model light curve is used to compare lensed and unlensed cases and we find that significant changes in shape can occur because of microlensing. We briefly discuss the probability of observing such effects as well. Throughout the paper we consider the redshift of the source to be 1.0 and we limit our work to spherically symmetric deflectors.

For a long time, gravitational instability in the disk of planetesimals has been suspected to be the main engine responsible for the beginning of dust growth, its advantage being that it provides for rapid growth. Its real importance in planetary formation is still debated, mainly because the potential presence of turbulence can prevent the settling of particles into a gravitationally unstable layer. However, several mechanisms could yield strongly inhomogeneous distributions of solids in the disk: radial drift, trapping in vortices, perturbations by other massive bodies, etc. In this paper we present a numerical study of a gravitationally unstable layer. This allows us to go beyond the classical analytical study of linear perturbations, exploring a highly non-linear regime. A hierarchical growth of structure in the presence of dissipation (gas drag) can yield large, virialized clusters of planetesimals, the first time such clusters have been observed in the context of planetesimal disks.

The Asiago Supernova Catalog is used to carry out a comparative study of supernova absolute-magnitude distributions. An overview of the absolute magnitudes of the supernovae in the current observational sample is presented, and the evidence for subluminous and overluminous events is examined. The fraction of supernovae that are underluminous (M_B > -15) appears to be higher (perhaps much higher) than one fifth but it remains very uncertain. The fraction that are overluminous (M_B < -20) is lower (probably much lower) than 0.01. The absolute-magnitude distributions for each supernova type, restricted to events within 1 Gpc, are compared. Although these distributions are affected by observational bias in favor of the more luminous events, they are useful for comparative studies. We find mean absolute blue magnitudes (for H_0=60) of -19.46 for normal Type Ia supernovae (SNe Ia), -18.04 for SNe Ibc, -17.61 and -20.26 for normal and bright SNe Ibc considered separately, -18.03 for SNe II-L, -17.56 and -19.27 for normal and bright SNe II-L considered separately, -17.00 for SNe II-P, and -19.15 for SNe IIn.

There is increasing evidence that many km-sized bodies in the Solar System are piles of rubble bound together by gravity. We present results from a project to map the parameter space of collisions between km-sized spherical rubble piles. The results will assist in parameterization of collision outcomes for Solar System formation models and give insight into fragmentation scaling laws. We use a direct numerical method to evolve the positions and velocities of the rubble pile particles under the constraints of gravity and physical collisions. We test the dependence of the collision outcomes on impact parameter and speed, impactor spin, mass ratio, and coefficient of restitution. Speeds are kept low (< 10 m/s, appropriate for dynamically cool systems such as the primordial disk during early planet formation) so that the maximum strain on the component material does not exceed the crushing strength. We compare our results with analytic estimates and hydrocode simulations. Off-axis collisions can result in fast-spinning elongated remnants or contact binaries while fast collisions result in smaller fragments overall. Clumping of debris escaping from the remnant can occur, leading to the formation of smaller rubble piles. In the cases we tested, less than 2% of the system mass ends up orbiting the remnant. Initial spin can reduce or enhance collision outcomes, depending on the relative orientation of the spin and orbital angular momenta. We derive a relationship between impact speed and angle for critical dispersal of mass in the system. We find that our rubble piles are relatively easy to disperse, even at low impact speed, suggesting that greater dissipation is required if rubble piles are the true progenitors of protoplanets.

Until recently, most asteroids were thought to be solid bodies whose shapes were determined largely by collisions with other asteroids. It now seems that many asteroids are little more than rubble piles, held together by self-gravity; this means that their shapes may be strongly distorted by tides during close encounters with planets. Here we report on numerical simulations of encounters between a ellipsoid-shaped rubble-pile asteroid and the Earth. After an encounter, many of the simulated asteroids develop the same rotation rate and distinctive shape (i.e., highly elongated with a single convex side, tapered ends, and small protuberances swept back against the rotation direction) as 1620 Geographos. Since our numerical studies show that these events occur with some frequency, we suggest that Geographos may be a tidally distorted object. In addition, our work shows that 433 Eros, which will be visited by the NEAR spacecraft in 1999, is much like Geographos, which suggests that it too may have been molded by tides in the past.