Issue 
A&A
Volume 643, November 2020



Article Number  A67  
Number of page(s)  6  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/202038812  
Published online  03 November 2020 
Ring dynamics around an oblate body with an inclined satellite: the case of Haumea
Department of Physics and Astronomy, University of Padova, Via Marzolo 8,
35131
Padova,
Italy
email: Francesco.Marzari@pd.infn.it
Received:
1
July
2020
Accepted:
31
July
2020
Context. The recent discovery of rings and massive satellites around minor bodies and dwarf planets suggests that they may often coexist, as for example around Haumea.
Aims. A ring perturbed by an oblate central body and by an inclined satellite may disperse on a short timescale. The conditions under which a ring may survive are explored both analytically and numerically.
Methods. The trajectories of ring particles are integrated under the influence of the gravitational field of a triaxial ellipsoid and (a) massive satellite(s), including the effects of collisions.
Results. A ring initially formed in the equatorial plane of the central body will be disrupted if the satellite has an inclination in the Kozai–Lidov regime (39.2° < i < 144.8). For lower inclinations, the ring may relax to the satellite orbital plane thanks to an intense collisional damping. On the other hand, a significant J2 term easily suppresses the perturbations of an inclined satellite within a critical semimajor axis, even in the case of Kozai–Lidov cycles. However, if the ring is initially inclined with respect to the equatorial plane, the same J2 perturbations are not a protective factor but instead disrupt the ring on a short timescale. The ring found around Haumea is stable despite the rise in the impact velocities that is due to the asymmetric shape of the body and the presence of a 3:1 resonance with the rotation of the central body.
Conclusions. A ring close to an oblate central body should be searched for in the proximity of the equatorial plane, where the J2 perturbations protect it against the perturbations of an external inclined satellites. In an inclined configuration, the J2 term is itself disruptive.
Key words: planets and satellites: rings / celestial mechanics
© ESO 2020
1 Introduction
In the Kuiper Belt, all six of the largest bodies with diameters greater than 1000 km have a satellite system with a secondary–to–primary mass ratio ranging from 10^{−1} to 10^{−3}. According to Arakawa et al. (2019), these satellites likely formed via giant impacts in the early stage of solar system formation. Another interesting feature, possibly related to the outcome of a less energetic collision, is the presence of a ring around some of these bodies. The two centaurs 10 199 Chariklo and 2060 Chiron are known to have ring systems detected by stellar occultation (BragaRibas et al. 2014; Ortiz et al. 2015; Ruprecht et al. 2015); recently, Ortiz et al. (2017) found a ring around the transNeptunian dwarf planet Haumea, which is also orbited by two satellites, Namaka and Hi’iaka. These rings may have originated from the ejection of cratering fragments from the parent body, the disruption of a satellite (Pan & Wu 2016), or the partial tidal disruption of the mantle during a close encounter with a giant planet (Hyodo et al. 2016).
The longterm survival of a ring around a body possessing one or more inclined satellites can be jeopardized by their secular perturbations, which would act on a much shorter timescale than the viscous spreading (Salmon et al. 2010). These perturbations would be particularly disruptive for the ring coherence in the presence of Kozai–Lidov configurations, which force wide oscillations in both eccentricity and inclination (Naoz 2016). However, even a mild inclination may lead to a misalignment of the ring particle orbits once the node longitudes are randomized due to the difference in semimajor axes.
There are two ways a ring can survive in the presence of a massive satellite on an inclined orbit. In the absence of strong Kozai–Lidov perturbations, the collisional damping may be intense enough to align the ring to the satellite orbital plane, where it will survive without losing coherence. Alternatively, even for high inclinations of the satellite, the J2 term can suppress the satellite perturbations, keeping the ring flat in the equatorial plane, if the central body is oblate and the ring lies in its equatorial plane.
I will explore these scenarios and, in particular, I will focus on the contribution from the central body oblateness to the stability of the ring. The dynamics of the ring around the dwarf planet Haumea, a rapidly rotating triaxial ellipsoid with an elongated shape, will also be investigated. I will focus on the interplay between different evolutionary mechanisms that act on the ring, such as collisions and the gravitational perturbations from the irregular shape of the central body and from its two satellites.
In Sect. 2, I describe in detail the numerical model exploited to perform the numerical simulations of the ring evolution. Section 3 is devoted to the rivalry between the Kozai–Lidov cycles and the J2 perturbations by the central body in the presence of a highly inclined satellite. The less dramatic but still potentially harmful effects of secular perturbations from a mildly inclined satellite are explored in Sect. 4. In Sect. 5, I illustrate the scenario where the ring is not initially placed in the equatorial plane of an oblate central body and the J2 perturbations lead to a quick disruption of the ring. Section 6 focuses on the dynamics of Haumea’s ring, while the results of this whole dynamical exploration are discussed and commented on Sect. 7.
2 The numerical model
The ring evolution was modeled by numerically integrating the trajectories of the ring particles, treated as massless bodies, that evolve in the gravity field of a central body of mass M and one or more external perturbers (satellites). The 15th order Radau integrator (Everhart 1985) was used to better deal with the fast changes in the gravitational force acting on the particles. The central body was assumed to be a triaxial ellipsoid with three independent semiaxes, r_{a}, r_{b}, and r_{c}. The potential of this body was computed with MacCullagh’s formula (Murray & Dermott 1999): (1)
The inertial moments A, B, and C, in the case of a triaxial ellipsoid, are related to the semiaxes by the following relations: (3)
where ρ is the body bulk density. In this case, J_{2} can be obtained as (Turcotte et al. 2002) (4)
In the code, all possible twobody encounters were tested at each time step. Since the number of particles is limited by the CPU load and the number of impacts would be too low if I took the particles’ true sizes, an “inflated radius” was adopted for each particle (Brahic 1976). This inflated radius was made large enough to provide reliable impact statistics. The evolutionary speed was the same if the number of bodies was scaled by the ratio of the real particle radius r and the inflated radius R according to .
Each collision was treated as an inelastic rebound modeled with the algorithm of Brahic (1977), which is based on the assumption that the selfattraction of the particles can be neglected and that all particles have the same radius and mass. During a collision, the radial relative velocity was damped with a rebound coefficient η < 1 while the grazing one was left unaltered (frictionless spheres). The sizes of the ring particles were scaled as suggested by Brahic (1976), Thébault & Brahic (1998), and Charnoz et al. (2001).
3 Kozai–Lidov cycles suppression by the J2 perturbations
If the body has a massive satellite on a highly inclined orbit, the Kozai–Lidov perturbations may be a relevant disruption mechanismfor a ring in the equatorial plane of the central body. In this scenario, the Kozai–Lidov perturbations on the ring particles lead to a fast disintegration of the ring coherence if the central body is spherical, and most particles are ejected or impact the central body. This occurs even if the mutual collisions between the ring particles tend to damp the eccentricity and inclination.
I first performed a test numerical simulation with a spherical central body (r_{a} = r_{b} = r_{c}) that has the same mass as the TNO (transNeptunian object) Haumea (M_{c} = 4.006 × 10^{21} kg) and a satellite similar to Hi’iaka with a mass equal to m_{s} = M_{c}∕223. The orbit of the satellite has an eccentricity of 0.05 and an inclination of 55°. A set of 10 000 particles simulating the ring were initially distributed in the equatorial plane of the central body between 5000 and 10 000 km, with an eccentricity lower than 10^{−3} and an inclination in the range of 0° to 1°. An inflated diameter of 10 km was adopted, translating to approximately 10^{9} decimetersized particles according tothe scaling. The rebound coefficient η was set to 0.3.
After 10000 days of evolution, and despite the damping effects of the collisions, all the ring particles end up either on hyperbolic orbits or impact the central body because of the wide Kozai–Lidov oscillations of eccentricity and inclination. In this scenario, the ring is disrupted and could not survive longer than a few thousand days.
When I included the effects of the J2 term in the potential of the central body, assuming a 2:1 ratio between the r_{c} and r_{a} = r_{b} principal axes, the evolution is totally different. The ring maintains its coherence, even in the absence of collisions, and both the eccentricities and inclinations of the ring particles remain low, except at the locations of mean motion resonances with the satellite. The protection against the Kozai–Lidov perturbations is due to the fast precession of the pericenter argument induced by the J2 term. The ring evolution is shown in Fig. 1, where I compare thecases with and without collisions. In both scenarios, the ring particles, after 250 000 days, are still on loweccentricity and lowinclination orbits. At the resonance locations with the satellite (3:1 at about 9600 km, 4:1 at about 8000 km, etc.), the particles are excited by the resonant perturbations and spikes in both orbital elements are observed.
The critical distance within which the J2 disturbing potential term wins over the Kozai–Lidov mechanism can be estimated analytically by comparing the precession rates induced by the two perturbations. The equation for the variation of the pericenter argument is (Bertotti et al. 2003): (5)
In these equations, n is the satellite mean motion, R_{c} is the mean radius of thecentral body, and a, e, and i are the semimajor axis, eccentricity, and inclination of the ring particle, respectively.
The timescale associated with the precession of the pericenter argument can be estimated as and it is (6)
If the orbital inclination is small, cos(i) ~ 1, the above equation can be approximated as: (7)
For the Kozai–Lidov perturbations, the timescale associated with the precession of the argument of pericenter, in the quadrupole approximation (small eccentricity of the outer perturber), is given by (Antognini 2015; Naoz 2016) (8)
where m_{s} is the satellite mass, n_{s} is its mean motion, and e_{s} is its eccentricity. By comparing t_{KL} with t_{J2}_{i = 0}, a critical semimajor axis for which the two timescales are equal can be derived: (9)
with a_{s} being the semi–major axis of the satellite orbit.
Within a_{c}, the J2 term will prevail in determining the precession of the argument of the ring particle pericenter, and the Kozai–Lidov oscillations will be suppressed. Beyond a_{c}, the Kozai–Lidov perturbations dominate over the J2 effects and the wide changes in eccentricity and inclination will be restored, causing a progressive disruption of a putative ring.
In the top panel of Fig. 2, the two timescales t_{KL} and t_{J2} _{i = 0} are compared for different values of the ratio between r_{c} and r_{a}, showing that a larger flattening of the central body shifts the critical semimajor axis, within which the ring is stable against the Kozai–Lidov perturbations, outward. In the bottom panel, the value of a_{c} is computed for decreasing values of the r_{c}∕r_{a}, showing the range in which the J2 perturbations can protect an equatorial ring. These calculations show that a significant oblateness of the central body is crucial to ensuring the survival of an equatorial ring around the body in the presence of a highly inclined satellite.
Fig. 1 Orbital distribution of the ring particles after 250 000 days from the beginning of the simulation in the case without collisions (green circles) and with collisions (magenta circles). The observed spikes in both e and i are due to mean motion resonances with the satellite. 
4 Secular perturbations of an inclined satellite and final ring inclination
If the inclination of the perturbing satellite is mild (lower than ~ 39°), the evolution of the inclination and eccentricity is less dramatic. I ran a model with the inclination of the perturber set to 15°, switching the contribution from the J2 term on and off; the results are shown in Fig. 3. In the case of a spherical central body, the combined effect of collisions and secular perturbations by the satellite forces the ring to relax to the same orbital plane as the satellite. As illustrated in Fig. 3, the inclination of the ring particles has small oscillations around the inclination of the satellite, which may be further damped at later times, while the nodal longitudes are aligned; this dynamical configuration maintains the coherence of the ring (Fig. 4). Collisions act as a viscous force that tends to align the ring to the plane of the perturber. The ring appears slightly warped close to the central body (see Fig. 3) while it is well aligned in the central part. The warping is probably a temporary feature and may disappear at later times.
In the J2 term case, the eccentricities and inclinations are on average always small, and the ring is coherent and does not depart from the equatorial plane of the central body. The outer regions of the disk, approximately beyond 8000 km, show some scattering around the middle plane where the protective effect of the J2 term weakens. An analytical estimate of the limiting semimajor axis, within which the J2 holds the ring in the equatorial plane, is more complex in this case. The second order Lagrange–Laplace secular theory is a good approximation for low values of eccentricity and inclination. In this particular case, the inclination is set to 15° and the theory may not be very accurate. However, this secular theory can give an idea of the timescales involved in the dynamical evolution of the ring. According to Murray & Dermott (1999), the proper oscillation timescale for the pericenter longitude is given by T_{sec} = 2π∕A, where A reads: (10)
This frequency has a complex dependence on the semi–major axis of the ring particle due to the Laplace coefficient . For this reason, numerical computations lead to an easier evaluation of the value of semimajor axis a_{c}, for which the two time derivatives of the pericenter longitude, the one due to the J2 perturbations and the secular one, match. In Fig. 5, the two timescales associated with the frequencies are compared. The critical semimajor axis is located at approximately 15 000 km and is in good agreement with the numerical modeling (Fig. 3).
Fig. 2 Top panel: t_{KL} shown as a function of the semimajor axis (dashed black line) and t_{J2}_{i = 0} plotted for different values of the ratio between the principal axes of the ellipsoid r_{c} ∕r_{a}. The blue line is for r_{c}∕r_{a} = 1∕2, the magenta line is for r_{c}∕r_{a} = 3∕4, and the green line is for r_{c}∕r_{a} = 9∕10. Bottom panel: critical semimajor axis drawn as a function of the r_{c}∕r_{a} ratio. 
Fig. 3 Orbital distribution of the ring particles after 150 000 days in the J2 term (magenta circles) and spherical cases (green circles). As in Fig. 1, the observed spikes in both e and i are due to mean motion resonances with the satellite. 
Fig. 4 Distribution of the node longitudes of the ring particles after 150 000 days in the case without the J2 term. The alignment, at least up to 8000 km, ensures the compactness of the ring. 
5 Disruption of a ring misaligned with respect to the equatorial plane of the central body
If the ring isformed on a plane that is tilted with respect to the equatorial plane of an oblate central body, the J2 term quickly disrupts the ring. The forced circulation of the nodal longitudes dephases the inclined orbits of the ring particles, and the initial ring is transformed into a torus. The fattening of the ring depends on its initial inclination with respect to the equatorial plane. In Fig. 6, the threedimensional distribution of the ring particles is shown. They quickly disperse on a timescale on the order of a few 10^{3} yr and do not form a flat coherent structure anymore. In this case, the J2 perturbation is not a protective factor but rather a disruptive influence. As a consequence, potential rings around an oblate body must be sought in the equatorial plane.
Fig. 5 Comparison of the timescale for the secular precession of the pericenter longitude (dashed black line) and the timescale due to the J2 term (blue continuous line). 
Fig. 6 Threedimensional distribution of the ring particles after 2000 days. The nodal precession forced by J2 transforms the ring into a torus. 
6 The case of Haumea: interplay among its complex shape, fast rotation, and two satellites
Haumea, a dwarf planet orbiting beyond Neptune, has a very elongated shape and a fast rotation around the principal axis of inertia (~ 3.92 h, Lellouch et al. 2010). It is modeled as a triaxial ellipsoid, and its semiaxes are r_{a} = 1161 ± 30 km, r_{b} = 852 ± 4 km, and r_{c} = 513 ± 16 km (Ortiz et al. 2017) with a ratio between r_{c} and r_{a} of about 1:2 (same as the test case adopted in the previous sections). Haumea is orbited by two satellites, Namaka and Hi’iaka, and by a ring coplanar with the equator of Haumea with a radius of ~ 2287 km and a width of ~70 km (Ortiz et al. 2017). According to the previous results, the coplanarity of the ring with the equator of Haumea is expected since the body is a strongly oblate body (see Sect. 5).
The more massive satellite Hi’iaka (its mass is approximately 1/223 times the mass of Haumea) has a semimajor axis of 49 880 km (Ragozzine & Brown 2009); it is almost on a circular orbit (e ~ 0.05) and, according to Ortiz et al. (2017), has an orbital plane coplanar with the equator of Haumea. One therefore expects that its perturbations on the inclination of the ring are negligible. Namaka has a mass about ten times smaller than that of Hi’iaka, but it is closer to Haumea since its semimajor axis is 25 657 km, about half that of Hi’iaka. It is on an eccentric (e ~ 0.25) and inclined (i ~ 13°) orbit with respect to the equatorial plane of Haumea, and it may cause longterm perturbations on the inclination of the ring particles, potentially dispersing them.
In Fig. 7, the evolution of a single particle initially set in the middle of the ring is shown in two different configurations. The blue dashed line shows the inclination evolution of the particle for a spherical central body (r_{a} = r_{b} = r_{c}) while the red line shows the evolution of the same particle when the real values of r_{a}, r_{b}, and r_{c} and the rotation of the body are included. In the first case, the inclination of the particle regularly oscillates between 1° and 4° due to the secular perturbations of Namaka. As expected from Sect. 4, when the J2 term is included and I consider the full model for Haumea, the inclination is almost constant with very small oscillations around the initial value. Even in the case of Haumea’s ring, the J2 term helps in keeping the ring coherent in the equatorial plane of the body.
To explore the evolution of the ring, also including the effects of mutual collisions, I generated a population of 10 000 particles with initial semimajor axes between 2235 and 2340 km, eccentricity lower than 10^{−4}, and inclination lower than 0.5° (both randomly selected). An inflated diameter of 3 km was adoptedin order to attain significant collisional activity. In Fig. 8, the eccentricity and inclination distribution is compared after 10^{4} days of evolution in the case where only the J2 term is included (r_{b} − r_{a} = 0) and in the real case where r_{b} − r_{a}≠0. When r_{b} = r_{a}, the eccentricity of all particles in the ring is very low and the magenta dots, which represent them in Fig. 8, are barely visible above zero. The inclination does not change significantly from the initial value and, as shown in Fig. 7, it is not expected to grow since the J2 perturbations suppress the secular perturbations of Namaka. When the perturbations due to r_{b} ≠ r_{a} are included (full model), the potential of the central body excites the eccentricities of the ring particles to higher values, as shown in the top panel of Fig. 8, despite the damping effect of collisions. This is potentially related to the 3:1 resonance between the rotation period of Haumea and the orbital period of the ring particles (Winter et al. 2019). On the other hand, the inclination is on average smaller than the case with J2only, as shown in the bottom panel of Fig. 8.
The high values of eccentricity induced by the r_{a}≠r_{b} term of the potential lead to an increase in the mutual velocities between the ring particles. This may speed up the erosion of the ring if they are too high and may possibly cause fragmentation of the impactors most of the time. InFig. 9, I show the normalized distribution of the impact velocities when the steady stateis reached (Fig. 8). The large majority of collisions occur at a speed lower than 1 m s^{−1}, which, according to Bukhari Syed et al. (2017), is below the erosion limit for cmsize particles. Therefore, the increase in the impact velocity due to the high eccentricity, possibly related to the 3:1 resonance, is not enough to cause a quick erosion of Haumea’s ring.
Fig. 7 Time evolution of the inclination of a ring particle for a spherical central body (blue dashed line) and for a rotating triaxial ellipsoid (red line), like the real Haumea. 
Fig. 8 Orbital distribution of the ring particles after 10^{4} days. Top panel: eccentricity vs. semimajor axis is shown in two different models. The magenta dots illustrate the case where B = A≠C (only the J2 term) while the green ones show the outcome of the full model. The inclination vs. semimajor axis is shown in the top panel. 
Fig. 9 Distribution of the mutual impact velocities between the ring particles of Haumea when a collisional steady state is reached. 
7 Conclusions
The dynamical evolution of a ring around a minor body or dwarf planet that also possesses one or more satellites is complex. Different dynamical mechanisms come into play and may lead to the quick disruption of the ring. If the satellite is inclined with respect to the ring plane, the oscillations in inclination, in particular in the Kozai–Lidov regime, can fully disperse the ring on a short timescale. In addition, if the central body is significantly oblate, the fast pericenter precession induced by the J2 term quickly disperses the putative rings that are significantly inclined with respect to the equatorial plane.
To determine under which conditions a ring can survive, I exploited a numerical model that includes the oblateness of the central body, the mutual collisions between the ring particles, and the gravitational perturbations of one or more satellites. Analytical calculations have also been performed to evaluate the relative strength and contribution of the different perturbing mechanisms.
I find that the J2 term is doublefaced. It is disruptive for inclined rings, but it protects equatorial rings from the perturbations of an inclined satellite. This occurs for the strong Kozai–Lidov oscillations and also for the milder secular perturbations that occur when the inclination is lower than 39°. The J2 term suppresses the inclination perturbations of the satellite within a critical semimajor axis, which depends on the oblateness of the central body, by forcing a fast pericenter precession that dominates over that induced by the secular terms of the satellite. This is a general result that can be applied to any system, and it allows us to predict that rings formed around oblate bodies must be preferentially searched for in the equatorial plane. Rings that are born inclined with respect to this plane can survive only if the central body has a small J2. The dynamical structure of the central body and satellite determines which kind of rings can survive after their formation.
The Haumea case was also investigated with numerical simulations, and I find that its ring is forced by the strong J2 perturbations of the oblate triaxial ellipsoid to lie in the equatorial plane. In addition, the r_{b} ≠ r_{a} gravitational perturbations excite large eccentricities among the ring particles that are not damped by collisions. However, the relative impact velocities are lower than the expected fragmentation velocity and the survival of the ring is guaranteed.
Acknowledgements
I would like to thank the referee, Anthony Dobrovolskis, for his comments and suggestions and Mariasole Maglione for her help with the analytical calculations.
References
 Antognini, J. M. O. 2015, MNRAS, 452, 3610 [NASA ADS] [CrossRef] [Google Scholar]
 Arakawa, S., Hyodo, R., & Genda, H. 2019, Nat. Astron., 3, 802 [CrossRef] [Google Scholar]
 Bertotti, B., Farinella, P., & Vokrouhlick, D. 2003, Physics of the Solar System – Dynamics and Evolution, Space Physics, and Spacetime Structure (Dordrecht: Kluwer Academic Publishers), 293 [Google Scholar]
 BragaRibas, F., Sicardy, B., Ortiz, J. L., et al. 2014, Nature, 508, 72 [NASA ADS] [CrossRef] [Google Scholar]
 Brahic, A. 1976, J. Comput. Phys., 22, 171 [NASA ADS] [CrossRef] [Google Scholar]
 Brahic, A. 1977, A&A, 54, 895 [NASA ADS] [Google Scholar]
 Bukhari Syed, M., Blum, J., Wahlberg Jansson, K., & Johansen, A. 2017, ApJ, 834, 145 [NASA ADS] [CrossRef] [Google Scholar]
 Charnoz, S., Thébault, P., & Brahic, A. 2001, A&A, 373, 683 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Everhart, E. 1985, Astrophys. Space Sci. Lib., 115, 185 [NASA ADS] [CrossRef] [Google Scholar]
 Hyodo, R., Charnoz, S., Genda, H., & Ohtsuki, K. 2016, ApJ, 828, L8 [CrossRef] [Google Scholar]
 Lellouch, E., Kiss, C., SantosSanz, P., et al. 2010, A&A, 518, L147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Murray, C. D., & Dermott, S. F. 1999, Solar System Dynamics (Cambridge: Cambridge University Press) [Google Scholar]
 Naoz, S. 2016, ARA&A, 54, 441 [Google Scholar]
 Ortiz, J. L., Duffard, R., PinillaAlonso, N., et al. 2015, in European Planetary Science Congress, EPSC2015–230 [Google Scholar]
 Ortiz, J. L., SantosSanz, P., Sicardy, B., et al. 2017, Nature, 550, 219 [NASA ADS] [CrossRef] [Google Scholar]
 Pan, M., & Wu, Y. 2016, ApJ, 821, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Ragozzine, D., & Brown, M. E. 2009, AJ, 137, 4766 [NASA ADS] [CrossRef] [Google Scholar]
 Ruprecht, J. D., Bosh, A. S., Person, M. J., et al. 2015, Icarus, 252, 271 [NASA ADS] [CrossRef] [Google Scholar]
 Salmon, J., Charnoz, S., Crida, A., & Brahic, A. 2010, Icarus, 209, 771 [CrossRef] [Google Scholar]
 Thébault, P., & Brahic, A. 1998, Planet. Space Sci., 47, 233 [NASA ADS] [CrossRef] [Google Scholar]
 Turcotte, D., Schubert, G., & Schubert, J. 2002, Geodynamics (Cambridge: Cambridge University Press) [CrossRef] [Google Scholar]
 Winter, O. C., BorderesMotta, G., & Ribeiro, T. 2019, MNRAS, 484, 3765 [CrossRef] [Google Scholar]
All Figures
Fig. 1 Orbital distribution of the ring particles after 250 000 days from the beginning of the simulation in the case without collisions (green circles) and with collisions (magenta circles). The observed spikes in both e and i are due to mean motion resonances with the satellite. 

In the text 
Fig. 2 Top panel: t_{KL} shown as a function of the semimajor axis (dashed black line) and t_{J2}_{i = 0} plotted for different values of the ratio between the principal axes of the ellipsoid r_{c} ∕r_{a}. The blue line is for r_{c}∕r_{a} = 1∕2, the magenta line is for r_{c}∕r_{a} = 3∕4, and the green line is for r_{c}∕r_{a} = 9∕10. Bottom panel: critical semimajor axis drawn as a function of the r_{c}∕r_{a} ratio. 

In the text 
Fig. 3 Orbital distribution of the ring particles after 150 000 days in the J2 term (magenta circles) and spherical cases (green circles). As in Fig. 1, the observed spikes in both e and i are due to mean motion resonances with the satellite. 

In the text 
Fig. 4 Distribution of the node longitudes of the ring particles after 150 000 days in the case without the J2 term. The alignment, at least up to 8000 km, ensures the compactness of the ring. 

In the text 
Fig. 5 Comparison of the timescale for the secular precession of the pericenter longitude (dashed black line) and the timescale due to the J2 term (blue continuous line). 

In the text 
Fig. 6 Threedimensional distribution of the ring particles after 2000 days. The nodal precession forced by J2 transforms the ring into a torus. 

In the text 
Fig. 7 Time evolution of the inclination of a ring particle for a spherical central body (blue dashed line) and for a rotating triaxial ellipsoid (red line), like the real Haumea. 

In the text 
Fig. 8 Orbital distribution of the ring particles after 10^{4} days. Top panel: eccentricity vs. semimajor axis is shown in two different models. The magenta dots illustrate the case where B = A≠C (only the J2 term) while the green ones show the outcome of the full model. The inclination vs. semimajor axis is shown in the top panel. 

In the text 
Fig. 9 Distribution of the mutual impact velocities between the ring particles of Haumea when a collisional steady state is reached. 

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.