{"id":"1f00449c-fe9e-442e-9220-16e405756468","arxiv_id":"2506.15022","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The orbit of hot Jupiter WASP-4 b is decaying at a rate corresponding to a tidal quality factor Q'_* ≈ 80,000, which only the most massive, oldest, and largest stellar models can explain.","lead":"Astronomers combined 216 transit timings, including 37 new observations, and confirmed that the hot Jupiter WASP-4 b is slowly spiraling into its host star, with its orbital period shrinking by a few milliseconds per year. The measured dissipation is faster than most stellar models predict, echoing the mystery of the similar planet WASP-12 b and pointing to gaps in our understanding of tides in Sun-like stars.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The unmodeled 11.8 s EXOFAST-vs-literature timing offset (Sec. 3.3.1) is the main threat: a ~10 s-level zero-point drift over 16 yr could bias the quadratic coefficient, whose total cumulative signal is only ~160 s.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: time-correlated systematics in the TTV dataset could mimic or bias the parabolic curvature. The paper's own numbers make this concrete: the 11.8 s mean offset between the authors' EXOFAST timings and published timings for 109 identical light curves, with 98/109 in the same direction, is direct evidence of a pipeline-level zero-point problem, and the preferred quadratic fit has reduced chi^2 = 3.61, indicating underestimated scatter. A constant 11.8 s offset would not produce curvature, but the root cause is explicitly unidentified, so an epoch-dependent component cannot be excluded. The cumulative signal is only ~160 s over 16 years, so a drift of order 10 s is not negligible. The own-only fit weakens the concern that the result comes from mixing two timing populations, but it does not eliminate the possibility of a time-dependent bias within the homogeneous set. Previous independent analyses by Bouma et al. (2019), Turner et al. (2022), Harre et al. (2023), and Ma et al. (2025) also find decay, which is independent support for the qualitative conclusion; nevertheless, the precise decay rate and the inferred Q'_* depend on systematics that the paper has not resolved. The CONDITIONAL verdict is therefore appropriate, and no verdict change is needed. The proposed test directly checks whether the zero-point is constant in epoch and whether the decay signal is present in instrument-by-instrument subsets, which would settle whether the concern actually lands.","tokens_in":27987,"tokens_out":8429,"duration_ms":98231,"concrete_test":"For the 109 light curves with both EXOFAST and published timings, compute Δt(E) = T_EXOFAST - T_published and fit constant, linear, and quadratic models in epoch. If the slope or quadratic term is significant at >2σ, the zero-point drifts with time and the fitted coefficient a must be corrected. Independently, fit the quadratic model separately to each homogeneous subset (Danish-only, TESS-only, ETD-only, other literature-only, and each individual telescope), requiring that the negative quadratic coefficient appears consistently in each subset and is mutually consistent within errors. Finally, recompute the model comparison after inflating all timing errors by sqrt(chi^2_nu) to confirm that the preference for the quadratic model and the significance of a survive a conservative treatment of the underestimated scatter.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the assumption, implicit in Sec. 3.3.2, that the 216-point TTV diagram is free of epoch-correlated zero-point systematics. The paper itself reports (Sec. 3.3.1, Fig. 3) that its homogeneous EXOFAST v1 timings are systematically later than published timings of the same 109 light curves by 11.8 s on average, with 98/109 positive residuals, and states that the root cause is unidentified. The preferred quadratic coefficient a = (-9.81 ± 1.21) × 10^-11 d orbit^-2 corresponds to only ~160 s of cumulative curvature over the 15.98 yr baseline; a zero-point drift of order 10 s correlated with the changing mix of instruments and literature-only timings could therefore shift a by tens of percent, and a more adversarial drift could in principle mimic the signal. The own-only fit (a = -8.52 ± 1.47 × 10^-11) shows the result is not an artifact of simply combining the two timing populations, but it does not rule out a time-dependent bias inside the homogeneous EXOFAST pipeline itself. A simple late-epoch own/literature mixture would oppose the observed sign, so the offset is not automatically damning; the danger is a non-constant drift not captured by the two-set comparison. This is the weakest link in the chain from timings to Q'_* = (8.08 ± 1.17) × 10^4.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents 37 new transit light curves of WASP-4 b from several ground-based telescopes plus TESS data, re-analyzes an extensive set of literature and ETD light curves with EXOFAST v1 to derive homogeneous mid-transit timings, and combines these with 58 literature timings into a 216-point TTV diagram spanning 15.98 yr. The authors fit linear, quadratic, and apsidal models; the quadratic model is strongly preferred (ΔBIC = -51.96, Bayes factor ~2×10^11), yielding a negative quadratic coefficient and an implied modified tidal quality factor Q'_* = (8.08±1.17)×10^4 under the orbital-decay interpretation. They then compute theoretical Q'_* from internal gravity wave dissipation in MESA stellar models, finding that main-sequence models give Q'_* ~ 2-5×10^5, while only more massive/evolved models (with radii larger than observed) can reach the observed value.","tokens_in":28253,"tokens_out":12883,"duration_ms":130199,"significance":"The TTV dataset is the largest assembled for WASP-4 b, and the model comparison is thorough. If robust, the inferred Q'_* ≈ 8×10^4 makes WASP-4 b one of the best examples of tidal orbital decay and challenges current tidal theory, which the paper addresses with independent MESA-based calculations. The theoretical analysis is a strength: it uses reproducible stellar modeling (MESA) and an established formalism (Barker 2020), and it transparently states the tension between the observed Q'_* and main-sequence models. However, the astrophysical conclusion rests on the quadratic coefficient being free of time-correlated systematic errors in the timing data.","major_comments":[{"comment":"The systematic offset between the authors' EXOFAST v1 timings and published timings for the same 109 light curves (average 11.8 s, 98/109 positive, including 43/46 for TESS) is unexplained and is the main threat to the central claim. Because the cumulative quadratic signal over the 16-yr baseline is only ~160 s, a time-dependent component of this offset could bias the fitted quadratic coefficient in sign or magnitude. The authors should explicitly test for a time trend in the offset (e.g., comparing residuals in early vs. late epochs), include a free constant offset between the re-measured and literature-only subsets in the TTV model, or demonstrate that the result persists when using only a single homogeneous timing set (e.g., TESS-only or new ground-based data). The own-only fit in Table 6 reduces but does not eliminate this concern, because the re-measured subset itself combines heterogeneous instruments and epochs.","section":"Sec. 3.3.1, Fig. 3"},{"comment":"The reduced chi-squared values of 3.61 (all data) and 4.14 (own-only) for the preferred quadratic model indicate that the reported timing uncertainties are underestimated by about a factor of two. The quoted 1σ errors on the quadratic coefficient (a = -9.81±1.21 × 10^-11 days/orbit^2) and on Q'_* must therefore be regarded as lower limits. A re-analysis including a jitter term (or an explicit rescaling of the errors) is needed to obtain realistic uncertainties on a and Q'_* and to confirm that the model comparison (ΔAIC, ΔBIC, Bayes factor) is not driven by the underestimated scatter. This is directly relevant to the theoretical comparison in Section 4, which depends on whether Q'_* is consistent with ~6×10^4 or could be as high as ~2×10^5.","section":"Sec. 3.3.2, Table 6"}],"minor_comments":[{"comment":"Many instances of '9P' and '9M' (e.g., in the Introduction and Section 4) are LaTeX rendering errors for \\dot{P} and \\dot{M}; they should be fixed before publication.","section":"Throughout"},{"comment":"The reported numbers of new observations are inconsistent: the abstract states 37 new observations, Table 1 lists 28 Danish + 5 Ckoirama + 2 TRAPPIST + 4 El Sauce = 39 new light curves, and the conclusion says '37 newly obtained observations ... of which three were excluded.' Please clarify the exact counts and which telescopes contributed.","section":"Sec. 2.6 / Sec. 5"},{"comment":"The phrase 'the average difference between the measurements from the same light curves is only 11.83 seconds' is misleading, since 11.8 s is comparable to the typical 25-27 s uncertainties and is a substantial systematic; it should be flagged as a concern rather than 'only'.","section":"Sec. 3.3.1"},{"comment":"The y-axis label reads 'TTV (mi utes)'; the 'n' in 'minutes' is missing.","section":"Fig. 4"},{"comment":"When dismissing the 387.96-day periodic signal, the authors state that a Keplerian fit 'ends up within much poorer fit statistics' but do not give the resulting AIC/BIC or amplitude; please quantify this statement.","section":"Sec. 3.3.2"}],"recommendation":"major_revision","confidential_remarks":"The qualitative conclusion (orbital decay with Q'_* ~ 8×10^4) is consistent with several independent previous studies, so the scientific direction is probably correct. The new analysis, however, must address the unexplained systematic timing offset and the underestimated uncertainties before the result can be considered quantitative. This is a fixable issue within the manuscript's scope; I would not reject."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The orbital decay of WASP-4 b is not new - Bouma, Turner, Harre, and Ma all reported it. What this paper adds is the largest homogeneous TTV dataset (216 points, 37 new transits) and an honest, detailed attempt to confront the inferred tidal quality factor with MESA-based internal gravity wave models. The model comparison is clean: the quadratic ephemeris beats linear and apsidal models by a wide margin (BIC ~ -52, Bayes factor ~2x10^11). The own-only fit also gives a negative quadratic coefficient, so the signal is not merely an artifact of splicing together heterogeneous literature timings.\n\nThe authors are transparent about their data selection, light curve quality cuts, and the fact that their EXOFAST timings run systematically ~12 s later than published values for the same 109 light curves. That level of candor is welcome. The soft spots are real but not fatal. The ~12 s offset, averaged over the whole baseline, corresponds to a nontrivial fraction of the ~160 s cumulative curvature signal. If that offset drifts with time or with the changing mix of instruments, it could bias the quadratic coefficient by tens of percent. The reduced chi-square of 3.6 for the quadratic fit also says the timing errors are underestimated; the quoted Q'_* = (8.08 +/- 1.17) x 10^4 should be read as precise but not accurate.\n\nAnd Q'_* is an algebraic transform of the fitted coefficient, so it is not an independent check - the paper knows this. The IGW calculation is genuinely independent, and that is the most valuable part of the discussion: ordinary main-sequence models give Q' ~2-5x10^5, and only older, more massive, larger-radius models reach the observed ~8x10^4, echoing the WASP-12 situation. The theoretical caveats are clearly stated.\n\nMy take: the orbital decay is almost certainly real - too many independent groups see it - but the exact decay rate and the tension with stellar models should be treated as provisional until the timing systematics are understood. This paper deserves a serious referee, not a desk rejection. The referee should ask the authors to either explain the 12 s offset or at least test for a time-varying zero-point, and to propagate the systematic uncertainty into Q'_*. A thoughtful revision would make this a solid reference for the WASP-4 system.","headline":"Confirms WASP-4 b orbital decay with the largest TTV dataset yet; the 12 s systematic timing offset is a real concern but not disqualifying.","tokens_in":679,"tokens_out":1661,"would_cite":true,"duration_ms":32093,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The orbit of hot Jupiter WASP-4 b is shrinking: a quadratic timing model beats linear and apsidal-motion alternatives by a Bayes factor near $2\\times10^{11}$, implying a stellar tidal quality factor $Q'_\\star \\approx 8\\times10^4$.","keywords":["transit timing variations","orbital decay","hot Jupiter","WASP-4 b","tidal quality factor","internal gravity waves","TTV diagram","tidal dissipation"],"falsifier":"A concrete test is to continue monitoring WASP-4 b transits for another 5 to 10 years: if the quadratic coefficient is real orbital decay, the timing residuals should continue along the same parabola, with the period decreasing at the predicted rate; if the curvature flattens, reverses, or jumps, the decay interpretation fails. A shorter test is to re-reduce the shared light curves with a single common pipeline and check whether the 11.8-second systematic offset varies with time or source; if it drifts coherently across epochs, the fitted quadratic coefficient could be an artifact.","tokens_in":27741,"feed_emoji":"🪐","tokens_out":10191,"duration_ms":97266,"temperature":0.7,"pith_summary":"This paper argues that the hot Jupiter WASP-4 b is undergoing true orbital decay, not apsidal precession or an unseen companion. Using the largest transit-timing dataset for the system, 216 mid-transit times spanning 16 years, the authors find that a quadratic ephemeris is overwhelmingly preferred over a linear one ($\\Delta\\mathrm{BIC}=-51.96$, Bayes factor $\\approx 2\\times10^{11}$) and over an apsidal-motion model. The fitted quadratic coefficient gives a period change of about $-2\\times10^{-10}$ days per orbit, corresponding to a modified tidal quality factor $Q'_\\star=(8.08\\pm1.17)\\times10^4$ if stellar tides drive the decay. The paper then tests this value against tidal theory and finds that internal gravity-wave damping in the star's radiative core can reach the required efficiency only if the host star is near the end of its main-sequence life, with a slightly larger mass and radius than the best-fit observed parameters. A sympathetic reader would care because WASP-4 b is a primary case for tidal orbital decay, and the result puts a concrete number on a longstanding theoretical mismatch.","feed_headline":"WASP-4 b's orbit is in decay, 216 transits show","feed_subtitle":"A quadratic timing model beats linear and precession alternatives, pointing to a star that dissipates tides fast.","key_machinery":"The central object is the transit-timing variation (TTV) diagram: 216 mid-transit times of WASP-4 b plotted against orbital epoch, assembled from 158 homogeneously modelled light curves (37 new from ground-based telescopes plus TESS, ETD, and literature data) and 58 literature timings. The load-bearing comparison is between three parametric models fitted to this diagram, linear ephemeris, quadratic ephemeris, and apsidal motion, with model selection via AIC, BIC, and the resulting Bayes factor. The quadratic coefficient $a$ is the mechanism: it directly measures the secular period change through the relation $dP/dE=2a$ and, via the formula of Goldreich and Soter and of Patra et al., maps to the modified tidal quality factor $Q'_\\star$. On the theory side, the paper computes $Q'_\\star$ for gravity waves launched at the radiative-convective interface and fully damped in the core, using stellar-structure models; the Brunt–Väisälä frequency profile near that interface controls the dissipation efficiency. The mismatch between the observed and predicted $Q'_\\star$ is the argument's engine: it forces the conclusion that only an evolved, slightly overluminous host can explain the decay.","core_discovery":"The central claim is that WASP-4 b's transit times deviate from a constant-period ephemeris with a curvature best described by a negative quadratic term in the epoch, $T(E)=T_0+P_0E+aE^2$ with $a=(-9.81\\pm1.21)\\times10^{-11}$ days per orbit squared, equivalently $(-0.98\\pm0.12)\\times10^{-10}$ days per orbit squared. The quadratic model is strongly preferred in AIC, BIC, and Bayes-factor comparisons, and no significant periodic signal survives in the residuals; the 387.96-day peak found in a periodogram does not phase-fold and fails as a Keplerian. The authors therefore adopt orbital decay as the only viable explanation, convert the curvature to a tidal quality factor, and confront it with stellar models. Their main-sequence models predict $Q'_\\star\\approx(2\\text{--}5)\\times10^5$ from internal gravity waves, about two to four times too large (too little dissipation); matching the observed $Q'_\\star\\approx8\\times10^4$ requires an older, more massive, slightly larger host near the end of the main sequence, which does not quite agree with the observationally inferred radius. On the paper's own terms, the discovery is that WASP-4 b's orbit is decaying and that the decay rate demands a more evolved host star than the nominal parameters suggest.","pith_inferences":["A decisive extension would be to apply the authors' own pipeline to every archival light curve and test whether the 11.8-second systematic offset depends on epoch, telescope, or filter; a drifting offset would mimic part of the quadratic curvature and should be subtracted before trusting $a$.","If the decay is real, WASP-4 b joins a small sample of hot Jupiters whose measured tidal dissipation exceeds standard main-sequence predictions; comparing decay rates across this sample could reveal whether the common resolution is an evolved host star or a missing dissipation mechanism in the theory.","The near-absence of persistent periodic TTV signals, despite a 387.96-day peak that fails to phase-fold, suggests that the remaining scatter is dominated by noise or spot-induced asymmetries; high-cadence, multi-year monitoring at the level of a few seconds would be needed to detect the expected small signal from the proposed 7001-day companion."],"forward_implications":["If the quadratic model is right, the observed TTVs are dominated by a steady secular shrinkage of the orbit, and mechanisms like apsidal precession, line-of-sight acceleration, or a wide companion cannot explain them.","The implied tidal quality factor $Q'_\\star\\approx8\\times10^4$ places WASP-4 b among the most efficiently tidally dissipating hot-Jupiter hosts, comparable to WASP-12 b.","Standard main-sequence stellar models dissipate only about a third to a quarter as much tidal energy via internal gravity waves, so explaining the decay requires the host to be near the end of the main sequence with a slightly larger radius than observed.","Future transit monitoring can distinguish the decay model from systematic errors: the 11.8-second systematic offset between the authors' homogeneous timings and published timings is smaller than the accumulated TTV but must not drift with time."],"supporting_citations":[{"why":"Established the orbital decay model for WASP-4 b's TTVs, ruled out spot-crossing, apsidal precession, and external perturbers, and supplies the prior system parameters.","marker":"Bouma et al. (2019)"},{"why":"Added 22 transit light curves and provided the reference linear ephemeris used to construct the TTV diagram.","marker":"Southworth et al. (2019)"},{"why":"Contributed TESS sectors 28 and 29 timing measurements and the RV-based companion model used to exclude a perturber explanation.","marker":"Turner et al. (2022)"},{"why":"Added CHEOPS light curves and independently found the orbital decay model superior to apsidal and Keplerian models.","marker":"Harre et al. (2023)"},{"why":"Re-analysis of the radial velocities that rules out the line-of-sight acceleration model and estimates $Q'_\\star\\approx60000$.","marker":"Baluev et al. (2020)"},{"why":"Found the quadratic TTV trend superior to the linear model with $\\Delta\\mathrm{BIC}=6.7$, motivating continued monitoring.","marker":"Maciejewski (2022)"},{"why":"Provides the Eq. 41 formalism for computing $Q'_\\star$ from internal gravity-wave dissipation in stellar interiors.","marker":"Barker (2020)"},{"why":"Supplies the standard formula relating the tidal quality factor to the period derivative.","marker":"Goldreich & Soter (1966)"},{"why":"Applies the tidal quality factor formula to hot Jupiter decay, giving the estimator used in Eq. 4.","marker":"Patra et al. (2017)"},{"why":"The modelling suite used to measure all mid-transit timings homogeneously from the selected light curves.","marker":"Eastman et al. (2013)"}],"fun_headline_variants":["WASP-4 b's orbit decays; star must be more evolved","216 transits show WASP-4 b's orbit is shrinking","WASP-4 b's decay rate points to an evolved host star","Orbital decay of WASP-4 b demands a more evolved star","WASP-4 b's decaying orbit hints at an older star"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole argument hinges on the assumption that the 216 measured transit times, assembled from many telescopes and sources over 16 years, carry no time-correlated systematic error that mimics the parabolic curvature, for example a drift in the 11.8-second systematic offset the authors find between their homogeneous timings and published timings of the same light curves.","fun_headline_variants_meta":{"raw":{"variants":["WASP-4 b's orbit decays; star must be more evolved","216 transits show WASP-4 b's orbit is shrinking","WASP-4 b's decay rate points to an evolved host star","Orbital decay of WASP-4 b demands a more evolved star","WASP-4 b's decaying orbit hints at an older star"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00067,"raw_usage":{"total_tokens":3136,"prompt_tokens":1108,"completion_tokens":2028,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":724,"completion_tokens_details":{"reasoning_tokens":1933}},"tokens_in":724,"tokens_out":2028,"duration_ms":14445,"temperature":1.0,"reasoning_tokens":1933,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:46:40.264370+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test is to continue monitoring WASP-4 b transits for another 5 to 10 years: if the quadratic coefficient is real orbital decay, the timing residuals should continue along the same parabola, with the period decreasing at the predicted rate; if the curvature flattens, reverses, or jumps, the decay interpretation fails. A shorter test is to re-reduce the shared light curves with a single common pipeline and check whether the 11.8-second systematic offset varies with time or source; if it drifts coherently across epochs, the fitted quadratic coefficient could be an artifact.","supporting_citations":[{"cited_title":"V., et al., 2020, @doi [ ] 10.1093/mnrasl/slaa069 , https://ui.adsabs.harvard.edu/abs/2020MNRAS.496L..11B 496, L11","cited_arxiv_id":null,"evidence_quote":"Re-analysis of the radial velocities that rules out the line-of-sight acceleration model and estimates $Q'_\\star\\approx60000$."}],"review_version":2}