DYNAMICS AND ORIGIN OF THE 2:1 ORBITAL RESONANCES OF THE GJ 876 PLANETS

摘要

The discovery by Marcy and coworkers of two planets in 2:1 orbital resonance about the star GJ 876 has been supplemented by a dynamical fit to the data by Laughlin & Chambers, which places the planets in coplanar orbits deep in three resonances at the 2:1 mean-motion commensurability. The selection of this almost singular state by the dynamical fit means that the resonances are almost certainly real, and with the small amplitudes of libration of the resonance variables, indefinitely stable. Several unusual properties of the 2:1 resonances are revealed by the GJ 876 system. The libration of both lowest order mean-motion resonance variables and the secular resonance variable, θ_1 = λ_1 - 2λ_2 + w_1. θ_2 = λ_1 - 2λ_2 + w_2, and θ_3 = w_1 - w_2, about 0° (where λ_(1, 2)( are the mean longitudes of the inner and outer planet and w_(1, 2) are the longitudes of periapse) differs from the familiar geometry of the Io-Europa pair, where θ_2 and θ_3 librate about 180°. By considering the condition that w_1 = w_2 for stable simultaneous librations of θ_1 and θ_2, we show that the GJ 876 geometry results from the large orbital eccentricities e_i, whereas the very small eccentricities in the lo-Europa system lead to the latter's geometry. Surprisingly, the GJ 876 configuration, with θ_1, θ_2, and θ_3 all librating, remains stable for e_1 up to 0.86 and for amplitude of libration of θ_1 approaching 45° with the current eccentricities―further supporting the indefinite stability of the existing system. Any process that drives originally widely separated orbits toward each other could result in capture into the observed resonances at the 2:1 commensurability. We find that forced inward migration of the outer planet of the GJ 876 system results in certain capture into the observed resonances if initially e_1 approx< 0.06 and e_2 approx< 0.03 and the migration rate |a_2/a_2| approx< 3 x 10~(-2)(a_2/AU)~(-3/2) yr~(-1). Larger eccentricities lead to likely capture into higher order resonances before the 2:1 commensurability is reached. The planets are sufficiently massive to open gaps in the nebular disk surrounding the young GJ 876 and to clear the disk material between them, and the resulting planet-nebular interaction typically forces the outer planet to migrate inward on the disk viscous timescale, whose inverse is about 3 orders of magnitude less than the above upper bound on | a_2/a_2| for certain capture. If there is no eccentricity damping, eccentricity growth is rapid with continued migration within the resonance, with et exceeding the observed values after a further reduction in the semimajor axes a_i of only 7%. With eccentricity damping e_i/e_i = -K|a_i/a_i|, the eccentricities reach equilibrium values that remain constant for arbitrarily long migration within the resonances. The equilibrium eccentricities are close to the observed eccentricities for K ≈ 100 if there is migration and damping of the outer planet only, but for K ≈ 10 if there is also migration and damping of the inner planet. This result is independent of the magnitude or functional form of the migration rate a_i as long as e_i/e_i = -K|a_i/a_i|. Although existing analytic estimates of the effects of planet-nebula interaction are consistent with this form of eccentricity damping for certain disk parameter values, it is as yet unclear that such interaction can produce the large value of K required to obtain the observed eccentricities. The alternative eccentricity damping by tidal dissipation within the star or the planets is completely negligible, so the observed dynamical properties of the GJ 876 system may require an unlikely fine-tuning of the time of resonance capture to be near the end of the nebula lifetime.