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.     

Ringlet–satellite interactions

We study the interaction of a satellite and a nearby ringlet on eccentric and inclined orbits. Secular torques originate from mean motion resonances and the secular interaction potential which represents the m  = 1 global modes of the ring. The torques act on the relative eccentricity and inclination. The resonances damp the relative eccentricity. The inclination instability owing to the resonances is turned off by a finite differential eccentricity of the order of 0.27 for nearly coplanar systems. The secular potential torque damps the eccentricity and inclination and does not affect the relative semi-major axis; also, it suppresses the inclination instability that persists at small differential eccentricities. The damping of the relative eccentricity and inclination forces an initially circular and planar small mass ringlet to reach the eccentricity and inclination of the satellite. When the planet is oblate, the interaction of the satellite damps the proper precession of a small mass ringlet so that it precesses at the satellite's rate independently of their relative distance. The oblateness of the primary modifies the long-term eccentricity and inclination magnitudes and introduces a constant shift in the apsidal and nodal lines of the ringlet with respect to those of the satellite. These results are applied to Saturn's F-ring, which orbits between the moons Prometheus and Pandora.     

The motion of a lunar satellite

Presented in this theory is a semianalytical solution for the problem of the motion of a satellite in orbit around the moon. The principal perturbations on such a body are due to the nonspherical gravity field of the moon, the attraction of the earth, and, to a lesser degree, the attraction of the sun. The major part of the problem is solved by means of the celebrated von Zeipel Method, first successfully applied to the motion of an artificial earth satellite by Brouwer in 1959. After eliminating from the Hamiltonian all terms with the period of the satellite and those with the period of the moon, it is suggested to solve the remaining problem with the aid of numerical integration of the modified equations of motion.This theory was written in 1964 and presented as a dissertation to Yale University in 1965. Since then a great deal has been learned about the gravity field of the moon. It seems that quite a number of recently determined gravity coefficients would qualify as small quantities of order two. Hence, according to the truncation criteria employed, they should be considered in the present theory. However, the author has not endeavored to update the work accordingly. The final results, therefore, are incomplete in the lunar gravitational perturbations. Nevertheless, the theory does give the largest such variations and it does present the methods by which perturbations may be derived for any gravity terms not actually developed.     

On the energy integral for first post-Newtonian approximation

The post-Newtonian approximation for general relativity is widely adopted by the geodesy and astronomy communities. It has been successfully exploited for the inclusion of relativistic effects in practically all geodetic applications and techniques such as satellite/lunar laser ranging and very long baseline interferometry. Presently, the levels of accuracy required in geodetic techniques require that reference frames, planetary and satellite orbits and signal propagation be treated within the post-Newtonian regime. For arbitrary scalar W and vector gravitational potentials \(W^j (j=1,2,3)\), we present a novel derivation of the energy associated with a test particle in the post-Newtonian regime. The integral so obtained appears not to have been given previously in the literature and is deduced through algebraic manipulation on seeking a Jacobi-like integral associated with the standard post-Newtonian equations of motion. The new integral is independently verified through a variational formulation using the post-Newtonian metric components and is subsequently verified by numerical integration of the post-Newtonian equations of motion.     

Gravitational interaction between a planet and an optically thin disc

We investigate the gravitational interaction between a planet and an optically thin protoplanetary disc, performing local three-dimensional hydrodynamical simulations. In the present study, we take account of radiative energy transfer in optically thin discs. Before the stage of planetary accretion, dust opacity is expected to decrease significantly because of grain growth and planetesimal formation. Thus, it would be reasonable to consider optically thin discs in the disc–planet interaction. Furthermore, we focus on small planets that can neither capture disc gas nor open a disc gap. The one-sided torque exerted on a planet by an optically thin disc is examined for various values of the disc optical thickness (<1). In optically thin discs, the temperature behind the density waves is lower than the unperturbed value because of radiative cooling. Heating due to shock dissipation is less effective than radiative cooling. Because of radiative cooling, the density distribution around the planet is not axisymmetric, which exerts an additional torque on the planet. The torque enhancement becomes maximum when the cooling time is comparable with the Keplerian period. The enhancement is significant for low-mass planets. For planets with  3 M_⊕  , the additional one-sided torque can be 40 per cent of the torque in the isothermal case. The radiative cooling is expected to change the differential torque and the migration speed of planets, too.     

Capture of grains into resonances through Poynting-Robertson drag

We review here some relevant problems connected to the evolution of circumstellar dust grains, subjected to Poynting-Robertson (PR) drag, and perturbed by first-order resonances with a planet on a circular orbit. We show that only outer mean motion resonances are able to counteract the damping effect of PR drag. However, the high orbital eccentricities reached by the particle lead to orbit crossings with the planet. This is a serious difficulty for a permanent trapping to be achieved. In any case, we show that the time spent in the resonance is long enough for statistical effects (accumulation at the resonant radius) to be significant. We underline some difficulties associated with this problem, namely, the non-adiabaticity of motion in the resonance phase space and the existence of close encounters with the planet at high eccentricities.     

Dynamics of a cloud of fast electrons travelling through the plasma

Numerical analysis of quasi-linear relaxation has been made for four models of electron beam with a finite length travelling through the plasma. In Model 4, a model atmosphere of the corona is adopted and also an increase in the cross-section of the electron beam is taken into account. The electron velocity distribution generally becomes a quasi-plateau form in limited velocity and time ranges. If, however, collisional decay of the fast electrons is too strong and the initial beam density is not high enough, the plateau does not appear. Collisional damping of plasma waves cannot be neglected, since the growth rate of the waves is strongly suppressed by the appearance of the quasi-plateau.An approximate formula for the velocity distribution of the solar electrons passing through the corona has been derived analytically taking into account not only the interaction with plasma waves, but also the collisional damping of the plasma waves and collisions with thermal particles. By the use of this formula, we can easily compute the time profile of the plasma waves caused by these solar electrons at any given place in the interplanetary space. The validity of this semi-analytical approach is checked by the numerical analysis of Model 4, showing a satisfactory fit between the numerical and semi-analytical results.The direct application of this method to the problems of type III radio bursts is left to a later paper.     

On the secular decrease in the semimajor axis of Lageos's orbit

The semimajor axis of the Lageos satellite's orbit is decreasing secularly at the rate of 1.1 mm day^–1. Ten possible mechanisms are investigated to discover which one (s), if any, might be causing the orbit to decay. Six of the mechanisms, resonance with the Earth's gravitational field, gravitational radiation, the Poynting-Robertson effect, transfer of spin angular momentum to the orbital angular momentum, drag from near-earth dust, and atmospheric drag by neutral hydrogen are ruled out because they are too small or require unacceptable assumptions to account for the observed rate of decay. Three other mechanisms, the Yarkovsky effect, the Schach effect, and terrestrial radiation pressure give perturbations whose characteristic signatures do not agree with the observed secular decrease (terrestrial radiation pressure appears to be too small in any case); hence they are also ruled out. Charged particle drag with the ions at Lageos's altitude is probably the principal cause of the orbital decay. An estimate of charged particle drag based upon laboratory experiments and satellite measurements of ion number densities accounts for 60 percent of the observed rate of decrease in the semimajor axis, assuming a satellite potential of –1V. This figure is in good agreement with other estimates based on charge drag theory. A satellite potential of –1.5V will explain the entire decay rate. Atmospheric drag from neutral hydrogen appears to be the next largest effect, explaining about 10 percent of the observed orbital decay rate.     

Co-Orbital Motion with Slowly Varying Parameters

We consider the dynamics of a test particle co-orbital with a satellite of mass m _s which revolves around a planet of mass M ₀ m _s with a mean motion n _s and semi-major axis a _s. We study the long term evolution of the particle motion under slow variations of (1) the mass of the primary, M ₀, (2) the mass of the satellite, m _s and (3) the specific angular momentum of the satellite J _s. The particle is not restricted to small harmonic oscillations near L ₄ or L ₅, and may have any libration amplitude on tadpole or horseshoe orbits. In a first step, no torque is applied to the particle, so that its motion is described by a Hamiltonian with slowly varying parameters. We show that the torque applied to the satellite, as measured by _s = j_s/(n _s J _s) induces an distortion of the phase space which is entirely described by an asymmetry coefficient = _s/, where = m _s/M. The adiabatic invariance of action implies furthermore that the long term evolution of the particle co-orbital motion depends only on the variation of m _s a _s with time. Applying a constant torque to the particle, as measured by _s = j_s/(n _s J _p) is then merely equivalent to replacing = _s/ by = (_s – _p)/. However, if the torque acting on the particle exhibits a radial gradient, then the action is no more conserved and the evolution of the particle orbit is no more controlled by m _s a _s only. We show that even mild torque gradients can dominate the orbital evolution of the particle, and eventually decide whether the latter will be pulled towards the stable equilibrium points L ₄ or L ₅, or driven away from them. Finally, we show that when the co-orbital bodies are two satellites with comparable masses m ₁ and m ₂, we can reduce the problem to that of a test particle co-orbital with a satellite of mass m ₁ + m ₂. This new problem has then parameters varying at rates which are combinations, with appropriate coefficients, of the changes suffered by each satellite.     

Motion of a satellite in the equatorial plane of a spheroid

An exact, closed-form solution of the problem of the motion of a satellite in the equatorial plane of an oblate body is obtained. It is shown that the classic formula for the motion of the perihelion is a first order approximation to the exact formula.     

Secular Motion in a 2nd Degree and Order-Gravity Field with no Rotation

The motion of a particle about a non-rotating 2nd degree and order-gravity field is investigated. Averaging conditions are applied to the particle motion and a qualitative analysis which reveals the general character of motion in this system is given. It is shown that the orbit plane will either be stationary or precess about the body's axis of minimum or maximum moment of inertia. It is also shown that the secular equations for this system can be integrated in terms of trigonometric, hyperbolic or elliptic functions. The explicit solutions are derived in all cases of interest.This revised version was published online in October 2005 with corrections to the Cover Date.     

Analytical description of charged particle motion in a reconnecting current sheet

We have obtained an analytical solution to the equation of motion in the guiding center approximation for nonrelativistic charged particles in a reconnecting current sheet with a three-component magnetic field. Given the electric field attributable to magnetic reconnection, the solution describes stable and unstable three-dimensional particle orbits. We have found the domain of input parameters at which the motion is stable. A physical interpretation of the processes affecting the stability of the motion is given. Charge separation is shown to take place in the sheet during the motion: oppositely charged particles are localized mostly in different regions of the current sheet. A formula is derived for the particle energy in stable and unstable orbits. The results obtained by numerical and analytical methods are compared.     

Radiation Pressure and the Poynting-Robertson Effect for Fluffy Dust Particles

A correct understanding of the dynamical effect of solar radiation exerted on fluffy dust particles can be achieved with assistance of a light scattering theory as well as the equation of motion. We reformulate the equation of motion so that the radiation pressure and the Poynting-Robertson effect on fluffy grains are given in both radial and nonradial directions from the center of the Sun. This allows numerical estimates of these radiation forces on fluffy dust aggregates in the framework of the discrete dipole approximation, in which the first term of the scattering coefficients in Mie theory determines the polarizability of homogeneous spheres forming the aggregates.The nonsphericity in shape turns out to play a key role in the dynamical evolution of dust particles, while its consequence depends on the rotation rate and axis of the grains. Unless a fluffy dust particle rapidly revolves on its randomly oriented axis, the nonradial radiation forces may prevent, apart from the orbital eccentricity and semimajor axis, the orbital inclination of the particle from being preserved in orbit around the Sun. However, a change in the inclination is most probably controlled by the Lorentz force as a consequence of the interaction between electric charges on the grains and the solar magnetic field. Although rapidly and randomly rotating grains spiral into the Sun under the Poynting-Robertson effect in spite of their shapes and structures, fluffy grains drift inward on time scales longer at submicrometer sizes and shorter at much larger sizes than spherical grains of the same sizes. Numerical calculations reveal that the dynamical lifetimes of fluffy particles are determined by the material composition of the grains rather than by their morphological structures and sizes. The Poynting-Robertson effect alone is nevertheless insufficient for giving a satisfactory estimate of lifetimes for fluffy dust grains since their large ratios of cross section to mass would reduce the lifetimes by enhancing the collisional probabilities. We also show that the radiation pressure on a dust particle varies with the orbital velocity of the particle but that this effect is negligibly small for dust grains in the Solar System.     

Formation of the regular satellites of giant planets in an extended gaseous nebula II: satellite migration and survival

For a satellite to survive in the disk the time scale of satellite migration must be longer than the time scale for gas dissipation. For large satellites (∼1000 km) migration is dominated by the gas tidal torque. We consider the possibility that the redistribution of gas in the disk due to the tidal torque of a satellite with mass larger than the inviscid critical mass causes the satellite to stall and open a gap (W.R. Ward, 1997, Icarus 26, 261-281). We adapt the inviscid critical mass criterion to include gas drag, and m-dependent nonlocal deposition of angular momentum. We find that such a model holds promise of explaining the survival of satellites in the subnebula, the mass versus distance relationship apparent in the saturnian and uranian satellite systems, the concentration of mass in Titan, and the observation that the satellites of Jupiter get rockier closer to the planet whereas those of Saturn become increasingly icy. It is also possible that either weak turbulence (close to the planet) or gap-opening satellite tidal torque removes gas on a similar time scale (10⁴-10⁵ years) as the orbital decay time of midsized (200-700 km) regular satellites forming in the inner disk (inside the centrifugal radius (I. Mosqueira and P.R. Estrada, 2003, Icarus, this issue)). We argue that Saturn’s satellite system bridges the gap between those of Jupiter and Uranus by combining the formation of a Galilean-sized satellite in a gas optically thick subnebula with a strong temperature gradient, and the formation of smaller satellites, closer to the planet, in a disk with gas optical depth ?1, and a weak temperature gradient.Using an optically thick inner disk (given gaseous opacity), and an extended, quiescent, optically thin outer disk, we show that there are regions of the disk of small net tidal torque (even zero) where satellites (Iapetus-sized or larger) may stall far from the planet. For our model these outer regions of small net tidal torque correspond roughly to the locations of Callisto and Iapetus. Though the precise location depends on the (unknown) size of the transition region between the inner and outer disks, the result that Saturn’s is found much farther out (at ∼3rcS, where rcS is Saturn’s centrifugal radius) than Jupiter’s (at ∼ 2rcJ, where rcJ is Jupiter’s centrifugal radius) is mostly due to Saturn’s less massive outer disk and larger Hill radius. However, despite the large separation between Ganymede and Callisto and Titan and Iapetus, the long formation and migration time scales for Callisto and Iapetus (I. Mosqueira and P.R. Estrada, 2003, Icarus, this issue) makes it possible (depending on the details of the damping of acoustic waves) that the tidal torque of Ganymede and Titan clears the gas disk out to their location, thus stranding Callisto and Iapetus far from the planet. Either way, our model provides an explanation for the presence of regular satellites outside the centrifugal radii of Jupiter and Saturn, and the absence of such a satellite for Uranus.     

A problem in resonance

Saturn's satellite Hyperion experiences large perturbations by Titan, the largest of Saturn's satellites, because of the closeness of Hyperion's orbital period to three-quarters of that of Titan.The motion of Hyperion is a superposition of periodic fluctuations (both free and forced) onto a motion which is periodic in a suitable uniformly rotating frame of reference, and in which Hyperion would be at a maximum distance from Saturn at each conjunction with Titan.Successive attempts to determine the mass of Titan from observed perturbations of Hyperion have suffered from omissions in the theory of terms subsequently found to be significant.An attempt is in progress which, it is believed, comprises all long-period changes in the osculating elements which are of second degree in the mass of Titan, and of third degree in the eccentricity of its orbit. Results so far obtained indicate that the period of the free motion of the orbit plane of Hyperion is better determined by Woltjer's reduction of the observations than by his theoretical calculations.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 27–September 2, 1972.     

Accretion of terrestrial planets from oligarchs in a turbulent disk

We have investigated the final accretion stage of terrestrial planets from Mars-mass protoplanets that formed through oligarchic growth in a disk comparable to the minimum mass solar nebula (MMSN), through N-body simulation including random torques exerted by disk turbulence due to Magneto-Rotational Instability. For the torques, we used the semi-analytical formula developed by Laughlin et al. [Laughlin, G., Steinacker, A., Adams, F.C., 2004. Astrophys. J. 608, 489-496]. The damping of orbital eccentricities (in all runs) and type-I migration (in some runs) due to the tidal interactions with disk gas is also included. Without any effect of disk gas, Earth-mass planets are formed in terrestrial planet regions in a disk comparable to MMSN but with too large orbital eccentricities to be consistent with the present eccentricities of Earth and Venus in our Solar System. With the eccentricity damping caused by the tidal interaction with a remnant gas disk, Earth-mass planets with eccentricities consistent with those of Earth and Venus are formed in a limited range of disk gas surface density (∼10⁻⁴ times MMSN). However, in this case, on average, too many (?6) planets remain in terrestrial planet regions, because the damping leads to isolation between the planets. We have carried out a series of N-body simulations including the random torques with different disk surface density and strength of turbulence. We found that the orbital eccentricities pumped up by the turbulent torques and associated random walks in semimajor axes tend to delay isolation of planets, resulting in more coagulation of planets. The eccentricities are still damped after planets become isolated. As a result, the number of final planets decreases with increase in strength of the turbulence, while Earth-mass planets with small eccentricities are still formed. In the case of relatively strong turbulence, the number of final planets are 4-5 at 0.5-2 AU, which is more consistent with Solar System, for relatively wide range of disk gas surface density (∼¹⁰⁻⁴-¹⁰⁻² times MMSN).     

A Critical review of theoretical models of negatively polarized light scattered by atmosphereless solar system bodies

About a dozen physical mechanisms and models aspire to explain the negative polarization of light scattered by atmosphereless celestial bodies. This is too large a number for the reliable interpretation of observational data. Through a comparative analysis of the models, our main goal is to answer the question: Does any one model have an advantage over the others? Our analysis is based on new laboratory polarimetric and photometric data as well as on theoretical results. We show that the widely used models due to Hopfield and Wolff cannot realistically explain the phase-angle dependence of the degree of polarization observed at small phase angles. The so-called interference or coherent backscattering mechanism is the most promising model. Models based on that mechanism use well-defined physical parameters to explain both negative polarization and the opposition effect. They are supported by laboratory experiments, particularly those showing enhancement of negative polarization with decreasing particle size down to the wavelength of light. According to the interference mechanism, pronounced negative branches of polarization, like those of C-class asteroids, may indicate a high degree of optical inhomogeneity of light-scattering surfaces at small scales. The mechanism also seems appropriate for treating the negative polarization and opposition effects of cometary dust comae, planetary rings, and the zodiacal light.     

Gravitational perturbations of equatorial orbits

Asymptotic solutions are developed for the motion of a geocentric satellite in the equatorial plane due to gravitational perturbations such as nonsphericity (especially oblateness) of the primary body. Axisymmetric potentials are considered. A class of transformations is developed and the equations of motion are solved by the method of generalized multiple scales. Further it is shown that the equations of motion can be transformed into the required form to within any specified degree of accuracy. The transformations form an Abelian group of infinite order which leaves the differential equations of motion invariant. Solutions are developed in terms of elementary functions instead of elliptic or other higher transcendental functions and are shown to agree with known results.This investigation was carried out under NASA Grant NGR-31-001-152 with the author as a consultant to Princeton University.