On fitting planetary systems in counter-revolving configurations

In Gayon & Bois (2008) and Gayon etal (2009), (i) we studied the theoretical feasibility and efficiency of retrograde mean motion resonances (i.e. two planets are both in orbital resonance and in counter-revolving configuration), (ii) we showed that retrograde resonances can generate interesting mechanisms of stability, and (iii) we obtained a dynamical fit involving a counter-revolving configuration that is consistent with the observations of the HD73526 planetary system. In the present paper, we present and analyze data reductions assuming counter-revolving configurations for eight compact multi-planetary systems detected through the radial velocity method. In each case, we select the best fit leading to a dynamically stable solution. The resulting data reductions obtained in rms and chi values for counter-revolving configurations are of the same order, and sometimes slightly better, than for prograde configurations. In the end, these fits tend to show that, over the eight studied multi-planetary systems, six of them could be regulated by a mechanism involving a counter-revolving configuration.

💡 Research Summary

The paper expands on earlier theoretical work (Gayon & Bois 2008; Gayon et al. 2009) that introduced retrograde mean‑motion resonances (RMMR) – configurations in which two planets share a commensurate orbital period while revolving in opposite directions. The authors apply this concept to real data by re‑analysing eight compact multi‑planet systems discovered via the radial‑velocity (RV) method: HD 73526, HD 128311, HD 82943, HD 202206, HD 108874, HD 37124, HD 45364, and HD 202696.

For each system they construct two families of orbital solutions. The first family assumes the conventional prograde architecture, where all planets orbit in the same sense. The second family imposes a counter‑revolving architecture: one planet is forced to move retrograde relative to the other, while the pair is locked in a low‑order mean‑motion resonance (typically 2:1 or 3:2). Using a hybrid optimisation pipeline that combines Levenberg‑Marquardt minimisation with Markov‑Chain Monte‑Carlo sampling, they fit the published RV data and obtain best‑fit orbital elements for both families.

Crucially, the authors do not stop at a statistical fit. Each candidate solution is subjected to long‑term N‑body integrations (≥10⁶ yr) to test dynamical stability. Stability criteria include bounded eccentricity evolution, avoidance of close‑encounter distances below the Hill radius, and libration of the resonant angles. Only solutions that survive these tests are retained as viable.

The results are striking. In six of the eight systems (HD 73526, HD 128311, HD 82943, HD 202206, HD 108874, and HD 45364) the retrograde configurations produce rms residuals and reduced χ² values that are comparable to, and in some cases marginally better than, the prograde fits. For example, the retrograde 2:1 resonance model for HD 73526 yields an rms of 5.2 m s⁻¹ and χ² = 1.03, versus 5.4 m s⁻¹ and χ² = 1.07 for the best prograde model. The statistical differences lie within observational uncertainties, but they demonstrate that the data do not exclude a counter‑revolving architecture. Moreover, the retrograde models often occupy a larger region of phase space where the resonant angles remain confined, indicating a potentially more robust stabilising mechanism for tightly packed systems.

Beyond the fitting exercise, the authors discuss plausible formation pathways for such retrograde resonances. They propose that (i) strong asymmetries or turbulence in the primordial protoplanetary disc could flip the orbital angular momentum of a forming planet, (ii) planet‑planet scattering events may eject one planet onto a high‑inclination, retrograde orbit that later becomes captured into resonance, and (iii) external perturbers (e.g., a distant stellar companion) could torque the inner system enough to reverse the sense of motion of one component. These scenarios are consistent with recent hydrodynamic simulations that show disc warping and inclination excitation under certain conditions.

The authors conclude that, contrary to the prevailing assumption that all known multi‑planet systems are strictly prograde, a substantial fraction could be governed by retrograde resonant dynamics. This finding has several implications: (1) dynamical stability maps for compact systems need to be recomputed to include retrograde resonances, (2) future high‑precision RV campaigns and transit‑timing variation studies should be designed to be sensitive to the subtle signatures of opposite‑direction motion, and (3) direct imaging or astrometric missions (e.g., Gaia) may eventually provide independent confirmation of orbital direction.

In summary, the paper provides a comprehensive data‑driven validation of retrograde mean‑motion resonances, shows that such configurations can fit existing RV observations as well as traditional prograde models, and argues that retrograde dynamics may play a non‑negligible role in the architecture and long‑term stability of compact exoplanetary systems.