Free Access
 Issue A&A Volume 582, October 2015 A69 8 Planets and planetary systems https://doi.org/10.1051/0004-6361/201525939 09 October 2015

## 1. Introduction

Following observations of the Moon, Cassini (1693) established three empirical laws on its rotational motion. The first stated that the rotation rate and the orbital mean motion are synchronous, the second that the angle between Moon’s equator and the ecliptic is constant, and the third that the Moon’s spin axis and the normals to its orbital plane and ecliptic remain coplanar. The observed physical librations are described as departures of the rotational motion from these three equilibrium laws.

Colombo (1966) has shown that the second and third laws are independent of the first one, in the sense that even if the rotation rate is not synchronous, the second and third laws can still be satisfied since they correspond to the minimum dissipation of energy for the spin axis. For a non-synchronous Moon, only the angle between its equator and the ecliptic would change. Indeed, while the first law requires an triaxial ellipsoid to work, the two other laws only require an oblate spheroid. Moreover, Colombo (1966) generalised his theory to any satellite or planet whose nodal line on the invariant plane shifts because of perturbations, which can have a different origin, such as the oblateness of the central body, perturbations from a third body, or both. Peale (1969) further generalised second and third laws to include the effects of an axial asymmetry and rotation rates commensurable with the orbital mean motion.

In the classical approach, the Hamiltonian of a slightly aspherical body is developed in a reference frame that precesses with the orbit. If the angular momentum and the energy are approximately conserved, the precession of the spin axis relative to the coordinate system fixed in the orbital plane is determined by the intersection of a sphere and a parabolic cylinder. The spin axis is fixed relative to the precessing orbit when the energy has an extreme value. Thus, these equilibria states for the spin axis can be the end point of dissipation, usually with a tidal origin (Ward 1975), and they received the name of Cassini states.

An important hypothesis required to apply the classical approximation is that the orbit precesses about an inertial plane with uniform angular velocity and constant inclination. For small satellites, this requirement is usually met; they keep a constant inclination either to the ecliptic (like the Moon) or to the equatorial bulge of the central planet. Two planets in the solar system are also expected to occupy Cassini states, namely Mercury (Peale 2006) and Venus (Correia & Laskar 2003). The classical hypothesis still works on these cases, because none of these planets have satellites. However, for a planet with a huge satellite (such as the Earth-Moon or the Pluto-Charon systems), the classical approximations fail for several orbital configurations (Boué & Laskar 2006). It also fails if one wants to inspect the spin of a star with high-mass close-in planetary companions.

In this paper we intend to generalise the theory of Cassini states further. For simplicity we consider the secular three-body hierarchical problem with spin, since it has been shown to be integrable (Touma & Wisdom 1994; Boué & Laskar 2006). It thus allows us to compute the level curves of the Hamiltonian. In Sect. 2 we revise the secular dynamics of the three-body problem and introduce a new set of variables that allows computing Cassini states straightforwardly. In Sect. 3 we present the general conditions to obtain the Cassini states and apply it to the case of a star with two planetary companions and to the Earth-Moon system. In Sect. 4 we derive some conclusions and explain how to extend the model to n-body systems.

 Fig. 1Jacobi coordinates, where r1 is the position of m1 relative to m0 (inner orbit), and r2 the position of m2 relative to the center of mass of m0 and m1 (outer orbit). The body with mass m0 is considered an oblate ellipsoid, where s is the spin axis. Open with DEXTER

## 2. Secular dynamics

We consider a three-body hierarchical system composed of a central pair with masses m0 and m1, together with an external companion with mass m2 (Fig. 1). The body with mass m0 is considered an oblate ellipsoid with gravity field coefficient J2 and rotation rate ω = ωs, where s is also the axis of maximal inertia (gyroscopic approximation) with moment C. The rotational angular momentum is then simply given by (1)For the orbits we use Jacobi canonical coordinates, with r1 the position of m1 relative to m0 (inner orbit), and r2 the position of m2 relative to the centre of mass of m0 and m1 (outer orbit). The Keplerian orbital angular momenta of each orbit are (2)where ki is the unit vector normal to the orbital plane (i = 1,2), ai is the semi-major axis, ei is the eccentricity, β1 = m0m1/m01 and β2 = m01m2/ (m01 + m2) are the reduced masses, , , with m01 = (m0 + m1), and is the gravitational constant.

### 2.1. Hamiltonian

For simplicity, we restrict our analysis to hierarchical systems for which r1r2, where ri = | | ri | |. Similarly, we assume that the ellipsoid mean radius is Rr1. Thus, the interaction potential between the two orbits and the body oblateness can be restricted to terms in (r1/r2)2 and (R/ri)2; that is to say, we adopt a quadrupolar approximation for the Hamiltonian (e.g., Touma & Wisdom 1994).

Since we are only concerned with the secular evolution of the spin axis, we average the Hamiltonian of the quadrupolar three-body problem over the mean anomalies of the two orbits and also over the argument of perihelion of the inner orbit. For the non-constant parts, we simply get (Goldreich 1966; Touma & Wisdom 1994; Boué & Laskar 2006) (3)where (4)

### 2.2. Equations of motion

In the secular conservative problem, all quantities appearing in the Hamiltonian (Eq. (3)) are constant, except for the unit vectors s and ki, which can be related to the angular moment components (Eqs. (1) and (2)). The evolution of the system can therefore be described by the evolution of the angular momentum components, which can be obtained from the Hamiltonian through Poisson brackets (e.g., Dullin 2004; Breiter et al. 2005) (5)which gives for the unit vectors (6)and (7)with ij = 1,2. From previous expressions we also see that the total angular momentum is conserved: (8)

 Fig. 2Reference planes for the definition of the direction cosines and the precession angles. Open with DEXTER

### 2.3. Reduced problem

Following Goldreich (1966) and Boué & Laskar (2006), the equations of motion can be simplified if we consider only the relative position in space of the unit vectors s, k1, and k2, given by the direction cosines1 (Fig. 2) (9)together with the “berlingot” shaped volume (10)The equations of motion (6) and (7) are then rewritten as (11)and (12)We can get w directly from expression (10), but this last equation can be useful for determining whether w is positive or negative. In addition, we still have two remaining integrals, one from the Hamiltonian (Eq. (3)), (13)and another from the total angular momentum (Eq. (8)), (14)so Eqs. (11) reduce to an integrable problem.

### 2.4. New reduction

The set of variables (x,y,z) allow us to find an integrable solution for the spin-orbit problem, but this solution is not straightforward (see Sect. 3.4 in Boué & Laskar 2006). Therefore, we introduce here a new set of variables (u,v) that are more intuitive and natural to the Cassini states’ problem.

#### 2.4.1. Projection on the inner orbit

We let (15)and (16)where φ is the angle measured along the inner orbit from the interception with the outer orbit to the interception with the equatorial plane (Fig. 2). Thus, when φ = 0 the unit vectors normal to these planes (s, k1, k2) lie in the same plane. With this choice, x only depends on the new variables2(17)while y and z still depend on each other (18)However, replacing y above in expression (14), we get (19)which can be explicitly solved for z as (20)with (21)(22)Therefore, y and z also depend only on (x,u), hence on the new variables (u,v), as well as the Hamiltonian (Eq. (13)) (23)The corresponding equations of motion are

#### 2.4.2. Projection on the outer orbit

A similar reduction (u,v) could be obtained for the projection of the spin on the outer orbit, by defining (26)(27)where ϕ is the angle measured along the outer orbit from the interception with the inner orbit to the interception with the equatorial plane (Fig. 2). Thus, when ϕ = 0 the unit vectors normal to these planes still lie in the same plane. With this choice, the Hamiltonian (13) can also be expressed in the new variables (28)using the transformations (29)(30)(31)with (32)(33)

### 2.5. Classical approximation

For a small satellite, we assume that m0m1<m2 and Ra1a2. As a consequence, with C ~ m0R2, it follows from expressions (1)(4) that LG1G2 and α2α1. The Hamiltonian (13) can then be simplified as (34)and from expression (14) (35)Replacing z in expression (34), we obtain the Hamiltonian that is often used to study Cassini states (e.g., Colombo 1966; Ward 1975; Henrard & Murigande 1987) (36)where (37)Here, g/L is the constant precession rate of k1 about k2, while I0 = cos-1(z0) is the constant inclination between these two vectors. Colombo (1966) has shown that the previous Hamiltonian represents a family of parabolas, whose interception with the unit sphere gives the possible trajectories for the spin axis, s. In the precessing frame, one can express k1 = (0,0,1), k2 = (sinI0,0,cosI0), and s = (sinθcosφ,sinθsinφ,cosθ), thus (38)If one adopts the variables (u,v) from the very beginning, we can obtain the same Hamiltonian without introducing the precessing frame. Indeed, with the assumptions done for the classical approximation, we have Z(x,u) ≈ z0(1−Lx/G1) (Eq. (21)) and G(x,u) ≈ G1(1 + Lx/G1) (Eq. (22)). We then rewrite (20) as (39)Replacing z above in the general Hamiltonian (34), we directly obtain expression (38) again. The huge advantage of the new description presented in Sect. 2.4 is that it still holds in more general situations when L ~ G1 ~ G2 and α1 ~ α2, for which the precession rate of k1 about k2 and the angle between these two vectors are no longer constant.

## 3. Cassini states

Cassini states correspond to equilibria of the spin axis. They can thus be given by the extrema of the Hamiltonian (Eq. (23)): (40)Since H0 = H0(x,u), we have for the derivative with respect to v(41)We then conclude that v = 0 is always a possible equilibrium solution (equivalent to φ = 0), where the unit vectors s, k1, and k2 remain coplanar. Replacing v = 0 in the derivative with respect to u (Eq. (40)) provides a generalised implicit condition for coplanar Cassini states: (42)with uc = sinθc, and (Eqs. (17)(20)) (43)The roots of (42) can be found in the interval uc ∈ [−1,1] using numerical methods or simply by plotting its graph.

Alternatively, coplanar Cassini states can be obtained as stationary solutions for the equations of motion (), for which v = 0. Therefore, they can be simply obtained by setting v = 0 and (Eq. (25)) or, equivalently, for w = 0 and = 0 (Eq. (12)). For a given value of the total angular momentum of the system, K0, the Cassini states then verify the condition (44)Since v = v (Eq. (27)), we can obtain an equivalent condition for the coplanar Cassini states in terms of (45)or (46)with (Eqs. (29)(33)) (47)

### 3.1. Classical states

In the classical approximation, the Hamiltonian can be simplified by (Eq. (38)), therefore (48)which is equivalent to (49)This expression corresponds to the commonly used condition for finding the equilibrium points for the spin axis (e.g., Colombo 1966; Peale 1969; Ward 1975). It is usually expressed in terms of the obliquity uc = sinθc, xc = cosθc, and inclination z0 = cosI0 as (e.g., Ward & Hamilton 2004) (50)This condition could also have been obtained from expression (44), performing the same approximations as in Sect. 2.5, (51)since (Eq. (18)), and zcz0 (Eq. (35)).

### 3.2. Stars with close-in companions

Unlike satellites, the rotational angular momentum of stars is often comparable to the orbital angular momentum of close-in companions. Therefore, the classical approximation from Sect. 2.5 is not valid, and we need to apply the more general method presented in Sect. 2.4.

We first consider the case of a single star with two Jupiter-like planetary companions, for instance the HAT-P-13 system, for which L ~ G1G2 (Table 1). This system is well constrained, since data was collected combining radial velocity and transit measurements (Winn et al. 2010). The inner planet is a transiting hot Jupiter in a 2.9 day quasi-circular orbit, while the outer body has a 1.2 yr eccentric orbit and a minimum mass of about 14 Jupiter masses. The true mass and orbital inclination of the outer companion are unknown, so it is the mutual inclination between the two orbits. To enhance Cassini states, we set the initial value of I = 30° when v = 0. By modelling the Rossiter-McLaughlin effect, Winn et al. (2010) also show that the inner orbit angular momentum vector and the stellar spin vector are nearly aligned on the sky, so we take the initial θ = 0°. The total angular momentum of the system (Eq. (14)) is then (52)The rotation of the star is unknown. This parameter is important for computing the rotational angular momentum of the star (Eq. (1)), but also for estimating its oblateness through (e.g., Correia & Rodríguez 2013) . If the rotation is synchronous with the orbital period of the inner orbit, we will get a rotation period of 2.9 days. However, Winn et al. (2010) estimate the projected stellar rotation rate to be 1.66 ± 0.37 km s-1, which gives an upper limit of 48 days for the rotation period. Adopting k2 = 0.028 (Mecheri et al. 2004) and assuming a rotation period of 10 days, we have J2 ~ 10-6.

Table 1

Observed parameters for the HAT-P-13 system (Winn et al. 2010) and the Earth-Moon system (Yoder 1995).

In Fig. 3 (bottom), we show the secular trajectories for the spin projected on the inner orbit plane, obtained by plotting the level curves H0(u,v,K0) = const. (Eq. (23)), that is, without integrating the equations of motion. We observe that there are two Cassini states, uc = −0.586 and uc = 0.743, which correspond to θc ≈ −36° and θc ≈ 48°, respectively. This is a striking result, because in the classical approximation, only one final state was expected. Moreover, Cassini states correspond to the final outcome of tidal evolution, but we presently observe θ ≈ 0°. Therefore, either the present state is still precessing around a Cassini state (i.e., it is not yet damped), and we got it close to zero by chance, or the mutual inclination of the system is lower than 30°. Indeed, adopting a smaller initial inclination I = 3° in expression (52), we get θc ≈ −3.4° and θc ≈ 4.9°, which are more compatible with the observational data.

In Fig. 3 we additionally show the secular trajectories for the spin projected on the inner orbit normal k1 (middle), obtained using x = x(u,v) (Eq. (17)), and for the mutual inclination (top), obtained using z = z(u,v,K0) (Eq. (20)). The (u,v) that we use for each trajectory are those that simultaneously verify a given level curve H0(u,v,K0) = const. Therefore, we obtain the variations in the direction cosines without integrating the equations of motion. We observe that the mutual inclination undergoes some oscillations around the Cassini equilibria states, which was not possible in the classical approximation. Furthermore, the two Cassini states have quite different values for Ic ≈ 6.2° and Ic ≈ 42.2°, respectively. Here, the conserved quantity for all trajectories is only the total angular momentum K0 and no longer z, contrarily to the classical description, for which zK0/ (G1G2) = const. (Eq. (35)).

We now consider a fictitious system with the exact same parameters as in HAT-P-13 (Table 1), but where the inner planet is ten times more massive. As before, we also adopt a rotation period of 10 days and initial I = 30°, θ = 0°, and v = 0 (Eq. (52)). In Fig. 4 we show the secular trajectories for the modified HAT-P-13 system as in Fig. 3 for the standard system.

In this case we have LG1 ~ G2. The classical approximation is still inappropriate, although the level curves of the Hamiltonian recall those plotted by Ward & Hamilton (2004) for the spin of Saturn. We count three different Cassini states, uc = −0.889, uc = −0.308, and uc = 0.910, which correspond to θc ≈ −63°, θc ≈ −18°, and θc ≈ 65°, respectively. The smaller one corresponds to a hyperbolic unstable point, but the spin can be stabilised in the other two states. However, in contrast to the classical case, here the mutual inclination undergoes significant variations. Moreover, the Cassini states also present different values for the equilibrium inclination, Ic ≈ 18.8°, Ic ≈ 27.9°, and Ic ≈ 28.1°, respectively.

When we modify the mass of the inner planet m1, we change the quantities α1 and γ (Eq. (4)) and the angular momentum of the inner orbit, G1 (hence K0). Since the ratios α1/G1 and γ/G1 remain almost unchanged, different m1 values only lead to different precession rates of the spin axis s (Eq. (5)). In the classical approximation (Sect. 2.5), only the variations in the precession of the spin are significant (through the term in α1/L). Therefore, modifications in the Cassini states equilibrium points are usually studied as a function of the ratio α1/g (Eq. (50)). In Fig. 5 we show the Cassini states equilibria as a function of m1. Since g is also constant for different m1 values (Eq. (37)), this figure is equivalent to the classical maps for the ratio α1/gm1 (see, for instance, Fig. 3 in Ward & Hamilton 2004). For m1MJup, the number of Cassini states is the same as in the classical case, but for lower masses, an additional Cassini state exists for uc < 0, while the state for uc > 0 can reach values very close to 1.

 Fig. 3Secular trajectories in the HAT-P-13 system (Table 1). We show the mutual inclination (top), the stellar spin projected on the inner orbit normal (middle), and its projection on the orbital plane (bottom). Cassini states are marked with a dot. Open with DEXTER

 Fig. 4Secular trajectories in the HAT-P-13 modified system with m1 = 8.51 MJup. We show the mutual inclination (top), the stellar spin projected on the inner orbit normal (middle), and on the orbital plane (bottom). Cassini states are marked with a dot. Open with DEXTER

 Fig. 5Cassini states as a function of m1 for a HAT-P-13-like system with initial I = 30°. These equilibria are obtained by solving Eq. (42). Vertical dotted lines correspond to the configurations shown in Figs. 3 and 4. Open with DEXTER

We thus see that the classical approximation from Sect. 2.5 is unable to correctly describe the secular motion of HAT-P-13-like systems. Some trajectories for the spin may present similarities with the classical case, but they can also present very different behaviours. In particular, Cassini states cannot be given by expression (50), since z is not constant, so we do need to find the roots when solving the more general Eqs. (42) or (44). Since Cassini states correspond to the end point of tidal evolution, if we are able to estimate the obliquity of tidally evolved stars we can put some constraints on the geometry of the orbits.

### 3.3. The Earth-Moon system

For planets possessing a massive satellite, as in the Earth-Moon system (or in the Pluto-Charon system), the rotational angular momentum of the planet can be comparable to the orbital angular momentum of the satellite. However, in the case of a planet, the outer orbit’s angular momentum is much larger than the inner orbit’s, so we have L ~ G1G2. As a consequence, the outer orbit is almost an inertial frame, and it is more common to express the spin of the planet with respect to this plane. We thus adopt the set of variables (u,v) from Sect. 2.4.2 here.

We first consider the present orbital configuration of the Earth-Moon system with a1 ≈ 60 R and I = 5.145° (Table 1). For the spin of the Earth we adopt the present rotation period with obliquity ε = 23.44° and J2 ≈ 10-3 (Yoder 1995). The total angular momentum of the system (Eq. (14)) is obtained setting v = 0 (which gives θ = 18.3°): (53)In Fig. 6 (right), we show the secular trajectories for the spin projected on the plane of the ecliptic (outer orbit). We plot the level curves H0(u,v,K0) = const. (Eq. (28)), together with the relative positions of the unit vectors s, k1, and k2 (Eq. (9)). These trajectories are obtained without integrating the equations of motion, unlike in previous studies on the Earth-Moon system. We observe that the spin axis of the Earth describes an almost perfect circle at constant obliquity around the Cassini state εc = 22′′, which is very close to the ecliptic pole. The mutual inclination between the orbit of the Moon and the ecliptic is also nearly constant, since the Laplacian plane of the Earth-Moon system almost coincides with the ecliptic (e.g., Tremaine et al. 2009).

 Fig. 6Secular trajectories for the Earth-Moon system for different values of the semi-major axis, a1 = 4 R⊕ (left), a1 = 11 R⊕ (middle) and a1 = 60 R⊕ (right). We show the angles between the unit vectors (I, θ and ε), together with the projection of the Earth’s spin axis on the ecliptic (u∗,v∗). Cassini states are marked with a dot. To be compared with Figs. 1 and 2 in Touma & Wisdom (1994) and Figs. 6 and 7 in Boué & Laskar (2006). Open with DEXTER

Actually, at the present Earth-Moon distance a1 ≈ 60 R, we have LG1, so we could have used the classical approximation (Sect. 2.5). However, the Moon probably formed very close to the Earth: a Mars-sized body hit the nearly formed proto-Earth, blasting material into orbit around it, which accreted to form the Moon (e.g., Canup & Asphaug 2001). In Fig. 6 (left), we show the trajectories of the spin short-time after this impact, more precisely for a1 = 4 R, which gives LG1 (the remaining parameters are those in Table 1). During the early stages of the system, the Laplacian plane is close to the equatorial plane of the Earth, so the angle between the inner orbit and the equator, θ, is nearly constant (e.g., Goldreich 1965). As a consequence, the obliquity ε and the mutual inclination I are no longer constant, except for the Cassini states, and , which correspond to εc ≈ −16° and εc ≈ 18°, respectively. The precession of the spin axis projected on the ecliptic is not circular, and in the second case, it does not even encircle the ecliptic pole. In the first Cassini state, the orbit of the satellite almost coincides with the Earth’s equator, θc ≈ 0.08°, while in the other it keeps a significant tilt with respect to the equator, θc ≈ 32°.

As a result of tidal dissipation, the Moon evolved from the primordial close-in orbit into the present one (e.g., Touma & Wisdom 1994). On its way, there is a critical distance around 11 R where the Laplacian plane progressively shifts from the equator to the ecliptic. In Fig. 6 (middle), we show the trajectories of the spin for a1 = 11 R. At this evolutionary stage, none of the direction cosines (Eq. (9)) are constant unless the spin is trapped in the Cassini states, and , which correspond to εc ≈ −20° and εc ≈ 14°, respectively. The precession of the spin axis projected on the ecliptic is similar to the previous case with a1 = 4 R, except that the circulation area around εc ≈ 14° is larger. Indeed, as the Moon moves away from the Earth, the Cassini state with approaches zero, and the area around it grows, while the Cassini state with is shifted to the left until it disappears (Fig. 6, right).

As for the HAT-P-13 system from previous section, when the Moon is closer to the Earth, the rotational and the orbital angular momenta have similar magnitudes. Therefore, the classical approximations to find the Cassini states do not work, and we need to solve the more general Eqs. (45) or (46) to determine them. Moreover, the general method presented here also provides an easy way of determining the limits for the variations in the relative positions of the unit vectors s, k1, and k2. This information is very useful for climatic models when we inspect the habitability of new worlds, and it does not require performing numerical simulations.

Finally, since this method only relies on the total angular momentum of the system (Eq. (14)), it can also be used to quickly determine constraints for the past history of our planet. For simplicity, the evolution shown in Fig. 6 only accounts for variations in the semi-major axis. A more rigorous analysis requires that the conserved quantity in expression (14) is K instead of K0, and that the amount of angular momentum lost in G1 is transferred to L by increasing the rotation rate of the Earth.

## 4. Conclusions

In this paper we have presented a simple method for determining Cassini states and the trajectories of the spin in the secular three-body problem. This method is more general than previous approaches because it does not require the rotational angular momentum to be much smaller than the orbital one. Therefore, it can be used to study stars with close-in companions or planets with heavy satellites, for which the precession rate and the mutual inclination of the orbits are not constant. Our method only depends on the geometry of the Hamiltonian and thus does not require an integration of the equations of motion.

We have shown that previously unknown Cassini states may exist at high obliquities. As a consequence, the spin of tidally evolved stars with close-in companions can be significantly misaligned, provided that the orbits of the companions are also not coplanar. Thus, if we are able to determine the obliquity of these stars, we can put some constraints on the relative inclination between the two orbits. Planets with large close-in satellites can also present unexpected equilibrium configurations.

Our method relies on the conservation of the total angular momentum of the system. Thus, it can be useful to easily track the dynamical evolution of the system when it is subject to dissipation. Indeed, if we assume adiabatic evolution, the spin axis will travel across the constant energy levels of the Hamiltonian H(u,v,K0) towards stable Cassini states. We can thus predict the final configuration without performing numerical simulations.

Our model has some limitations. When m0 ~ m1 and r1 ~ R (for instance, a system of close binary stars), the spin and the J2 of the companion mass m1 should also be taken into account in our analysis. The problem is no longer integrable, but close solutions can still be found. The Hamiltonian (Eq. (3)) was obtained in the frame of the quadrupolar non-restricted problem, i.e., we assumed r1r2. Therefore, when octupole or resonant perturbations become important (close semi-major axis and/or very eccentric orbits), our method also only gives approximate results. Finally, since we averaged the Hamiltonian over the argument of the perihelion of the inner orbit, our method is not valid when the perihelion is in libration. This can be the case for extremely high values of the mutual inclination, where exchanges between e1 and I may occur (e.g., Farago & Laskar 2010).

In our model, we considered only the three-body problem. Although this represents many situations observed in nature, planetary systems usually contain more bodies. The n-body problem with spin is very complex, and it has a large number of degrees of freedom (e.g., Boué & Fabrycky 2014). However, we can generalise our method to those situations in the same way as done for the classical studies on Cassini states (e.g., Ward & Hamilton 2004; Peale 2006). The secular perturbations of a n-body system on the inner orbit can be decomposed in quasi-periodic series of the secular forcing frequencies in the system (e.g., Laskar 1988). For an isolated term associated

with a specific frequency, the Hamiltonian can be simplified and made integrable as shown here. Cassini states can therefore be found for the dominating perturbations in the Hamiltonian.

1

The variables x and y here are switched with respect to those in Boué & Laskar (2006), because they define the index 2 for the inner orbit and 1 for the outer orbit. They also differ from the notations in Goldreich (1966), where x = s·k2, y = k1·k2, and z = s·k1.

2

We consider here that x> 0, but this method is still valid for x< 0 adopting .

## Acknowledgments

A.C. thanks G. Boué and J. Laskar for discussions, and acknowledges support from CIDMA strategic project UID/MAT/04106/2013.

## All Tables

Table 1

Observed parameters for the HAT-P-13 system (Winn et al. 2010) and the Earth-Moon system (Yoder 1995).

## All Figures

 Fig. 1Jacobi coordinates, where r1 is the position of m1 relative to m0 (inner orbit), and r2 the position of m2 relative to the center of mass of m0 and m1 (outer orbit). The body with mass m0 is considered an oblate ellipsoid, where s is the spin axis. Open with DEXTER In the text
 Fig. 2Reference planes for the definition of the direction cosines and the precession angles. Open with DEXTER In the text
 Fig. 3Secular trajectories in the HAT-P-13 system (Table 1). We show the mutual inclination (top), the stellar spin projected on the inner orbit normal (middle), and its projection on the orbital plane (bottom). Cassini states are marked with a dot. Open with DEXTER In the text
 Fig. 4Secular trajectories in the HAT-P-13 modified system with m1 = 8.51 MJup. We show the mutual inclination (top), the stellar spin projected on the inner orbit normal (middle), and on the orbital plane (bottom). Cassini states are marked with a dot. Open with DEXTER In the text
 Fig. 5Cassini states as a function of m1 for a HAT-P-13-like system with initial I = 30°. These equilibria are obtained by solving Eq. (42). Vertical dotted lines correspond to the configurations shown in Figs. 3 and 4. Open with DEXTER In the text
 Fig. 6Secular trajectories for the Earth-Moon system for different values of the semi-major axis, a1 = 4 R⊕ (left), a1 = 11 R⊕ (middle) and a1 = 60 R⊕ (right). We show the angles between the unit vectors (I, θ and ε), together with the projection of the Earth’s spin axis on the ecliptic (u∗,v∗). Cassini states are marked with a dot. To be compared with Figs. 1 and 2 in Touma & Wisdom (1994) and Figs. 6 and 7 in Boué & Laskar (2006). Open with DEXTER In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.