The Astrometric Resoeccentric Degeneracy: Eccentric Single Planets Mimic 2:1 Resonant Planet Pairs in Astrometry

TL;DR

The paper identifies an astrometric analogue of the radial-velocity resoeccentric degeneracy: a single mildly eccentric planet can mimic two coplanar circular planets in a 2:1 resonance in Gaia-like astrometry. It derives a harmonic decomposition showing fundamental and first-harmonic signals, maps these to an effective eccentricity $e_{\rm eff}$, and connects $e_{\rm eff}$ to planet mass and period via $e_{\rm eff} = 2\, \frac{M_{p,2}}{M_{p,1}} \left(\frac{P_2}{P_1}\right)^{2/3}$, with a strict 2:1 case giving $\frac{M_{p,2}}{M_{p,1}} = 2^{-1/3} e_{\rm eff}$. Simulations with Gaia DR4/DR5-like data show the degeneracy can bias eccentricity inferences for coplanar 2:1 pairs, while mutual inclination disrupts the degeneracy, revealing the true multi-planet architecture. The work discusses implications for Gaia-era exoplanet demographics, resonant planet occurrence, and dynamical histories, and emphasizes the need for robust multi-modal inference and combined astrometry–RV follow-up to accurately characterize long-period giant planets. It also highlights how Hill stability bounds and resonant dynamics shape the plausible parameter space for detectable systems and suggests strategies to break the degeneracy in future analyses.

Abstract

Detections of long-period giant exoplanets will expand dramatically with Gaia Data Release 4 (DR4), but interpreting these signals will require care. We derive the astrometric resoeccentric degeneracy: an astrometric analogue of the well-known radial velocity degeneracy in which a single eccentric planet can mimic two circular planets near a 2:1 period ratio. To first order in eccentricity, the sky-projected motion of a single eccentric orbit decomposes into a fundamental mode and first harmonic with an amplitude proportional to that eccentricity. A pair of coplanar, circular planets in a 2:1 orbital resonance produces the same harmonic structure: the outer planet sets the fundamental mode, while the inner planet supplies an apparent first harmonic. We present a mapping between the harmonic amplitudes and effective eccentricity ($e_\mathrm{eff}$) of a single planet that mimics a 2:1 configuration, demonstrating that $e_\mathrm{eff} = \, 2^{1/3}(M_{p,2}/M_{p,1})$, the masses of the inner and outer planets, respectively. Using simulated Gaia data we show that (1) coplanar 2:1 systems are statistically indistinguishable from a single eccentric planet and (2) mutual inclination can break this degeneracy. This bias favors detecting mutually inclined systems, often fingerprints of a dynamically hot history -- traces for processes such as planet-planet scattering or secular chaos. Determining the planetary architectures in which this degeneracy holds will be essential for measuring cool-giant occurrence rates with Gaia and for inferring their dynamical evolution histories.

Figures (2)

• Figure 1: Gaia DR5-like simulated orbits (vertical line shows Gaia DR4 baseline), coplanar system ($T_\mathrm{span}=5.5\,\mathrm{yr}$, $i_\mathrm{mut}=0^\circ$): 2:1 system with $P_1=5.2\,\mathrm{yr}$, $M_{p,1}=12\,M_\mathrm{Jup}$, $e_1=0.0$, $P_2=2.6\,\mathrm{yr}$, $M_{p,2}=2.84\,M_\mathrm{Jup}$, $e_2=0.0$, compared to a single-planet model with $P=5.2\,\mathrm{yr}$, $e=0.3$, $M_p=12\,M_\mathrm{Jup}$. Star assumed to be solar-mass at 50 pc. Here, $\Delta\alpha\cos\delta = \Delta\alpha^\ast$ is the RA and $\Delta\delta$ is the declination of the stellar reflex motion projected onto the sky plane of a star; $\Delta\eta$ is the Gaia along–scan (AL) coordinate; Res are the residuals between the simulated Gaia-like data with the eccentric model and the 2:1 model, respectively.
• Figure 2: Gaia DR5-like simulated orbits (vertical line shows Gaia DR4 baseline), mutually inclined system ($T_\mathrm{span}=5.5\,\mathrm{yr}$, $i_\mathrm{mut}\approx45^\circ$): 2:1 system with $P_1=5.2\,\mathrm{yr}$, $M_{p,1}=12\,M_\mathrm{Jup}$, $e_1=0.0$, $P_2=2.6\,\mathrm{yr}$, $M_{p,2}=2.84\,M_\mathrm{Jup}$, $e_2=0.0$, compared to a single-planet model with $P=5.2\,\mathrm{yr}$, $e=0.3$, $M_p=12\,M_\mathrm{Jup}$. Star assumed to be solar-mass at 50 $pc$. Here, $\Delta\alpha\cos\delta = \Delta\alpha^\ast$ is the RA and $\Delta\delta$ is the declination of the stellar reflex motion projected onto the sky plane of a star; $\Delta\eta$ is the Gaia along–scan (AL) coordinate; Res are the residuals between the simulated Gaia-like data with the eccentric model and the 2:1 model, respectively.