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.     

Rotation rate and velocity dispersion of planetary ring particles with size distribution II. Numerical simulation for gravitating particles

We examine rotation rates of gravitating particles in low optical depth rings, on the basis of the evolution equation of particle rotational energy derived by Ohtsuki [Ohtsuki, K., 2006. Rotation rate and velocity dispersion of planetary ring particles with size distribution. I. Formulation and analytic calculation. Icarus 183, 373-383]. We obtain the rates of evolution of particle rotation rate and velocity dispersion, using three-body orbital integration that takes into account distribution of random velocities and rotation rates. The obtained stirring and friction rates are used to calculate the evolution of velocity dispersion and rotation rate for particles in one- and two-size component rings as well as those with a narrow size distribution, and agreement with N-body simulation is confirmed. Then, we perform calculations to examine equilibrium rotation rates and velocity dispersion of gravitating ring particles with a broad size distribution, from 1 cm up to 10 m. We find that small particles spin rapidly with 〈ω²〉^1/2/Ω?¹⁰²-¹⁰³, where ω and Ω are the particle rotation rate and its orbital angular frequency, respectively, while the largest particles spin slowly, with 〈ω²〉^1/2/Ω?1. The vertical scale height of rapidly rotating small particles is much larger than that of slowly rotating large particles. Thus, rotational states of ring particles have vertical heterogeneity, which should be taken into account in modeling thermal infrared emission from Saturn's rings.     

Collision mechanics and the structure of planetary ring edges

A numerical simulation of collisional evolution, originally developed to model planetary accretion processes, is applied to a hypothetical ring with parameters modeled after Saturn's rings in order to study changes in radial structure near ring edges. The tendency of rings to spread so as to conserve angular momentum while energy is dissipated in collisions is confirmed if random motion is in equilibrium. Even with no energy loss (coefficient of restitution in velocity ε = 1), spreading occurs becase random motion is increasing. With a moderately side-scattering collisional model, characteristic of collisions of nonrotating spheres (the slippery “billiard-ball” model), random motion increases for ε > 0.63, in agreement with analytical models. For isotropic scattering, which may be more realistic given particle rotation, damping dominates for ε up to 0.83. As long as random motion is damped, ring edges may contract rather than spread, producing concentrations of material just inside the ring edges reminiscent of results of earlier stimulation which did not precisely conserve angular momentum.     

Simulations of dense planetary rings: IV. Spinning self-gravitating particles with size distribution

Previous self-gravitating simulations of dense planetary rings are extended to include particle spins. Both identical particles as well as systems with a modest range of particle sizes are examined. For a ring of identical particles, we find that mutual impact velocity is always close to the escape velocity of the particles, even if the total rms velocity dispersion of the system is much larger, due to collective motions associated to wakes induced by near-gravitational instability or by viscous overstability. As a result, the spin velocity (i.e., the product of the particle radius and the spin frequency) maintained by mutual impacts is also of the order of the escape velocity, provided that friction is significant. For the size distribution case, smaller particles have larger impact velocities and thus larger spin velocities, particularly in optically thick rings, since small particles move rather freely between wakes. Nevertheless, the maximum ratio of spin velocities between the smallest and largest particles, as well as the ratio for translational velocities, stays below about 5 regardless of the width of the size distribution. Particle spin state is one of the important factors affecting the temperature difference between the lit and unlit face of Saturn's rings. Our results suggest that, to good accuracy, the spin frequency is inversely proportional to the particle size. Therefore, the mixing ratio of fast rotators to slow rotators on the scale of the thermal relaxation time increases with the width of the particle size distribution. This will offer means to constrain the particle size distribution with the systematic thermal infrared observations carried by the Cassini probe.     

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.     

A new formulation of the viscosity in planetary rings

We present a new formulation of the viscosity in planetary rings, where particles interact through their gravitational forces and direct collisions. In the previous studies on the viscosity in self-gravitating rings, the viscosity consists of three components, which are defined separately in different ways. The complex definitions make it difficult to evaluate the viscosity in N-body simulation of rings. In our new formulation, the viscosity is expressed in terms of changes in orbital elements of particles due to particle interactions. This makes the expression of the viscosity simple. The new formulation gives a simple way to evaluate the viscosity in N-body simulation. We find that for practical evaluation of the viscosity of planetary rings, only energy dissipation at direct inelastic collisions is needed.For tenuous particle disks (i.e., optically thin disks), we further derive a formula of the viscosity. The formula requires only a numerical coefficient that can be obtained from three-body calculation. Since planetesimal disks are also tenuous, the viscosity in planetesimal disks can be also obtained from this formula. In a subsequent paper, we will evaluate this coefficient through three-body calculation and obtain the viscosity for a wide range of parameters such as the restitution coefficient and the radial location in rings.     

On the rotation of a moonlet embedded in planetary rings

We examine the rotation of a small moonlet embedded in planetary rings caused by impacts of ring particles, using analytic calculation and numerical orbital integration for the three-body problem. Taking into account the Rayleigh distribution of particles' orbital eccentricities and inclinations, we evaluate both systematic and random components of rotation, where the former arises from an average of a large number of small impacts and the latter is contribution from large impacts. Calculations for parameter values corresponding to inner parts of Saturn's rings show that a moonlet would spin slowly in the prograde direction if most impactors are small particles whose velocity dispersion is comparable to or smaller than the moonlet's escape velocity. However, we also find that the effect of the random component can be significant, if the velocity dispersion of particles is larger and/or impacts of large particles comparable to the moonlet's size are common: in this case, both prograde and retrograde rotations can be expected. In the case of a small moonlet embedded in planetary rings of equal-sized particles, we find that the systematic component dominates the moonlet rotation when m/M?0.1 (m and M are the mass of a particle and a moonlet, respectively), while the random component is dominant when m/M?0.3. We derive the condition for the random component to dominate moonlet rotation on the basis of our results of three-body orbital integration, and confirm agreement with N-body simulation.     

Saturn's rings: Particle composition and size distribution as constrained by observations at microwave wavelengths: II. Radio interferometric observations

The sizes, composition, and number of particles comprising the rings of Saturn may be meaningfully constrained by a combination of radar- and radio-astronomical observations. In a previous paper, we have discussed constraints obtained from radar observations. In this paper, we discuss the constraints imposed by complementary “passive” radio observations at similar wavelengths. First, we present theoretical models of the brightness of Saturn's rings at microwave wavelengths (0.34–21.0 cm), including both intrinsic ring emission and diffuse scattering by the rings of the planetary emission. The models are accurate simulations of the behavior of realistic ring particles and are parameterized only by particle composition and size distribution, and ring optical depth. Second, we have reanalyzed several previously existing sets of interferometric observations of the Saturn system at 0.83-, 3.71-, 6.0-, 11.1-, and 21.0-cm wavelengths. These observations all have spatial resolution sufficient to resolve the rings and planetary disk, and most have resolution sufficient to resolve the ring-occulted region of the disk as well. Using our ring models and a realistic model of the planetary brightness distribution, we are able to establish improved constraints on the properties of the rings. In particular, we find that: (a) the maximum optical depth in the rings is ～ 1.5 ± 0.3 referred to visible wavelengths; (b) a significant decrease in ring optical depth from λ3.7 to λ21.0 cm allows us to rule out the possibility that more than ～30% of the cross section of the rings is composed of particles larger than a meter or so; this assertion is essentially independent of uncertainties in particle adsorption coefficient; and (c) the ring particles cannot be primarily of silicate composition, independently of particle size, and the particles cannot be primarily smaller than ～0.1 cm, independently of composition.     

HST observations of azimuthal asymmetry in Saturn's rings

From 378 Hubble Space Telescope WFPC2 images obtained between 1996-2004, we have measured the detailed nature of azimuthal brightness variations in Saturn's rings. The extensive geometric coverage, high spatial resolution (), and photometric precision of the UBVRI images have enabled us to determine the dependence of the asymmetry amplitude and longitude of minimum brightness on orbital radius, ring elevation, wavelength, solar phase angle, and solar longitude. We explore a suite of dynamical models of self-gravity wakes for two particle size distributions: a single size and a power law distribution spanning a decade in particle radius. From these N-body simulations, we calculate the resultant wake-driven brightness asymmetry for any given illumination and viewing geometry. The models reproduce many of the observed properties of the asymmetry, including the shape and location of the brightness minimum and the trends with ring elevation and solar longitude. They also account for the “tilt effect” in the A and B rings: the change in mean ring brightness with effective ring opening angle, |Beff|. The predicted asymmetry depends sensitively on dynamical ring particle properties such as the coefficient of restitution and internal mass density, and relatively weakly on photometric parameters such as albedo and scattering phase function. The asymmetry is strongest in the A ring, reaching a maximum amplitude A∼25% near a=128,000 km. Here, the observations are well-matched by an internal particle density near 450 kg m⁻³ and a narrow particle size distribution. The B ring shows significant asymmetry (∼5%) in regions of relatively low optical depth (τ∼0.7). In the middle and outer B ring, where τ?1, the asymmetry is much weaker (∼1%), and in the C ring, A<0.5%. The asymmetry diminishes near opposition and at shorter wavelengths, where the albedo of the ring particles is lower and multiple-scattering effects are diminished. The asymmetry amplitude varies strongly with ring elevation angle, reaching a peak near |Beff|=10° in the A ring and at |Beff|=15-20° in the B ring. These trends provide an estimate of the thickness of the self-gravity wakes responsible for the asymmetry. Local radial variations in the amplitude of the asymmetry within both the A and B rings are probably caused by regional differences in the particle size distribution.     

The Christiansen Effect in Saturn’s narrow dusty rings and the spectral identification of clumps in the F ring

Stellar occultations by Saturn’s rings observed with the Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft reveal that dusty features such as the F ring and the ringlets in the Encke and the Laplace Gaps have distinctive infrared transmission spectra. These spectra show a narrow optical depth minimum at wavelengths around 2.87 μm. This minimum is likely due to the Christiansen Effect, a reduction in the extinction of small particles when their (complex) refractive index is close to that of the surrounding medium. Simple Mie-scattering models demonstrate that the strength of this opacity dip is sensitive to the size distribution of particles between 1 and 100 μm across. Furthermore, the spatial resolution of the occultation data is sufficient to reveal variations in the transmission spectra within and among these rings. In both the Encke Gap ringlets and F ring, the opacity dip weakens with increasing local optical depth, which is consistent with the larger particles being concentrated near the cores of these rings. The Encke Gap ringlets also show systematically weaker opacity dips than the F ring and Laplace Gap ringlet, implying that the former has a smaller fraction of grains less than ∼30 μm across. However, the strength of the opacity dip varies most dramatically within the F ring; certain compact regions of enhanced optical depth lack an opacity dip and therefore appear to have a greatly reduced fraction of grains in the few-micron size range. Such spectrally-identifiable structures probably represent a subset of the compact optically-thick clumps observed by other Cassini instruments. These variations in the ring’s particle size distribution can provide new insights into the processes of grain aggregation, disruption and transport within dusty rings. For example, the unusual spectral properties of the F-ring clumps could perhaps be ascribed to small grains adhering onto the surface of larger particles in regions of anomalously low velocity dispersion.     

A numerical model of cohesion in planetary rings

We present a numerical method that incorporates particle sticking in simulations using the N-body code pkdgrav to study motions in a local rotating frame, such as a patch of a planetary ring. Particles stick to form non-deformable but breakable aggregates that obey the (Eulerian) equations of rigid-body motion. Applications include local simulations of planetary ring dynamics and planet formation, which typically feature hundreds of thousands or more colliding bodies. Bonding and breaking thresholds are tunable parameters that can approximately mimic, for example, van der Waals forces or interlocking of surface frost layers. The bonding and breaking model does not incorporate a rigorous treatment of internal fracture; rather the method serves as motivation for first-order investigation of how semi-rigid bonding affects the evolution of particle assemblies in high-density environments.We apply the method to Saturn’s A ring, for which laboratory experiments suggest that interpenetration of thin, frost-coated surface layers may lead to weak cohesive bonding. These experiments show that frost-coated icy bodies can bond at the low impact speeds characteristic of the rings. Our investigation is further motivated by recent simulations that suggest a very low coefficient of restitution is needed to explain the amplitude of the azimuthal brightness asymmetry in Saturn’s A ring, and the hypothesis that fine structure in Saturn’s B ring may in part be caused by large-scale cohesion.This work presents the full implementation of our model in pkdgrav, as well as results from initial tests with a limited set of parameters explored. We find a combination of parameters that yields aggregate size distribution and maximum radius values in agreement with Voyager data for ring particles in Saturn’s outer A ring. We also find that the bonding and breaking parameters define two strength regimes in which fragmentation is dominated either by collisions or other stresses, such as tides. We conclude our study with a discussion of future applications of and refinements to our model.     

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.     

Features around embedded moonlets in Saturn's rings: The role of self-gravity and particle size distributions

This paper presents the results of N-body simulations of moonlets embedded in broad rings, focusing specifically on the saturnian A ring. This work adds to previous efforts by including particle self-gravity and particle size distributions. The discussion here focuses primarily on the features that form in the background particles as a result of the moonlet. Particle self-gravity tends to damp out features produced by embedded moonlets and this damping is enhanced if the moonlet is simply the largest member of a continuous size distribution. Observable features around an embedded moonlet appear to require that the largest ring particles be no more massive than 1/30 the mass of the moonlet. These results, compared with current and future Cassini observations, will provide insight into the nature of the particle population in the saturnian rings. Some time is also spent analyzing the way in which the background particles cluster around the moonlet. The accretion of small particles onto the moonlet can be limited by disruptive collisions with the largest ring particles in the particle size distribution.     

A numerical algorithm for dissipative Keplerian particle discs

We present a numerical algorithm designed to study the evolution of a distribution of solid particles orbiting around the Sun that could be applied to similar systems in which f _c ≤ Ω/2π, where f _c is the frequency of collisions, and Ω is the orbital angular speed. A number of sample particles are used to represent the whole system. Binary collisions are treated in a novel way using a Monte Carlo method that works as follows. Orbits are locally sampled to compute the velocity dispersion. Then the velocity vectors of the sample particles are modified according to random collisions with virtual particles which have velocities taken from a normal distribution computed using the previously found local velocity dispersion. The energy and momentum taken up by the virtual particles are redistributed among the neighbours of the sample particle undergoing the collision, so that conservation laws are satisfied. Simulations using this model give an estimation of the final distribution of inclinations and the associated evolutionary time-scale.     

The vertical structure and thickness of Saturn's rings

We have considered the steady state vertical structure of Saturn's rings with regard to whether collapse to a monolayer due to collisions between particles, the end state predicted by Jeffreys (1947a), may be prevented by any of a variety of mechanisms. Given a broad distribution of particle sizes such as a typical power law n(R) = n₀R^?3, it is found that gravitational scattering of small particles by large particles maintains a true ring thickness of several times the radius of the largest particles, or many times the radius of the smallest particles. Thus the “many-particle-thick” condition which best satisfies optical observations, such as the opposition effect, may be reconciled with ongoing particle collisions. If we consider the obvious sources of energy available for such a process, we find that a ring thickness of only tens of meters may be sustained over the lifetime of the solar system. This implies a maximum particle size on the order of a few meters.     

Size distribution of particles in planetary rings

Harris (Icarus24, 190–192) has suggested that the maximum size of particles in a planetary ring is controlled by collisional fragmentation rather than by tidal stress. While this conclusion is probably true, estimated radius limits must be revised upward from Harris' values of a few kilometers by at least an order of magnitude. Accretion of particles within Roche's limit is also possible. These considerations affect theories concerning the evolution of Saturn's rings, of the Moon, and of possible former satellites of Mercury and Venus. In the case of Saturn's rings, comparison of various theoretical scenarios with available observational evidence suggests that the rings formed from the breakup of larger particles rather than from original condensation as small particles. This process implies a distribution of particle sizes in Saturn's rings possibly ranging up to ～100 km but with most cross-section in cm-scale particles.     

Thermal emission from a multiple scattering model of Saturn's rings

Models of Saturn's B ring have been investigated which include the shadowing mechanism, realistic phase functions for the ring particles, and the effects of multiple scattering and a particle size dispersion. These models are based on the assumption that the rings form a layer many particles thick. A power law relation dn ∝ ?^?s is used for the size dispersion law of the ring particles, where dn is the number of particles with radii between ? and ? + d?. In the calculation of the infrared brightness temperature of the rings, the effect of mutual heating among the ring particles is considered quantitatively for the first time. The parameters of the polydisperse s = 2 model can be chosen to satisfy both optical (λ ? 1.1 μ) and infrared data, but the situation could be much clarified if a good phase curve for the rings were available in the red, if the ring brightness were known accurately for λ > 1 μ, and if it could be established whether the ring particles are rotating synchronously.     

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.     

Saturn's rings: Particle size distributions for thin layer models

The small physical thickness of Saturn's rings requires that radio occultation observations be interpreted using scattering models with limited amounts of multiple scatter. A new model in which the possible order of near-forward scatter is strictly limited allows for the small physical thickness, and can be used to relate Voyager 1 observations of 3.6-and 13-cm wavelength microwave scatter from Saturn's rings to the ring particle size distribution function n(a), for particles with radius 0.001 ≤ a ≤ 20 m. This limited-scatter model yields solutions for particle size distribution functions for eight regions in Saturn's rings, which exhibit approximately inverse-cubic power-law behavior, with large-size cutoffs in particle radius ranging from about 5 m in ring C to about 10 m in parts of ring A. The power-law index is about 3.1 in ring C, about 2.8 in the Cassini division, and increases systematically with radial location in ring A from 2.7 at 2.10R_s to slightly more than 3.0 at 2.24R_s. Corresponding mass densities are 32–43 kg/m² in ring C, 188 kg/m² in the Cassini division, and 244–344 kg/m² in ring A, under the assumption that the material density of the particles is 0.9 g/cm³. These values are a factor of 1 to 2 lower than first-order mass loading estimates derived from resonance phenomena. In view of the uncertainties in the measurements and in the linear density wave model, and the strong arguments for icy particles with specific gravity not greater than about 1, we interpret this discrepancy as being indicative of possible differences in the regions studied, or systematic errors in the interpretation of the scattering results, the density wave phenomena, or some combination of the above.     

Saturn's rings in the thermal infrared

This paper reviews our current knowledge of Saturn's rings’ physical properties as derived from thermal infrared observations. Ring particle composition, surface structure and spin as well as the vertical structure of the main rings can be determined. These properties are the key to understand the origin and evolution of Saturn's rings. Ring composition is mainly constrained by observations in the near-infrared but the signature of some probable contaminants present in water ice may also be found at mid-infrared wavelengths. The absence of the silicate signature limits nowadays their mass fraction to 10^−7±1. Recent measurements on the thermal inertia of the ring particle surface show it is very low, of the order of 5±2 Jm⁻² K⁻¹ s^−1/2. New models and observations of the complete crossing of the planetary shadow are needed to attribute this low value either to compact regoliths covered by cracks due to collisions and thermal stresses or to large fluffy and irregular surfaces. Studies of the energy balance of ring particles show a preference for slowly spinning particles in the main rings. Supplementary observations at different phase angles, showing the temperature contrast between night and day sides of particles, and new models including finite spin and thermal inertia, are needed to constrain the actual spin distribution of ring particles. These results can then be compared to numerical simulations of ring dynamics. Many thermal models have been proposed to reproduce observations of the main rings, including alternative mono- or many-particles-thick layers or vertical heterogeneity, with no definitive answer. Observations on the lit and dark faces of rings as a function of longitude, at many incidence and emission angles, would provide prime information on the vertical thermal gradient due to interparticle shadowing from which constraints on the local vertical structure and dynamics can be produced. Future missions such as Cassini will provide new information to further constrain the ring thermal models.