{"id":"54cad6c7-8310-4735-aba8-defd62be3753","arxiv_id":"2505.18941","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Transit-timing data for TOI-2109b rule out fast orbital decay and favor a slow decay rate of a few milliseconds per year consistent with a young host star.","lead":"TOI-2109b, an ultra-hot Jupiter with a 16-hour orbit, is decaying very slowly, if at all, according to new transit-timing measurements. The data disfavor rapid orbital decay and support a young host star, but the signal is only marginal.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Old/young discrimination rests on an unquantified M_crit threshold from Barker (2020); if M_crit exceeds 5 M_J at 2.65 Gyr, the TTV constraint no longer favors a young host.","rationale":"The reader's weakest_assumption correctly identifies the IGW wave-breaking threshold as the linchpin of the age interpretation. I considered two alternative concerns: (1) the TTV fit depending on the single SuperWASP-N point, and (2) unmodeled dynamical TTV sources (companion, oblateness, GR) biasing Pdot. The first is not load-bearing: with a 4-year TESS + CHEOPS baseline and per-transit precisions of tens of seconds, a 10 ms/yr decay would accumulate ~43 s of shift, giving a Pdot precision of order 1–2 ms/yr even without the 2006–2011 SuperWASP-N point; the old-star 1107 ms/yr would be ruled out by TESS sectors 52 and 79 alone, as the paper notes. The second is mitigated because a 0.2 M_J companion producing 1–2 min TTVs would show up as large residuals in the quadratic fit; the small reported Pdot uncertainty implies such signals are not present at that amplitude, though a joint fit would still be a useful check. The IGW threshold, by contrast, is explicitly used to split the theoretical predictions into two scenarios differing by more than two orders of magnitude, and the paper does not quantify the uncertainty in M_crit or the age at which the transition occurs. The authors themselves flag the old-star results as needing caution. The paper's strengths are real: the new 20-s cadence TESS sector 79 data, the use of the established PdotQuest code, the reproducibility of the timing data on GitHub, and the consistency with Harre et al. (2024). Those support the observational Pdot constraint. But the leap from 'Pdot is small' to 'the host star is young' requires the critical-mass threshold, so that is the single most load-bearing assumption. The reader's CONDITIONAL verdict is appropriate; specifying the M_crit scaling and exploring continuous IGW efficiency would directly address it.","tokens_in":27760,"tokens_out":13088,"duration_ms":119662,"concrete_test":"Recompute M_crit for a 1.45 M_sun star (with TOI-2109's metallicity, if available) at ages 1.0–3.0 Gyr using the Barker (2020) nonlinear wave-breaking criterion, e.g., via MESA stellar models. If M_crit(2.65 Gyr) > 5 M_J, the old-star fast-decay scenario is invalidated and the TTV data do not single out a young host; if M_crit(1.77 Gyr) < 5 M_J, the young scenario may also enter the breaking regime, further weakening the age claim. As a minimal check, request the authors to state the exact M_crit value and uncertainty they read from fig. 9 of Barker (2020) for this star.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the TTV data 'support a young host star' depends on the old-star scenario predicting Pdot_TIDE ≈ −1107 ms/yr (Section 4.2). That prediction is activated only if IGWs are fully damped, which requires the planet mass (5.02 ± 0.75 M_J) to exceed the critical mass M_crit for wave breaking. The authors assert, citing fig. 9 of Barker (2020), that M_crit at t_age = 2.65 Gyr (the 1σ upper age limit) lies between 1 and 10 M_J, so a ~5 M_J planet would trigger full damping. If M_crit at this age is actually >5 M_J, or if the true age is below the threshold, the old-star scenario reverts to inefficient IGW dissipation (Q'_IGW ≈ 10^8.5), and the predicted decay rate becomes comparable to the young-star value (~4 ms/yr). In that case the measured Pdot_OBSV = −2.6 ± 1.3 ms/yr would no longer discriminate between young and old host stars; the conclusion would reduce to a Pdot upper limit. The paper neither quotes M_crit explicitly nor quantifies its dependence on stellar age/mass, and the authors themselves caution that the old-star results 'should be taken with a grain of salt' (Section 6). The observational TTV constraint itself is valuable and likely robust, but the abstract's age interpretation is only as secure as this binary threshold.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the orbital evolution of the ultra-hot Jupiter TOI-2109b using a two-layer tidal model that includes inertial waves in the convective envelope and internal gravity waves in the radiative region, with stellar age as the key parameter. For a 'young' host star the model predicts a slow period decay of approximately -4 ms/yr, while for an 'old' host star (t_age = 2.65 Gyr) it predicts a much faster decay of about -1100 ms/yr. The authors combine TESS data from three sectors (including new 20-s cadence data from Sector 79), CHEOPS transit times from the literature, and ground-based light curves to perform a TTV analysis with the PdotQuest code. They find a best-fit period decay of Pdot_OBSV = (-2.616 ± 1.285) ms/yr under a 3-sigma rejection scheme, which rules out the high end of the 10-740 ms/yr range suggested by Wong et al. (2021) and yields a 95% lower limit of Q'_* > 3.7 × 10^7. The authors interpret this as support for a young host star and a constant-period orbit, and they also simulate TTV contributions from an outer companion, stellar/planetary oblateness, and general relativity to aid interpretation.","tokens_in":28095,"tokens_out":7006,"duration_ms":60335,"significance":"The observational TTV result is a solid contribution: it incorporates new TESS 20-s cadence data, re-reduces earlier TESS and ground-based transits, and reports a consistent small period decay under two sigma-clipping schemes. This is a genuine, model-independent timing constraint that tightens the lower limit on Q'_* and rules out the fast-decay end of previous predictions. The paper also makes its timing data and reduction code publicly available, which supports reproducibility. If the young-star interpretation holds, the result implies that TOI-2109b is not currently spiraling into its host star at an observable rate and that stellar age is a controlling factor in tidal dissipation efficiency. The theoretical tidal modeling, however, rests on two under-specified inputs—the Barker (2020) critical-mass threshold and the scaling of the two-layer inertial-wave model—so the central age interpretation is less secure than the timing measurement itself.","major_comments":[{"comment":"The division into 'young' and 'old' host-star scenarios relies on an unquantified critical-mass threshold from Barker (2020). The authors state that at t_age = 2.65 Gyr (the 1-sigma upper age limit in Table 1) M_crit lies between 1 and 10 M_J, so a ~5 M_J planet fully damps IGWs; they do not quote M_crit itself nor its dependence on stellar mass and age. If M_crit exceeds 5 M_J at this age, or if the true age is below 2.65 Gyr, the 'old' star reverts to Q'_IGW ≈ 10^8.5 and Pdot_TIDE ≈ 4 ms/yr, making it observationally indistinguishable from the 'young' case. Since the abstract's claim that TTVs 'support a young host star' depends entirely on this binary threshold, the authors should compute or quote M_crit from Barker (2020) and show how it varies across the allowed age range, and present the old-star prediction as a conditional scenario rather than a fixed alternative.","section":"§4.2 and Table 1"},{"comment":"The two-layer inertial-wave model is said to underpredict dissipation in F-type stars by 1-2 orders of magnitude, and the authors account for the difference by using Barker (2020) to scale Q'_*, but no numerical scale factor, equation, or methodology is provided. As a result, the quoted Pdot_TIDE = (-4.186 ± 0.797) ms/yr for the young-star scenario is not reproducible, and its uncertainty almost certainly understates the model error from the unspecified scaling. This matters because the central comparison in Section 5.3 between Pdot_TIDE and Pdot_OBSV = (-2.616 ± 1.285) ms/yr is used to support the young-star interpretation; a factor of a few in the scaling would change the theoretical prediction by an order of magnitude and could spoil the agreement.","section":"§3.1 and §4.1"},{"comment":"The conversion from the fitted Pdot_OBSV to the headline lower limit Q'_*,OBSV > 3.7 × 10^7 is not shown. Unlike Section 4.1, where the Goldreich & Soter relation and adopted stellar/planetary parameters are described, Section 5.3 simply states the result. The authors should present the formula used to convert a measured period derivative into Q'_* and list the assumed values of M_p, R_*, a, and the stellar moment of inertia, so that the reader can verify the one-sided 95% lower limit. This is a central quantitative claim of the paper.","section":"§5.3"}],"minor_comments":[{"comment":"There are a few typographical errors: 'With a a cadence' in Section 5.3 and 'evolutionaryS timescales' in Section 6.","section":"§5.3 and §6"},{"comment":"The statement that the result 'rules out any Pdot > 10 ms/yr' is imprecise. The 3-sigma interval on Pdot is approximately [-6.5, 1.3] ms/yr, so a decay rate of 10 ms/yr is excluded at about 5.7 sigma; the formal significance should be stated instead.","section":"§5.3"},{"comment":"When combining TESS, CHEOPS, and heterogeneous ground-based data, the analysis does not quantify correlated red noise or possible systematic offsets between instruments. The reported uncertainties on Pdot would be more robust if a jitter term or per-instrument offset were included; at least a brief discussion of this limitation should be added.","section":"§5.1 and §5.2"},{"comment":"The two sigma-rejection fits use different samples (53 vs 102 transits); it would aid interpretation to plot both data sets and the best quadratic curves with residuals in a single panel, in addition to the TTV plots shown.","section":"Figure 9"},{"comment":"The simulated TTV amplitudes from the companion, oblateness, and general relativity are presented without a quantitative comparison to the observed TTV residuals. Even a simple chi-square or amplitude comparison would clarify whether these effects are consistent with the data or are merely illustrative.","section":"§5.4"}],"recommendation":"major_revision","confidential_remarks":"The observational core of the paper is strong and fits the journal's scope; the new TESS data and the reproducible TTV analysis are valuable. The main risk is the overinterpretation of the Barker (2020) critical-mass threshold in the age discrimination, which is fixable with a more transparent presentation and a sensitivity test. I would invite a revision. There is also a pattern of self-citation to the authors' earlier tidal evolution work, which is not inappropriate but could be trimmed if it is not directly necessary."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the observational TTV analysis is solid, and the theoretical wrapper around it is shaky. The genuinely new piece is TESS Sector 79's 20-second cadence data, which adds 22 mid-transit times to a four-year baseline. The timing analysis is careful: two sigma-clipping schemes agree, the inclusion of CHEOPS and ground-based points is reasonable, and the resulting Pdot = -2.6 ± 1.3 ms/yr robustly excludes decay rates above roughly 10 ms/yr. The lower limit Q'_star > 3.7e7 is defensible, and the paper ships its data and code, which makes the result reproducible.\n\nThe soft spots are all in the tidal interpretation. The two-layer inertial-wave calculation is scaled to Barker's F-star results by hand, and the scaling factor is never quantified. More importantly, the 'old star' scenario that produces ~1100 ms/yr only activates if IGWs are fully damped, which requires the planet to exceed a critical mass M_crit taken from fig. 9 of Barker (2020). The authors never quote M_crit, and if M_crit at 2.65 Gyr is actually above ~5 M_J, the old-star decay rate drops to something comparable to the young-star value. In that case the TTV data would just give an upper limit, not a young-host discrimination. The authors themselves wave this away in Section 6 ('should be taken with a grain of salt'), so the abstract overstates the age conclusion.\n\nNone of this hurts the measured Pdot. The observed constraint is independent of the tidal machinery and will be useful for anyone modeling ultra-hot Jupiter evolution. The hypothetical companion TOI-2109c is clearly labeled as a test case, not a detection, so that is fine. The paper deserves a serious referee; the referee should ask for the Q' scaling to be written out and for the M_crit threshold to be either quantified or removed from the headline claim.","headline":"New TESS timing data give a solid upper limit on TOI-2109b's orbital decay, but the paper's young-versus-old star conclusion leans on an unquantified wave-breaking threshold.","tokens_in":28719,"tokens_out":2188,"would_cite":true,"duration_ms":21357,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"TOI-2109b, the shortest-period ultra-hot Jupiter known, shows transit-timing variations consistent with a constant orbital period rather than rapid tidal decay.","keywords":["ultra-hot Jupiter","orbital decay","transit timing variations","tidal dissipation","inertial waves","internal gravity waves","stellar tidal quality factor","TOI-2109b"],"falsifier":"Continue high-cadence transit timing of TOI-2109b for several more years. If the true decay were in the 10–740 ms/yr range, the cumulative mid-transit-time shift over the 2020–2028 baseline would reach tens of minutes to hours, vastly exceeding current uncertainties of roughly two minutes, so an absence of such drift would confirm the constant-period conclusion. Independently, a direct stellar age measurement (for example by asteroseismology or gyrochronology) would settle whether the star is young enough for the slow-decay branch to apply.","tokens_in":2210,"feed_emoji":"🪐","tokens_out":4606,"duration_ms":87349,"temperature":0.7,"pith_summary":"This paper argues that the ultra-hot Jupiter TOI-2109b—the shortest-period hot Jupiter known, with a 16-hour orbit—is not currently spiraling into its host star at a detectable rate. Re-analyzing transit timings from TESS, CHEOPS, and ground-based telescopes spanning 2020–2023, the authors find a best-fit period decay of −2.616 ± 1.285 ms/yr, consistent with a constant-period orbit and ruling out the 10–740 ms/yr range predicted when the planet was discovered. The key insight is that the efficiency of tidal dissipation in the F-type host star depends sharply on the star's age: inertial waves in the convective envelope dominate for a young star, while internal gravity waves in radiative regions would break and dissipate strongly only if the star is near its upper age limit. The data favor the young-star branch, implying a stellar tidal quality factor Q'_* > 3.7×$10^{7}$ and a decay slow enough that mid-transit times shift by less than about 10 seconds over three years. If correct, TOI-2109b becomes a testbed for quiet tidal evolution rather than a candidate for imminent orbital death.","feed_headline":"Shortest-period hot Jupiter isn't spiraling inward after all","feed_subtitle":"Transit-timing data cap the orbital decay near 3 ms/yr, ruling out the fast-decay prediction.","key_machinery":"The engine of the argument is the coupled set of tidal evolution equations (constant time-lag model) with frequency-averaged tidal quality factors computed from two-layer interior models: inertial waves (IWs) in the convective envelope, internal gravity waves (IGWs) in radiative regions, and a viscoelastic solid core for the planet. The decisive element is the stellar-age dependence: internal gravity waves are efficiently damped only if the planet exceeds a critical mass $M_\\mathrm{crit}$ that decreases with stellar age, so the same $\\sim5\\,M_\\mathrm{J}$ planet produces either slow decay (young star, $Q'_{\\star,\\mathrm{IGW}} \\approx 10^{8.5}$) or fast decay (old star, $Q'_{\\star,\\mathrm{IGW}} \\propto (P_\\mathrm{tide}/0.5\\,\\mathrm{d})^{8/3}$). The transit-timing measurements themselves are carried by a quadratic-ephemeris fit to 102 mid-transit times from TESS, CHEOPS, and ground-based follow-up, using a dedicated period-derivative fitting procedure.","core_discovery":"The central claim is that transit-timing variations of TOI-2109b across four years of data favor a rather constant-period orbit, with a best-fit $\\dot{P} = (-2.616 \\pm 1.285)$ ms/yr at $3\\sigma$, which excludes the fast-decay range of 10–740 ms/yr proposed in the discovery paper and supports a young host star with $Q'_\\star > 3.7\\times10^7$. In the authors' tidal evolution model, a young star dissipates tides mainly through inertial waves in the convective envelope, giving $\\dot{P} \\approx -4.2$ ms/yr, whereas an old star whose internal gravity waves reach the wave-breaking regime would decay at about $-1107$ ms/yr—a rate that would have produced easily visible shifts in the 2022 and 2024 TESS transits and is rejected by the data. The paper also shows that gravitational perturbations from a possible outer companion, stellar and planetary oblateness, and general-relativistic precession can each generate TTV signals that mimic or mask orbital decay, so those must be subtracted before attributing any trend to tides.","pith_inferences":["The paper leaves stellar age as the dominant unknown; a precise age determination would break the young/old degeneracy and either confirm the slow decay or reopen the fast-decay scenario.","The assumption of a solid planetary core with rigidity derived from Jupiter-Io tidal interaction influences the planetary dissipation term; updated Jupiter interior models with a diluted core could shift the modeled $Q'_p$ and slightly alter the predicted $\\dot{P}$.","If the suggested $\\sim 0.2\\,M_\\mathrm{J}$ outer companion exists, its ~1–2 minute periodic TTV signal should be detectable in a few more CHEOPS or TESS sectors, offering a direct test independent of the decay question."],"forward_implications":["If the period is truly nearly constant, TOI-2109b is not on the verge of tidal disruption, and its 16-hour orbit can be used to probe stellar tidal quality factors near $Q'_\\star \\sim 10^7$ rather than a fast-decay endpoint.","The young-star scenario implies $Q'_\\star > 2.3\\times10^7$ from tidal modeling and $>3.7\\times10^7$ from timing, values at the high end for F stars, which constrains how efficiently inertial waves dissipate in rapidly rotating convective envelopes.","Under the accepted decay rate, mid-transit times advance by less than about 10 seconds over three years, a signal detectable with high-cadence space photometry and a concrete target for future campaigns.","The demonstration that a possible outer companion, oblateness, and relativistic precession can produce TTV signals resembling decay means future analyses of ultra-short-period planets must model those effects before claiming orbital decay."],"supporting_citations":[{"why":"Discovery paper that provides the system parameters, initial ephemeris, and the 10–740 ms/yr decay prediction that this work's timing analysis rules out.","marker":"I. Wong et al. (2021)"},{"why":"Supplies the inertial-wave and internal-gravity-wave dissipation formalisms, the $Q'_{\\star,\\mathrm{IGW}} \\propto (P/0.5\\,\\mathrm{d})^{8/3}$ scaling, and the critical-mass threshold that separates the young and old host-star scenarios.","marker":"A. J. Barker (2020)"},{"why":"Provides independent CHEOPS mid-transit times, a $\\dot{P} \\sim -1$ ms/yr measurement, and the outer-companion hypothesis used for comparison and TTV decomposition.","marker":"J. V. Harre et al. (2024)"},{"why":"The transit-timing analysis code (IV22) used to extract TESS mid-transit times and detrend the light curves.","marker":"E. S. Ivshina & J. N. Winn (2022)"},{"why":"Provides PdotQuest, the package used to fit the quadratic ephemeris and derive the best-fit period derivative.","marker":"W. Wang et al. (2024)"},{"why":"Establishes the frequency-averaged inertial-wave dissipation formalism in convective envelopes that the two-layer model applies to the star.","marker":"G. I. Ogilvie (2013)"},{"why":"Gives the planetary two-layer model with a solid core, including the viscoelastic dissipation term and aspect-ratio parameters adopted for TOI-2109b.","marker":"M. Guenel et al. (2014)"},{"why":"Defines the tidal quality factor $Q'$ and the relation used to convert measured decay rates into lower limits on $Q'_\\star$.","marker":"P. Goldreich & S. Soter (1966)"}],"fun_headline_variants":["Data rule out fast decay for shortest-period hot Jupiter","TOI-2109b's orbit decays at just a few ms per year","Fast orbital decay ruled out for 16-hour exoplanet","Ultra-hot Jupiter's orbit is steadier than predicted"],"cache_read_input_tokens":30592,"weakest_assumption_plain":"The argument that the data single out a young host star rests on the wave-breaking threshold: the authors assume a 5-Jupiter-mass planet fully damps internal gravity waves only when the star is near its 2.65 Gyr upper age limit, because the critical mass for wave breaking from the adopted tidal model falls between 1 and 10 Jupiter masses at that age; if that threshold were different for a 1.45-solar-mass F star, the old-star scenario would not necessarily decay at the fast rate, and the timing data would no longer favor a young host.","fun_headline_variants_meta":{"raw":{"variants":["Data rule out fast decay for shortest-period hot Jupiter","TOI-2109b's orbit decays at just a few ms per year","Fast orbital decay ruled out for 16-hour exoplanet","Ultra-hot Jupiter's orbit is steadier than predicted"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001637,"raw_usage":{"total_tokens":6625,"prompt_tokens":1178,"completion_tokens":5447,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":794,"completion_tokens_details":{"reasoning_tokens":5375}},"tokens_in":794,"tokens_out":5447,"duration_ms":29428,"temperature":1.0,"reasoning_tokens":5375,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:23:28.792325+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Continue high-cadence transit timing of TOI-2109b for several more years. If the true decay were in the 10–740 ms/yr range, the cumulative mid-transit-time shift over the 2020–2028 baseline would reach tens of minutes to hours, vastly exceeding current uncertainties of roughly two minutes, so an absence of such drift would confirm the constant-period conclusion. Independently, a direct stellar age measurement (for example by asteroseismology or gyrochronology) would settle whether the star is young enough for the slow-decay branch to apply.","supporting_citations":[{"cited_title":"V., Smith, A","cited_arxiv_id":null,"evidence_quote":"Provides independent CHEOPS mid-transit times, a $\\dot{P} \\sim -1$ ms/yr measurement, and the outer-companion hypothesis used for comparison and TTV decomposition."},{"cited_title":"2014, A&A, 566, L9, doi: 10.1051/0004-6361/201424010","cited_arxiv_id":null,"evidence_quote":"Gives the planetary two-layer model with a solid core, including the viscoelastic dissipation term and aspect-ratio parameters adopted for TOI-2109b."}],"review_version":1}