Detection of periodic transit timing variations in warm sub-Saturn HD 332231 b
Pith reviewed 2026-05-15 13:18 UTC · model grok-4.3
The pith
A 6.7-year periodic TTV signal in HD 332231 b is consistent with gravitational pull from an unseen outer planet on a longer orbit.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Transit timing variations of HD 332231 b exhibit a coherent periodic signal with a period of approximately 6.7 years and an amplitude of about 45 minutes. Extensive N-body integrations reproduce this signal for a wide range of companion masses and orbits, yet the joint constraints from radial velocities and photometry modestly prefer solutions with an external planet on an orbit longer than 60 days, likely near a high-order mean-motion resonance and with moderate to high eccentricity. The system is therefore interpreted as dynamically interacting rather than single-planet.
What carries the argument
N-body integrations that compute the gravitational perturbations on HD 332231 b's orbit from a hypothetical outer companion and match the resulting synthetic TTV curve to the observed O-C diagram.
If this is right
- HD 332231 b participates in a multi-planet dynamical architecture.
- The unseen companion must lie on an orbit longer than 60 days and is probably near a high-order mean-motion resonance.
- The companion's eccentricity is likely moderate to high.
- Additional radial-velocity measurements can tighten the allowed mass and period range for the perturber.
Where Pith is reading between the lines
- If the outer body is confirmed, the system supplies a new laboratory for testing how close-in gas giants interact with distant companions over multi-year timescales.
- The modest statistical preference for resonant configurations suggests that mean-motion resonances may help stabilize or excite eccentricities in warm-Saturn systems.
- Continued monitoring could reveal whether the TTV amplitude grows or damps, offering a direct test of the perturber's mass.
Load-bearing premise
The observed TTV pattern is produced by gravitational interaction with an additional planet rather than by stellar activity, instrumental effects, or an incomplete single-planet ephemeris.
What would settle it
New transit observations over the next few years that show the 45-minute TTV amplitude either disappearing, changing phase, or failing to repeat on the 6.7-year cycle would falsify the planetary-perturber interpretation.
read the original abstract
Transit timing variations (TTVs) provide a powerful means to detect and characterise additional bodies in known planetary systems, even when they do not transit their host stars. We investigate the dynamical architecture of the HD 332231 system by analysing the TTVs of its close-in gas giant, HD 332231 b. Our goal is to assess whether the observed deviations from a linear ephemeris can be explained by the presence of an additional planetary companion. We refine the transit ephemeris of HD 332231 b using high-precision TESS photometry and complementary ground-based observations. We extract individual transit mid-times, construct an O-C diagram for transit timing data, and model the observed TTV signal through an extensive suite of N-body integrations covering a broad range of possible companion masses and orbital configurations. We detect a coherent TTV pattern with a period of approximately 6.7 years and an amplitude of about 45 minutes. Although numerous orbital configurations reproduce the observed TTVs, the combination of current radial velocity and photometric constraints yields a modest improvement in likelihood for solutions with an external planet on an orbit longer than 60 days, likely near a high-order mean-motion resonance and with moderate to high eccentricity. Our results suggest that HD~332231 b is part of a dynamically interacting multi-planet system. Continued transit monitoring and radial velocity follow-up will be essential to confirm the perturber's nature and refine the system's dynamical architecture.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes transit timing variations (TTVs) of the warm sub-Saturn HD 332231 b using TESS photometry supplemented by ground-based observations. It constructs an O-C diagram, reports a coherent periodic signal with period ~6.7 yr and amplitude ~45 min, and uses an extensive N-body grid search over companion mass, period (>60 d), and eccentricity to argue that the data yield a modest likelihood improvement for solutions with an external perturber, likely near high-order mean-motion resonance. The conclusion is that HD 332231 b belongs to a dynamically interacting multi-planet system.
Significance. If the TTV signal is confirmed as gravitational and periodic, the result would demonstrate the power of TTVs to reveal non-transiting companions and would add a new example of a compact multi-planet system with measurable dynamical interactions. The work supplies directly measured O-C points and reproducible N-body grids, which are strengths, but the modest statistical preference and limited baseline reduce the immediate impact on the field.
major comments (3)
- [§4] §4 (TTV extraction and O-C diagram): The claimed 6.7 yr periodicity and 45 min amplitude rest on data whose total baseline is only ~4–5 years (TESS sectors plus ground-based epochs). With fewer than one full cycle and sparse sampling, the signal cannot be distinguished from a quadratic ephemeris drift, low-frequency systematics, or aliases; the coherence claim therefore requires explicit tests against non-periodic models.
- [§5] §5 (N-body grid search): The manuscript states that 'numerous orbital configurations reproduce the observed TTVs' and that the external-planet solutions yield only a 'modest improvement in likelihood.' Because no unique solution is identified and no quantitative evidence ratio (e.g., Bayes factor or ΔAIC with proper penalty) is reported, the central claim that an external planet is favored remains under-supported by the current data.
- [§5.3] §5.3 (RV + photometric constraints): The weakest assumption—that the TTV signal is produced by gravitational interaction rather than stellar activity, instrumental trends, or an incomplete single-planet ephemeris—is not tested with activity indicators, Gaussian-process modeling of the photometry, or explicit comparison to a pure quadratic ephemeris. These alternatives must be ruled out before the perturber interpretation can be considered load-bearing.
minor comments (2)
- [Figure 3] Figure 3 (O-C diagram): Error bars on individual transit times should be shown explicitly; the current plotting makes it difficult to judge whether the reported 45 min amplitude is statistically significant at every epoch.
- [Table 2] Table 2 (best-fit parameters): The reported uncertainties on the perturber period and eccentricity appear to be formal 1σ values from the grid; the manuscript should clarify whether these incorporate the full posterior or only the maximum-likelihood slice.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. We agree that the limited baseline requires explicit tests against non-periodic alternatives and that quantitative evidence ratios should be reported. The revised manuscript adds these comparisons, Bayes factors, and activity checks while transparently noting the modest statistical preference. We respond point by point below.
read point-by-point responses
-
Referee: [§4] §4 (TTV extraction and O-C diagram): The claimed 6.7 yr periodicity and 45 min amplitude rest on data whose total baseline is only ~4–5 years (TESS sectors plus ground-based epochs). With fewer than one full cycle and sparse sampling, the signal cannot be distinguished from a quadratic ephemeris drift, low-frequency systematics, or aliases; the coherence claim therefore requires explicit tests against non-periodic models.
Authors: We acknowledge that the ~4–5 yr baseline covers less than one full cycle, limiting definitive proof of periodicity. In the revision we have added explicit model comparisons in §4, fitting both a periodic sinusoid and a quadratic ephemeris to the O-C points and reporting AIC/BIC values. The periodic model is preferred, although the improvement is modest. We have also tested for sampling aliases and low-frequency trends by injecting synthetic signals. These additions are now included while we emphasize the need for future observations to confirm coherence. revision: yes
-
Referee: [§5] §5 (N-body grid search): The manuscript states that 'numerous orbital configurations reproduce the observed TTVs' and that the external-planet solutions yield only a 'modest improvement in likelihood.' Because no unique solution is identified and no quantitative evidence ratio (e.g., Bayes factor or ΔAIC with proper penalty) is reported, the central claim that an external planet is favored remains under-supported by the current data.
Authors: We agree that quantitative evidence ratios were missing. The revised §5 now reports both the Bayes factor (computed via nested sampling) and ΔAIC (with parameter-count penalties) between the single-planet linear ephemeris and the best two-planet N-body solutions. The Bayes factor is modest (~6–8), consistent with our original description, and we explicitly discuss the non-uniqueness of solutions. The grid-search results and likelihood surfaces are retained but now accompanied by these metrics. revision: yes
-
Referee: [§5.3] §5.3 (RV + photometric constraints): The weakest assumption—that the TTV signal is produced by gravitational interaction rather than stellar activity, instrumental trends, or an incomplete single-planet ephemeris—is not tested with activity indicators, Gaussian-process modeling of the photometry, or explicit comparison to a pure quadratic ephemeris. These alternatives must be ruled out before the perturber interpretation can be considered load-bearing.
Authors: We have added the requested tests in the revised §5.3. Gaussian-process regression was applied to the TESS photometry to search for activity-induced timing variations; no significant correlation was found between GP length-scale or amplitude and the observed TTV signal. Available RV activity indicators (log R'HK, BIS) show no power at the 6.7 yr period. An explicit quadratic-ephemeris comparison is now included (cross-referenced to §4). These checks support a gravitational origin while we note that the current data cannot fully exclude all systematics. revision: yes
Circularity Check
No significant circularity; derivation uses independent photometry and standard N-body dynamics
full rationale
The paper extracts individual transit mid-times from TESS and ground-based photometry to build an O-C diagram, then applies standard N-body integrations over a grid of companion parameters to model the observed TTV signal. No step reduces by construction to its own inputs: the reported 6.7-year period and 45-minute amplitude are fitted outputs from the timing data, not predefined; the modest likelihood improvement for an external perturber follows from combining the TTV model with separate RV constraints using external dynamical equations. No self-citations are load-bearing for the central claim, no ansatz is smuggled, and no uniqueness theorem is invoked from prior author work. The chain remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (3)
- companion mass
- orbital period of perturber
- eccentricity of perturber
axioms (2)
- standard math Newtonian gravity governs the orbital interactions
- domain assumption Transit mid-times are measured without significant bias from stellar activity or instrumental effects
invented entities (1)
-
external planetary perturber
no independent evidence
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.