The RoPES project with HARPS and HARPS-N III. Two candidate planets orbiting the G-type star HD 161098
Pith reviewed 2026-05-10 15:43 UTC · model grok-4.3
The pith
Radial velocity data reveal two candidate planets orbiting the nearby G-type star HD 161098, one potentially in the optimistic habitable zone.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Analysis of radial velocity data from multiple high-precision spectrographs reveals two candidate planets orbiting the G-type star HD 161098. The inner candidate, HD 161098 b, has an orbital period of 72.578 days and a minimum mass of 3.63 Earth masses. The outer candidate, HD 161098 c, orbits every 682.5 days with a minimum mass of 7.8 Earth masses and, if real, would sit in the optimistic habitable zone. The study also determines the star's rotation period as 28.22 days and a magnetic cycle of 4090 days. These findings are obtained through careful separation of planetary signals from stellar activity and instrumental effects, paving the way for confirmation observations.
What carries the argument
Joint modeling of radial velocities with activity indicators to isolate Keplerian planetary signals from stellar activity.
If this is right
- If confirmed, the outer candidate would represent a super-Earth in the optimistic habitable zone of a nearby G-type star.
- The work achieves sub-meter-per-second radial velocity precision at orbital periods of hundreds of days.
- The analysis characterizes the star's rotation and long magnetic cycle, providing context for similar planet searches.
- Future targeted campaigns can collect more data to confirm or refute the planetary interpretation.
Where Pith is reading between the lines
- Transit photometry follow-up could test whether either candidate crosses the stellar disk and yield density constraints.
- Confirmation would add a benchmark system for studying super-Earth occurrence rates around G-type stars.
- The joint modeling approach may be applied to other long-baseline radial velocity datasets to uncover additional small planets.
Load-bearing premise
The two periodic radial velocity signals represent true planetary orbits rather than residual stellar activity or instrumental systematics.
What would settle it
Additional radial velocity observations spanning several years that fail to recover consistent signals at 72.6 or 682 days or show those signals correlating with activity indicators instead.
Figures
read the original abstract
Context. The development of refined instruments and techniques for a detailed analysis of the radial velocities (RVs) of stars other than the Sun allows the presence of planetary signals of amplitude below 1 ms$^{-1}$ to be investigated. Long-term RV surveys allow the detection of Earth-like and super-Earth-like planets in the habitable zones of Sun-like stars, prime targets for future missions for the atmospheric characterization of exoplanets. Aims. We present the analysis of the nearby G8 V-type star HD 161098 (V = 7.68 mag, d = 29.75 pc). We searched for terrestrial planets in the habitable zone. Methods. We combined historical datasets with new data collected in an ongoing blind search program with HARPS, HARPS-N, and ESPRESSO. We utilized recently developed tools to extract RVs and to deal with the analysis of stellar activity. We performed a joint analysis of RVs and activity indicators to separate the planetary signals from those related to activity. Results. We detected two sub-ms$^{-1}$ signals that we claim as candidate planets. We are not able to confirm their nature with certainty. Candidate HD 161098 b has an orbital period of 72.578$_{-0.060}^{+0.059}$ d and a minimum mass of 3.63 $\pm$ 0.59 MEarth. HD 161098 c has an orbital period of 682.5$_{-9.9}^{+9.5}$ d and a minimum mass of 7.8$_{-1.4}^{+1.5}$ MEarth.If confirmed, candidate HD 161098 c would reside in the optimistic habitable zone of the star. We find a magnetic cycle of 4090$_{-130}^{+140}$ d period and a rotation period of 28.22$_{-0.35}^{+0.30}$ d. Our analysis sets the stage for future observing campaigns of the star, finalized for the confirmation of our results. Conclusions. We are entering the sub-ms$^{-1}$ era at long orbital periods with a combination of stellar activity treatment and long-term campaigns.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes combined radial velocity time series for the G8 V star HD 161098 obtained with HARPS, HARPS-N, and ESPRESSO. A joint fit of the RVs with activity indicators is used to isolate two periodic signals interpreted as candidate planets: HD 161098 b (P = 72.578 d, m sin i = 3.63 M⊕) and HD 161098 c (P = 682.5 d, m sin i = 7.8 M⊕, in the optimistic habitable zone). The analysis also recovers a stellar rotation period of 28.22 d and a magnetic cycle of 4090 d. The authors present the signals as candidates only and call for additional observations to confirm their planetary nature.
Significance. If the two signals survive more stringent tests for residual activity, the result would demonstrate that sub-m s^{-1} amplitudes at periods of hundreds of days can be extracted from heterogeneous long-baseline RV data sets once activity indicators are modeled jointly. This would be a useful benchmark for the RoPES program and for future HZ searches around bright G dwarfs, even though the present work stops short of claiming confirmed planets.
major comments (2)
- [§4 (Results)] The central claim that the 72.578 d and 682.5 d signals are planetary rather than residual activity rests on the joint RV + activity-indicator modeling. The manuscript should provide quantitative diagnostics (e.g., posterior distributions or residual periodograms after subtraction of the best-fit activity model) showing that no significant power remains at the rotation period, its harmonics, or the 4090 d cycle; without these, the separation cannot be assessed at the reported sub-m s^{-1} amplitudes.
- [§3.3 (Joint modeling) and Table 2] The long-period candidate (682.5 d) lies near the optimistic HZ and has an amplitude comparable to plausible residual instrumental offsets across the three instruments. The manuscript should report the covariance between the planetary semi-amplitude and the instrument-specific zero-point offsets (or any additional linear/quadratic trends) to demonstrate that the 7.8 M⊕ mass is not driven by those systematics.
minor comments (3)
- [Abstract] The abstract states that the signals are 'sub-ms^{-1}' but does not quote the actual semi-amplitudes; these values should be added for immediate context.
- [Figure 5] Figure captions for the phase-folded RVs should explicitly state whether the plotted points have had the activity model subtracted or whether the full model is over-plotted.
- [§4.3] The uncertainty on the magnetic-cycle period (4090_{-130}^{+140} d) is reported, but the manuscript does not discuss whether this period is consistent with the length of the observational baseline or whether it could be an artifact of the window function.
Simulated Author's Rebuttal
We thank the referee for the positive assessment and constructive feedback on our manuscript. We agree that additional quantitative diagnostics will strengthen the presentation of our results. We have revised the manuscript accordingly and provide point-by-point responses below.
read point-by-point responses
-
Referee: [§4 (Results)] The central claim that the 72.578 d and 682.5 d signals are planetary rather than residual activity rests on the joint RV + activity-indicator modeling. The manuscript should provide quantitative diagnostics (e.g., posterior distributions or residual periodograms after subtraction of the best-fit activity model) showing that no significant power remains at the rotation period, its harmonics, or the 4090 d cycle; without these, the separation cannot be assessed at the reported sub-m s^{-1} amplitudes.
Authors: We agree with the referee that explicit diagnostics are important to demonstrate the robustness of the signal separation. In the revised manuscript, we have added a new figure showing the residual periodograms of the RVs after subtraction of the best-fit joint model (including activity indicators). These periodograms exhibit no significant peaks at the stellar rotation period (28.22 d), its first harmonic, or the magnetic cycle (4090 d). Furthermore, we include posterior distribution corner plots for the key activity parameters in an appendix to illustrate the constraints and lack of degeneracy with the planetary signals. These additions confirm that the planetary candidates remain the dominant signals at their respective periods with amplitudes well above the residual noise. revision: yes
-
Referee: [§3.3 (Joint modeling) and Table 2] The long-period candidate (682.5 d) lies near the optimistic HZ and has an amplitude comparable to plausible residual instrumental offsets across the three instruments. The manuscript should report the covariance between the planetary semi-amplitude and the instrument-specific zero-point offsets (or any additional linear/quadratic trends) to demonstrate that the 7.8 M⊕ mass is not driven by those systematics.
Authors: We appreciate this point regarding potential covariances with instrumental parameters. In the revised manuscript, we have expanded Table 2 to include the correlation coefficients between the semi-amplitude of HD 161098 c and the zero-point offsets for HARPS, HARPS-N, and ESPRESSO. The reported correlations are all below 0.25 in absolute value, indicating minimal covariance. We have also added a brief discussion in §3.3 noting that the planetary signal persists when the fit is performed with and without additional linear trends, and we provide the full covariance matrix excerpt in a new supplementary table. This supports that the derived minimum mass of 7.8 M⊕ is not an artifact of the instrumental modeling. revision: yes
Circularity Check
No circularity: planetary parameters obtained by direct least-squares/MCMC fitting after separate activity-indicator modeling
full rationale
The paper reports orbital periods and minimum masses as the direct output of a joint RV + activity-indicator fit (Keplerian signals plus GP or similar activity model). The reported periods (72.578 d, 682.5 d) are numerically distinct from the fitted stellar rotation (28.22 d) and cycle (4090 d), and no equation defines one in terms of the other. No self-citation is invoked to forbid alternative models or to rename a fitted parameter as an independent prediction. The derivation chain therefore remains self-contained against external benchmarks and does not reduce to any of the enumerated circular patterns.
Axiom & Free-Parameter Ledger
free parameters (5)
- orbital period of candidate b =
72.578 d
- minimum mass of candidate b =
3.63 MEarth
- orbital period of candidate c =
682.5 d
- minimum mass of candidate c =
7.8 MEarth
- magnetic cycle period =
4090 d
axioms (1)
- domain assumption Radial velocity time series can be decomposed into planetary Keplerian orbits plus stellar activity and noise using joint modeling of RVs and activity indicators
Reference graph
Works this paper leans on
-
[1]
A., Dumusque , X., & Halverson , S
Aigrain, S. & Foreman-Mackey, D. 2023, ARA&A, 61, 329 Aigrain, S., Pont, F., & Zucker, S. 2012, MNRAS, 419, 3147 Anderson, D. R., Collier Cameron, A., Hellier, C., et al. 2011, ApJ, 726, L19 Barragán, O., Gillen, E., Aigrain, S., et al. 2023, MNRAS, 522, 3458 Burt, J. A., Dumusque, X., & Halverson, S. 2025, arXiv e-prints, arXiv:2511.01954 Chaplin, W. J.,...
-
[2]
To mitigate the tendency of the GPs to overfit (Nari et al
toO(n), with n the number of measurements. To mitigate the tendency of the GPs to overfit (Nari et al. 2026), we can train them together with activity indicators in a multidimensional GP fashion (Rajpaul et al. 2015; Barragán et al. 2023). This framework was first inspired by the FF’ method (Aigrain et al. 2012). TheFF’method models the RV variation with ...
work page 2026
-
[3]
C.3: Zoomed-in image of the phase-folded plot of candidate planets of HD 161098
2 0 2 RV (m/s) 72.6 d planet HARPS 03 HARPS 15 HARPS-N ESPRESSO 682.5 d planet 0.0 0.5 1.0 Phase 4 2 0 2 4 Residuals RV (m/s) 0.0 0.5 1.0 Phase Fig. C.3: Zoomed-in image of the phase-folded plot of candidate planets of HD 161098. Table C.1: Significance of stellar activity-only mod- els when searching for signals with periods compat- ible with the periods...
work page 2000
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.