Formation of the JanusEpimetheus system through collisions^{⋆}
^{1}
Univ. Estadual Paulista – UNESP, Grupo de Dinâmica Orbital &
Planetologia, Guaratinguetá, CEP, 12516410
SP, Brazil
^{2}
Theoretische Physik II, Physikalisches Institut, Universität
Bayreuth, 95440
Bayreuth,
Germany
email: lucas.t@treffenstaedt.de
Received: 18 December 2014
Accepted: 5 September 2015
Context. Coorbital systems are bodies that share the same mean orbit. They can be divided into different families according to the relative mass of the coorbital partners and the particularities of their movement. Janus and Epimetheus are unique in that they are the only known coorbital pair of comparable masses and thus the only known system in mutual horseshoe orbit.
Aims. We aim to establish whether the JanusEpimetheus system might have formed by disruption of an object in the current orbit of Epimetheus.
Methods. We assumed that four large main fragments were formed and neglected smaller fragments. We used numerical integration of the full Nbody problem to study the evolution of different fragment arrangements. Collisions were assumed to result in perfectly inelastic merging of bodies. We statistically analysed the outcome of these simulations to infer whether coorbital systems might have formed from the chosen initial conditions.
Results. Depending on the range of initial conditions, up to 9% of the simulations evolve into coorbital systems. Initial velocities around the escape velocity of Janus yield the highest formation probability. Analysis of the evolution shows that all coorbital systems are produced via secondary collisions. The velocity of these collisions needs to be low enough that the fragments can merge and not be destroyed. Generally, collisions are found to be faster than an approximate cutoff velocity threshold. However, given a sufficiently low initial velocity, up to 15% of collisions is expected to result in merging. Hence, the results of this study show that the considered formation scenario is viable.
Key words: planets and satellites: formation / planets and satellites: dynamical evolution and stability
Movies associated to Figs. 7 and 8 are available in electronic form at http://www.aanda.org
© ESO, 2015
1. Introduction
In a coorbital system, two or more bodies share a mean orbit around a central body in a 1:1 orbital resonance. The possible orbits have been described by Lagrange in his Analytical Mechanics in 1788. We speak of tadpole orbits when one body librates around the Lagrange points L_{4} or L_{5} of another. Libration around L_{3}, L_{4}, and L_{5} constitutes a horseshoe orbit (Brown 1911). Many trojan asteroids are in coorbital trajectories with planets, but only in the Saturnian system are satellites known to be in coorbital motion. While there are many examples of bodies with much smaller coorbital partners, the Saturn moons Janus and Epimetheus are a unique system of two bodies with similar mass in horseshoe orbit. This coorbital motion was confirmed in 1980 by the Voyager I probe (Synnott et al. 1981). Harrington & Seidelmann (1981) showed the stability of this system through numerical integrations including the oblateness of Saturn. Dermott & Murray (1981) developed an analytical study estimating the radial and angular amplitudes of each horseshoe orbit. A theory of motion was proposed two years later by Yoder et al. (1983), who explored the phase space of coorbital systems with comparable masses, identifying the equilibrium points and showing the possible trajectory types (tadpole and horseshoe). Considering a system with up to nine equalmass coorbital bodies, Salo & Yoder (1988) found the location of the equilibrium points and their stability. Following a Hamiltonian approach, Sicardy & Dubois (2003) demonstrated the connection between coorbital systems of comparable masses with that of the restricted case.
Parameters of Janus and Epimetheus (Jacobson et al. 2008; Thomas et al. 2013).
The formation of coorbital systems has been the subject of intensive studies. In 1979, Yoder (1979) proposed three possible mechanisms for the formation of the Jovian trojans:

1.
early capture during the accretion of Jupiter;

2.
the continuous capture of asteroids; and

3.
transfer of Jovian satellites to the Lagrange points.
In the first scenario, gas drag plays an important role. Peale (1993) and Murray (1994) showed that, in the presence of a gas, L_{4} and L_{5} can still be stable for protoJupiter. In 1995, Kary & Lissauer (1995) showed that a protoplanet on an eccentric orbit can indeed trap planetesimals in a 1:1 resonance, decaying from a largelibration horseshoe orbit to a tadpole orbit around L_{5}. Fleming & Hamilton (2000) explored the effect due to mass accretion and orbital migration on the dynamics of coorbital systems. They showed how these perturbations affect the radial and angular amplitudes of the trajectories. Kortenkamp (2005) later proposed a mechanism to efficiently capture planetesimals on a satellite orbit after gravitational scattering by the protoplanet, making this a more likely fate than trapping in a coorbital. Chanut et al. (2008) determined the size limit of trapping planetesimals via nebular gas drag to be 0.2 m per AU of orbit of the protoplanet. They (Chanut et al. 2013) later extended their approach to a nonuniform gas with a cavity in the vicinity of a heavy planet, but found little to no coorbitals for eccentricities comparable to Janus and Epimetheus. All these scenarios result in coorbital configurations with a large secondary body and a small partner, that is, with a low mass ratio. This almost excludes trapping as a mechanism for the formation of the JanusEpimetheus system. Furthermore, Hadjidemetriou & Voyatzis (2011) showed that a system of two equalmass bodies in a 1:1 resonance tends to evolve under gas drag into a closely bound binary satellite system. Another take on the subject by Cresswell & Nelson (2006) was the study of the evolution of a swarm of similarly sized planetesimals in a gas disk. They showed that the formation of meanmotion resonances, including 1:1 resonances, is commonplace, but that the bodies are still prone to be lost to accretion by the central body. Giuppone et al. (2012) followed a similar approach to Chanut et al. (2013), but studying a hypothetical extrasolar system. They showed that a pair of inwardly migrating planets can enter into a 1:1 resonance at the edge of a cavity in a gas disk. Whether their result is applicable to Janus and Epimetheus is not known to us.
Laughlin & Chambers (2002) showed that stable 1:1 resonances of equalsized planets around a central body exist. They hypothesized that at the L_{5} point of a planet in the accretion disk of a protostar, another planet could accrete. Beaugé et al. (2007) investigated formation by accretion of a swarm of planetesimals in a tadpole orbit around L_{4} or L_{5} of a giant extrasolar planet. They showed that an Earthtype body can form, but with a very low efficiency and thus a low final mass. Chiang & Lithwick (2005) considered this a likely process for the formation of the Neptune trojans. Izidoro et al. (2010) studied the mechanism for the Saturnian system and showed that formation of coorbital satellite pairs with a relative mass of 10^{13} to 10^{9} is possible. This again rules out the formation of a coorbital pair with similar masses.
While much research has been done on the subject of coorbital systems, there is as yet no satisfactory explanation for the origin of the JanusEpimetheus system in particular. Our approach here is based on the collision of two bodies at or close to the current orbit of Epimetheus. This is consistent with the equal densities (see Table 1) and compositions of Janus and Epimetheus. Moreover, results by Charnoz et al. (2010) show a recent formation of the Saturn moonlets from its rings through series of subsequent collisions. Crida & Charnoz (2012) later extended this scenario to all regular moons of Saturn. In this context, lowvelocity collisions are frequent and are almost inevitable, especially for the predecessors of Janus and Epimemetheus. We assume the disruption into four main fragments that later collide again constructively to form a pair of bodies in a coorbital configuration. We establish the statistical likelihood of this outcome in 8000 numerical simulations. We used Chambers’ Mercury (Chambers 1999) NBody integrator with a BurlischStoer algorithm.
In Sect. 2, we explain our model in detail and introduce the methods used to obtain our results. Section 3 shows the setup and results of our simulations. Finally, we discuss the viability of our model considering cases with six and eight fragments, and propose opportunities for future studies in Sect. 4.
2. Model and method
Fig. 1 Sketch of typical coorbital formation scenario from our simulations. 

Open with DEXTER 
We investigated whether the Janus/Epimetheus coorbital system might have formed by disruption of a large body in the mean orbit of Epimetheus. Figures 1a, b illustrates the scenario that we believe can lead to coorbital configurations. We assumed that four large fragments are formed by a collision event. While the initial collision is not subject of this study, we note that a low collision velocity and low density (compare Table 1) of the colliding bodies is a requirement. The collision is neither completely destructive nor is it constructive. After the collision event, fragments spread out on orbits around Saturn. Fragments can be lost by ejection from the system, collision with Saturn, or collision with each other. In the last case, fragments can unite to form a larger body, given a sufficiently low collision velocity (see Kortenkamp & Wetherill 2000).
Our simulations start right after the initial collision event, with the four fragments α, β, γ, and δ on the vertices of a square, moving in opposite directions with equal speed (see Fig. 2 for details). We always chose equal masses for bodies opposite on the square. The length a_{0} of the sides of the square was chosen to be a_{0} = 2R_{1} + d_{0} , where R_{1} is the radius of the largest body and d_{0} = 100m. The parameters varied in this study include the masses m_{i} of the fragments (given as fractions of the mass of Janus and Epimetheus, see Tables 1 and 3), the orientation of the square in the orbital plane denoted by the angle ϕ (going from 0 to π because of symmetry), and the initial speed v_{init} (see Fig. 2) of the fragments in the coordinate system moving with v_{0} relative to Saturn. The position and velocity { r_{0},v_{0} } of the centre of the square are equal to the state vector of Epimetheus at an arbitrarily chosen point of its orbit. This is the same for all simulations.
Properties of Saturn as used in the simulations.
Parameters for our simulations.
We used the package Mercury to numerically integrate the system using a BurlischStoer algorithm (Chambers 1999). We included the moons Mimas and Tethys into our simulations, since they exert the strongest influence on the JanusEpimetheus system, but neglected other objects. We also considered the oblateness of Saturn through J_{2}, J_{4}, and J_{6} (Table 2). To obtain the state vectors of Mimas, Tethys, and Epimetheus (centre of the square from which the fragments start), we used the data given by Jacobson et al. (2008) and followed the approach presented by Renner & Sicardy (2006) to transform the geometric elements. We took the mean distance between the system JanusEpimetheus and Saturn as the unit of length R_{0}. The simulation code was set to record the state vectors of all bodies every 1 / 100 day. We need this high time resolution to study the velocities and angles of collisional events between fragments. Mercury assumes a perfectly inelastic collision where the colliding bodies aggregate into a new body of the added masses of the colliding bodies, and the momentum is preserved. We therefore needed to analyse the collisions to ensure that this assumption is realistic. The outcome of a collision is strongly dependent on the collision velocity and the physical characteristics (density, strength, etc.) of each body. Kortenkamp & Wetherill (2000) established an escape velocity term (1)for the encounter of rocky bodies, where m_{i} is the mass of body i, R_{i} is its radius, and G is the gravitational constant. For collision velocities below v_{max}, a collision can be assumed to be constructive. For a collision of Janus and Epimetheus, v_{max} ≈ 58 m/s. However, their study involved bodies with a density about five times as high as that of Janus and Epimetheus. Since in a collision event between fluffy bodies (such as Janus and Epimetheus, whose density is lower than 700 kg / m^{3}) kinetic energy is lost to deformation, higher velocity collision events can still be constructive. Highvelocity collisions can also result in the catastrophic disruption of one or both colliding bodies. Multiple studies on the disruption threshold (Colwell 1994; Benz & Asphaug 1999; Leinhardt & Stewart 2009) have been conducted. Using Fig. 11 from Leinhardt & Stewart (2009) and a target radius of 6 × 10^{4} m (the rough size of Epimetheus), we obtain a disruption threshold of J/kg. This gives a cutoff velocity (2)(where m_{T} is the mass of the target and m_{P} the mass of the projectile) of approximately 60 m/s for a collision with Janus. Therefore we used the first estimate of 58 m/s as a cutoff, keeping in mind that the actual threshold might still vary and is in fact dependent on which of the fragments collide.
Fig. 2 Initial configuration of the simulations. are given in the reference frame moving with v_{0}. The initial velocity of each fragment in the simulation frame of reference is given by v_{i} = . The direction of is given by the vector from the centre of the square to the respective vertex. The collision fragments α, β, γ, and δ are positioned on the vertices of a square centred on r_{0}, the position vector with respect to Saturn. The angle ϕ gives the rotation of the square around the vector r_{0} × v_{0} (the normal of the orbital plane). For ϕ = 0α, γ, and Saturn are connected by a straight line, with γ between α and Saturn. The fragments are always set on the vertices in a counterclockwise sense, as shown in the sketch. The dimensions here are not to scale. 

Open with DEXTER 
Because of the large number of simulations, we needed an automatic filtering algorithm to find the systems that evolved into a coorbital. We used a fourstep filtering process.

1.
We selected the systems in which at least two fragments havesurvived until the end of the simulation (20–30 yr). Fragments canbe lost through ejection, collision with Saturn, or collision witheach other.

2.
From the remaining systems, we selected the systems with exactly one pair of fragments that survived until the end of the simulation and overlap in their range of semimajor axes for the time after the last collision in the simulation.

3.
From these, we selected the systems in which the mean semimajor axis of the surviving fragments is smaller than 1.2 R_{0} and the minima and maxima of the semimajor axes are no more than 0.02 R_{0} apart.

4.
We plot the semimajor axes of the remaining fragments against time and inspected them visually to check for periodic crossings. We thus determined false positives of the automatic filtering process and removed them.
3. Simulations and results
Fig. 3 Cumulative distribution of collision velocities for coorbital systems in batch 3. 

Open with DEXTER 
In each simulation, we sampled the parameters v_{init} and ϕ in equally spaced steps. We started our investigation with a batch of 4000 simulations, covering 20 different initial speeds v_{init} and 200 different angles ϕ (see Table 3–1). The initial speeds were between 10 and 100 times as high as the escape velocity of Janus. We found that around 1.5% of the simulations evolved into coorbital systems. From the analysis of these systems, we were able to establish a common formation scenario as depicted in Fig. 1: after the initial collision event, the fragments first spread out on orbits around Saturn. Then, they recollide pairwise, most commonly α with γ and β with δ. The velocities of these collision events are very high (1000 m/s and above), far beyond where a perfectly inelastic joining, as is assumed by Mercury, is realistic. But the moons created in this way quickly (in just a few orbits) enter into coorbital horseshoe motion after the second moon has been formed. No coorbital systems are the result of ejection of fragments or collision with Saturn.
Then, we investigated whether lower initial velocities (Table 3–2) result in lower collision velocities. We reduced the parameter range of v_{init} by a full order of magnitude. This indeed reduced the collision velocities by nearly an order of magnitude, but not to the point where simple merging of bodies was a realistic outcome for any of the collisions. We note, however, that coorbital formation increased to 5.7%.
Our third batch of simulations was divided into three parts: first, we simulated with the same masses as before, but at even lower initial velocities (Table 3–3a). As expected, the collision velocities decrease even more (see Fig. 3 for a comparison of collision velocities). The portion of coorbitals increases to 8.4%. Another result is that the frequency of coorbital formation decreases for very low initial velocities (below ~0.015 R_{0}/ d or around half of the Janus escape velocity). The highest observed frequency of coorbitals is around 0.025 R_{0}/ d.
Fig. 4 Frequency of coorbitals as a function of the initial speed v_{init} of the fragments for batch 3. 

Open with DEXTER 
Fig. 5 Period of the horseshoe orbits. Note the similar distribution for all masses versus selected masses. The masses of the bodies have almost no influence on the period. 

Open with DEXTER 
For the second and third part of this batch, we decreased the masses of the fragments by a half. One part was ran at initial velocities like in the previous batches (Table 3–3b) and the other at lower ones (Table 3–3c). By comparing simulations 3a and 3c, we see that the mass of the fragments has little effect on the distribution of generated coorbital systems. Only the decrease in occurrence of coorbitals for the lowest initial velocities (see Fig. 4) is less prominent for lower mass. The decrease at these initial velocities indicates a lower limit for the formation of coorbitals. Supplementary simulations indicate that for even lower velocities v_{init}< 2 × 10^{5}R_{0}/ d, formation of coorbitals is quite unlikely (~0.7%) compared to the other simulations. A cutoff would probably depend on the mass of the fragments, assuming it is related to the escape velocity of the fragments. This could explain the higher formation probability for lowermass fragments at low initial velocities.
Fig. 6 Four classes of coorbitals by mass ratio and (see Table 4) The variables are inversely proportional for horseshoe orbits (Dermott & Murray 1981, Eq. (1)). The red circle indicates the only pair in a tadpole orbit we observed in this batch of simulations. 

Open with DEXTER 
We thus have two constraints on the initial conditions if we wish to determine a range likely to spawn coorbitals: the collision velocities are highly dependent on the initial velocities. Low collision velocities are required for the fragments to collide constructively, thus low initial velocities. At the same time, the formation frequency of coorbitals decreases after a certain point, so the initial velocity cannot be arbitrarily low. The existence and width of a sweet spot for the formation of coorbitals is key to decide the viability of this model.
Fig. 7 Example orbits from all four classes. Radial axis is the semimajor axis relative to the mean semimajor axis. Angular is the mean longitude in a frame rotating with the mean motion associated with the mean semimajor axis. To view the evolution of each of these cases, see the associated movies available online. 

Open with DEXTER 
We examined the results of simulation batch 3c more closely. In many but not all cases, the fragments collided in such a way that the final mass of the bodies was equal to Janus and Epimetheus. This is the case for 36 of the 87 coorbital systems we obtained in this batch, or 3.6% of all simulations in the batch (see Table 4). The collision velocities reach up to about 450 m/s. The collisions occur at angles smaller than 5° for nearly all cases. Assuming a limit of 58 m/s for constructive collision, we find that 15.8% of all collision events in this batch remain (see Fig. 3). By multiplying the original success rate of 3.6% by 15.8% squared (since both collision events need to be below the threshold), we obtain a corrected success rate of 0.09% to generate Janus and Epimetheus or 0.2% to generate any coorbital.
We compared the libration period of the horseshoe orbits to the period of the horseshoe orbit of Janus and Epimetheus (see Fig. 5). Janus and Epimetheus cross orbits roughly every four years. This period is highly variable for different simulations and shows no particular connection to the mass of the moons. Instead, there seems to be a general preference for short periods.
We see more influence of the mass of the fragments on the trajectories by examining the relative range of the semimajor axes. With our setup of four fragments, two each of identical mass, there are four possible combinations of a final coorbital pair, listed in Table 4.
Classes of coorbital systems by mass.
We now consider the ratio of radial amplitudes (3)where a_{1} is the set of semimajor axes of the heavier body of the final coorbital pair after the last collision event in the simulation, and a_{2} is the same set for the lighter body. For the JanusEpimetheus system, this is (4)By plotting against the ratio of masses (see Fig. 6), we recover the expected inversely linear correlation between the two. Each class has similar ratio of radial amplitudes. There is one simulation with significant deviance (see Fig. 6, red circle) from this law. Closer investigation reveals that in this case, the surviving bodies enter into a tadpole orbit (see Fig. 8). We have found no other case of tadpole orbits in our simulations.
Fig. 8 Single case of tadpole orbit. The fragments have the mass of Janus and Epimetheus. The plot shows the semimajor axis relative to the mean semimajor axis plotted against the mean longitude in a frame rotating with the mean motion of the mean semimajor axis. (Online movie.) 

Open with DEXTER 
4. Conclusion and outlook
We here investigated the possibility of a formation process of the JanusEpimetheus system starting from a collision event that spawns four large fragments. Additional debris was neglected. We used numerical integration of the Nbody problem to obtain the evolution of 8000 initial configurations. We selected the configurations that led to the formation of a coorbital system. By analysing these, we established a common scenario that leads to the formation of coorbital systems: after the fragments spread out on orbit, they recollide and form a pair of coorbiting bodies.
We investigated the circumstances of these secondary collisions and showed that the initial velocity has a high influence on the collision velocities. At the same time, we found that lower initial velocities yield a higher probability of coorbital formation. This reinforces our assertion that for our assumed initial condition to be achieved, the prior collision event must have occurred at a low velocity. It is unclear, however, whether the collision velocities observed in these simulations are low enough that a perfectly inelastic collision and joining of the bodies is a good approximation. We assume that relative velocities below 58 m/s allow the fragments to collide constructively. Previous studies on the collisions of celestial bodies are not fully applicable to our system. Most works dealt either with collisions of rocky bodies and/or of a small projectile hitting a large target. Therefore it is necessary to study the collision behaviour of large lowdensity fluffy objects such as Janus and Epimetheus in more detail. Studies in this area would help to ultimately decide the viability of our formation model. We also found that there is a lower cutoff below which the formation of coorbitals decreases significantly. We attribute this to the inability of the fragments to escape their combined gravity when the kinetic energy is too low. This imposes two constraints on the initial velocity: it must be low enough to obtain reasonably low collision velocities. At the same time, it must be high enough to allow the fragments to escape their mutual gravitational pull. Our results indicate that this sweet spot is quite narrow, but indeed exists.
Another limiting feature of our physical model is the number of initial fragments adopted in the simulations. To strengthen the scenario proposed here, we performed two sets of simulations, one considering six, the other eight initial fragments. Considering initial velocities as those of simulation 3c, 100 simulations were performed for each case (six and eight fragments). The results showed an even higher rate of coorbital generation (18–20%). Therefore, the scenario proposed here is somewhat robust.
Furthermore, we are aware that our initial conditions are not very natural. Instead, we chose a simple configuration to constrain the otherwise quite large number of free parameters. The main result of this study is therefore that the proposed formation model is possible. We are confident to find coorbital formation for a wider range of initial conditions building upon this result. Ideally, a future work would include the initial collision as part of the simulation and use a better approximation for the secondary collision events.
Acknowledgments
We thank Rafael Sfair for providing file storage for our results on short notice. Othon C. Winter thanks FABESP (proc. 2011/081713) and CNPq. The constructive comments and suggestions of the referee significantly improved this paper.
References
 Beaugé, C., Sándor, Z., Érdi, B., & Suli, A. 2007, A&A, 463, 359 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Benz, W., & Asphaug, E. 1999, Icarus, 142, 5 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, E. W. 1911, MNRAS 71, 438 [Google Scholar]
 Chambers, J. E. 1999, MNRAS, 304, 793 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Chanut, T., Winter, O., & Tsuchida, M. 2008, A&A, 481, 519 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chanut, T., Winter, O., & Tsuchida, M. 2013, A&A, 552, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Charnoz, S., Salmon, J., & Crida, A. 2010, Nature, 465, 752 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Chiang, E., & Lithwick, Y. 2005, ApJ, 628, 520 [NASA ADS] [CrossRef] [Google Scholar]
 Colwell, J. 1994, Planet. Space Sci., 42, 1139 [NASA ADS] [CrossRef] [Google Scholar]
 Cresswell, P., & Nelson, R. 2006, A&A, 450, 833 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Crida, A., & Charnoz, S. 2012, Science, 338, 1196 [NASA ADS] [CrossRef] [Google Scholar]
 Dermott, S. F., & Murray, C. D. 1981, Icarus, 48, 12 [NASA ADS] [CrossRef] [Google Scholar]
 Fleming, H. J., & Hamilton, D. P. 2000, Icarus, 148, 479 [NASA ADS] [CrossRef] [Google Scholar]
 Giuppone, C., BenítezLlambay, P., & Beaugé, C. 2012, MNRAS, 421, 356 [NASA ADS] [Google Scholar]
 Hadjidemetriou, J. D., & Voyatzis, G. 2011, Celest. Mech. Dyn. Astron., 111, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Harrington, R., & Seidelmann, P. 1981, Icarus, 47, 97 [NASA ADS] [CrossRef] [Google Scholar]
 Izidoro, A., Winter, O., & Tsuchida, M. 2010, MNRAS, 405, 2132 [NASA ADS] [Google Scholar]
 Jacobson, R., Spitale, J., Porco, C., et al. 2008, AJ, 135, 261 [NASA ADS] [CrossRef] [Google Scholar]
 Kary, D. M., & Lissauer, J. J. 1995, Icarus, 117, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Kortenkamp, S. J. 2005, Icarus, 175, 409 [NASA ADS] [CrossRef] [Google Scholar]
 Kortenkamp, S. J., & Wetherill, G. W. 2000, Icarus, 143, 60 [NASA ADS] [CrossRef] [Google Scholar]
 Laughlin, G., & Chambers, J. E. 2002, AJ, 124, 592 [NASA ADS] [CrossRef] [Google Scholar]
 Leinhardt, Z. M., & Stewart, S. T. 2009, Icarus, 199, 542 [NASA ADS] [CrossRef] [Google Scholar]
 Murray, C. D. 1994, Icarus, 112, 465 [NASA ADS] [CrossRef] [Google Scholar]
 Peale, S. 1993, Icarus, 106, 308 [NASA ADS] [CrossRef] [Google Scholar]
 Renner, S., & Sicardy, B. 2006, Celest. Mech. Dyn. Astron., 94, 237 [NASA ADS] [CrossRef] [Google Scholar]
 Salo, H., & Yoder, C. 1988, A&A, 205, 309 [NASA ADS] [Google Scholar]
 Sicardy, B., & Dubois, V. 2003, Celest. Mech. Dyn. Astron., 86, 321 [NASA ADS] [CrossRef] [Google Scholar]
 Synnott, S. P., Peters, C. F., Smith, B. A., & Morabito, L. A. 1981, Science, 212, 191 [NASA ADS] [CrossRef] [Google Scholar]
 Thomas, P., Burns, J., Hedman, M., et al. 2013, Icarus, 226, 999 [NASA ADS] [CrossRef] [Google Scholar]
 Yoder, C. F. 1979, Icarus, 40, 341 [NASA ADS] [CrossRef] [Google Scholar]
 Yoder, C. F., Colombo, G., Synnott, S., & Yoder, K. 1983, Icarus, 53, 431 [NASA ADS] [CrossRef] [Google Scholar]
Online material
Movie of Fig. 7 (Access here)
Movie of Fig. 7 (Access here)
Movie of Fig. 7 (Access here)
Movie of Fig. 7 (Access here)
Movie of Fig. 8 (Access here)
All Tables
All Figures
Fig. 1 Sketch of typical coorbital formation scenario from our simulations. 

Open with DEXTER  
In the text 
Fig. 2 Initial configuration of the simulations. are given in the reference frame moving with v_{0}. The initial velocity of each fragment in the simulation frame of reference is given by v_{i} = . The direction of is given by the vector from the centre of the square to the respective vertex. The collision fragments α, β, γ, and δ are positioned on the vertices of a square centred on r_{0}, the position vector with respect to Saturn. The angle ϕ gives the rotation of the square around the vector r_{0} × v_{0} (the normal of the orbital plane). For ϕ = 0α, γ, and Saturn are connected by a straight line, with γ between α and Saturn. The fragments are always set on the vertices in a counterclockwise sense, as shown in the sketch. The dimensions here are not to scale. 

Open with DEXTER  
In the text 
Fig. 3 Cumulative distribution of collision velocities for coorbital systems in batch 3. 

Open with DEXTER  
In the text 
Fig. 4 Frequency of coorbitals as a function of the initial speed v_{init} of the fragments for batch 3. 

Open with DEXTER  
In the text 
Fig. 5 Period of the horseshoe orbits. Note the similar distribution for all masses versus selected masses. The masses of the bodies have almost no influence on the period. 

Open with DEXTER  
In the text 
Fig. 6 Four classes of coorbitals by mass ratio and (see Table 4) The variables are inversely proportional for horseshoe orbits (Dermott & Murray 1981, Eq. (1)). The red circle indicates the only pair in a tadpole orbit we observed in this batch of simulations. 

Open with DEXTER  
In the text 
Fig. 7 Example orbits from all four classes. Radial axis is the semimajor axis relative to the mean semimajor axis. Angular is the mean longitude in a frame rotating with the mean motion associated with the mean semimajor axis. To view the evolution of each of these cases, see the associated movies available online. 

Open with DEXTER  
In the text 
Fig. 8 Single case of tadpole orbit. The fragments have the mass of Janus and Epimetheus. The plot shows the semimajor axis relative to the mean semimajor axis plotted against the mean longitude in a frame rotating with the mean motion of the mean semimajor axis. (Online movie.) 

Open with DEXTER  
In the text 