The hierarchical stability of the seven known large size ratio triple asteroids using the empirical stability parameters

In this study, the hierarchical stability of the seven known large size ratio triple asteroids is investigated. The effect of the solar gravity and primary’s J2 are considered. The force function is expanded in terms of mass ratios based on the Hill’s approximation and the large size ratio property. The empirical stability parameters are used to examine the hierarchical stability of the triple asteroids. It is found that the all the known large size ratio triple asteroid systems are hierarchically stable. This study provides useful information for future evolutions of the triple asteroids.


Introduction
The triple asteroids are newly found systems in the Solar System. The system (87) Sylvia is the first known triple asteroid system to be confirmed in 2004, for which the outer satellite Romulus was discovered in 2001 using Keck Margot and Brown 2001) and was also detected in HST images (Storrs et al. 2001), and the inner satellite Remus was discovered in 2004 using adaptive optics system on the Very Large Telescope (Marchis et al. 2005). Recently, Fang et al. (2012) derived the masses, orbits and primary's oblateness of (87) Sylvia system, and examined the short-term and long-term stability of the orbits of two satellites numerically. At the time of writing this paper, nine triple asteroids (including two Kuiper Belt objects) have been identified, namely, (45) Eugenia (Marchis et al. 2010;Beauvalet et al. 2011), (87) Sylvia (Marchis et al. 2005;Fang et al. 2012), (93) Minerva (Marchis et al. , 2013, (216) Kleopatra (Descamps et al. 2011), (136617) 1994 CC (Fang et al. 2011;Brozović et al. 2011), (153591) 2001SN263 (Fang et al. 2011 Haumea (Ragozzine andBrown 2009), (47171) 1999 TC 36 (Benecchi et al. 2010), and (3749) Balam (Marchis et al. 2008(Marchis et al. , 2012. The first seven are large size ratio triple asteroids, while the last two triple asteroids have their central components of a comparable size and orbiting very close to each other. In connection with these recently discovered triple asteroids, the hierarchical stability of the large size ratio triple asteroid system including the effect of the solar gravity and primary's J 2 is analyzed in this study. The system is said to be hierarchically stable if its hierarchical configuration is not altered over a time scale that is much longer than the basic periods of the system (Milani and Nobili 1983;Walker and Roy 1983). Thus, either orbit crossing or escape of the bodies from the system is not possible. Although the fourbody that consists of the triple asteroids and the Sun is nonintegrable, the hierarchical stability can be examined using the Walker-Roy empirical stability criteria (Walker 1983a;Walker and Roy 1983). In this study, only the seven large size ratio triple asteroids are concerned because it is convenient to expand the force function in terms of mass ratios. The effect of the primary's J 2 is added to the traditional empirical stability parameters. In previous studies, numerical simulations were performed to investigate orbital stability of (87) Sylvia (Winter et al. 2009;Frouard and Compère 2012;Fang et al. 2012) and stability regions around (153591) 2001 SN263 (Araujo et al. 2012).
For clarity, this paper is structured as follows: in Sect. 2, the model of the large size ratio triple asteroid system involving the Sun in the four-body problem is described, and the force function is expanded in terms of mass ratios. In Sect. 3, the empirical stability parameters are used to examine the hierarchical stability of the triple asteroids. Finally, our conclusions are presented in Sect. 4.
2 Model of the large size ratio triple asteroid system involving the solar gravity and primary's J 2 in the four-body problem The large size ratio triple asteroids involving the Sun in the four-body problem consists of a nonspherical primary with mass m 1 , two small satellites (also known as moonlets) with masses m 2 and m 3 , and the Sun with mass m 4 . The simple four-body hierarchy instead of the double hierarchy is adopted here (Milani and Nobili 1983;Walker 1983a;Walker and Roy 1983). In other words, these four bodies are arranged that the body m i (i = 2, 3, 4) is in orbit about the barycenter of all of its inner bodies, i.e., m 1 , . . . , m i−1 .
Hill's approximation is used here since the triple asteroids are close to each other, their center of mass is far away from the Sun, and the mass of the Sun dominates the four-body system (Hamilton and Krivov 1997;Scheeres 1998). Additionally, for the large size ratio triple asteroids, the mass of the primary is much larger than its moonlets, i.e. m 1 m 2 , m 3 . Based on the previous study (cf. Walker 1983a, the equations before Eq. (17)) and adding the effect of J 2 , the equations of motion in the expansion of the mass ratios are written as where ; P r is the Legendre function of degree r; J 2 is the primary's oblateness coefficient; R e is the primary's reference radius; and ϕ i is the latitude of m i in the primary-fixed coordinate system. Taking the (216) Kleopatra system as an example, the primary is Kleopatra with the oblateness coefficient J 2 = 0.6, and two moonlets are S/2008 ((216)) 1 and S/2008 ((216)) 2 (Descamps et al. 2011). Therefore, the four-body system consists of the primary Kleopatra, S/2008 ((216)) 1, S/2008 ((216)) 1, and the Sun. Because the semi-major axis of the moonlet is usually several times larger than the primary's reference radius (Richardson and Walsh 2006;Johnston 2013), the effects of higher harmonic expansions are negligible. Therefore, only J 2 term is considered in this paper.
Note that Eq. (1) is different from Eq. (17) in Walker (1983a). The reason is that m 1 dominates the whole system in Walker (1983a) while m 4 dominates our system and m 1 only dominates the large size ratio triple asteroids in the present paper.
Take m 1 as the unit of the mass, the scaled masses m 2 , m 3 , and m 4 are denoted as μ 2 , μ 3 , and μ 4 , respectively. Reminding of the large size ratio property (m 1 m 2 , m 3 ) and the Hill's approximation (ρ 4 ρ 2 , ρ 3 and m 4 m 1 , m 2 , m 3 ), the equations of motion take a compact form Note that the expressions are expanded in terms of mass ratios due to the properties of the large size ratio triple asteroid systems, which is different from Walker and Roy (1983) where the expressions were expanded in terms of the ratios of the orbital radii.

Hierarchical stability of the large size ratio triple asteroids
According to Walker and Roy (1983), in our problem at the collinear configuration m 1 -m 2 -m 3 -m 4 where both the mutual perturbations of the bodies and the perturbation due to primary's J 2 are greatest. Note that this collinear configuration is only an extreme scenario in the mathematical sense. In this scenario, C ij = 1 and ϕ i = 0, so that P r (C ij ) = 1 and P 2 (sin ϕ i ) = −1/2 ∀r, i, j ; there is also the equation S ij = 1/(1 − α ij ). Walker and Roy (1983) introduced the empirical stability parameters that are the sum of the disturbing terms at this configuration, when they expanded the force function in terms of the ratios of the orbital radii. In the  Beauvalet et al. (2011). c Fang et al. (2012). d Marchis et al. (2013). e Descamps et al. (2011). f Ragozzine and Brown (2009). g Fang et al. (2011) four-body problem, the physical significance of each empirical stability parameter Σ i (i = 2, 3, 4) can be considered as the measure of the disturbance of the system's other components on the orbit of m i relative to M i−1 ; in order to ensure the hierarchical stability, the values of Σ i must be less than unity and generally less than 10 −2 (Walker and Roy 1983). Walker and Roy (1983) also pointed out that the eccentricity's effect on the Σ i parameters is small. Thus, the eccentricity's effect is not considered in this study, just as Walker et al. (1980), Walker (1983b), and Veras and Armitage (2004) when they investigated the Hill stability of the triple systems using analytic methods. In our problem for the large size ratio triple asteroid systems, the empirical stability parameters Σ 2 , Σ 3 , and Σ 4 at the collinear configuration m 1 -m 2 -m 3 -m 4 can be expressed as Σ 4 = J 2 R 2 e /2/ρ 2 4 + μ 2 α 2 24 + μ 3 α 2 34 .
3.1 Without consideration of the solar gravity effect and primary's J 2 If the effect of the solar gravity and primary's J 2 is not taken into account, i.e. μ 4 = 0, ρ 4 = 0, and Σ 4 = 0, the system only consists of three masses. Table 1 lists ratios of masses and orbital radii for the seven known large size ratio triple asteroids. Based on Eqs. (5)-(6), the values of Σ 2 and Σ 3 for the seven known large size ratio triple asteroids are shown in Fig. 1 using the parameters listed in Table 1. It is easy to see that the values of Σ 2 and Σ 3 are all much less than 10 −2 . Therefore, these seven large size ratio triple asteroids are all hierarchically stable when without considering the solar gravity effect and primary's J 2 . Fig. 1 The values of empirical stability parameters Σ 2 and Σ 3 for the seven known large size ratio triple asteroids without considering the solar gravity effect and primary's J 2 . The actual large size ratio triple asteroids cases are marked with asterisk 3.2 With consideration of the solar gravity effect and primary's J 2 When the solar gravity effect and primary's J 2 is taken into consideration, the case is a four-body problem, which is a bit more complicated. According to Eqs. (5)- (7), the values of Σ 2 , Σ 3 and Σ 4 for the seven known large size ratio triple asteroids are shown in Fig. 2. It can be seen that the values of Σ 2 , Σ 3 and Σ 4 are all much less than 10 −2 for all the known high size ratio asteroids. It means that all the known high size ratio asteroids are hierarchically stable. This result is consistent with numerical investigations of (87) Sylvia (Winter et al. 2009;Frouard and Compère 2012;Fang et al. 2012), which showed that the triple system (87) Sylvia is orbitally stable.

Conclusions
In the present paper, the hierarchical stability of the seven known large size ratio triple asteroids is analyzed. The solar gravity and primary's J 2 are considered. The force function Fig. 2 The values of empirical stability parameters Σ 2 , Σ 3 and Σ 4 for the seven known large size ratio triple asteroids when considering the solar gravity effect and primary's J 2 . The actual large size ratio triple asteroids cases are marked with asterisk for the four-body problem is expanded in terms of mass ratios based on the Hill's approximation and the large size ratio property. The hierarchical stability of the triple asteroids is examined using the empirical stability parameters. Without considering the solar gravity and primary's J 2 , it is found that all the seven known high size ratio asteroids are hierarchically stable. While considering the solar gravity and primary's J 2 , all the known high size ratio asteroids are still hierarchically stable. We concluded that the hierarchical arrangement of the system will not be altered for (45)