Spinorbit coupling and chaotic rotation for circumbinary bodies^{⋆}
Application to the small satellites of the PlutoCharon system
^{1}
CIDMA, Departamento de FísicaUniversidade de Aveiro,
Campus de Santiago,
3810193
Aveiro, Portugal
email:
correia@ua.pt
^{2}
ASD, IMCCECNRS UMR8028, Observatoire de Paris,
UPMC, 77 Av.
DenfertRochereau, 75014
Paris,
France
^{3}
Université Pierre et Marie Curie, UPMC −
Paris 06,
France
Received: 21 June 2015
Accepted: 30 July 2015
We investigate the resonant rotation of circumbinary bodies in planar quasicircular orbits. Denoting n_{b} and n the orbital mean motion of the inner binary and of the circumbinary body, respectively, we show that spinorbit resonances exist at the frequencies n ± kν/2, where ν = n_{b} − n, and k is an integer. Moreover, when the libration at natural frequency has the same magnitude as ν, the resonances overlap and the rotation becomes chaotic. We apply these results to the small satellites in the PlutoCharon system, and conclude that their rotations are likely chaotic. However, the rotation can also be stable and not synchronous for small axial asymmetries.
Key words: planetary systems / planets and satellites: dynamical evolution and stability / planets and satellites: individual: Pluto / planets and satellites: individual: Charon / celestial mechanics
Figure 4 and Appendices are available in electronic form at http://www.aanda.org
© ESO, 2015
1. Introduction
Circumbinary bodies are objects that orbit around a more massive binary system. In the solar system, the small satellites of the PlutoCharon system are the best example (e.g., Brozović et al. 2015). Planets orbiting two stars, often called circumbinary planets, have also been reported (e.g., Correia et al. 2005; Doyle et al. 2011; Welsh et al. 2012). Most of these bodies are close enough to the central binary to undergo tidal dissipation, which slowly modifies the rotation rate until it becomes close to the mean motion (e.g., MacDonald 1964; Correia et al. 2014).
For an object that moves in an eccentric orbit around a single massive body, the rotation rate can be captured in a halfinteger commensurability with the mean motion, usually called a spinorbit resonance (Colombo 1965; Goldreich & Peale 1966). Moreover, for very eccentric orbits or large axial asymmetries, the rotational libration width of the individual resonances may overlap, and the rotation becomes chaotic (Wisdom et al. 1984; Wisdom 1987). However, for nearly circular orbits, all these equilibria disappear and the only possibility for the spin is synchronous rotation (e.g., Correia & Laskar 2009).
When a third body is added to the classic twobody problem, the mutual gravitational perturbations introduce additional spinorbit resonances at the perturbing frequency (Goldreich & Peale 1967; Correia & Robutel 2013). In the case of circumbinary bodies there is, in addition, a permanent misalignment of the long inertia axis of the rotating body from the radius vector pointing to each inner body (Fig. 1). The resulting torque on the rotating body’s figure induces some rotational libration that may give rise to some unexpected behaviors for the rotation rate (Showalter & Hamilton 2015). In this Letter we investigate all possibilities for the rotation of circumbinary bodies, and apply them to the spin evolution of the small satellites of the PlutoCharon system.
2. Spinorbit coupling
We consider a planar threebody hierarchical system composed of a central binary with masses m_{0} and m_{1}, together with an external circumbinary companion with mass m (Fig. 1), where m ≪ m_{1} ≤ m_{0}. For the orbits we use Jacobi canonical coordinates, with r_{b} being the position of m_{1} relative to m_{0} (inner orbit), and r the position of m relative to the center of mass of m_{0} and m_{1} (outer orbit). The body with mass m has principal moments of inertia A ≤ B ≤ C and angular velocity given by ω. In Appendix A we provide the full equations of motion for the spin. Since tidal dissipation usually damps the obliquity to zero (e.g., Hut 1980; Correia 2009), for simplicity, here we describe the motion for ω normal to the orbit. In absence of tides, the equation for the rotation angle, θ, is then given by (Eq. (A.5)) (1)where a and n are the semimajor axis and the mean motion of the outer orbit, respectively, (r_{i},f_{i}) are the radial and the angular coordinates of r_{i}, The parameter δ is the mass ratio of the inner binary, while σ is approximately the frequency of smallamplitude rotational librations in an unperturbed synchronous resonance.
Fig. 1 Reference angles and coordinates. The system is coplanar and is composed of an inner binary with masses m_{0} and m_{1}, and an outer body m that rotates with angular velocity ω. 

Open with DEXTER 
We can express , sinf_{0}, and cosf_{0} as (4)and (5)Similar expressions can be found for , sinf_{1}, and cosf_{1}, if we replace δ by (δ − 1). Therefore, since   r_{b}   <   r  , we can develop Eq. (1) in power series of r_{b}/r. In addition, for the hierarchical problem the orbits are approximately ellipses. Thus, r and r_{b} can also be expanded in terms of Hansen coefficients (e.g., Laskar & Boué 2010). We then rewrite Eq. (1) as (6)where e is the eccentricity; M is the mean anomaly; β_{l,p,q} ∝ α^{l}, with α = a_{b}/a; and p, q are halfintegers. Spinorbit resonances occur whenever . In Appendix B we provide the complete expression for Eq. (6) to the first order in the eccentricities and to the second order in α.
3. Near circular orbits
For near circular orbits (e ≈ 0 and e_{b} ≈ 0), we can neglect both eccentricities, thus r ≈ a and r_{b} ≈ a_{b}. Retaining terms in α^{3} in Eq. (6) gives
(7)where γ = θ − M, φ = M_{b} − M, (8)We see that we have several islands of rotational libration, (9)Each term individually behaves like a pendulum where the rotation can be trapped. Therefore, together with the classic synchronous equilibrium at , there are additional possibilities for the spin at the super and subsynchronous resonances , with k ∈Z.
Terms with k ≠ 0 depend on σ^{2}, δ, and α^{2}, that is, on the axial asymmetry of the body, on the mass ratio of the central binary, and on the semimajor axis ratio. Thus, for circumbinary bodies far from the inner binary (α ≪ 1), the amplitude of these higher order terms decreases very quickly to zero. For δ ≈ 0 (i.e., m_{1} ≪ m_{0}) we have ρ_{2} = ρ_{3} = 0, thus only the synchronous resonance subsists, as in the classic circular twobody problem (e.g., Goldreich & Peale 1966). The effect of the binary mass ratio is maximized for δ = 1/2 (i.e., m_{0} = m_{1}), but only for the second order resonances (terms in ρ_{2}), since third order resonances also vanish for equal masses (ρ_{3} = 0). Interestingly, the size of the subsynchronous resonance with k< 0 is always larger than the symmetrical supersynchronous resonance with k> 0.
3.1. Chaotic motion
When the libration amplitudes of some individual resonant islands overlap, the rotational motion can be chaotic (Chirikov 1979; Morbidelli 2002). This means that the rotation exhibits random variations in short periods of time. This phenomenon was first described for the rotation of Hyperion, which is chaotic because its orbit is eccentric (Wisdom et al. 1984).
In the circular approximation, each individual resonance is placed at , with libration width σ, where β_{k} = ∑ _{l}β_{l,(1 − k)/2,k/2}. Overlap between two resonances with k_{1}<k_{2} then occurs whenever (10)The synchronous resonance (k = 0) has the largest width, so overlap is more likely for this resonance. The third order resonance with k = −1 is the nearest resonance with larger amplitude, thus chaos sets in for (11)Rotation in the chaotic zone is also attitude unstable (Wisdom et al. 1984; Wisdom 1987). A nonzero obliquity introduces additional resonant terms to the rotation (Eq. (7)), increasing the chances of overlapping (e.g., Correia & Laskar 2010). Therefore, in order to correctly account for the chaotic behavior, we need to integrate the full equations of motion (Appendix A).
3.2. Global dynamics
Fig. 2 Stability analysis of the rotation rate close to the synchronization for δ = 0.1, α ∈ [ 0.25:0.5 ], and (B − A) /C ∈ [ 0:0.9 ]. We adopt three different initial obliquity values ψ = 1°, 30°, and 60°. The color index indicates the proportion of chaotic orbits inside the studied domain: from dark blue for fully regular to red for entirely chaotic. The dashed line (Eq. (11)) roughly delimits the zone where chaotic motion is possible. The dots identify two representative cases that we detail in Fig. 3. 

Open with DEXTER 
The general problem of the spinorbit coupling for near circular circumbinary bodies can be reduced to the analysis of , which depends on the libration frequency σ/n, the mass ratio δ, and the semimajor axis ratio α (Eq. (7)). However, since δ ∈ [ 0,1/2 ] and it is a wellknown parameter for a specific binary system, the global dynamics is approximately controlled by the only two free parameters: α (or n_{b}/n) and σ/n (or (B − A) /C).
Rescaling the time with n = 1, we can perform a stability analysis of in the plane (α,σ) to quickly identify the rotational regime for any system of quasicircular circumbinay bodies that is near the synchronous equilibrium. If isolated, the halfwidth of the synchronous resonant island in the direction of is equal to σ. When perturbed, chaos is expected for σ ~ ν (Eq. (11)). Thus, for a given (α,σ) pair, we fix the initial value of γ at the libration center of the synchronous resonance, and select 400 equispaced initial values of . The corresponding solutions are integrated using the equations for the threebody problem together with Eqs. (A.1) and (A.2) for the rotational motion (without tidal dissipation). The dynamical nature (stable/unstable) is deduced from frequency analysis (Laskar 1990, 1993), which gives the fraction of chaotic trajectories.
In Fig. 2 we show the general results for δ = 0.1 using three different initial obliquity values ψ = 1°, 30°, and 60°. All the remaining variables are initially set to zero. The color index indicates the proportion of chaotic orbits inside the studied domain: from dark blue for fully regular to red for entirely chaotic. Depending on the value of α and σ, the rotation can present different behaviors, ranging from nonsynchronous equilibria to chaotic motion. Figure 2 remains almost unchanged for δ> 0.1 since the global dynamics is not very sensitive to δ. The only exception would be for δ ≪ 0.1 since in that case the system would behave like the classic circular twobody problem where the rotation can only be synchronous.
A common way of visualizing and understanding the different regimes is to use frequency map analysis (Laskar 1993). We selected two representative pairs (α,σ) and plotted the corresponding crosssection maps for the three initial obliquities in Fig. 3. We are able to identify the stable and the chaotic regions very easily. We observe that stable islands can exist in the middle of chaos, implying a large diversity of possibilities for the final rotation. In the chaotic regions for the rotation, the obliquity is also chaotic and can vary between 0° and 180°. For high obliquities, the size of the stable synchronous island shrinks and the chaotic zone is extended. Therefore, depending on the obliquity, the rotation alternates between more or less chaotic. In times of small obliquity, the rotation can be trapped in a stable resonant island, preventing the obliquity from increasing again, and the spin can become stable (see Fig. 4).
Fig. 3 Cross sections representing the spin dynamics of Hydra (top) and Nix (bottom) at different initial obliquities ψ = 1°, 30°, and 60°. The dots show the main frequency of γ(t), denoted by ξ, for different initial values of . It is obtained by numerical integration of the threebody problem together with Eqs. (A.1) and (A.2) for the rotational motion (without dissipation). θ is initially set equal to f so the sections go through the middle of the synchronous resonance (). All the remaining variables are initially set to zero. The dots are blue for stable rotation (libration or circulation), and red for chaotic motion. The distinction between stable and unstable trajectories is based on the estimate of the second derivative of ξ with respect to the initial condition (see Laskar 1993). The plateaus correspond to the resonance crossings (ξ/ν = k/2). The red dots randomly distributed between the plateaus indicate the overlapping of the associated spinorbit resonances, which generates chaotic motion. 

Open with DEXTER 
For small values of α and/or σ, the individual resonances (Eq. (7)), associated with the plateaus in Fig. 3, are well separated. Thus, for dissipative systems, the rotation can be captured in individual spinorbit resonances and stay there. For rotation rates increasing from lower values, the subsynchronous resonance with has the largest amplitude, so this is the most likely spinorbit resonance to occur. For rotation rates decreasing from higher values, the supersynchronous resonance with becomes the most likely possibility. This case is very interesting as it corresponds to a synchronous rotation with the inner orbit period. However, since the amplitude of this resonance is small, the rotation can easily escape it and subsequently evolve into the classic synchronous rotation with the outer orbit ().
The amplitude of the resonant terms increases with α and σ. At some point, the individual libration islands merge and chaotic motion can be expected (see Sect. 3.1). The transition between the two regimes is roughly given by the dashed curved, obtained with expression (11). Near the transition regime, chaotic rotation can be observed around the separatrix of the synchronous resonance, but the motion is still regular nearer the center (Fig. 3, left). Therefore, when the rotation of a body is evolving by tidal effect, the rotation becomes chaotic when approaching the synchronous resonance. However, depending on the strength and on the geometry of the tidal torque, it is still possible to find a pathway to the synchronous rotation. The wandering in the chaotic region may also provide a path into another stable super or subsynchronous resonance (see Fig. 4).
For α> 0.4 and σ ~ 1, there is a large overlap between several resonant terms, in particular for the negative ones (k< 0). As a consequence, there is a large chaotic zone for the spin (Fig. 3, bottom). Some of the individual resonances may still subsist, including the synchronous one, but they have small stable widths and they can only be reached at periods of small obliquity. The chaotic motion is maximized for σ ≈ 1, which corresponds to (B − A) /C ≈ 0.35, because secondary resonances generate strong instabilities inside the synchronous island (Robutel et al. 2012). Finally, for α ≈ 0.49 the system reaches the 3/1 mean motion resonance, which introduces additional forcing to the rotation, and all trajectories become chaotic.
4. Application to the PlutoCharon system
In 1978, regular series of observations of Pluto showed that the images were consistently elongated revealing the presence of Charon (Christy & Harrington 1978). The orbital parameters hence determined have shown that the two bodies evolved in an almost circular orbit with a 6.4day period (e.g., Tholen & Buie 1997). Charon has an important fraction of the total mass (about 11%), and therefore the system is considered a binary planet rather than a planet and a moon. Later, it was found that four tiny satellites move around the barycenter of the PlutoCharon system (Weaver et al. 2006; Brozović et al. 2015), also in nearly circular and coplanar orbits (see Table 1).
Mean orbital and physical parameters for Pluto’s moons (Brozović et al. 2015).
Pluto’s brightness also varies with a period of 6.4 days (e.g., Tholen & Tedesco 1994). Since Charon is too dim to account for the amplitude of the variation, this period has been identified as the rotation of Pluto. Therefore, at present, the spin of Pluto is synchronous with the orbit of Charon, a configuration acquired from the action of tidal torques raised on Pluto by Charon (e.g., Farinella et al. 1979; Cheng et al. 2014). Tidal torques raised on Charon and on the remaining satellites by Pluto are even stronger, so all these bodies are presumably also tidally evolved.
Showalter & Hamilton (2015) have measured the brightness variations of the small satellites. They concluded that Hydra and Nix show no obvious pattern, suggesting that their rotation is chaotic. Assuming that these two satellites are uniform triaxial ellipsoids, they additionally estimate their figure semiaxis ratios, which allows us to compute B/C ≈ 0.92 and A/C ≈ 0.44 for Hydra, and B/C ≈ 0.94 and A/C ≈ 0.31 for Nix.
In Fig. 2 we already show the global dynamics for the PlutoCharon system, for which δ = 0.1085 ≈ 0.1 (Brozović et al. 2015). Therefore, we plot vertical lines at the α values corresponding to the small satellites in the system (Table 1). Moreover, since we have an estimation for the axial asymmetries of Hydra and Nix, we also show a dot for their (B − A) /C values in Fig. 2. These two satellites are in a different dynamical regime: the rotation of Hydra can present a stable spinorbit coupling (Fig. 3, top), while the rotation of Nix is most likely chaotic (Fig. 3, bottom), as discussed in Sect. 3.2. The spin dynamics of Kerberos is probably similar to that of Hydra, while the rotation of Styx can be even more chaotic than that of Nix, depending on their (B − A) /C values.
To test the reliability of the dynamical picture described in the previous section, we now add the effect from tides to our model, so that we can follow the longterm evolution of the spin of these satellites. For simplicity, we adopt a constant timelag linear model, whose contribution to the spin is given by expression (A.13). In Fig. 4 we show some examples for the final evolution of Hydra’s spin with tidal dissipation, starting with slightly different initial values of θ. We adopt an initial retrograde rotation of 4.4 days (), and 30° for the initial obliquity, to force the rotation to cross the large amplitude subsynchronous resonances. An initial retrograde rotation for such small bodies is as likely as a prograde one (e.g., Dones & Tremaine 1993). For the uncertain parameters we used R = 45 km, k_{2}/Q = 10^{4}, and C/(mR^{2}) = 0.4.
In one example (Fig. 4a) the rotation of Hydra is trapped in the spinorbit resonance with , while in the others the rotation reaches the chaotic zone. However, after some wandering in this zone some simulations are able to find a path into a stable spinorbit resonance with k> 0 (Fig. 4b) or with k< 0 (Fig. 4c). Interestingly, while for resonances with k> 0 the spin axis stabilizes near zero degrees, for resonances with k< 0, the spin axis stabilizes with a high obliquity value. The last example (Fig. 4d) remained chaotic for the length of the integration and therefore can represent the observed present state.
Over 10 Myr, most of our simulations remained chaotic, but several captures in stable nonsynchronous resonances also occurred. The final scenario depends on the initial conditions, and also on the tidal model. A constantQ model would prevent any capture in resonance (Goldreich & Peale 1966), while a viscoelastic model would increase the chances of capture (e.g., Makarov 2012; Correia et al. 2014). Strange attractors can also exist in the chaotic zone, which may prevent the spin from stabilizing (e.g., Batygin & Morbidelli 2011). However, the global dynamics described in Sect. 3.2 is very robust and reliable, since it does not depend on dissipative forces (Fig. 2).
5. Discussion
This work was motivated by the recent observations on the rotation of Hydra and Nix (Showalter & Hamilton 2015), which appear to be chaotic. Our model confirms that chaotic rotation is a likely scenario for both satellites, but stable spinorbit coupling could also be possible, in particular for Hydra. This model assumes a hierarchical threebody system with coplanar orbits. However, we also integrated the spins of Hydra and Nix using the ephemerides for the full system provided by Brozović et al. (2015) and SPICE routines (Acton 1996), and we have found no evidence of any substantial differences.
In our study, we conclude that stable spinorbit coupling is also a plausible scenario for near circular circumbinary bodies with small α and/or small σ (Fig. 2). Equilibrium rotation occurs for , with ν = n_{b} − n and k ∈Z. The largest amplitude nonsynchronous resonance corresponds to the subsynchronous resonance at n − ν = 2n − n_{b}. Bodies captured in this resonance present retrograde rotation if n_{b}/n> 2. This condition is verified for the small satellites of the PlutoCharon system (Table 1), and likely for any circumbinary system, since large orbital
instabilities are expected for period ratios below the 2/1 mean motion resonance. Therefore, bodies trapped in a nonsynchronous resonance are also likely to present retrograde rotation.
Lately, many planets have been detected around binary stars (e.g., Welsh et al. 2012). So far, all these planets are gaseous giants, for which the axial asymmetry is very low, for instance (B − A) /C ~ 10^{7} for Jupiter (Jacobson 2001). Therefore, although many of these planets are close enough to their stars to undergo tidal dissipation, spinorbit coupling is very unlikely. However, for smaller mass Earthlike circumbinary planets these states are possible since (B − A) /C ~ 10^{5} (e.g., Yoder 1995). In particular, nonsynchronous rotation is possible, which is an important point to take into account in future habitability studies (e.g., Selsis et al. 2007). Unlike the PlutoCharon system, circumbinary exoplanets usually present eccentric orbits. As a consequence, the number of spinorbit resonances drastically increases (Eq. (B.3)). Overlap of the different contributions is then easier, so chaotic rotation can also be more likely in this case.
Appendix A: Equations of motion
We let () be a noninertial frame attached to the body’s principal inertial axes, with inertia tensor I = diag(A,B,C). In this frame, ω = (ω_{a},ω_{b},ω_{c}) and the angular momentum L = (Aω_{a},Bω_{b},Cω_{c}). The equations of motion are given by (A.1)where T_{i} is the gravitational torque on the body’s figure. For second order interactions, we have (e.g., Goldstein 1950) (A.2)where r_{i} has coordinates (x_{i},y_{i},z_{i}) in the body’s frame. Thus, projecting Eq. (A.1) over each axis () gives To solve these equations, a set of generalized coordinates to specify the orientation of the satellite must be chosen. We adopt the modified Euler angles (θ,ϕ,ψ) as defined in Wisdom et al. (1984). Starting the ĉ axis coincident with the normal to the orbital plane, and â along an inertial direction, we first rotate the body about the ĉ axis by an angle θ, then we rotate about â by an angle ϕ, and finally we rotate about the by an angle ψ. Then, and The full equations of motion for the spin are thus described by the set of variables (ω_{a},ω_{b},ω_{c},θ,ϕ,ψ), whose derivatives are given by Eqs. (A.3)−(A.5) together with (A.9)−(A.11). The angle ε between the ĉ axis and the normal to the orbit, usually called obliquity, can be obtained from (A.12)For the dissipative tidal torque, we adopt a constant timelag linear model, whose contribution to the spin (Eq. (A.1)) is given by (e.g., Mignard 1979) (A.13)
where k_{2} is the Love number, Δt is the time lag, Q^{1} ≡ nΔt is the dissipation factor, and R is the radius of the rotating body. Also note that is the derivative of r_{i} in an inertial reference frame.
Appendix B: Series with eccentricity
Equation (1) for the rotational motion can be expanded in power series of the eccentricities and semimajor axis ratio α as given by the general expression (6) (B.1)with (B.2)When we truncate the series to the first order in the eccentricities and to the second order in α, we get (B.3)
Acknowledgments
We thank Doug Hamilton and Mark Showalter for helpful suggestions. We acknowledge support from the “conseil scientifique” of the Observatory of Paris and CIDMA strategic project UID/MAT/04106/2013.
References
 Acton, C. H. 1996, Planet. Space Sci., 44, 65 [NASA ADS] [CrossRef] [Google Scholar]
 Batygin, K., & Morbidelli, A. 2011, Cel. Mech. Dyn. Astron., 111, 219 [NASA ADS] [CrossRef] [Google Scholar]
 Brozović, M., Showalter, M. R., Jacobson, R. A., & Buie, M. W. 2015, Icarus, 246, 317 [NASA ADS] [CrossRef] [Google Scholar]
 Cheng, W. H., Lee, M. H., & Peale, S. J. 2014, Icarus, 233, 242 [NASA ADS] [CrossRef] [Google Scholar]
 Chirikov, B. V. 1979, Phys. Rep., 52, 263 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Christy, J. W., & Harrington, R. S. 1978, AJ, 83, 1005 [NASA ADS] [CrossRef] [Google Scholar]
 Colombo, G. 1965, Nature, 208, 575 [NASA ADS] [CrossRef] [Google Scholar]
 Correia, A. C. M. 2009, ApJ, 704, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Correia, A. C. M., & Laskar, J. 2009, Icarus, 201, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Correia, A. C. M., & Laskar, J. 2010, Icarus, 205, 338 [NASA ADS] [CrossRef] [Google Scholar]
 Correia, A. C. M., & Robutel, P. 2013, ApJ, 779, 20 [NASA ADS] [CrossRef] [Google Scholar]
 Correia, A. C. M., Udry, S., Mayor, M., et al. 2005, A&A, 440, 751 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Correia, A. C. M., Boué, G., Laskar, J., & Rodríguez, A. 2014, A&A, 571, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dones, L., & Tremaine, S. 1993, Icarus, 103, 67 [NASA ADS] [CrossRef] [Google Scholar]
 Doyle, L. R., Carter, J. A., Fabrycky, D. C., et al. 2011, Science, 333, 1602 [NASA ADS] [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 Farinella, P., Milani, A., Nobili, A. M., & Valsecchi, G. B. 1979, Moon Planets, 20, 415 [NASA ADS] [CrossRef] [Google Scholar]
 Goldreich, P., & Peale, S. 1966, AJ, 71, 425 [NASA ADS] [CrossRef] [Google Scholar]
 Goldreich, P., & Peale, S. 1967, AJ, 72, 662 [NASA ADS] [CrossRef] [Google Scholar]
 Goldstein, H. 1950, Classical mechanics (Reading: AddisonWesley) [Google Scholar]
 Hut, P. 1980, A&A, 92, 167 [NASA ADS] [Google Scholar]
 Jacobson, R. A. 2001, BAAS, 33, 1039 [NASA ADS] [Google Scholar]
 Laskar, J. 1990, Icarus, 88, 266 [NASA ADS] [CrossRef] [Google Scholar]
 Laskar, J. 1993, Physica D Nonlinear Phenomena, 67, 257 [NASA ADS] [CrossRef] [Google Scholar]
 Laskar, J., & Boué, G. 2010, A&A, 522, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 MacDonald, G. J. F. 1964, Rev. Geophys., 2, 467 [NASA ADS] [CrossRef] [Google Scholar]
 Makarov, V. V. 2012, ApJ, 752, 73 [NASA ADS] [CrossRef] [Google Scholar]
 Mignard, F. 1979, Moon Planets, 20, 30 [NASA ADS] [CrossRef] [Google Scholar]
 Morbidelli, A. 2002, Modern celestial mechanics: aspects of solar system dynamics (London: Taylor & Francis) [Google Scholar]
 Robutel, P., Rambaux, N., & El Moutamid, M. 2012, Cel. Mech. Dyn. Astron., 113, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Selsis, F., Kasting, J., Levrard, B., et al. 2007, A&A, 476, 1373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Showalter, M. R., & Hamilton, D. P. 2015, Nature, 522, 45 [NASA ADS] [CrossRef] [Google Scholar]
 Tholen, D. J., & Buie, M. W. 1997, Icarus, 125, 245 [NASA ADS] [CrossRef] [Google Scholar]
 Tholen, D. J., & Tedesco, E. F. 1994, Icarus, 108, 200 [NASA ADS] [CrossRef] [Google Scholar]
 Weaver, H. A., Stern, S. A., Mutchler, M. J., et al. 2006, Nature, 439, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Welsh, W. F., Orosz, J. A., Carter, J. A., et al. 2012, Nature, 481, 475 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Wisdom, J. 1987, AJ, 94, 1350 [NASA ADS] [CrossRef] [Google Scholar]
 Wisdom, J., Peale, S. J., & Mignard, F. 1984, Icarus, 58, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Yoder, C. F. 1995, in Global Earth Physics: A Handbook of Physical Constants (Washington D.C: American Geophysical Union), 1 [Google Scholar]
Online material
Fig. 4 Possible evolutions for the spin of Hydra. We show the rotation (in red) and the cosine of the obliquity (in blue). We numerically integrated the equations of the threebody problem together with Eqs. (A.1), (A.2), and (A.13) for the rotational motion. We adopt an initial retrograde rotation of 4.4 days with different initial values of θ, and an initial obliquity of 30°. 

Open with DEXTER 
All Tables
All Figures
Fig. 1 Reference angles and coordinates. The system is coplanar and is composed of an inner binary with masses m_{0} and m_{1}, and an outer body m that rotates with angular velocity ω. 

Open with DEXTER  
In the text 
Fig. 2 Stability analysis of the rotation rate close to the synchronization for δ = 0.1, α ∈ [ 0.25:0.5 ], and (B − A) /C ∈ [ 0:0.9 ]. We adopt three different initial obliquity values ψ = 1°, 30°, and 60°. The color index indicates the proportion of chaotic orbits inside the studied domain: from dark blue for fully regular to red for entirely chaotic. The dashed line (Eq. (11)) roughly delimits the zone where chaotic motion is possible. The dots identify two representative cases that we detail in Fig. 3. 

Open with DEXTER  
In the text 
Fig. 3 Cross sections representing the spin dynamics of Hydra (top) and Nix (bottom) at different initial obliquities ψ = 1°, 30°, and 60°. The dots show the main frequency of γ(t), denoted by ξ, for different initial values of . It is obtained by numerical integration of the threebody problem together with Eqs. (A.1) and (A.2) for the rotational motion (without dissipation). θ is initially set equal to f so the sections go through the middle of the synchronous resonance (). All the remaining variables are initially set to zero. The dots are blue for stable rotation (libration or circulation), and red for chaotic motion. The distinction between stable and unstable trajectories is based on the estimate of the second derivative of ξ with respect to the initial condition (see Laskar 1993). The plateaus correspond to the resonance crossings (ξ/ν = k/2). The red dots randomly distributed between the plateaus indicate the overlapping of the associated spinorbit resonances, which generates chaotic motion. 

Open with DEXTER  
In the text 
Fig. 4 Possible evolutions for the spin of Hydra. We show the rotation (in red) and the cosine of the obliquity (in blue). We numerically integrated the equations of the threebody problem together with Eqs. (A.1), (A.2), and (A.13) for the rotational motion. We adopt an initial retrograde rotation of 4.4 days with different initial values of θ, and an initial obliquity of 30°. 

Open with DEXTER  
In the text 