Report

Forced Resonant Migration of Pluto's Outer Satellites by Charon

See allHide authors and affiliations

Science  25 Aug 2006:
Vol. 313, Issue 5790, pp. 1107-1109
DOI: 10.1126/science.1127293

Abstract

Two small moons of Pluto have been discovered in low-eccentricity orbits exterior to Pluto's large satellite, Charon. All three satellite orbits are nearly coplanar, implying a common origin. It has been argued that Charon formed as a result of a giant impact with primordial Pluto. The orbital periods of the two new moons are nearly integer multiples of Charon's period, suggesting that they were driven outward by resonant interactions with Charon during its tidal orbital expansion. This could have been accomplished if Charon's orbit was eccentric during most of this orbital evolution, with the small moons originating as debris from the collision that produced Charon.

Hubble Space Telescope observations in May 2005 revealed two previously undetected satellites of Pluto, S/2005 P1 and S/2005 P2 (P1 and P2) (13). The diameters of P1 and P2 are respectively ∼160 km and ∼120 km if they are as dark as cometary nuclei, but only ∼35 km and ∼30 km for a Charon-like reflectivity, with each object containing ≤0.2% of Charon's mass. The orbital period, P, of P1 is 38.2065 ± 0.0014 days, whereas that of P2 is 24.8562 ± 0.00013 days. This implies that P1 and P2 have periods of about six and four times the 6.38723-day orbital period of Charon (13), and a nearly 3:2 mutual ratio. Such orbital period ratios are referred to as mean motion commensurabilities, where the mean motion is n = 2π /P. The moons orbit in the same sense and plane as Charon (1), consistent with all three having formed by a common process rather than, for example, through sequential capture events (4, 5).

The favored explanation for Charon's origin is a large impact (6, 7) similar to that believed to have produced Earth's Moon. Charon appears most likely to have formed intact, on an eccentric orbit (7). It seems plausible that the impact generated some accompanying orbiting debris as source material for the outer small satellites. However, the present orbital distances of P2 at 41.9 RP and P1 at 54.8 RP, where RP = 1180 ± 15 km is Pluto's radius, are far outside where debris from the type of collision favored for Charon's origin would be expected [e.g., typically within 15 RP (7, 8)].

Charon occupies a tidal end-state with an orbital period synchronous with Pluto'sdayand has likely undergone considerable outward tidal migration, perhaps by a factor of ∼4, to arrive at its current distance of 16.6 RP (9). The near 6:1 and 4:1 commensurabilities of P1 and P2 naturally lead to the supposition that they were driven outward by resonant torques from Charon (2). However, the required degree of expansion is large, and most observed solar system resonant configurations (e.g., that between Neptune and Pluto itself) act to excite the eccentricity of the resonantly perturbed object as the orbital radius of the perturbing object approaches it. This would have led to large, destabilizing eccentricities for P2 and P1. We propose instead that Pluto's small moons were captured into exterior corotation resonances with an eccentrically orbiting Charon. This is the same general type of resonance thought to confine the Neptune ring arcs (10, 11) and does not excite eccentricities, thus providing an explanation for how Pluto's distant new moons could have originated from the same impact event as Charon.

When the period ratio of two satellites is nearly that of two integers, resonant forcing of their orbits acts to maintain the locations of their orbital conjunctions relative to the pericenter of one or both of their orbits. Their mutual gravitational potential can be expanded as a Fourier series (12), but near a commensurability, the motion will be dominated by only a few resonant terms. An exterior moon with a period relative to Charon's of P/PC ∼ (m + 1), with m = 3 for P2 or m = 5 for P1, will be subject to multiple resonant terms of the form Φml (a, ac, e, ec) cos φml, with Embedded Image(1) and 0 ≤ lm (13). Here, λ and w̃ are the longitudes of the satellite and its orbital pericenter. Subscripts C refer to Charon, whereas unsubscripted quantities refer to a moon. Resonance capture causes the argument, φml, of the resonance to librate about some fixed value. The amplitude of a resonant term has the form Embedded Image(2) where a and e are semi-major axis and eccentricity, the quantity fml(α) is a combination of Laplace coefficients and their derivatives (12), and MC is Charon's mass. For any single resonance, the equations of motion admit two constants (1416): the Brouwer integral Embedded Image(3) and the Jacobi integral Embedded Image(4)

If an exterior moon were forced to migrate as Charon tidally evolved, Eq. 3 gives the moon's eccentricity, e(t), as a function of a(t). For resonances with lm, the resonant argument in Eq. 1 contains the moon's longitude of pericenter, w̃, and evolution in resonance increases the moon's eccentricity as a increases. Even for an initially zero eccentricity, a four-fold expansion of semi-major axis excites an eccentricity Embedded Image(5)

For m = 3 (i.e., the 4:1), the l μ m resonances would produce e30 = 0.781, e31 = 0.661, and e32 = 0.484. Eccentricities for m = 5 (the 6:1) are also substantial (e ≥ 0.4) for all l μ m terms. Because objects in the 6:1 and 4:1 resonances are separated in semi-major axis by only a1a2 ≈ 0.31a2 ≈ 0.24a1, crossing orbits and instability between such objects would be expected for e ≥ 0.3. However, for the special case of l = m, the resonance argument becomes φmm = (m + 1)λ – λCmC, which does not contain the outer moon's pericenter longitude nor excite its eccentricity. Because the resonant amplitude, Φmm, is proportional to e mC (Eq. 2), corotation resonances with m ≥ 1 require an eccentric perturber.

To describe the character of libration in a corotation resonance, we define the corotation distance as where n(acr) = Ωps with Ωps = nC/(m + 1) + m(dC/dt)/(m + 1) being the so-called pattern speed, and consider a perturbation in a, i.e., a = acr + Δa, Δaacr. The trajectory is found from the Jacobi integral, which can be recast as (16) Embedded Image(6) Embedded Image and describes an ellipse with major axis aφmax/m and minor axis Embedded Image, as shown in Fig. 1. The maximum radial width of the librating region is Embedded Image, where μCMC/MP = 0.116 and MP is Pluto's mass (3, 17). Thus, for example, with a Charon eccentricity of eC ∼0.2, moons librating in the 4:1 and 6:1 corotation resonances would have a Δa/acrw/acr ∼0.03 and 0.002, respectively. The libration period 2π/ωlib, with Embedded Image (16), is much longer than the + orbital period, so that over multiple orbits in inertial space a moon completes a single libration, as seen in the frame rotating with the resonant term (Fig. 2) (18).

Fig. 1.

Schematic of a corotation resonance island, shown in a frame rotating with the resonant pattern speed, Ωps. Dashed curve indicates the separatrix, dotted line is the corotation distance, and the primary is in the negative y-axis direction. The separatrix separates resonantly librating orbits from non-resonant circulating orbits. Particle trajectories displaced from a stable equilibrium point (filled circles) librate about the point, with each island of libration separated by unstable equilibrium points (open circles).

Fig. 2.

The m + 1 = 4 and m + 1 = 6 stability islands for the corotation + resonances near P2 and P1, shown when Charon is at pericenter. For each resonance, the (m + 1) islands are confined by a separatrix. The island widths are arbitrarily set to 0.1 for clarity; the actual widths are a function of Charon's eccentricity. The solid point contained in each island is a stable equilibrium about which trajectories librate. Each pattern rotates counter-clockwise at the pattern speed ΩpsnC/(m + 1), so for each Charon orbit, a new stable + equilibrium point is brought to conjunction at pericenter.

Now consider the response of a trapped moon to a slow tidal expansion of Charon's orbit according to (9) Embedded Image(7) Embedded Image where QP ∼102 is Pluto's tidal dissipation parameter, and kP ∼0.055 is its tidal Love number (9). Objects remain trapped in a moving resonance so long as the time, Embedded Image, that it takes for Charon to migrate a distance comparable to the libration zone half-width w, is much longer than the libration period (13, 16). This yields an adiabatic condition, emC f(α) δ (3π)(kP/QP)(RP/aC)5. Adopting aC(0) ≈ 4RP leads to a threshold value for Charon's initial eccentricity Embedded Image (19). This requires eC ≈ 0.2 for capture of P1, the more difficult to capture moon. Simulations of Charon's intact formation by impact (7) find cases with initial eccentricities and semi-major axes comparable to these values. Early resonant trapping of P1 and P2 would allow for them to have originated from impact debris extending only to ∼13 RP, which is reasonable given results of impact simulations (7, 8).

We are interested in the evolution of ec as Charon's orbit tidally expands, which is given by (20) Embedded Image(8) where ȧC is the expansion rate due only to Pluto tides, σ ≡ 2ωP – 3nC, with ωP being the spin frequency of Pluto, and A ≈ (kC/kP)(Qp/QC)(RP/RC) where similar densities have been assumed for Pluto and Charon. The first term on the right side of Eq. 8 is due to tides raised on Pluto, whereas the second term is due to tidal dissipation in Charon. Assuming that the k's are dominated by similar rigidities, then kC /kP ∼ (RC /RP)2, ARCQP/RPQC, and ec increases if Embedded Image. Given the uncertainties, Charon's eccentricity could have either grown or decayed during most of its orbital expansion. However, once ωP ≤ 3n/2 toward the end of its expansion, σ reversed sign and eC decayed.

Because the current eccentricity of Charon is very low as a result of tidal damping, P1 and P2 are probably not now in corotation resonances. On the other hand, it is unlikely they escaped simultaneously. At present 4n2nC = 3.14 × 10–7s–1 (3), so that P2 lies slightly inside the nominal positions of the m = 3 terms given by n2nC/4 + (3 – lC2/4, where ωC2 = ≈3.87 × 10–8s–1 is Charon's contribution to P2's apse precession (16). This suggests that P2 escaped from corotation first (perhaps because of a larger libration amplitude), just before Charon's orbital migration halted. Continued occupancy of the 3:2 resonance could have protected P2 from the remaining m = 3 resonances by controlling the rate of apse precession, d2/dt ≈ 3n1 – 2n2. This would have caused the resonance variable φ33 to circulate, and as Charon continued to migrate outward, the other φml resonances variables would circulate as well (16). Finding P2 in 3:2 libration currently (21) would be strong support for this notion, although it is also possible that a small free eccentricity could have caused a transition from libration to circulation as P1 and P2 separated somewhat.

Because most resonances excite large eccentricities, much of the original impact-produced debris may have been destabilized by mutual collisions or scattering into Charon or Pluto. Debris captured into some corotation islands could also have been dislodged through encounters with other high-eccentricity material. Nevertheless, it seems plausible that a fraction comparable to the tiny masses of P1 and P2 might have survived such stochastic removal processes.

What about corotation resonances other than the 4:1 and 6:1? For an eccentric Charon, the 3:1 corotation resonance is nearly overlapped by its Hill sphere at apocenter and was likely not a stable niche. Corotation resonances also occur when (m + p)npnC for p > 1, but these fall at distances Embedded Image and are shifted inward. Thus, those that fall in the vicinity of P1 and P2 have amplitudes that are dependent on a higher power of Charon's eccentricity (i.e., ∼Embedded Image, Embedded Image) and are weaker. Although transient forced eccentricities may interfere with the stability of adjacent p = 1 resonances, it remains an intriguing possibility that smaller, yet undetected moons may orbit Pluto near the 5:1.

Because the corotation resonances we invoke no longer exist, direct diagnostic evidence of this mechanism is elusive. However, a circumstantial case can be made by considering the alternatives. Although there are capture mechanisms (4, 5) to create well-separated secondaries such as some Kuiper belt binaries, they do not select a common orbital plane or direction. In addition, the subsequent hardening of these configurations tends to produce large eccentricities that could not be damped by tidal forces given the small masses of P1 and P2 (2). Alternatively, if a protosatellite disk were to extend to sufficient distance to allow the accretion of P1 and P2 in situ, there is no obvious reason why they should be found in near-resonant orbits, because tidal torques are also too weak to migrate them into such configurations. A final unanswered question is how the moons were initially trapped in corotation resonances. One possibility is that a small amount of vapor and/or small particles extended past the location of the 6:1 resonance (∼3.3 aC) and their free eccentricities damped by collisional viscosity. This could have initially populated many resonance sites, but most would be later cleared as eccentricities were excited by resonant migration. Indeed, material comprising P1 and P2 may have begun as ring arcs, except that by lying external to Pluto's Roche limit, single moons were able to accumulate.

Supporting Online Material

www.sciencemag.org/cgi/content/full/1127293/DC1

SOM Text

References and Notes

View Abstract

Stay Connected to Science

Navigate This Article