Stability of particulate discs

Collisionally-induced amplification of density fluctuations can also produce non-axisymmetric local condensations in particulate discs if the optical thickness is between definite values. Gravitational instability occurs above this interval. The theory of both phenomena is derived from collisional equations. The conventional criterion for gravitational instability in a gaseous medium cannot be used for particulate discs, in which the equilibrium depends on the collisional energy loss. These instabilities can produce an unbounded growth in density or a gravitational coagulation of particles, but the typical consequence is the formation of highly elongated clouds which are denser than the background matter and have a relatively long lifetime before decay. The third type of instability, the thermal one, appears at low values of velocity dispersion. It only affects the random motion of particles without producing condensations.     

Numerical simulations of collisions and gravitational encounters in systems of non-identical particles

Numerical simulations of 200 mutually colliding non-identical particles indicate that the equipartition of random kinetic energy is possible only in systems having a narrow distribution of particle masses. Otherwise the random energy is concentrated on heavy particles. The form of the velocity distribution versus particle mass depends also on the elastic properties of the particles, and on the relative importance of the particle size. If the coefficient of restitution is a weakly decreasing function of impact velocity, a large difference in the equilibrium velocities of largest and smallest particles is possible. On the other hand, if the elasticity drops to a low level even in the small velocity regime, the dispersion of velocities is maintained by finite size and differential rotation, and the velocities of smallest particles are, at most, slightly larger than those of the largest ones. The results of simulations are consistent with the predictions of the collisional theory of non-identical particles (Hämeen-Anttila, 1984). The application to Saturn's rings indicates that the geometric thickness of cm-sized particles is of the order of 50 m in the rarefied regions of the rings. Without the gravitational encounters a thickness of about 30 m is derived. These estimations are made by using the latest measurements (Bridges et al., 1984) for the restitution coefficient of icy particles.     

Generalized theory of impacts in particulate systems

The theory of collisional systems is generalized for an arbitrary geometry and forces acting in the system, mixtures of different particle types, friction, small deviations from the ideal spherical form, axial rotation, finite size of particles and gravitational interactions. Terms for the formation of new particles and destruction of old ones are also included, and other unspecified parameters can be introduced. Although some approximations are made to simplify the basic equations and to avoid excessive numerical interactions, a comparison with computer simulations shows a good agreement. The tests were continued up to the optical thickness = 5.     

Negative diffusion in planetary rings with a nearby moon

This paper analyzes a process that has been observed in simulations of numerous systems where ring material is strongly perturbed by a nearby moon. If the ring particles can be imparted with a forced eccentricity on the order of 10⁻⁵ in a single pass by the moon, particle orbits are observed to move towards regions of higher density as a result of the organized collisions that occur in the dense peaks of the satellite wake. The width of the ring can decrease by as much as 90% if the forced eccentricity is greater than 3 × 10⁻⁵ and the unperturbed geometric optical depth is greater than 0.03. The fractional change in ring width is relatively insensitive to the particle size so long as the particle radius is much less than the product of the semimajor axis and the forced eccentricity. Including a power law particle size distribution with slope of −2.8 spanning a decade in particle radius reduces the fractional width change by about 10% compared to the uniform particle-size case. Adding gravitational interactions between ring particles only has a significant effect on ring confinement if the unperturbed geometric optical depth exceeds .03, but a 40% reduction in ring width is still achieved in a self-gravitating ring of geometric optical depth 0.3 if the forced eccentricity exceeds 3 × 10⁻⁵. This process does not require the material to be in resonance with the moon, nor does it have any minimum mass constraints because particle self-gravity is not required. The collisional damping of satellite wakes therefore provides a simple mechanism by which a single moon can reduce the radial extent of any ringlet that is close to it and has sufficient optical depth for collisions to be significant.     

Gravitational accretion of particles in Saturn's rings

Gravitational accretion in the rings of Saturn is studied with local N-body simulations, taking into account the dissipative impacts and gravitational forces between particles. Common estimates of accretion assume that gravitational sticking takes place beyond a certain distance (Roche distance) where the self-gravity between a pair of ring particles exceeds the disrupting tidal force of the central object, the exact value of this distance depending on the ring particles' internal density. However, the actual physical situation in the rings is more complicated, the growth and stability of the particle groups being affected also by the elasticity and friction in particle impacts, both directly via sticking probabilities and indirectly via velocity dispersion, as well as by the shape, rotational state and the internal packing density of the forming particle groups. These factors are most conveniently taken into account via N-body simulations. In our standard simulation case of identical 1 m particles with internal density of solid ice, ρ=900 kg m⁻³, following the Bridges et al., 1984 elasticity law, we find accretion beyond a=137,000-146,000 km, the smaller value referring to a distance where transient aggregates are first obtained, and the larger value to the distance where stable aggregates eventually form in every experiment lasting 50 orbital periods. Practically the same result is obtained for a constant coefficient of restitution εn=0.5. In terms of rp parameter, the sum of particle radii normalized by their mutual Hill radius, the above limit for perfect accretion corresponds to rp<0.84. Increased dissipation (εn=0.1), or inclusion of friction (tangential force 10% of normal force) shifts the accretion region inward by about 5000 km. Accretion is also more efficient in the case of size distribution: with a q=3 power law extending over a mass range of 1000, accretion shifts inward by almost 10,000 km. The aggregates forming in simulations via gradual accumulation of particles are synchronously rotating.     

Formation of Asteroid Families by Catastrophic Disruption: Simulations with Fragmentation and Gravitational Reaccumulation

This paper builds on preliminary work in which numerical simulations of the collisional disruption of large asteroids (represented by the Eunomia and Koronis family parent bodies) were performed and which accounted not only for the fragmentation of the solid body through crack propagation, but also for the mutual gravitational interaction of the resulting fragments. It was found that the parent body is first completely shattered at the end of the fragmentation phase, and then subsequent gravitational reaccumulations lead to the formation of an entire family of large and small objects with dynamical properties similar to those of the parent body. In this work, we present new and improved numerical simulations in detail. As before, we use the same numerical procedure, i.e., a 3D SPH hydrocode to compute the fragmentation phase and the parallel N-body code pkdgrav to compute the subsequent gravitational reaccumulation phase. However, this reaccumulation phase is now treated more realistically by using a merging criterion based on energy and angular momentum and by allowing dissipation to occur during fragment collisions. We also extend our previous studies to the as yet unexplored intermediate impact energy regime (represented by the Flora family formation) for which the largest fragment's mass is about half that of the parent body. Finally, we examine the robustness of the results by changing various assumptions, the numerical resolution, and different numerical parameters. We find that in the lowest impact energy regime the more realistic physical approach of reaccumulation leads to results that are statistically identical to those obtained with our previous simplistic approach. Some quantitative changes arise only as the impact energy increases such that higher relative velocities are reached during fragment collisions, but they do not modify the global outcome qualitatively. As a consequence, these new simulations confirm previous main results and still lead to the conclusion that: (1) all large family members must be made of gravitationally reaccumulated fragments; (2) the original fragment size distribution and their orbital dispersion are respectively steeper and smaller than currently observed for the real families, supporting recent studies on subsequent evolution and diffusion of family members; and (3) the formation of satellites around family members is a frequent and natural outcome of collisional processes.     

An improved and generalized theory for the collisional evolution of Keplerian systems

A calculation of collisional integrals with a higher accuracy yields excellent agreement between computer simulations and the collisional theory of Keplerian systems. Inclusion of axial rotation of particles modifies the evolution but does not introduce qualitatively new phenomena. Friction between the particles has a stabilizing influence, while deviations from an exactly spherical shape produce an opposite effect. The rotation of spherical or irregular bodies cannot prevent a final flattening of the system into a monolayer without also causing its disintegration. Computer simulations with a small number of particles do not represent the typical collisional evolution. They provide a test for the theory, but may sometimes lead to a misinterpretation of astronomical phenomena.     

Stochastic scattering and other contributions to the Sun's wake in the local hydrogen cloud

Gas streaming through the solar system experiences both destructive and scattering processes, the latter primarily in collisional interactions with the solar wind protons. The scattering interactions can be important in filling the downstream wake. They may effectively increase the velocity dispersion and also cause discrete orbit changes.The downstream intensity moment is here evaluated analytically for particles suffering a single, discrete collision, and compared with the moment from a thermal velocity dispersion (both in the absence of a central force field). The elastic scattering collisions of protons in H-gas lead to a contribution to theL backscatter from the wake equivalent to an initial thermal velocity of about 1 km s^–1, giving an intensity for cool gas of the order of 10R. This exceeds the contribution due to focussing in the solar gravitational field if the radiation pressure is not less than 0.8 of the gravitational attraction.     

Statistical mechanics of Keplerian orbits

Theoretical predictions agree with computer simulations at least for those collisional systems in which the restitution coefficient is independent of impact velocity. An uncertainty principle for the orbits restricts the validity of the theory and its predictions. Discussion of the whole theory and of computer simulations shows that a velocity-dependent restitution coefficient provides the only astronomically interesting applications of the collisional processes. The Saturnian and Uranian ring systems correspond very well to theoretical expectations if the restitution coefficient is of this type.     

Empirical formulae to describe some physical properties of small groups of protogalaxies with multiplicity

By means of identical cubic elements, we generate a partition of a volume in which a particle-based cosmological simulation is carried out. In each cubic element, we determine the gas particles with a normalized density greater than an arbitrarily chosen density threshold. By using a proximity parameter, we calculate the neighboring cubic elements and generate a list of neighbors. By imposing dynamic conditions on the gas particles, we identify gas clumps and their neighbors, so that we calculate and fit some properties of the groups so identified, including the mass, size and velocity dispersion, in terms of their multiplicity(here defined simply as the number of member galaxies). Finally, we report the value of the ratio of kinetic energy to gravitational energy of such dense gas clumps, which will be useful as initial conditions in simulations of gravitational collapse of gas clouds and clusters of gas clouds.     

Hoag’s object,remnant of a vanished bar?

The beautiful ringed Hoag’s object, named after its discoverer, is an interesting galaxy. Because of the roundness of its ring-like structure, it has been proposed to be a collisional ring galaxy; however, there is no obvious nearby culprit galaxy that could have collided with it. Considering an alternative, much gentler hypothesis, we study the development of the observed structure via a turning, bar perturbation in the disk potential. However, there is currently no obvious bar present, and rings produced by bars are typically oval. On the basis of much recent work improving our understanding of bar evolution, we assume the bar grows and then vanishes. In simulations of a disk of particles, under such a bar turning in the disk plane, we obtain a bulge core, empty void, and circular ring in the disk that mimic the observations of Hoag’s object. We conclude the inner edge of the ring is just beyond the outer Lindblad resonance (OLR) with the bar pattern speed. We estimate the amount of gas mass in the bulge core to be twice that of the ring. Our simulations indicate that the Hoag Object ring could survive at least 6 billion years after the bar vanishes.     

Velocity dispersion in particulate disks with size distribution

Quasi-equilibrium solutions for the pre-planetary disk are studied in terms of Hämeen-Anttila's theory (1984) of collisional, self-gravitating systems. The distribution of particle sizes is assumed to follow simple power-law distributions, with a power index in the range of 1.5–5.0. The treatment includes mutual impacts with a velocity dependent coefficient of restitution, as well as gravitational encounters with dynamical friction. The mean gravitational field of the disk is also taken into account. The results indicate that the energy(equi)-partition depends mainly on the index of size distribution, but is also affected by the optical thickness of the system, as well as on the vertical thickness as compared to the particle size. The vertical component of the gravitational field is found to be important, especially when the mass of the system is concentrated on the large particles.     

Stochastic diffusion in inverse square fields: General formulation for interplanetary gas

The Fokker-Planck equation for small stochastic changes to particles in Kepler orbits has to be formulated in terms of the integrals of motion. We generalize the modelling of proton and electron collisional perturbations to gas particles on trajectories through the solar system in order to include both spatial and velocity diffusion. The general solution is obtained in terms of a 4-dimensional normal distribution. Treatment of the singularity in the Fokker-Planck operator reduces the dimensionality by one. In addition to extending earlier results for anisotropic collisional heating in the thermal approximation, the present formulation gives the changes in density due to the mean repulsive force and to perturbations of trajectories (spatial diffusion). The net diffusion is almost everywhere towards the sun and the density increase is significant in the downstream hydrogen wake, particularly where destructive depletion is strong and gravitational focussing weak.     

Dust growth in protoplanetary disks-a comprehensive experimental / theoretical approach

More than a decade of dedicated experimental work on the collisional physics of protoplanetary dust has brought us to a point at which the growth of dust aggregates can-for the first time-be self-consistently and reliably modeled. In this article, the emergent collision model for protoplanetery dust aggregates, as well as the numerical model for the evolution of dust aggregates in protoplanetary disks, is reviewed. It turns out that, after a brief period of rapid collisional growth of fluffy dust aggregates to sizes of a few centimeters, the protoplanetary dust particles are subject to bouncing collisions, in which their porosity is considerably decreased. The model results also show that low-velocity fragmentation can reduce the final mass of the dust aggregates but that it does not trigger a new growth mode as discussed previously. According to the current stage of our model, the direct formation of kilometer-sized planetesimals by collisional sticking seems unlikely, implying that collective effects, such as the streaming instability and the gravitational instability in dust-enhanced regions of the protoplanetary disk, are the best candidates for the processes leading to planetesimals.     

Apse Alignment of the Uranian Rings

An explanation of the dynamical mechanism for apse alignment of the eccentric uranian rings is necessary before observations can be used to determine properties such as ring masses, particle sizes, and elasticities. The leading model (P. Goldreich and S. Tremaine 1979, Astron J.84, 1638-1641) relies on the ring self-gravity to accomplish this task, yet it yields equilibrium masses which are not in accord with Voyager radio measurements. We explore possible solutions such that the self-gravity and the collisional terms are both involved in the process of apse alignment. We consider limits that correspond to a hot and a cold ring, and we show that pressure terms may play a significant role in the equilibrium conditions for the narrow uranian rings. In the cold ring case, where the scale height of the ring near periapse is comparable to the ring particle size, we introduce a new pressure correction pertaining to a region of the ring where the particles are locked in their relative positions and jammed against their neighbors and the velocity dispersion is so low that the collisions are nearly elastic. In this case, we find a solution such that the ring self-gravity maintains apse alignment against both differential precession (m=1 mode) and the fluid pressure. We apply this model to the uranian α ring and show that, compared to the previous self-gravity model, the mass estimate for this ring increases by an order of magnitude. In the case of a hot ring, where the scale height can reach a value as much as 50 times the particle size, we find velocity dispersion profiles that result in pressure forces which act in such a way as to alter the ring equilibrium conditions, again leading to a ring mass increase of an order of magnitude. We find that such a velocity dispersion profile would require a different mechanism than is currently envisioned for establishing a heating/cooling balance in a finite-sized, inelastic particle ring. Finally, we introduce an important correction to the model of E. I. Chiang and P. Goldreich (2000, Astrophys. J.540, 1084-1090.). These authors relied on collisional forces in the last ∼100 m of an ∼10 km wide ring to increase ring equilibrium masses by up to a factor of ∼100. However, their treatment of ring edges as one-sided surface density drops leads to a strong dependence of the ring mass on the adjustable parameter λ (the length scale over which the ring's optical depth drops from order unity to zero at the edge). A treatment of the ring edges that takes into account their ridgelike structure retains the increase of ring mass of the order of ∼100 for a 10 km wide ring, while exhibiting weak dependence on λ. We conclude that a modified Chiang-Goldreich model can likely account for the masses of narrow, eccentric planetary rings; however, the role of shepherd satellites both in forming ring edges and in altering the streamline precession conditions near them needs to be explored further. It is also unclear whether a fully self-consistent ring model allows for the possibility of rings with negative eccentricity gradients.     

Expectations for Cassini observations of ring material with nearby moons

This paper describes N-body simulations of two regions of the saturnian ring system and examines what we might expect the Cassini orbiter to see in those areas. The first region is the edge of the Encke gap in the A ring that is perturbed by the satellite, Pan. Our previous simulations of this region neglected particle self-gravity [Lewis and Stewart, 2000a, Bull. Am. Astron. Soc. 34, 883]. Here we examine the interactions of the wakes caused by Pan with the wakes that form from local gravitational instabilities. We find that the two phenomena do not normally coexist and predict that measurements of particle sizes between the moon wakes should reflect the true particle size distribution of the region and not what is caused by gravitational aggregation. The region between the Encke gap edge and the first wake peak is an exception to this rule because our simulations exhibit the formation of exceptionally large gravity-induced wakes in this region. We also describe simulations of the F ring and explain the nature of braid-like structures that form naturally when the ring is perturbed by a single moon on an eccentric orbit. Finally, we discuss the very dynamic nature of the F ring system and how this should be taken into account when interpreting observations and even when planning future observations of this system.     

Planetary growth with collisional fragmentation and gas drag

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct N-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about ∼0.1 Earth mass within 10⁷ years. This brings another difficulty in formation of gas giant planets, which requires cores with ∼10 Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos’ collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the Solar System.     

Aggregate impacts in Saturn's rings

Ring particle aggregates are formed in the outer parts of Saturn's main rings. We study how collisions between aggregates can lead to destruction or coalescence of these aggregates, with local N-body simulations taking into account the dissipative impacts and gravitational forces between particles. Impacts of aggregates with different mass ratios are studied, as well as aggregates that consist of particles with different physical properties. We find that the outcome of the collision is very sensitive to the shape of the aggregate, in the sense that more elongated aggregates are more prone to be destroyed. We were interested in testing the accretion criterion Barbara and Esposito [Barbara, J.M., Esposito, L.W., 2002. Icarus 160, 161-171] used in their F ring simulations, according to which accretion requires that the masses of the colliding bodies differ at least by a factor of 100. We confirm that such a critical mass ratio exists. In particular, simulations indicate that the exact critical mass ratio depends on the internal density and elasticity of particles, and the distance from the planet. The zone of transition, defined by the distance where individual particles or small aggregates first start to stick on the larger aggregate, and by the distance where two similar sized aggregates on the average eventually coalesce is only about 5000 km wide, if fixed particle properties are used. The rotational state of the aggregates that form via aggregate collision rapidly reaches synchronous rotation, similarly to the aggregates that form via gradual growth.     

Effect of rotation on magnetogravitational instability of finite conducting gas in presence of suspended particles

The effect of rotation on the self-gravitational instability of an infinite homogeneous magnetised gas-particle medium in the presence of suspended particles is investigated. The conductivity of the medium is assumed to be finite. The equations of the problem are linearized and the general dispersion relation is obtained. The rotation is assumed along two different directions separately and separate dispersion relation for each case is obtained. The dispersion relation for propagation parallel and perpendicular to the uniform magnetic field along with rotation is derived. It is found that in presence of suspended particles, magnetic field, finite conductivity, rotation and viscosity, Jeans's criterion determines the condition of gravitational instability of gas-particle medium.