Issue 
A&A
Volume 576, April 2015



Article Number  A92  
Number of page(s)  5  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201425102  
Published online  06 April 2015 
Research Note
Stability of rings around a triaxial primary
^{1}
École Normale Supérieure de Cachan, 61 Avenue du Président
Wilson,
94230
Cachan,
France
email:
antoine.lehebel@enscachan.fr
^{2}
Cornell University, Center for Radiophysics and Space Research,
Ithaca, NY
14853,
USA
Received: 2 October 2014
Accepted: 11 February 2015
Aims. Generally, the oblateness of a planet or moon is what causes rings to settle into its equatorial plane. However, the recent suggestion that a ring system might exist (or have existed) about Rhea, a moon whose shape includes a strong prolate component pointed toward Saturn, raises the question of whether rings around a triaxial primary can be stable. We study the role of prolateness in the behavior of rings around Rhea and extend our results to similar problems such as possible rings around exoplanets.
Methods. Using a Hamiltonian approach, we point out that the dynamical behavior of ring particles is governed by three different time scales: the orbital period of the particles, the rotation period of the primary, and the precession period of the particles’ orbital plane. In the case of Rhea, two of these are well separated from the third, allowing us to average the Hamiltonian twice. To study the case of slow rotation of the primary, we also carry out numerical simulations of a thin disk of particles undergoing secular effects and damping.
Results. In the case of Rhea, the averaging reduces the Hamiltonian to an oblate potential, under which rings would be stable only in the equatorial plane. This is not the case for Iapetus; rather, it is the lack of a prolate component to its shape that allows Iapetus to host rings. Plausible exoplanets should mostly be in the same regime as Rhea, though other outcomes are possible. The numerical simulations indicate that, even when the double averaging is irrelevant, rings settle in the equatorial plane on an approximately constant time scale.
Key words: celestial mechanics / minor planets, asteroids: general / planets and satellites: rings
© ESO, 2015
1. Introduction
The oblateness of a planet (i.e., the slight flattening along its spin axis), characterized by the coefficient J_{2}, defines the Laplace plane in which potential rings may be stable (Tremaine et al. 2009). However, some bodies, like moons orbiting their planet, are also elongated by tidal effects in the direction toward the object they are orbiting. This elongation is described by the coefficient C_{22}. Its effect on the stability of planetary rings is the subject of this work.
In 2008, the Cassini spacecraft observed some sharp symmetric drops in the flow of charged particles around Rhea (Jones et al. 2008). This observation was interpreted as being due to the presence of narrow rings orbiting the moon. Furthermore, Schenk & McKinnon (2009) found a narrow equatorial color signature that was ascribed to the fallout of ring material and was cited as further evidence of Rhean rings. However, Tiscareno et al. (2010) ruled out the current existence of rings with an indepth imaging study of the equatorial region, after which Schenk et al. (2011) suggested that the narrow equatorial color signature may still indicate the existence of rings in the past.
Rhea would have been the first triaxial object known to possess a ring system. As for most moons close to their planet, both Rhea’s J_{2} and C_{22} are significant. This work is motivated by the question of whether it is even possible for rings to be stable in such an environment.
The calculations we develop can be adapted to Iapetus, for which rings have been invoked to explain the origin of its equatorial mountain range (Ip 2006; Levison et al. 2011; Dombard et al. 2012). Iapetus is different from Rhea in that it is simply oblate, but a comparison with our model remains interesting.
Some exoplanets, such as hot Jupiters, are likely to be drawn out by tidal effects in the direction of their star. This deformation, added to the fast rotation of the planet, must have complex effects on the presence of putative rings.
We will carry out a Hamiltonian analysis, following the model developed by Dobrovolskis et al. (1989) for Neptune polar rings. SteimanCameron & Durisen (1984) used a similar model to describe gas rings surrounding triaxial galaxies, which is better adapted to the current case of slowly rotating bodies.
In the next section, we show that the inertial frame is better adapted to the description of the motion than the rotating frame fixed on Rhea. In Sect. 3, we derive the perturbative Hamiltonian and show that the prolateness has no average effect. In Sect. 4, we investigate the case of Iapetus and exoplanets. Section 5 contains a computational simulation that allows us to study slower rotations.
2. The choice of a frame
We consider a test particle moving around Rhea on a nearlyKeplerian orbit perturbed by the shape of Rhea and the presence of Saturn. In this work, the particle’s orbit is assumed to be circular. This is probably a good approximation for any putative rings.
2.1. Definition of the frame and time scales
We will first define a system of axes in which the expressions for the perturbations are simple (Fig. 1). In general, I is inclination, Ω is longitude of ascending node, and θ and φ are the particle’s colatitude and longitude, respectively. The x_{R} axis always points toward Saturn during the motion of Rhea; therefore, it corresponds to Rhea’s long axis and to the C_{22} term. The z_{R} axis is taken parallel to Rhea’s angular momentum, i.e., perpendicular to its orbital plane; it is the short axis. The y_{R} axis completes the other two to form a direct basis. The angles in this frame will be labeled with a subscript R: I_{R}, Ω_{R}, etc.
Fig. 1
Angles in the rotating frame. 
This frame is well adapted to express Rhea’s gravitational potential, but it will not allow a simple description of the behavior of putative rings. In order to better understand the importance of the frame, we have to focus on the three relevant time scales of our problem: τ, the period of one orbit of the testparticle around the moon; τ_{R}, the period of one orbit of the moon around the planet (identical to the rotation period of the moon); τ_{p}, the precession period of the test particle’s orbit due to the nonsphericity of the moon. This precession is responsible for the settlement of the rings in preferred planes, so we can take this time as a low estimation of the settling time.
If τ_{R} ≫ τ and τ_{R} ≫ τ_{p}, the particle has time to settle in the preferred planes relative to the moon, and the rotating frame coordinates is the most convenient. The rings turn with the moon during its motion.
If τ_{p} ≫ τ_{R}, the rings remain steady in the inertial frame of Saturn during one period of Rhea. Therefore, on the time scale on which the particles can settle into rings, they only see an average potential: the rotating C_{22} term averages to give an azimuthally symmetric shape to the satellite. In this case, the frame that rotates with Rhea is not the most convenient.
2.2. The case of Rhea
We estimate the different time scales as (1)where a is the semimajor axis and is the mean motion for the test particle on its orbit around Rhea, a_{R} and n_{R} are the same for Rhea on its orbit around Saturn, M_{R} and M_{S} are the respective masses of Rhea and Saturn, R_{R} is the radius and J_{2} is the oblateness parameter for Rhea, and G is Newton’s constant.
In τ_{p}, we estimated as from Eq. (6.250) in Murray & Dermott (1999). This is valid when the only perturbation is oblateness, but it should give the correct precession rate to the first order in our situation, as the J_{2} and C_{22} are almost equally significant.
In the case of Rhea, J_{2} = 889 × 10^{6}, M_{R} = 2.308 × 10^{21} kg, R_{R} = 764.3 km, and τ_{R} = 4.518d (Anderson & Schubert 2007). If we take a ≃ 1.5 R_{R}, as in Jones et al. (2008), the time scales are as given in Table 1.
Time scales of the problem for Rhea at a = 1.5 R_{R}.
Clearly, τ_{p} ≫ τ_{R}. Thus, the most convenient frame is the inertial frame centered on Rhea that does not rotate. The transformation between the inertial and the rotational frames is just a rotation of angle n_{R}t around the zaxis.
There is a factor of at least 50 between each time scale. This invites us to make a double averaging: the first on the orbit of the particle around Rhea, and the second on the orbit of Rhea around Saturn.
3. The perturbative Hamiltonian
3.1. Derivation of the Hamiltonian
The system is conservative, so the semimajor axis a is constant. Therefore, an orbit is described by its inclination I_{R} and the position of the ascending node Ω_{R}. Another set of coordinates is relevant here, the particle’s colatitude and longitude (θ_{R},φ_{R}), because the perturbative Hamiltonian is expressed in their terms. The timeevolution of I_{R} and Ω_{R} is governed by (Kaula 1968) (2)All Hamiltonians will be written per unit of the test particle’s mass. The most important perturbation for us is the one that results from the shape of Rhea (SteimanCameron & Durisen 1984): (3)The particle also feels the effect of the gravitational field from Saturn. The quadrupole term (i.e., order 2 in the development in a/a_{R}) is a good approximation for this perturbation. It is given by (Murray & Dermott 1999) (4)We will now take the average of the perturbation over the duration τ = 2π/n of one orbit of the particle around Rhea, , where H = H_{R} + H_{S}. Dropping the parts of the Hamiltonian that depend neither on I nor on Ω, we obtain (5)We employ a similar method to average the Hamiltonian once more, over a time τ_{R}, which corresponds to one orbit of Rhea. Now the C_{22} term vanishes, and the doubleaveraged Hamiltonian is (6)
3.2. Interpretation for Rhea
The sin^{2}I factor is characteristic of a J_{2} perturbation. The presence of Saturn is equivalent to an increase in the oblateness of Rhea, and can be taken into account with an effective J_{2}. From Eq. (6), this effective J_{2} is given by (7)Close to the surface of Rhea, i.e., for a = R_{R}, the second part is already equal to 375 × 10^{6}. This means that the effect of Saturn quickly becomes dominating when a increases.
4. Extension of the model
4.1. Validity of our approximations
The averaging over several orbits of the particle around Rhea requires that τ_{R} be much larger than τ. The other averaging, over several revolutions of Rhea necessitates that τ_{p} be much larger than τ_{R}. We will study the most favorable possible case of a particle close to the moon’s surface: a ≃ R_{R}. With the definitions given in Eq. (1) we can rewrite the two previous conditions as (8)where q is the rotation parameter defined by (9)Figure 2 represents the space . The gray region corresponds to the region where Eqs. (8) are satisfied, and the dark gray zone to the sets of parameters for which each pair of time scales is separated by at least 10x. For Rhea, J_{2} = 889 × 10^{6} and (Anderson & Schubert 2007), so it lies in the dark region.
Fig. 2
Phase space in terms of Eqs. (8). Rhea lies in the dark gray zone, which means that we can safely apply our doubleaveraging method. Iapetus is not in that region. The dashed curve gives the position of homogenous, synchronous ellipsoidal objects. 
4.2. Iapetus
Iapetus’ J_{2} has never been measured, but Thomas et al. (2007) give a = 747.4 km, b = c = 712.4 km, and mean radius R_{I} = 735.6 km, yielding J_{2} = 1890 × 10^{6} (Yoder 1995), under the simplifying assumption that Iapetus is homogenous. With the definitions of the time scales given in Eq. (1), we calculate the values of Table 2 at a = 1.5R_{I}.
Different time scales for Iapetus at a = 1.5 R_{I}.
The condition required to make the averaging over several Iapetus rotations is not well satisfied: τ_{I}/τ_{p} = 0.4 (see Fig. 2). This would be an important effect if Iapetus did not have an azimuthal asymmetry. Thus, equatorial rings can be stable around Rhea and Iapetus, but for different reasons: for Iapetus it is the absence of a significant C_{22} term, while for Rhea the C_{22} term is neutralized by the fast revolution. If Iapetus’ C_{22} had the same magnitude as its J_{2}, polar rings would also be stable (Dobrovolskis et al. 1989), because the moon’s rotation is slower than the precession of a particle in orbit around it.
4.3. Exoplanets
We now consider whether the nonspherical shape of some exoplanets (such as hot Jupiters) can destabilize a putative ring system. In such an environment, any ring system would have to be composed of refractory material, not ice (Schlichting & Chang 2011). We choose a simple model of a homogenous synchronous planet at hydrostatic equilibrium. Therefore, we have (Murray & Dermott 1999) (10)This relation further simplifies the problem: the condition for averaging the exoplanet’s shape over its rotation (τ_{p}/τ_{R} ≫ 1) then becomes , which is automatically satisfied with the first condition for averaging the particle along its orbit (τ_{R}/τ ≫ 1 or equivalently ). We represent the hydrostatic relation graphically in Fig. 2, with the dashed curve .
We express q in terms of more physical parameters: the mean density of the exoplanet , the mass of the central star M_{∗}, and the semimajor axis of the exoplanet a_{e}. If we want a factor of 10 at least between each time scale, we can look for the limit parameter when the two others are fixed. The restrictive condition can then be written as (11)If (M_{∗},a_{e}) is set, there is a minimal mean density for the exoplanet below which we cannot apply our averagings: . Figure 3 shows the minimum density of the exoplanet as a function of the semimajor axis for several given masses of the star from 0.1M_{⊙} to 10M_{⊙}.
Fig. 3
The minimum density for which doubleaveraging is appropriate, as a function of the semimajor axis. Contours are for different masses of the central star. 
Figure 3 indicates that above approximately 0.1 AU, all exoplanets are dense enough to satisfy the conditions. For M_{∗} = M_{⊙} and a between 0.012 AU and 0.022 AU, only the planets with a density greater than Jupiter’s are dense enough. For even smaller values of the semimajor axis, only telluric planets have the required density. Very few exoplanets orbit so close to their star that no reasonable density can fit the criteria. For hot Jupiters, the typical values are about a_{e} = 0.05 AU and M_{∗} = M_{⊙}, so the minimum density equals 113 kg/m^{3}. This value is even lower than that of HD 209458 b, which is one of the least dense of the known exoplanets, with kg/m^{3} (Brown et al. 2001). For these exoplanets the conclusions about the role of the prolateness are the same as for Rhea.
4.4. Hill’s sphere limit
In the case of exoplanets, the star also directly attracts the test particle. Hill’s sphere radius can be estimated by R_{H} = a_{e}(M_{e}/ 3M_{∗})^{1 / 3}. It would have no meaning to consider rings outside Hill’s sphere. Therefore, in a system that is able to shelter rings, Hill’s radius must be greater than a, the semimajor axis of the orbiting particle. Setting a = αR_{e}, this condition can also be written in terms of q: (12)Even for very close rings, the maximum allowed value of q is on the order of 0.1. This is the most unfavorable case, and yet it nearly already satisfies our hypotheses for the double averaging . To sum up, we should be careful when q increases too much, but in these cases the attraction of the central body is likely to prevent any rings from forming.
5. Computational simulation for slower rotations
5.1. Model and results
Not all planets are homogenous or at hydrostatic equilibrium. Iapetus is a typical example of a planet whose J_{2} is too large to be caused by rotation at its distance to Saturn. Thomas et al. (2007) called this a “fossil bulge”, based on the likely idea that it is left over from a time when Iapetus rotated faster. Therefore we consider what happens when the rotational period of a moon (e.g., Rhea) τ_{R} is not short compared with the precession time τ_{p}. The Hamiltonian to consider is then the one given by Eq. (5). Given its complexity, we will use a simulation.
Our simulation will describe the behavior of a thin disk of interacting particles. To take into account these interactions, we will simply give our particle disk a kinematic viscosity ν, as in Dobrovolskis et al. (1989) and SteimanCameron & Durisen (1988). We modify the system of Eqs. (2) by adding a dissipative term (see Eqs. (9) and (10) in Dobrovolskis et al. 1989): (13)We assume that a is constant. Indeed, the settling times are much shorter than the amount of time where the inflow is significant (Dobrovolskis et al. 1989). The ring is then modeled as a set of 100 annular elements. Forty kilometers separate the innermost annulus from the outermost, and all annuli have the same initial inclination. Then, we take a time step of τ_{R}/ 100. The kinematic viscosity is estimated to be about 1 m^{2}/s, using the model of SteimanCameron & Durisen (1988), with ice particles of millimeter size. The settling time is denoted τ_{s}.
Testing our simulation with an object having the characteristics of Rhea, we find satisfactory results. No matter what the initial inclination is, the rings always settle into the equatorial plane, as expected. The settling time we admitted for our theoretical study was equal to the precession period τ_{p}. The simulation gives settling times τ_{s} which are closer to ten times τ_{p}. However, our goal was only to obtain a theoretical behavior for the rings, not to predict the detailed settling process. Such a difference is therefore acceptable.
We can now try to vary the parameters of the problem, and especially to study what happens when τ_{R} ≃ τ_{p}. We will decrease the rotational frequency n_{R}. This is equivalent to moving Rhea outward from Saturn while retaining its current shape (J_{2} and C_{22}) as a fossil bulge. We will call this new frequency to avoid confusion with Rhea’s proper characteristics.
We find that the behavior of the rings is roughly independent of . The rings still settle in the equatorial plane during a nearly constant duration τ_{s}. Table 3 shows that τ_{s} is not dependent on the rotational speed.
Settling time at several rotational frequencies.
5.2. Caveat
The simulation results should be interpreted with some caution. While the simulation robustly reproduces secular effects, it may not represent any effects that are not azimuthally averaged, such as any imparting of eccentricity to the ring or any resonant effects due to the rotating character of the C_{22} term.
6. Discussion
The separation between the different time scales in the case of Rhea allows us to carry out a double averaging and obtain a more understandable Hamiltonian. We have shown that rings can settle only in the equatorial plane. This rules out the possibility of rings in a region around Rhea that would not have been carefully investigated. The results of Tiscareno et al. (2010) then seem to definitively exclude the present existence of these rings.
Our model turned out to be interesting for other celestial objects. The most likely objects to have a shape close to Rhea’s are exoplanets orbiting close to their stars. Nearly all of them should be suitable for our model, except for those orbiting at less than 0.02 AU from their star. A weak density could put such rare planets into a more complex regime. On the other hand, such close exoplanets are unlikely to shelter rings owing to their small Hill spheres. Therefore, results obtained for Rhea can be extended to almost all synchronous exoplanets.
Finally, we carried out a computational simulation focused on what happens when the rotational period becomes as significant as the precession time. This happens when the J_{2} is larger than it should be, as it is in the case of Iapetus. Though it neglects all effects that are not azimuthally averaged, our simulation indicates that, over a wide range of rotational periods, a thin particle disk will become rings in the equatorial plane. The duration for the settling is roughly independent of the rotational period.
Thus, in all cases that were studied, the prolate shape seems to have no determining effect on the shape of the rings, keeping
in mind that classical shape models for the exoplanets were assumed. Nothing prevents the discovery of more original shapes, in particular higher prolate deformations allowing more imaginative ring systems.
Acknowledgments
We thank Joe Burns, Phil Nicholson, Matt Hedman, and Tom SteimanCameron for helpful conversations. We acknowledge funding from the École Normale Supérieure de Cachan and from NASA Outer Planets Research (NNX10AP94G).
References
 Anderson, J. D., & Schubert, G. 2007, Geophys. Res. Lett., 34, 2202 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, T. M., Charbonneau, D., Gilliland, R. L., et al. 2001, ApJ, 552, 699 [NASA ADS] [CrossRef] [Google Scholar]
 Dobrovolskis, A. R., Borderies, N. J., & SteimanCameron, T. Y. 1989, Icarus, 81, 132 [NASA ADS] [CrossRef] [Google Scholar]
 Dombard, A. J., Cheng, A. F., McKinnon, W. B., & Kay, J. P. 2012, J. Geophys. Res., 117, 3002 [CrossRef] [Google Scholar]
 Ip, W.H. 2006, Geophys. Res. Lett., 33, 16203 [NASA ADS] [CrossRef] [Google Scholar]
 Jones, G. H., Roussos, E., Krupp, N., et al. 2008, Science, 319, 1380 [NASA ADS] [CrossRef] [Google Scholar]
 Kaula, W. M. 1968, An Introduction to Planetary Physics: The Terrestrial Planets (New York: Wiley) [Google Scholar]
 Levison, H. F., Walsh, K. J., Barr, A. C., & Dones, L. 2011, Icarus, 214, 773 [NASA ADS] [CrossRef] [Google Scholar]
 Murray, C. D., & Dermott, S. F. 1999, Solar System Dynamics (Cambridge University Press) [Google Scholar]
 Schenk, P. M., & McKinnon, W. B. 2009, AAS Division for Planetary Sciences Meeting Abstracts, 41, 3.03 [Google Scholar]
 Schenk, P., Hamilton, D. P., Johnson, R. E., et al. 2011, Icarus, 211, 740 [NASA ADS] [CrossRef] [Google Scholar]
 Schlichting, H. E. & Chang, P. 2011, ApJ, 734, 117 [NASA ADS] [CrossRef] [Google Scholar]
 SteimanCameron, T. Y., & Durisen, R. H. 1984, ApJ, 276, 101 [NASA ADS] [CrossRef] [Google Scholar]
 SteimanCameron, T. Y., & Durisen, R. H. 1988, ApJ, 325, 26 [NASA ADS] [CrossRef] [Google Scholar]
 Thomas, P. C., Burns, J. A., Helfenstein, P., et al. 2007, Icarus, 190, 573 [NASA ADS] [CrossRef] [Google Scholar]
 Tiscareno, M. S., Burns, J. A., Cuzzi, J. N., & Hedman, M. M. 2010, Geophys. Res. Lett., 37, 14205 [NASA ADS] [CrossRef] [Google Scholar]
 Tremaine, S., Touma, J., & Namouni, F. 2009, AJ, 137, 3706 [NASA ADS] [CrossRef] [Google Scholar]
 Yoder, C. F. 1995, in Global Earth Physics: A Handbook of Physical Constants, ed. T. J. Ahrens, 1 [Google Scholar]
All Tables
All Figures
Fig. 1
Angles in the rotating frame. 

In the text 
Fig. 2
Phase space in terms of Eqs. (8). Rhea lies in the dark gray zone, which means that we can safely apply our doubleaveraging method. Iapetus is not in that region. The dashed curve gives the position of homogenous, synchronous ellipsoidal objects. 

In the text 
Fig. 3
The minimum density for which doubleaveraging is appropriate, as a function of the semimajor axis. Contours are for different masses of the central star. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.