Hic Sunt Dracones: Uncovering Dynamical Perturbers Within the Habitable Zone

The continuing exploration of neighboring planetary systems is providing deeper insights into the relative prevalence of various system architectures, particularly with respect to the solar system. However, a full assessment of the dynamical feasibility of possible terrestrial planets within the Habitable Zones (HZ) of nearby stars requires detailed knowledge of the masses and orbital solutions of any known planets within these systems. Moreover, the presence of as yet undetected planets in or near the HZ will be crucial for providing a robust target list for future direct imaging surveys. In this work, we quantify the distribution of uncertainties on planetary masses and semi-major axes for 1062 confirmed planets, finding median uncertainties of 11.1% and 2.2%, respectively. We show the dependence of these uncertainties on stellar mass and orbital period, and discuss the effects of these uncertainties on dynamical analyses and the locations of mean motion resonance. We also calculate the expected radial velocity (RV) semi-amplitude for a Neptune-mass planet in the middle of the HZ for each of the proposed Habitable Worlds Observatory target stars. We find that for more than half of these stars, the RV semi-amplitude is less than 1.5 m/s, rendering them unlikely to be detected in archival RV data sets and highlighting the need for further observations to understand the dynamical viability of the HZ for these systems. We provide specific recommendations regarding stellar characterization and RV survey strategies that work toward the detection of presently unseen perturbers within the HZ.

💡 Research Summary

The paper presents a comprehensive statistical analysis of the uncertainties associated with the masses and semi‑major axes of 1,062 confirmed exoplanets drawn from the NASA Exoplanet Archive (as of 22 August 2024). By extracting the reported uncertainties for planetary mass (or M sin i), semi‑major axis, orbital period, eccentricity, and host‑star mass, the authors quantify the typical precision of the two orbital parameters most critical for dynamical modeling. They find that planetary mass uncertainties have a mean of 16 % and a median of 11 % (range 0.2–622 %), while semi‑major‑axis uncertainties average 4.2 % with a median of 2.2 % (range 0.005–60 %). Host‑star mass uncertainties average 8.4 % (median 5.7 %).

The analysis links these uncertainties to their physical origins. For radial‑velocity (RV) detections, the mass error budget is dominated by the measurement error on the RV semi‑amplitude K and by the uncertainty on the stellar mass. When the stellar mass error exceeds roughly 20 %, it becomes the primary source of planetary‑mass uncertainty, eclipsing the contribution from K. Orbital‑period uncertainties play only a minor role in the mass error budget. In contrast, the semi‑major‑axis error is governed mainly by the period uncertainty once the period error surpasses about 10 %; beyond ≈5 AU, period errors dominate the a‑error, while stellar‑mass errors are less influential.

Figure 1 illustrates these trends, separating planets detected solely by RV (cyan circles) from those with both RV and transit data (maroon diamonds). The well‑known “mass gap” between 0.1–0.3 M_J appears, and the outlier with a 622 % mass error is the circumbinary planet Kepler‑47 b, whose mass is only loosely constrained. Figure 2 shows that host‑star visual magnitude (V) does not strongly correlate with either mass or a‑uncertainty, indicating that the additional constraints provided by transits (precise period and epoch) offset the lower signal‑to‑noise of fainter stars.

The authors then discuss the implications for dynamical simulations. Large uncertainties in Mₚ and a translate directly into ambiguous locations of mean‑motion resonances (MMRs) and uncertain stability boundaries. For low‑mass planets (Mₚ ≤ 0.1 M_J), the median mass uncertainty rises to ~15 %, meaning that even modest perturbations could be mischaracterized in long‑term integrations. Consequently, improving host‑star mass determinations (e.g., via asteroseismology or Gaia astrometry) and refining orbital periods (through extended RV baselines) are essential for reliable dynamical assessments of habitable‑zone (HZ) stability.

A key part of the study focuses on the upcoming Habitable Worlds Observatory (HWO) mission. Using the proposed list of 164 HWO target stars (Mamajek & Stapelfeldt 2024), the authors compute the expected RV semi‑amplitude K for a hypothetical Neptune‑mass planet (≈17 M⊕) placed at the midpoint of each star’s HZ. The resulting K distribution shows that more than half of the targets would produce signals below 1.5 m s⁻¹. Since most archival RV data achieve precisions of ~1 m s⁻¹, such planets would remain undetected, especially around active or faint stars where stellar jitter dominates.

To address this detection gap, the paper proposes three complementary strategies:

1. Stellar Characterization: Reduce host‑star mass uncertainties to ≤5 % using asteroseismic scaling relations, high‑resolution spectroscopy, and Gaia DR4 luminosity constraints.

2. High‑Precision, Long‑Baseline RV Surveys: Deploy state‑of‑the‑art spectrographs (e.g., ESPRESSO, NEID, EXPRES) to reach ≤0.5 m s⁻¹ precision over multi‑year baselines, thereby uncovering low‑amplitude signals in the HZ.

3. Synergistic Astrometry: Leverage Gaia’s micro‑arcsecond astrometry to detect or constrain long‑period companions, providing independent orbital constraints that tighten RV fits and improve dynamical models.

The authors conclude that while current exoplanet catalogs already provide a solid statistical foundation, the uncertainties in planetary masses and semi‑major axes—especially those driven by host‑star properties—still limit our ability to assess the dynamical habitability of nearby systems. Targeted improvements in stellar mass determination, extended high‑precision RV monitoring, and the integration of astrometric data are essential steps before the HWO can reliably select truly habitable worlds for direct imaging.