Issue 
A&A
Volume 620, December 2018



Article Number  A90  
Number of page(s)  13  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201833735  
Published online  30 November 2018 
Modeling the evection resonance for Trojan satellites: application to the Saturn system
^{1}
Universidad Nacional de Córdoba, Observatorio Astronómico, IATE,
Laprida 854, 5000 Córdoba, Argentina
email: cristian@oac.unc.edu.ar
^{2}
Observatório Nacional,
Rio de Janeiro,
20921400, RJ, Brazil
Received:
27
June
2018
Accepted:
6
September
2018
Context. The stability of satellites in the solar system is affected by the socalled evection resonance. The moons of Saturn, in particular, exhibit a complex dynamical architecture in which coorbital configurations occur, especially close to the planet where this resonance is present.
Aims. We address the dynamics of the evection resonance, with particular focus on the Saturn system, and compare the known behavior of the resonance for a single moon with that of a pair of moons in coorbital Trojan configuration.
Methods. We developed an analytic expansion of the averaged Hamiltonian of a Trojan pair of bodies, including the perturbation from a distant massive body. The analysis of the corresponding equilibrium points was restricted to the asymmetric apsidal corotation solution of the coorbital dynamics. We also performed numerical Nbody simulations to construct dynamical maps of the stability of the evection resonance in the Saturn system, and to study the effects of this resonance under the migration of Trojan moons caused by tidal dissipation.
Results. The structure of the phase space of the evection resonance for Trojan satellites is similar to that of a single satellite, differing in that the libration centers are displaced from their standard positions by an angle that depends on the periastron difference ϖ_{2} −ϖ_{1} and on the mass ratio m_{2}∕m_{1} of the Trojan pair. In the Saturn system, the inner evection resonance, located at ~8 R_{S}, may capture a pair of Trojan moons by migration; the stability of the captured system depends on the assumed values of the dissipation factor Q of the moons. On the other hand, the outer evection resonance, located at >0.4 R_{Hill}, cannot exist at all for Trojan moons, because Trojan configurations are strongly unstable at distances from Saturn longer than ~0.15 R_{Hill}.
Conclusions. The interaction with the inner evection resonance may have been relevant during the early evolution of the Saturn moons Tethys, Dione, and Rhea. In particular, Rhea may have had Trojan companions in the past that were lost when it crossed the evection resonance, while Tethys and Dione may either have retained their Trojans or have never crossed the evection. This may help to constrain the dynamical processes that led to the migration of these satellites and to the evection itself.
Key words: celestial mechanics / methods: numerical / methods: analytical / planets and satellites: dynamical evolution and stability / planets and satellites: individual: Saturn
© ESO 2018
1 Introduction
In studying the dynamics of satellites in the solar system, it is important to establish whether the bodies are bound to the planets or may escape to heliocentric orbits. Therefore, the perturbation by the Sun is a key issue when the stability of these satellites is to be addressed. In the case of the Moon, the development of Hill’s lunar theory allowed identifying the socalled evection term (Brouwer & Clemence 1961) as the largest periodic correction to the mean longitude of the Moon. For any satellite, the evection term is associated with the harmonic cos(2λ_{⊙}− 2ϖ), where λ_{⊙} is the mean longitude of the Sun and ϖ is the longitude of the pericenter of the satellite. The evection resonance arises when the precession rate of ϖ equals the solar mean motion, causing the angle λ_{⊙}−ϖ to librate around an equilibrium point. In the classical evection, the pericenter precession is driven by the solar perturbation itself, but in principle, any other perturbation inducing a pericenter precession may originate an evection resonance.
In the restricted threebody problem, Henon (1969, 1970) showed that the classical evection resonance appears as a bifurcation of a family of simple periodic orbits at a value of the semimajor axis a = 0.45 R_{Hill}.
Hamilton & Krivov (1997) showed that in the case of prograde orbits, the evection resonance appears at a = 0.53 R_{Hill} and is characterized by the resonant angle λ_{⊙}−ϖ librating either around 0° or 180°. This alignment or antialignment of the satellite pericenter with the Sun direction induces cumulative perturbations that may cause the escape of the satellite (Nesvorný et al. 2003), which might be relevant for the stability of bodies that are migrating as a result of tidal evolution or gas drag.
The evection resonance plays an important role in sculpting the architecture of the satellite systems. Nesvorný et al. (2003) studied the orbital and collisional evolution of satellites and found that prograde satellite orbits with large semimajor axes are unstable because of the effect of the evection resonance. Ćuk & Gladman (2009) explored the fate of fictitious objects trapped in the lunar Trojan points after the formation of the Moon, and found that these bodies can survive the tidal migration of the Moon until they reach 38 Earth radii, where the evection resonance ejects them from the system. Ćuk et al. (2016) also studied the past evolution of the Tethys–Dione system and found that the evection resonance perturbing a pair of mediumsize moons is the most likely mechanism for triggering instability and massive collisions in the Saturn system. Finally, Spalding et al. (2016) studied the inward migration of giant exoplanets immersed in protoplanetary gas disks, and showed that dynamical interactions may naturally destroy the hypothetical moons of such planets through capture into the evection resonance. Several analytical approaches to the study of the evection resonance have been published. Yokoyama et al. (2008) analytically determined that the semimajor axis for the appearance of the evection resonance lies at a = 0.529R_{Hill}, while Frouard et al. (2010) used dynamical maps to show that this resonance constitutes the outermost region of stability for the prograde satellites of Jupiter, being surrounded by chaotic orbits. The authors concluded that the evection can harbor stable orbits at distances from Jupiter of ~0.43 au. Frouard et al. (2010) also showed that the inclusion of thirdorder terms in the expansion of the solar disturbing function in Legendre polynomial causes the evection resonance to display a nonsymmetric topology with respect to λ _{⊙}−ϖ = 0° and 180°. AndradeInes & Robutel (2018) constructed a secular theory for planetary motion around compact binaries and analyzed the occurrence of the evection resonance that is driven by the stellar companion. The study of the evection resonance using analytical models is not straightforward, however, because usually large perturbations act on the regions of interest. In this work, we study the occurrence and stability of the evection resonance in the satellite system of Saturn. This system is particularly interesting because it harbors the only known examples of coorbital configurations in satellites in the solar system. These are
 1.
the pair Janus and Epimetheus, which evolve in a relative horseshoe orbit (Murray & Dermott 1999), and may have originated from a former Trojan configuration (Gott 2005);
 2.
the Trojan system of Tethys, accompanied by the moons Telesto (at L_{4}) and Calypso (at L_{5}) (see Niederman et al. 2018, and references therein);
 3.
the Trojan system of Dione, accompanied by the moons Helene (at L_{4}) and Polydeuces (at L_{5}).
Our main goal is to provide some clues about the possible past dynamics of both single and coorbital moons. In Sect. 2 we introduce an analytical expansion of the Hamiltonian for two satellites in Trojan coorbital configuration, and we apply it to describe the basic topology of the evection resonance in such a system. In Sect. 3 we introduce our numerical analysis of the stability of the evection resonance in the Saturn system, and we apply it first to the case of a single satellite (Sect. 3.1), and then to the case of a pair of Trojan moons (Sect. 3.2). Section 4 is devoted to studying the effect of the evection resonance on satellites that migrate as a result of tidal evolution. Finally, our conclusions are presented in Sect. 5.
2 Analytical model
In this section, we set up a model to describe the evection resonance in the case of two satellites in coorbital Trojan motion. We then apply this model to the Saturn system.
2.1 Trojan pair under the influence of a distant perturber
We consider a hierarchical system of four bodies of masses m_{i} (i = 0, ..., 3), in which m_{0} is a central oblate body of radius R_{0} (planet), m_{1}, m_{2} ≪ m_{0} are two coorbital bodies around m_{0} (satellites), and m_{3} is a distant perturber (the Sun). u_{i}, w_{i} are the barycentric positions and momenta of the bodies, and we introduce the following set of canonical coordinates: (1)
where and M = m + m_{3}. The positions r _{1}, r _{2} are measured from m_{0}, while r _{3} is measured from the barycenter of the subsystem m_{0}, m_{1}, m_{2}. Their conjugated momenta are (2)
the Hamiltonian of the system reads (4)
where is the oblateness perturbation felt by m_{i}, (5)
is the coorbital disturbing function, and (6)
is the distant perturbation, with the auxiliary variable (7)
which is the position of m_{3} relative to m_{0} (Chambers 2010).
Assuming that R_{0} ≪ r_{i} (i = 1, 2), we may expand in powers of R_{0}∕r_{i} using Legendre polynomials. Up to the fourth order, we obtain (8)
where J_{k} are the zonal harmonic coefficients of m_{0}, and φ_{i} is the latitude of m_{i} over the m_{0} equator (Murray & Dermott 1999).
Similarly, assuming that r_{i} ≪r_{3} (i = 1, 2), we may expand in powers of r_{i}∕r_{3}. Taking into account Eq. (7) and introducing the small quantities (9)
we obtain up to the third order^{1} (10)
where γ_{i} is the angle between r_{3} and r _{i} and ϕ is the angle between r_{1} and r _{2}. The terms arise from the motion of the coorbital pair. When either m_{1} or m_{2} is zero, then and Eq. (10) reduces to the classical expansion of the lunar theory (Brouwer & Clemence 1961).
No expansion is made of the coorbital perturbation, since r_{1} ~ r_{2}, but taking into account that the osculating velocity of m_{i} is (15)
where ϕ′ is the angle between and .
Equations (8), (10), and (16) can be further expanded in terms of the orbital elements of the Keplerian motion. We restrict the motion to coplanar orbits and let the equator of m_{0} coincide with the reference plane (φ_{i} = 0). We introduce the set of canonical elements: (17)
where a_{i}, e_{i}, λ_{i}, and ϖ_{i} are the semimajor axis, eccentricity, mean longitude, and pericenter longitude, respectively, and we expand the Hamiltonian in powers of the eccentricities, around e_{i} = 0, up to .
The coorbital motion is characterized by the libration of the resonant angle θ _{1}, while the synodic angle θ_{2} is a fast angle that can be eliminated by an averaging: (18)
thus Z_{2} becomes an integral of motion.
In this way, the averaged coorbital perturbation, , reduces to an expansion similar to that obtained by Robutel & Pousse (2013). This expansion has a singularity when a_{1} = a_{2} and θ _{1} = 0, even if e_{i} ≠ 0, and only contains terms of degree 1, e^{2}, e^{4}, since the odddegree terms are eliminated by the average. The averaged expansion of the oblateness perturbation, , is straightforward and also contains terms of degree 1, e^{2}, e^{4}. Only the averaged expansion of the distant perturbation, , contains both even and odddegree terms in e. In particular, only contains terms of degree 1, e^{2}, while contains terms of degree 1, e^{2}, e^{4}. Both and contain terms of degree e, e^{3} only.
After rearranging terms, the averaged Hamiltonian may be written as follows: (19)
where is a one degree of freedom Hamiltonian that contains all the terms that do not depend on the eccentricities, is a perturbation that contains the terms that depend solely on e_{1}, e_{2}, and is a perturbation that contains all the termsdepending on e_{3}. Elimination of two degrees of freedom is immediate by setting e_{3} = 0. In that case, vanishes and the resulting Hamiltoniandoes not depend on W_{3}, ψ_{3}, nor on θ_{3}, thus Z_{3} becomes an integral of motion. This integral reflects the perturbation of the coorbital system onto the distant perturber, which in this approach has a fixed circular orbit with a varying semimajor axis.
The literal expansion of in noncanonical elements contains more than 50 terms, but we can restrict it to those of the lowest degree and to the most relevant terms. For the distant perturbation, we restrict the terms arising from the coorbital motion, and , to degree e^{2}, and keep the main terms, and , up to degree e^{3}. This is a reasonable assumption when we take into account that for most satellite systems. Therefore, the term is smaller than the term , while the term is on the same order as . For the coorbital perturbation, , we follow Robutel & Pousse (2013) and keep only terms of degree e^{2}. Finally, for the oblateness perturbation, we assume J_{4} = 0. The averaged Hamiltonian then reduces to the following relevant terms: (20)
where A_{ki}, B_{k}, C_{k}, and D_{i} are functions of θ_{1}, Z_{1}, and Z_{2}. Expressions for these functions in terms of noncanonical elements are given in the appendix.
We note that when m_{3} = 0, all the A_{ki} and C_{k} functions are zero and theexpansion reduces to the same expansion as reported by Robutel & Pousse (2013), with the addition of the oblateness terms D_{i}. On the other hand, when m_{2} = 0, all the B_{k} and C_{k} functions are zero, and the expansion reduces to the classical expansion of the lunar theory, again with the addition of the oblateness terms (Frouard et al. 2010).
For e_{1} = e_{2} = 0, the Hamiltonian (21)
reflects the zerothorder topology of the coorbital motion. This is shown in Fig. 1 (left panel) for the Dione–Helene system. We note that at variance with the model of Robutel & Pousse (2013), the inclusion in of the zerothorder terms of the solar perturbation (A_{0i}, C_{0}) and the planet oblateness (D_{i}) shift the libration center to values lower than 60° at increasing distances from the planet, as shown in Fig. 1 (right panel).
The Hamiltonian Eq. (20) represents a system with three degrees of freedom that might be reduced in principle, provided that the degrees of freedom evolve on very different timescales. This would allow the application of the adiabatic invariance principle to eliminate degrees of freedom, as in Beaugé & Roig (2001) or Saillenfest et al. (2016, 2017). This rigorous semianalytical treatment is beyond the scope of this work. Here, we instead use a very simple zerothorder approach to provide an overview of the Hamiltonian topology. We know that in the absence of other perturbations, the stable equilibrium of the Trojan coorbital motion arises when Δ λ = λ_{2} −λ_{1} ≃ 60°, Δ ϖ= ϖ_{2} −ϖ_{1} ≃ 60°, and e_{2} ≃ e_{1} (Giuppone et al. 2010). This corresponds to the socalled asymmetric apsidal corotation resonance solution, or ACR, which is schematically represented in Fig. 2. Therefore, we may freeze two degrees of freedom in Eq. (20) by fixing θ _{1}, Z_{1} at one of the stable equilibrium points of , and also setting ψ_{1} −ψ_{2} = 60° and e_{1} = e_{2}. In this way, the functions A_{ki}, B_{k}, C_{k}, and D_{i} become constant coefficients. The resulting level curves of the remaining Hamiltonian of one degree of freedom are shown in Figs. 3–6.
It is worth recalling that in freezing the system at the ACR solution, we avoid the occurrence of other possible equilibria of the Hamiltonian with three degrees of freedom in Eq. (20). Moreover, we cannot be sure a priori that the ACR solution represents an equilibrium solution when the solar perturbation and the oblateness perturbations are taken into account. Nevertheless, the numerical experiments presented in Sects. 3.2 and 4 demonstrate that the ACR assumption is still valid in this case, and the ACR is an effective equilibrium state of the system. Exploring other equilibria would require searching for the zeroamplitude solutions given by the stationary conditions: (22)
where J_{i} are the actions (Z_{1}, W_{1}, W_{2}) and Θ_{i} are the angles (θ_{1}, ψ_{1}, ψ_{2}). These solutions represent the extrema of Eq. (20) in a sixdimensional space and could be assessed using numerical (e.g., Giuppone & Leiva 2016), semianalytical (e.g., Giuppone et al. 2010), or analytical (e.g., Robutel & Pousse 2013) techniques. Such an analysis, however, is beyond the scope of this work.
Fig. 1 Left panel: phase space of the Hamiltonian for the Dione–Helene system. Right panel: shift of the trailing libration center (dashed line) in our model with respect to the simple coorbital motion (full line). 

Open with DEXTER 
Fig. 2 Schematic orientation of the orbits for the Trojan satellites in the ACR condition, assuming λ _{3} = 0° and ψ _{2} = 0°. 

Open with DEXTER 
Fig. 3 Topology of the inner evection resonance for different mass ratios. Panels a and b approximately correspond to the Dione–Helene pair. In all cases, a_{1} = a_{2} = 0.00326 au (488 000 km) and Δϖ = ϖ_{2} −ϖ_{1} = 60°. 

Open with DEXTER 
Fig. 4 Topology of the outer evection resonance for m_{2} ≪ m_{1}. Panels a andb: behavior of m_{1} at two different distances from the planet (in all cases a_{1} = a_{2}). In panel a, two libration islands at ψ_{1} = 0° and 180° are connected by an unstable point that is shifted to positive values of e_{1} cos ψ _{1}. In panel b only a kinematic libration around ψ_{1} = 180° is observed. Panel c: behavior of m_{2} at the same distance from the planet as in panel a. In this case, the topology is the same, but the libration islands are rotated by −Δϖ. In all cases, Δϖ = 60°. 

Open with DEXTER 
2.2 Application to the Saturn system
We applied our model to a physical system consisting of Saturn with m_{0} = 2.858 × 10^{−4} M_{⊙}, R_{0} = 60 268 km, J_{2} = 1.6298 × 10^{−2}, J_{4} = 0, a distance to the Sun a_{3} = 9.537 au, and assuming e_{3} = 0 (we recall that the eccentricity of Saturn is 0.054). We considered satellite masses on the order of those of Tethys and Dione or lower (≤ 6–10 × 10^{20} kg), and varied the mass ratio m_{2}∕m_{1} to describe the dynamics of the model. Hereafter, we refer to the angles ψ _{i} as the evectionangles.
Two different resonant regimes can be distinguished. The first occurs at very close distances from the central body and is driven by the perturbation of the planet oblateness. This regime corresponds to the inner evection resonance, and its topology displays twosymmetric libration islands that resemble a secondorder Andoyer Hamiltonian. For low mass ratios m_{2} ∕m_{1}, the libration center of the trailing Trojan, m_{1}, is located atψ_{1} = ±90° (Fig. 3a). This topology is similar to the case of a single satellite (i.e., m_{2} = 0). On the other hand, the behavior of the leading Trojan m_{2} is exactly the same, but the libration centers are shifted to ± 90° −Δϖ (Fig. 3b). When m_{2} = m_{1}, the topology is preserved, but the libration centers are shifted to ± 90° + Δϖ∕2 for the trailing and to to ± 90° −Δϖ∕2 for the leading Trojan (Figs. 3c and d). It is worth noting that the level curves in were computed for the same semimajor axes (i.e., the same values of Z_{1}), implying that the distance at which the inner evection resonance occurs is roughly independent of the mass ratio.
The second resonant regime occurs at large distances from the planet and is driven by the solar perturbation. This regime corresponds to the outer (or classical) evection resonance. The basic topology shows two asymmetric libration islands enclosed by a separatrix generated by a single saddle point close to the origin (Fig. 4a). The asymmetry is caused by the presence in the disturbing function of the odddegree terms in eccentricity, as previously shown by Frouard et al. (2010), in particular the term e_{i} cosψ_{i}, which is associated with the variation harmonic in the classical lunar theory. For low m_{2} ∕m_{1} ratios, the libration center of the trailing Trojan is located at either ψ_{1} = 180° or 0°, and this is again similar to the evection topology in the case of a single satellite. It is worth noting that for closer distances to the planet, the topology is not strictly resonant because there is no separatrix (Fig. 4b). In this case, we observe a kinematic libration around ψ _{1} = 180° forced by the odddegree terms in e. On the other hand, since the unstable equilibrium that generates the separatrix is always shifted from the origin, a resonant orbit around ψ _{1} = 180° but very close to the separatrix may eventually display a circulating angle ψ_{1}. Figure 4c shows the topology of the leading Trojan, which is the same as in Fig. 4a, but rotated by −Δϖ.
We recall that in the case of a single satellite, the precession rate caused by the oblateness of Saturn is (Roy 1978), while the precession rate caused by solar perturbation is (Innanen et al. 1997). Then, according to Li & Christou (2016), the critical semimajor axis at which one precession rate dominates is (23)
which defines the transition limit between the inner and outer evections in the Saturn system.
The locations of the libration centers of the inner and outer evection resonances, computed for the parameters of the Dione–Helene system (m_{2}∕m_{1} ≃ 2 × 10^{−5}), are shown in Fig. 5. These locations were obtained by searching for the conditions in the a, e plane that make . The outer evection shows two branches, one corresponding to ψ_{1} = 180°, which is always stable (blue line), and the other corresponding to ψ _{1} = 0°, which splits into the stable (blue line) and the unstable (red line) equilibria.
Figure 5 also shows the position of the prograde moons of Saturn. We note that all the moons are located inside of the location of the outer evection resonance, including the group of retrograde moons (not shown), whose outermost limit is indicated by the vertical dashed line. We also note that the regular moons, that is, those with e≲0.02 and i≲1°, are clustered inside of the location of the inner evection resonance (a≲0.003 au). In particular, all the moons involved in coorbital configurations (JanusEpimetheus, Tethys, and Dione) are in this group. The only exception among the regular moons is Rhea, which is located slightly outside of the inner evection, but has no coorbitals.
In the case of equal masses, the topology of the outer evection resonance is preserved, but the libration centers of both the trailing and leading Trojans are rotated by ±Δϖ∕2, respectively (e.g., Fig. 6a). At variance with the inner evection resonance, the existence existence of the libration islands in this case is very sensitive to the value of Δϖ, and they quickly disappear if the Trojan pair moves away from the ACR solution toward increasing Δ ϖ (Fig. 6b). Another major difference with respect to the inner evection is that the distance from the planet at which the resonant islands appear (more precisely, where the unstable equilibrium appears) is strongly dependent on the mass ratio, as shown in Fig. 6c (full line). This figure also shows how the libration center moves away from ψ _{1} = 0 as the mass ratio tends to 1 (dashed line).
Fig. 5 Location of the inner and outer evections (full lines) for the masses and eccentricities of the Dione–Helene system. The stable branches (related to the libration centers) are shown in blue, and the unstable branch (related to the saddle point of the outer evection) is plotted in red. The full dots indicate the current position of the prograde satellites of Saturn, and the vertical dashed line indicates the outermost limit of the irregular satellites. The letters specify the locations of Dione (D), Tethys (T), and Rhea (R), respectively. For the outer evection, the stable branch at ψ _{1} = 180° and e_{1} ≲ 0.14 is related to a kinematic libration (Fig. 4b). 

Open with DEXTER 
3 Numerical study
In this section we present a numerical survey of the stability of the Saturn satellites, with particular focus on the occurrence of the evection resonance in both the cases of a single satellite and a Trojan pair. Our numerical model takes into account the perturbations from the Sun as well as from the oblateness of Saturn by including the J_{2} zonal harmonic. The perturbations from other satellites in the Saturn system are not taken into account. The satellite orbits are described by the orbital elements a, e, λ, and ϖ. Hereafter, we maintain the subscripts 1,2 for the Trojan pair (trailing and leading, respectively), and we use the subscripts ⊙ and S for the Sun and Saturn, respectively.
The Newton equations of motion are numerically solved using a Burlisch–Stoer integrator with adaptive step size, which was modified to independently monitor the error in each variable; this imposes a relative precision better than 10^{−13}. We stopped the integrations when the mutual distance between Saturn and the satellite was smaller than the sum of their radii. The total integration time was set to encompass several periods of the orbital secular variations.
The stability was evaluated through the use of some chaos and resonance indicators, and we explored different grids of initial conditions in the space of orbital elements. In particular, for each initial condition of the grid, we computed the value of the indicator for the mean exponential growth of nearby orbits (MEGNO), , i = 1, 2 (Cincotta & Simó 2000; Maffione et al. 2011). This indicator is very efficient to quickly distinguish regular () from chaotic orbits (), but it does not provide a detailed representation of the structure of a resonance.
To better address the structure of the resonances, we also computed the amplitude of maximum variation of the orbital eccentricity of the satellites during the integrations, that is, , i = 1, 2. This indicator has proven to be an extremely useful tool for mapping the resonant structure in Nbody problems (e.g., Ramos et al. 2015). We recall that δe_{i} is not a direct measure of chaotic motion, but abrupt changes in δe_{i} are often tracers for the presence of resonances. Therefore, regions with large variations in δ e_{i} are more sensitive to perturbations and are very likely chaotic. Frouard et al. (2010) used a similar indicator in their study of the evection resonance in Jupiter satellites, defined as , where the total integration time, τ, was divided into two consecutive samples τ_{a} and τ_{b}. Nevertheless, we realized that this indicator is quite sensitive to the value of τ and may artificially create some structures that are not related with actual resonances.
A third indicator that we used is the amplitude of oscillation of the evection angle, defined as . This indicator has a minimum at the center of the resonance, or more generally, at the condition .
Fig. 6 Behavior of the outer evection resonance for equal masses of the Trojan pair. Panel a: topology of m_{1} at a distance of a_{1} = a_{2} = 0.39 au (5.88 × 10^{7} km), and assuming a periastron difference Δϖ = 60° (i.e., atthe ACR solution). Panel b: same as panel a, but for a periastron difference Δ ϖ = 70°. Panel c: distance from the planet at which the evection resonance appears (full line), and the angle of the resonance center (dashed line) as a function of the mass ratio. 

Open with DEXTER 
Fig. 7 Dynamical maps in the (a, e)plane for the regions of the inner evection (top row) and the outer evection (bottom row). The MEGNO, δ e, and δ ψ indicators are displayed in the left, middle, and right columns, respectively. Angles are in degrees. White areas correspond to unstable orbits. The gray areas in the δ ψ maps indicate nonlibrating orbits. The grids correspond to osculating orbital elements with initial ψ = 90° for the inner evection, and initial ψ = 0° for the outer evection. To allow comparison with previous works, the upper horizontal axes give the distance to the planet in convenient units: planetary radii (R_{S}) for the inner evection, and Hill radii (R_{Hill}) for the outer evection. 

Open with DEXTER 
3.1 Evection for a single moon
Figure 7 shows a set of dynamical maps in the (a, e)plane for the inner (top row, ψ = 90°) and the outer (bottom row, ψ = 0°) evections. The maps display the values of the three indicators described in Sect. 3: (left column), δ e (middle column), and δψ (right column). We drop here the subindex i for simplicity.
The inner evection region is quite regular, and both δe and δψ allow identifying the exact location of the resonance. For eccentricities lower than ~ 0.6, this location coincides with the eccentricity predicted by the analytical model (cf. Fig. 5). At higher eccentricities, the resonance location moves toward slightly larger semimajor axes, a behavior that the analytical model does not reproduce because the series is truncated early. The gray areas in the map of δψ correspond to conditions where the angle ψ circulates. White areas correspond to orbits that either collide with Saturn (a ≤ R_{S}) or escape from the system (a > a_{S}) during the integration timespan.
The outer evection resonance, on the other hand, occurs inside a stable island surrounded by a thin layer of chaotic motion (red in the map of ). The map of δ ψ shows that the libration amplitude around ψ = 0° is never larger than ~60°, in good agreement with the analytical model (cf. Fig. 4a). The map of δ e also shows a stable narrow region at low eccentricities (blue) that is related to the kinematic librations induced by the forced eccentricity mode. Nevertheless, we note that the resonant island occurs closer to the planet (a ~ 0.19 au) than predicted by the analytical model (a ~ 0.26 au). This is due in part to the fact that the grids in Fig. 7 are given in terms of the osculating initial semimajor axis, while in terms of the mean semimajor axis, the island shifts to higher values, as shown in Fig. 8. Another reason might be the early truncation of the disturbing function Eq. (10) and that the eccentricity of the Sun is not considered in the model. A similar behavior has been identified by Frouard et al. (2010).
Figure 8 shows the dynamical maps of δ ψ in the (a, ψ)plane for the outer evection. In this case, we have initially e = 0.4 and λ_{⊙} = 0°, which means that varying the initial ψ is equivalent to varying the initial ϖ. We note that the libration regions around 0° or 180° occur at different osculating semimajor axes (Fig. 8, left), but this dependence disappears when the mean semimajor axis is considered (Fig. 8, right).
Since the inner evection is driven by the planet oblateness, it is expected that the planet obliquity plays some role in the behavior of the resonance. In Fig. 9 we address the effect of the obliquity of Saturn (ϵ _{S} = 26.73°) and its orbital eccentricity (e_{S} = 0.054) on the inner evection. We note that a nonzero obliquity shifts the location of the resonance to slightly smaller distances (from ~ 8 to ~ 7 R_{S}), although the structure of the equilibrium family is preserved. On the other hand, the eccentricity of the Saturn orbit (or equivalently, the Sun) does not produce any significant changes.
Fig. 8 Dynamical maps of the δψ indicator in the (a, ψ)plane for the outer evection, with initial e = 0.4. The leftpanel corresponds to a grid of osculating elements, while the right panel is the same grid in terms of mean elements. Blue regions represent the conditions for which δ ψ ≲ 30°. Gray regionscorrespond to larger libration amplitudes or circulation of ψ. White regions are unstable orbits. The resonance occurs at a mean semimajor axis ā = 0.205 au. The upper horizontal axis is the distance in terms of R_{Hill}. 

Open with DEXTER 
Fig. 9 Value of δe as a function of the distance from the planet in the region of the inner evection and for e = 0. The curve in violet corresponds to the case of an obliquity ϵ_{S} = 0° and orbital eccentricity e_{S} = 0° for Saturn. The green dots correspond to the case of ϵ_{S} = 26.73°, either with e_{S} = 0° or e_{S} = 0.054°. The minimum inside the peak of the curves marks the center of the evection resonance. 

Open with DEXTER 
Fig. 10 Dynamical map of the δψ_{1} indicator in the osculating (a, e)plane for the inner evection, setting ψ_{1} = 120° initially. Case of Trojan satellites with equal masses. The resonance occurs at the same location as in the case of a single satellite (compare with Fig. 7). The upper horizontal axis gives the distance in R_{S}. Angles are in degrees. 

Open with DEXTER 
3.2 Evection for Trojan moons
In this section, we study the evolution of two coorbital moons in Trojan configuration initially near the equilibrium of the ACR condition, that is, σ ≡Δλ = 60°, Δ ϖ= 60°, e_{2} = e_{1} (Giuppone et al. 2010; Robutel & Pousse 2013). We consider two cases: m_{2} ∕m_{2} = 1 and m_{2} ∕m_{1} = 10^{−4}. All the simulations initially have λ _{⊙} = 0°. Figure 2 shows a schematic orbital diagram when ψ_{2} = 0°.
Figure 10 shows the dynamical map of the δ ψ _{i} indicator for a system of two Trojan moons with equal masses in the region of the inner evection. The indicator corresponds to those of the trailing Trojan (i = 1), but the situation is identical for the leading Trojan. A comparison with Fig. 7 shows that the evection resonance occurs at almost the same distance from the planet as for a single moon (a ~ 0.003 au at e~ 0), implying that the location of the inner evection is independent of the mass ratio of the Trojans. This agrees well with the predictions of the analytical model (Fig. 3).
On the other hand, the outer evection in the case of Trojan moons cannot be studied because the Trojan configuration does not survive at large distances from Saturn. This is shown in Fig. 11 for moons of equal masses. The white regions correspond to the conditions where an initial Trojan ACR configuration is destroyed, and the coorbital pair does not survive as such. Stable Trojan configurations exist up to distances of ~ 0.06 au only for almost circular orbits. At distances ≲0.02 au (blue regions), the Trojan ACR configuration is stable even for high eccentricities, except at the locations of the evection resonance and its higher harmonics (; ; ). The same qualitative results are obtained for a Trojan pair with In the region where the outer evection resonance occurs for single moons (a ~ 0.2 au; see Fig. 7), Trojan configurations are not sustainable for more than 100 yr. As a consequence, a pair of moons initially inside the evection island and in a Trojan ACR configuration quickly ends up as a nonTrojan pair with only one of the components trapped in the evection resonance. This is shown in Fig. 12 for a pair of equalmass moons. The top panel shows the case where ψ _{1} librates while ψ_{2} circulates, and the bottom panel shows the opposite case. Neither of these cases show σ librating around ± 60°, so Trojan motion no longer exists. Simultaneous libration of ψ_{1}, ψ_{2} is never found in the outer evection either. Figure 13 shows the same evolutions as Fig. 12, but in the space k_{i} = e_{i} cosψ_{i}, h_{i} = e_{i} sinψ_{i}.
The libration of either ψ_{1} or ψ _{2} is related to the initial longitudes of the periastrons relative to the longitude of the Sun (e.g., Fig. 2). The behavior of the former Trojan companion that becomes capture in the evection resonance is very similar to that of a single moon. In the (a, e)plane, the resonance occurs inside an island surrounded by a layer of chaotic motion, as shown in Fig. 14. We note, however, that the island corresponding to the libration of ψ _{1} (left panel) is smaller than that corresponding to the libration of ψ_{2} (right panel), indicating that even in the case of equal masses, the configuration may not be symmetric. Neither configuration is longterm stable, and at some point, the moon where ψ _{i} circulates is ejected or hits the other moon or the planet. Nevertheless, configurations that lead to the libration of ψ _{1} appear to be more stable (up to 10^{6} yr) than those leading to the libration of ψ_{2} (up to ~10^{3} yr; see Fig. 13).
Fig. 11 Dynamical maps in the osculating (a, e)plane displaying the amplitude δe_{1} of the eccentricity variation (top panel), the amplitude of the σ libration (middle panel), and the amplitude of the Δϖ libration (bottom panel) for a system of two equalmass moons initially in a Trojan ACR configuration. Angles are in degrees. Trojan motion is unstable in the white regions in less than 5 × 10^{4} yr, and it is most stable in the blue regions. The stripes of less stable Trojan motion observed for a < 0.03 au correspond to the locations of the evection resonance and its higher harmonics. 

Open with DEXTER 
4 Tidal evolution
The dynamical evolution of the Saturn system under the effects of tidal forces is relevant for constrainingthe past evolution of the moons. In particular, the proximity of the inner evection to the planet and to the few coorbital satellites of the system raises the question about the possible effect of this resonance on a pair of Trojan moons that evolve by tidal effects.
Short and long tidal migration rates of the Saturn moons have been explored in a wide variety of scenarios by Ćuk et al. (2016). The authors concluded that either the moons are significantly younger than the planet (formed only about 100 Myr ago), or their tidal evolution must be extremely slow (implying a dissipation factor for Saturn of Q > 80 000).
They also concluded that the inner evection resonance is located in a chaotic region, possibly because several harmonics of their own resonance overlap, although we showed here that this is not the case for isolated moons or Trojan pairs (see Figs. 7 and 10). Ćuk et al. (2016) proposed a mechanism to form a disk of protomoons that could produce a favorable environment for moons to be created, which increases the possibility that some of them might harbor Trojan companions.
The dissipation factor of Saturn has traditionally been estimated to be Q ~ 18 000, which gives a very low tidal migration rate. With this value, Enceladus, Tethys, Dione, and Rhea are not expected to have migrated significantly over the age of the solar system; they were probably formed at their current locations. However, based on an analysis of astrometric positions from old photographic plates and images of the Hubble Space Telescope, Lainey et al. (2012) revealed that the Saturn moons did migrate over the last century. They estimated a dissipation factor ten times smaller, Q ~ 1 682 ± 540. Based on this result, Crida & Charnoz (2014) suggested that the inner moons, up to Rhea, could have been formed from the primordial spreading ring of Saturn, and in this scenario, it may be expected that many primordial Trojan pairs have crossed the evection resonance.
Here, we study the effect of tidal forces in the evolution of hypothetical Trojan moons by incorporating a weakfriction equilibrium tidal model into our Nbody code. The tidal forces were modeled using the closed formulas of Mignard (1979) and following the approach of Rodríguez et al. (2013) for coorbital companions. For Saturn, we assumed a Love number k_{2} = 0.39 and Q = 1 682 (Lainey et al. 2012), and for the satellites, we assumed k_{2} = 0.05 and 20 ≤ Q ≤ 200 (Ćuk et al. 2016).
Numerical integrations using the nominal values of the system parameters are impracticable. Therefore, we scaled the total timespan and the tidal parameters down by the same factor α > 1 to reduce the total integration time, as suggested by Zoppetti et al. (2018). We considered values of α low enough to keep the accelerated tidal evolution adiabatic with respect to the resonant and secular gravitational timescales. The results shown here correspond to α = 200, but we performed some other numerical test using lower values to check the reliability of the results.
We know that in the adiabatic regime, the system evolves along the family of periodic orbits, thus we ran some simulations with a single satellite. We considered a fictitious Dionelike moon that was captured in the evection resonance at a = 0.00326 au, with ψ librating around ± 90°, as predicted by the analytical model. When we increased the value of α, the libration amplitude increased and the satellite was able to cross the evection resonance without any effect on it.
Then we ran experiments considering a pair of Trojan moons with mass ratios of m_{2} ∕m_{1} = 2 × 10^{−5}, m_{2} ∕m_{1} = 1, and m_{2} ∕m_{1} = 2 × 10^{5}. We recall that the first ratio approximately corresponds to the Dione–Helene system. In all cases, the initial semimajor axis was set to a = 0.00321 au, slightly inside of the evection resonance, and the initial eccentricity was set to the current values of Dione and Helene (0.0022). Typical results are shown in Fig. 15, where the evolution of m_{1} is identify in magenta and that of m_{2} in green. The resonant capture is observed in the bottom panels, where both evection angles ψ_{1}, ψ_{2} start to librate. For very low or very high mass ratios, the most massive component is always captured around 90° or 270°, while the less massive componentis captured around ± 60° of those values. For equal masses, the components are captured around 60° (or 140°) and 120° (or 300°), respectively,as predicted by the analytical model. The evolution of the Trojan pair inside the evection resonance implies that the eccentricities start to increase, moving along the family of periodic orbits with oscillations of higher amplitude for the less massive component. We also note that the ACR equilibrium configuration is preserved during the evolution inside the evection resonance. For low mass ratios of the Trojan pair, excitation of the eccentricities contributes to destabilize the system and may lead to the loss of the less massive component, either by hitting the planet or the other component.
The results are similar over the whole range of values of Q we tested for the satellites. On the other hand, when we considered lower values of Q for Saturn (Q < 1 000), which implies abandoning the adiabatic regime, the evection resonance causes a small jump in the eccentricities of the Trojan satellites that breaks the libration of Δϖ, causing the system to depart from the ACR solution. After this, the less massive satellite is captured with with ψ _{i} librating around ± 90°, the angle Δϖ starts to circulate, and σ → 180°, leading the pair into a horseshoe regime. The system survives under this condition until a collision occurs after ~ 10^{6} yr or sooner.
Fig. 12 Evolution of the resonant angles σ (cyan lines) and ψ_{i} (red/blue lines) for a pair of moons with equal masses, initially placed in a Trojan ACR configuration (σ = 60°, Δ ϖ = 60°) in the region of the outer evection (a_{1} = a_{2} ~ 0.2 au, e_{1} = e_{2}= 0.4). Top panel: case ψ_{1} = 0° and a_{1} = a_{2} = 0.1825 au initially. Bottom panel: case ψ_{2} = 0° and a_{1} = a_{2} = 0.1893 au initially. 

Open with DEXTER 
Fig. 13 Same as Fig. 12, but in the k_{i}, h_{i} space (see text) and for the full time span of the simulation. The top panels correspond to the top panel of Fig. 12; the system survives for at least 10^{6} yr. The bottom panels correspond to the bottom panel of Fig. 12; the system survives for ~ 10^{3} yr. In both cases, λ_{⊙} = 0° initially. 

Open with DEXTER 
Fig. 14 Dynamical maps of δψ_{i} in the osculating (a, e)plane for a Trojan pair of equal masses, initially placed in the ACR configuration. Initial conditions for the angles are the same as in Fig. 13. Blue islands correspond to libration of ψ _{1} (left panel) or ψ_{2} (right panel). Gray regions correspond to circulation of ψ_{i}, and white regions correspond to unstable orbits (compare to Fig. 7). The upper horizontal axis is the distance in terms of R_{Hill}. 

Open with DEXTER 
Fig. 15 Adiabatic evolution of pair of Trojan satellites that migrate as a result of tidal effects, using Q = 100. Each column of panels corresponds to a different mass ratio m_{2}∕m_{1}: 2 × 10^{−5} (left panels),1 (middle panels), and 2 × 10^{5} (right panels). Eccentricity is shown in the top row, σ and Δ ϖ in the middle row, and ψ_{1}, ψ_{2} in the bottom row. Magenta and green identify identify the trailing and leading Trojan, respectively, except in the middle panels, where they identify σ and Δ ϖ, 

Open with DEXTER 
5 Conclusions
We have addressed the dynamics of the evection resonance in the case of coorbital motion, and focused on Trojan pairs in the ACR configuration (Δλ = 60°, Δϖ = 60°). We applied an expansion of the averaged Hamiltonian describing the motion of two bodies in Trojan orbits, perturbed by a third distant body. The expansion combines the model developed by Frouard et al. (2010) for the evection of single bodies and the model described by Robutel & Pousse (2013) for the motion of Trojan bodies. This model can be applied for Trojan satellites in the solar system as well as for circumstellar Trojan planets in binary systems. Here we consider the Saturn system. We also performed a series of numerical simulations, both to validate the predictions of the model and to better assess the actual dynamics of Trojan satellites around Saturn. Our conclusions are summarized below.
Our model addressed the two types of evection resonance that may appear: the inner evection resonance, close to the planet and related to the precession of the periastron forced by the oblateness of the central body, and the outer evection resonance, far away from the planet and related to the precession of the periastron forced by the distant perturber.
The global structure of the phase space of the evection resonances for Trojan satellites is similar to that of a single satellite, differing in that the libration centers are displaced from their standard positions (± 90° for the inner and 0°, 180° for the outer evection) by an angle that depends on the periastron difference ϖ_{2} −ϖ_{1} and on the mass ratio m_{2}∕m_{1} of the Trojan pair.
The distance from the planet at which the inner evection resonance occurs is independent of the mass ratio of the Trojanpair. The resonance only moves toward slightly smaller distances when the obliquity of the planet is taken into account.
In the Saturn system, the evolution of Trojan moons into the inner evection resonance, located at ~8 R_{S}, is possible. Capture in this resonance can be driven by adiabatic tidal migration, exciting the eccentricities of the two Trojan components. This in turn may destabilize the resonant system, leading to the loss of one of the components.
The location of the outer evection resonance depends, in principle, on the mass ratio of the Trojan pair, moving to longer distances as m_{2}∕m_{1} → 1. For the Saturn system, this would occur at distances > 0.4 R_{Hill}.
Nevertheless, the outer evection in the Saturn system cannot exist at all for Trojan moons because Trojan configurations are strongly unstable at distances from Saturn longer than ~0.15 R_{Hill}. A pair of moons forced into a Trojan configuration in the region of the outer evection always ends in a noncoorbital configuration, with one of the moons captured in a stable orbit inside the evection resonance, while the other moon becomes unstable and is eventually lost.
In the Saturn system, the moons that are closer to the current location of the inner evection are Tethys, Dione, and Rhea. Of these, only the first two have Trojan companions and are located at smaller distances from Saturn than the evection. Rhea has no Trojans andis located slightly beyond the evection. If Rhea had Trojan companions in the past, it may have lost them when it crossed the evection resonance. Thus, it would be interesting to study the possibility that the craters on the surface of Rhea were produced by collisions with former Trojan companions that were destabilized by the evection. This study might be extended even for other regular satellites such as Titan, Hyperion, and Iapetus.
The high eccentricity of Titan, Hyperion, and Iapetus might be related to the excitation caused by the crossing of the inner evection resonance. The eccentricity of Rhea, on the other hand, may have been excited by the evection and may have later been damped to the current value by tidal dissipation.
On the other hand, the inner evection resonance moved in the past from its current location as a result of changes in the shape or obliquity of Saturn and because Saturn was closer to the Sun than it is today (e.g., Nesvorný & Morbidelli 2012). Then, Tethys and Dione might have crossed the resonance inthe past and preserved their Trojan companions. This might impose constraints on the dynamical processes that led to this type of evolution.
As for the outer evection, none of the current Saturn moons (not even the irregular moons) reach the current domain of this resonance. However, the resonance may have been relevant for the early evolution of the moons, when Saturn was closer to the Sun.
Acknowledgements
The authors acknowledge useful discussions with C. Beaugé that helped us to improve the speed of the numerical integrations with the Bulirsch–Stoer method. We also wish to thank the anonymous referee for their comments. This work has been supported by the National Council of Research of Argentina (CONICET) and the Rio de Janeiro State Science Foundation (FAPERJ), in the framework of a FAPERJ/CONICET bilateral project, and by the Brazilian National Council of Research (CNPq). Computations were performed at the BlaFis cluster of the University of Aveiro, Portugal, and at the SDumont cluster of the Brazilian System of HPC (SINAPAD).
Appendix A Coefficients of
We provide here the explicit expressions of the coefficient functions of the Hamiltonian Eq. (20) in terms of noncanonical orbital elements.
References
 AndradeInes, E., & Robutel, P. 2018, Celest. Mech. Dyn. Astron., 130, 6 [NASA ADS] [CrossRef] [Google Scholar]
 Beaugé, C., & Roig, F. 2001, Icarus, 153, 391 [NASA ADS] [CrossRef] [Google Scholar]
 Brouwer, D., & Clemence, G. M. 1961, Methods of Celestial Mechanics (New York: Academic Press) [Google Scholar]
 Chambers, J. E. 2010, NBody Integrators for Planets in Binary Star Systems, ed. N. Haghighipour (Dordrecht, The Netherlands: Springer), 239 [Google Scholar]
 Cincotta, P. M., & Simó, C. 2000, A&AS, 147, 205 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Crida, A., & Charnoz, S. 2014, in Complex Planetary Systems, Proceedings of the International Astronomical Union, IAU Symp., 310, 182 [NASA ADS] [CrossRef] [Google Scholar]
 Ćuk, M., & Gladman, B. J. 2009, Icarus, 199, 237 [NASA ADS] [CrossRef] [Google Scholar]
 Ćuk, M., Dones, L., & Nesvorný, D. 2016, ApJ, 820, 97 [NASA ADS] [CrossRef] [Google Scholar]
 Frouard, J., Fouchard, M., & Vienne, A. 2010, A&A, 515, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Giuppone, C. A., & Leiva, A. M. 2016, MNRAS, 460, 966 [NASA ADS] [CrossRef] [Google Scholar]
 Giuppone, C. A., Beaugé, C., Michtchenko, T. A., & FerrazMello, S. 2010, MNRAS, 407, 390 [NASA ADS] [CrossRef] [Google Scholar]
 Gott, J., & Richard, I. 2005, Ann. N. Y. Acad. Sci., 1065, 325 [NASA ADS] [CrossRef] [Google Scholar]
 Hamilton, D. P., & Krivov, A. V. 1997, Icarus, 128, 241 [NASA ADS] [CrossRef] [Google Scholar]
 Henon, M. 1969, A&A, 1, 223 [NASA ADS] [Google Scholar]
 Henon, M. 1970, A&A, 9, 24 [NASA ADS] [Google Scholar]
 Innanen, K. A., Zheng, J. Q., Mikkola, S., & Valtonen, M. J. 1997, AJ, 113, 1915 [NASA ADS] [CrossRef] [Google Scholar]
 Lainey, V., Karatekin, Ö., Desmars, J., et al. 2012, ApJ, 752, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Li, D., & Christou, A. A. 2016, Celest. Mech. Dyn. Astron., 125, 133 [NASA ADS] [CrossRef] [Google Scholar]
 Maffione, N. P., Darriba, L. A., Cincotta, P. M., & Giordano, C. M. 2011, Celest. Mech. Dyn. Astron., 111, 285 [NASA ADS] [CrossRef] [Google Scholar]
 Mignard, F. 1979, Moon Planets, 20, 301 [NASA ADS] [CrossRef] [Google Scholar]
 Murray, C. D., & Dermott, S. F. 1999, Solar System Dynamics (Cambridge, UK: Cambridge University Press) [Google Scholar]
 Nesvorný, D., & Morbidelli, A. 2012, AJ, 144, 117 [NASA ADS] [CrossRef] [Google Scholar]
 Nesvorný, D., Alvarellos, J. L. A., Dones, L., & Levison, H. F. 2003, AJ, 126, 398 [NASA ADS] [CrossRef] [Google Scholar]
 Niederman, L., Pousse, A., & Robutel, P. 2018, ArXiv eprints [arXiv:1806.07262] [Google Scholar]
 Ramos, X. S., CorreaOtto, J. A., & Beaugé, C. 2015, Celest. Mech. Dyn. Astron., 123, 453 [NASA ADS] [CrossRef] [Google Scholar]
 Robutel, P., & Pousse, A. 2013, Celest. Mech. Dyn. Astron., 117, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Rodríguez, A., Giuppone, C. A., & Michtchenko, T. A. 2013, Celest. Mech. Dyn. Astron., 117, 59 [NASA ADS] [CrossRef] [Google Scholar]
 Roy, A. E. 1978, Orbital Motion (Bristol, UK: Institute of Physics Publishing) [Google Scholar]
 Saillenfest, M., Fouchard, M., Tommei, G., & Valsecchi, G. B. 2016, Celest. Mech. Dyn. Astron., 126, 369 [NASA ADS] [CrossRef] [Google Scholar]
 Saillenfest, M., Fouchard, M., Tommei, G., & Valsecchi, G. B. 2017, Celest. Mech. Dyn. Astron., 127, 477 [NASA ADS] [CrossRef] [Google Scholar]
 Spalding, C., Batygin, K., & Adams, F. C. 2016, ApJ, 817, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Yokoyama, T., Vieira Neto, E., Winter, O. C., Sanchez, D. M., & de Oliveira, P. I. 2008, Math. Probl. Eng., 2008, 1 [CrossRef] [Google Scholar]
 Zoppetti, F. A., Beaugé, C., & Leiva, A. M. 2018, MNRAS, 477, 5301 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Left panel: phase space of the Hamiltonian for the Dione–Helene system. Right panel: shift of the trailing libration center (dashed line) in our model with respect to the simple coorbital motion (full line). 

Open with DEXTER  
In the text 
Fig. 2 Schematic orientation of the orbits for the Trojan satellites in the ACR condition, assuming λ _{3} = 0° and ψ _{2} = 0°. 

Open with DEXTER  
In the text 
Fig. 3 Topology of the inner evection resonance for different mass ratios. Panels a and b approximately correspond to the Dione–Helene pair. In all cases, a_{1} = a_{2} = 0.00326 au (488 000 km) and Δϖ = ϖ_{2} −ϖ_{1} = 60°. 

Open with DEXTER  
In the text 
Fig. 4 Topology of the outer evection resonance for m_{2} ≪ m_{1}. Panels a andb: behavior of m_{1} at two different distances from the planet (in all cases a_{1} = a_{2}). In panel a, two libration islands at ψ_{1} = 0° and 180° are connected by an unstable point that is shifted to positive values of e_{1} cos ψ _{1}. In panel b only a kinematic libration around ψ_{1} = 180° is observed. Panel c: behavior of m_{2} at the same distance from the planet as in panel a. In this case, the topology is the same, but the libration islands are rotated by −Δϖ. In all cases, Δϖ = 60°. 

Open with DEXTER  
In the text 
Fig. 5 Location of the inner and outer evections (full lines) for the masses and eccentricities of the Dione–Helene system. The stable branches (related to the libration centers) are shown in blue, and the unstable branch (related to the saddle point of the outer evection) is plotted in red. The full dots indicate the current position of the prograde satellites of Saturn, and the vertical dashed line indicates the outermost limit of the irregular satellites. The letters specify the locations of Dione (D), Tethys (T), and Rhea (R), respectively. For the outer evection, the stable branch at ψ _{1} = 180° and e_{1} ≲ 0.14 is related to a kinematic libration (Fig. 4b). 

Open with DEXTER  
In the text 
Fig. 6 Behavior of the outer evection resonance for equal masses of the Trojan pair. Panel a: topology of m_{1} at a distance of a_{1} = a_{2} = 0.39 au (5.88 × 10^{7} km), and assuming a periastron difference Δϖ = 60° (i.e., atthe ACR solution). Panel b: same as panel a, but for a periastron difference Δ ϖ = 70°. Panel c: distance from the planet at which the evection resonance appears (full line), and the angle of the resonance center (dashed line) as a function of the mass ratio. 

Open with DEXTER  
In the text 
Fig. 7 Dynamical maps in the (a, e)plane for the regions of the inner evection (top row) and the outer evection (bottom row). The MEGNO, δ e, and δ ψ indicators are displayed in the left, middle, and right columns, respectively. Angles are in degrees. White areas correspond to unstable orbits. The gray areas in the δ ψ maps indicate nonlibrating orbits. The grids correspond to osculating orbital elements with initial ψ = 90° for the inner evection, and initial ψ = 0° for the outer evection. To allow comparison with previous works, the upper horizontal axes give the distance to the planet in convenient units: planetary radii (R_{S}) for the inner evection, and Hill radii (R_{Hill}) for the outer evection. 

Open with DEXTER  
In the text 
Fig. 8 Dynamical maps of the δψ indicator in the (a, ψ)plane for the outer evection, with initial e = 0.4. The leftpanel corresponds to a grid of osculating elements, while the right panel is the same grid in terms of mean elements. Blue regions represent the conditions for which δ ψ ≲ 30°. Gray regionscorrespond to larger libration amplitudes or circulation of ψ. White regions are unstable orbits. The resonance occurs at a mean semimajor axis ā = 0.205 au. The upper horizontal axis is the distance in terms of R_{Hill}. 

Open with DEXTER  
In the text 
Fig. 9 Value of δe as a function of the distance from the planet in the region of the inner evection and for e = 0. The curve in violet corresponds to the case of an obliquity ϵ_{S} = 0° and orbital eccentricity e_{S} = 0° for Saturn. The green dots correspond to the case of ϵ_{S} = 26.73°, either with e_{S} = 0° or e_{S} = 0.054°. The minimum inside the peak of the curves marks the center of the evection resonance. 

Open with DEXTER  
In the text 
Fig. 10 Dynamical map of the δψ_{1} indicator in the osculating (a, e)plane for the inner evection, setting ψ_{1} = 120° initially. Case of Trojan satellites with equal masses. The resonance occurs at the same location as in the case of a single satellite (compare with Fig. 7). The upper horizontal axis gives the distance in R_{S}. Angles are in degrees. 

Open with DEXTER  
In the text 
Fig. 11 Dynamical maps in the osculating (a, e)plane displaying the amplitude δe_{1} of the eccentricity variation (top panel), the amplitude of the σ libration (middle panel), and the amplitude of the Δϖ libration (bottom panel) for a system of two equalmass moons initially in a Trojan ACR configuration. Angles are in degrees. Trojan motion is unstable in the white regions in less than 5 × 10^{4} yr, and it is most stable in the blue regions. The stripes of less stable Trojan motion observed for a < 0.03 au correspond to the locations of the evection resonance and its higher harmonics. 

Open with DEXTER  
In the text 
Fig. 12 Evolution of the resonant angles σ (cyan lines) and ψ_{i} (red/blue lines) for a pair of moons with equal masses, initially placed in a Trojan ACR configuration (σ = 60°, Δ ϖ = 60°) in the region of the outer evection (a_{1} = a_{2} ~ 0.2 au, e_{1} = e_{2}= 0.4). Top panel: case ψ_{1} = 0° and a_{1} = a_{2} = 0.1825 au initially. Bottom panel: case ψ_{2} = 0° and a_{1} = a_{2} = 0.1893 au initially. 

Open with DEXTER  
In the text 
Fig. 13 Same as Fig. 12, but in the k_{i}, h_{i} space (see text) and for the full time span of the simulation. The top panels correspond to the top panel of Fig. 12; the system survives for at least 10^{6} yr. The bottom panels correspond to the bottom panel of Fig. 12; the system survives for ~ 10^{3} yr. In both cases, λ_{⊙} = 0° initially. 

Open with DEXTER  
In the text 
Fig. 14 Dynamical maps of δψ_{i} in the osculating (a, e)plane for a Trojan pair of equal masses, initially placed in the ACR configuration. Initial conditions for the angles are the same as in Fig. 13. Blue islands correspond to libration of ψ _{1} (left panel) or ψ_{2} (right panel). Gray regions correspond to circulation of ψ_{i}, and white regions correspond to unstable orbits (compare to Fig. 7). The upper horizontal axis is the distance in terms of R_{Hill}. 

Open with DEXTER  
In the text 
Fig. 15 Adiabatic evolution of pair of Trojan satellites that migrate as a result of tidal effects, using Q = 100. Each column of panels corresponds to a different mass ratio m_{2}∕m_{1}: 2 × 10^{−5} (left panels),1 (middle panels), and 2 × 10^{5} (right panels). Eccentricity is shown in the top row, σ and Δ ϖ in the middle row, and ψ_{1}, ψ_{2} in the bottom row. Magenta and green identify identify the trailing and leading Trojan, respectively, except in the middle panels, where they identify σ and Δ ϖ, 

Open with DEXTER  
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.