A Quasi-Stationary Solution to Gliese 436bs Eccentricity

We investigate the possibility that the large orbital eccentricity of the transiting Neptune-mass planet Gliese 436b is maintained in the face of tidal dissipation by a second planet in the system. We find that the currently observed configuration can be understood if Gliese 436b and a putative companion have evolved to a quasi-stationary fixed point in which the planets’ orbital apses are co-linear and in which secular variations in the orbital eccentricities of the two planets have been almost entirely damped out. We adopt an octopole-order secular theory based on a Legendre expansion in the semi-major axis ratio to delineate well-defined regions of (P_c, M_c, e_c) space that can be occupied by a perturbing companion. We incorporate the evolutionary effect of tidal dissipation into our secular model of the system, and solve the resulting initial value problems for a large sample of the allowed configurations. We then polish the stationary configurations derived from secular theory with full numerical integrations. We present our results in the form of candidate companion planets to Gliese 436b. For these candidates, radial velocity half-amplitudes, K_c, are of order 3 m/s, and the maximum amplitude of orbit-to-orbit transit timing variations are of order Delta t=1 s to Delta t=5s. For the particular example case of a perturber with orbital period, P_c=40 d, mass, M_c=8.5 M_Earth, and eccentricity, e_c=0.58, we confirm our semi-analytic calculations with a full numerical 3-body integration of the orbital decay that includes tidal damping and spin evolution. Additionally, we discuss the possibility of many-perturber stationary configurations, utilizing modified Laplace-Lagrange secular theory.

💡 Research Summary

Gliese 436b is a transiting Neptune‑mass planet (≈23 M⊕) on a 2.64‑day orbit with an unusually high eccentricity (e≈0.15). Standard tidal theory predicts that such an eccentricity should be damped to near‑circular values on timescales far shorter than the system’s age, yet the observed eccentricity persists. This paper investigates whether a second, unseen planet can maintain the eccentricity by forcing the system into a quasi‑stationary secular fixed point. In this configuration the apsidal lines of both planets are aligned (co‑linear) and the secular exchange of eccentricity is essentially halted, so that tidal dissipation can no longer efficiently circularize the inner orbit.

The authors adopt an octupole‑order secular theory, expanding the disturbing function in powers of the semi‑major‑axis ratio a_b/a_c. While the classic Laplace‑Lagrange (quadrupole) approximation suffices for low eccentricities and very hierarchical systems, the octupole term becomes dominant when the outer companion has a substantial eccentricity (e_c ≳ 0.4) and a moderate period ratio (P_c/P_b ≈ 10–30). The octupole contribution introduces a non‑linear coupling that can produce a stable apsidal alignment and a fixed point in the (e_b, e_c) phase space.

To incorporate tidal effects, the authors embed a constant‑time‑lag tidal model with a planetary quality factor Q′≈10⁵–10⁶ and allow the inner planet’s spin to evolve toward synchronous rotation. They then integrate the coupled secular‑tidal equations for a large ensemble of initial conditions, spanning a wide range of possible companion periods (P_c), masses (M_c), and eccentricities (e_c). The integrations are carried out over 10⁸–10⁹ yr, long enough to assess whether the system can settle into the quasi‑stationary state and remain there for Gyr timescales.

The numerical experiments delineate a well‑defined region of parameter space where a companion can sustain Gliese 436b’s eccentricity. Viable companions have orbital periods between roughly 30 and 70 days, masses of 5–12 M⊕, and eccentricities of 0.4–0.7. Within this region the secular fixed point is robust: the inner planet’s eccentricity settles at e_b≈0.12–0.15 and varies only by a few percent thereafter. The predicted radial‑velocity semi‑amplitude of such companions is modest, K_c≈2–4 m s⁻¹, placing them at the detection threshold of current high‑precision spectrographs (e.g., HARPS, ESPRESSO). The associated transit‑timing variations (TTVs) are small, Δt≈1–5 s per orbit, but could be measured with space‑based photometry given sufficient cadence and baseline.

A concrete example is explored in detail: a companion with P_c = 40 days, M_c = 8.5 M⊕, and e_c = 0.58. The authors first locate the fixed point using the octupole secular model, then verify it with a full three‑body N‑body integration that includes tidal dissipation and spin evolution. The full integration reproduces the semi‑analytic predictions, confirming that the system remains near the fixed point for >10⁸ yr and that the inner planet’s eccentricity is effectively frozen.

Beyond the single‑companion scenario, the paper discusses the possibility of multi‑planet stationary configurations. By extending the Laplace‑Lagrange formalism to include higher‑order terms and mutual secular interactions among several low‑mass companions, the authors show that the collective mode spectrum can also produce long‑lived apsidal alignments. In such cases the parameter space for maintaining Gliese 436b’s eccentricity widens, suggesting that a more complex planetary architecture could be responsible.

In summary, the study provides a self‑consistent dynamical framework that couples octupole secular interactions with tidal dissipation, demonstrating that a modest‑mass, moderately eccentric outer planet can lock Gliese 436b into a quasi‑stationary state and preserve its observed eccentricity over gigayear timescales. The work delivers concrete observational predictions—radial‑velocity amplitudes of a few meters per second and sub‑second TTVs—that are within reach of present‑day instrumentation, offering a clear pathway to test the hypothesis and to potentially discover the hidden companion(s) shaping this intriguing exoplanet system.