{"id":"88896e88-e13e-4731-b309-159cab348f7a","arxiv_id":"2501.08992","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The slope of the hot Jupiter decay-time distribution implies the stellar tidal quality factor scales as Q'* proportional to P^alpha with alpha between -4.33 and -2, favoring resonance locking.","lead":"Hot Jupiters slowly spiral into their host stars, and the speed depends on how efficiently the star dissipates tidal energy. A new population-level analysis measures this efficiency and shows it drops steeply with orbital period, pointing to resonance locking.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inferred alpha range is not uniquely determined by the model; it depends on matching the observed power-law slope at a single logarithmic midpoint (tau_d = 0.03 Gyr) where the analytic steady-state slope is known to vary across the fitted interval, and the quoted error bars exclude this…","rationale":"The reader identified the steady-state/source-distribution assumption as the weakest assumption. My concern is closely related but more specific: the analytic comparison in Fig. 6 evaluates the steady-state slope at one log-space midpoint, while the observed and simulated slopes are power-law fits over the entire interval. Because the analytic slope is not constant over that interval, equating the two involves a systematic that is not propagated. This is a concrete, load-bearing technical issue with the alpha inference, distinct from (though related to) the source-distribution uncertainty. The Q0 normalization also rests on an assumed break age. The paper is honest about caveats and the qualitative result is supported by independent P18 consistency, so the verdict should remain CONDITIONAL rather than ACCEPT or REJECT. The reader's weakest_assumption partially overlaps with my concern, so agreement is partial.","tokens_in":28839,"tokens_out":1168,"duration_ms":12851,"concrete_test":"Recompute the inference using a consistent slope definition on both sides of the comparison: fit the analytic steady-state solution (Eq. 16) over the full fitting interval 0.001 < tau_d < 0.3 Gyr for each alpha, exactly as done for the simulated histograms, and compare that interval-averaged p_analytic(alpha) with the observed p. If the resulting alpha range shifts by more than the quoted bootstrap width, the quoted [-4.33, -2] range is not robust to the slope-definition choice. Also repeat the Q0 calibration with break ages of 1 and 3 Gyr instead of 2 Gyr and record the shift in the Q0 range; if the shift is comparable to a decade in Q0, the absolute normalization is not tightly constrained.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that the observed decay-time slope in 0.001 < tau_d < 0.3 Gyr equals the analytic steady-state slope p(alpha) from Eq. 17, giving alpha in [-4.33, -2]. That identification is not uniquely determined by the model. Eq. 17 gives p as a function of tau_d/tau_d,p, and the paper evaluates it at the midpoint tau_d = 0.03 Gyr with tau_d,p = 1 Gyr. But the analytic slope is not constant over the fitted interval: for alpha = 0 it varies from p = 1 at tau_d -> 0 to p = 0.84 at tau_d/tau_d,p = 0.03, and the simulations in Section 2 fit the whole interval with a single power law of p = 0.86. Thus the observed 'slope' is a fit over the full interval while the analytic comparison uses one point; the mismatch between these two slope definitions is not propagated into the quoted alpha range. The Q0 constraint is no more secure: it is normalized by assuming the distribution break is at approximately 2 Gyr, and the quoted range 10^5.5-10^7 includes only bootstrap scatter of the break position, not the uncertainty in the assumed 2 Gyr break age or in tau_break,106. Residual selection effects beyond the inverse transit probability weighting, which the paper does not model, could further bias the observed slope. These systematics are larger than the bootstrap-only uncertainties quoted for alpha and Q0. This does not invalidate the qualitative conclusion of strongly negative alpha, since consistency with the independent P18 result provides support, but the precise ranges should be treated as conditional rather than robust.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a population-level method for constraining the stellar tidal quality factor Q'* of short-period exoplanet host stars. The authors define a tidal decay timescale tau_d and use toy simulations plus a continuity-equation model to show that the low-tau_d portion of the distribution reaches a steady state with an approximately power-law slope p that depends on the frequency exponent alpha in the parameterization Q'* = Q0 (P/2 days)^alpha. They measure the observed slope in a sample of 655 transiting planets (and a hot-Jupiter subset) from the NASA Exoplanet Archive, weighting by inverse transit probability, and compare with analytic, numerical PDE, and simulated predictions. This yields alpha in [-4.33,-2] and Q0 in [10^5.5,10^7], implying strong frequency dependence, and the authors argue that the constraints are most consistent with resonance locking of g-modes. They also find that only the cool-star subsample shows the expected steady-state broken power law, and they give predictions for |P/Pdot| and discuss WASP-12 b and Kepler-1658 b.","tokens_in":29217,"tokens_out":12682,"duration_ms":124184,"significance":"If the central constraints hold, this is a valuable new empirical probe: it uses the population-level shape rather than individual decay detections, and the analytic steady-state framework is elegant and clearly presented. The numerical PDE solutions and Monte Carlo simulations are internally consistent and reproduce the analytic derivation, and the independent agreement with the Penev et al. (2018) spin-up constraints is encouraging. The authors are transparent about many caveats. However, the headline ranges for alpha and Q0 are less robust than the abstract suggests because of a slope-definition mismatch between the observed and theoretical comparisons, a conditional age calibration for Q0, and selection effects that are not quantified. The qualitative conclusion of strongly negative alpha and Q0 around 10^6-10^7 is plausible, but the precise quoted ranges need revision before they should be taken as standalone empirical constraints.","major_comments":[{"comment":"The comparison between observed and theoretical slopes uses inconsistent slope definitions. Observed slopes are power-law fits over 0.001 < tau_d < 0.3 Gyr (Section 5.2, Figure 5) or 0.001-0.4 Gyr (Figure 7), while the simulated slopes in Figure 6 are fits over tau_d < 1 Gyr per the caption, and the analytic and PDE curves are evaluated at a single point tau_d = 0.03 Gyr with tau_d,p = 1 Gyr. Equation 17 is not constant over the fitted interval: for alpha = -3 it drops from about 0.99 at tau_d = 0.001 Gyr to about 0.71 at tau_d = 0.3 Gyr. The effective slope of a power-law fit over the interval therefore depends on the interval choice, and this dependence is not propagated into the quoted alpha range. The fitting protocol should be identical for observed and simulated samples, and the uncertainty in tau_d,p (or the break location) should be marginalized over when deriving alpha.","section":"§5.3, Figure 6"},{"comment":"The Q0 constraint is calibrated by requiring the break in the observed tau_d distribution to be at 2 Gyr, and the text explicitly states that other choices are reasonable. Since Q0 scales linearly with the assumed break age, the quoted range 10^5.5 - 10^7 is conditional on a 2 Gyr population age. The bootstrap percentiles in Figure 11 do not include uncertainty in this age calibration. The abstract should either quote Q0 with the age dependence made explicit, for example Q0 approximately (tau_pop/2 Gyr) times [10^5.5, 10^7], or the analysis should marginalize over a realistic stellar age distribution for the sample.","section":"§5.6, Figure 11"},{"comment":"The only selection correction applied is inverse transit probability weighting. The sample is further restricted to transiting planets with measured masses and e < 0.02, and the detection and measurement completeness as a function of orbital period is not modeled. Because the observed slope p is the key inferential quantity, any selection function that correlates with tau_d (for example, RV mass-measurement feasibility, transit signal-to-noise, or the relationship between eccentricity filtering and circularization age) can bias the inferred alpha and Q0. The authors should quantify this effect, for instance by injecting a realistic selection function into the synthetic populations, or state explicitly that the quoted ranges assume such effects are negligible.","section":"§5.1"},{"comment":"The reported alpha range is internally inconsistent: the abstract and Section 5.3 state alpha in [-4.33,-2], while Section 5.5 gives alpha in [-4.33,-1.5] for the cool-star sample and Sections 5.6 and 5.7 use [-4.33,-1.5] for the Q0 constraints and for Figure 12. The paper should adopt a single agreed range or clearly explain why the cool-star sample is preferred in the summary statements. In addition, the lower bound -4.33 is the formal gamma > 0 boundary from Eq. 24 rather than an independently data-driven detection, so the wording should distinguish the parameterization's allowed range from the observationally favored range.","section":"Abstract vs. §5.5–5.7"}],"minor_comments":[{"comment":"The text reads 'Planets are are born'; this should be corrected to 'Planets are born'.","section":"§2, after Eq. 3"},{"comment":"The caption states that the analytic, PDE, and simulated curves agree for all alpha, but the curves are not identical; please specify the fitting intervals and the tolerance used for 'agreement'.","section":"Figure 6 caption"},{"comment":"It would be helpful to state explicitly that alpha = -13/3 is excluded because tau_d diverges there, so the formal range is alpha > -4.33 rather than alpha >= -4.33.","section":"§4.1, Eq. 24"},{"comment":"The code link and data statement say the code reproduces 'some' of the results; please provide a precise list of which figures and analyses are reproducible, and confirm the GitHub URL is correct (the typeset version contains a space).","section":"§7 and GitHub link"},{"comment":"The gray constraint lines in Figure 12 are described as spanning the cool-star alpha range [-4.33,-1.5], while the abstract cites [-4.33,-2]; the caption and text should use the same range so readers are not misled.","section":"§5.7, Figure 12"}],"recommendation":"major_revision","confidential_remarks":"The paper is well-written and the method is genuinely novel; the core qualitative conclusion is likely correct. The main issues are systematic uncertainties in how the observed slope is mapped to alpha, the conditional age calibration underlying Q0, and the internal inconsistency in the reported alpha range. These are fixable with additional analysis and clearer presentation, so I do not recommend rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's what you should know about arXiv:2501.08992: it introduces a genuinely new population-level way to measure stellar tidal dissipation, and the result is a strongly frequency-dependent Q' that matches resonance-locking predictions. The method is simple in the best way: map each observed hot Jupiter to its tidal decay timescale tau_d, look at the steady-state low-tau_d tail, and read off the exponent alpha from its power-law slope. The analytic continuity-equation model in Section 3 is clean and matches both simulations and the observed distributions. The constraints (10^5.5 <= Q0 <= 10^7, -4.33 <= alpha <= -2) are new, and the fact that they overlap with Penev et al. (2018)'s independent spin-up constraints is a real point in their favor.\n\nThe paper is also honest about its limits: Section 6.3 lists the big assumptions, Appendix B tests source-distribution sensitivity, and the resonance-locking interpretation is presented with the caveat that the mechanism is disputed. That counts.\n\nThe soft spots are real but not fatal. The observed slope is a power-law fit over a finite interval, while the analytic curve is evaluated at a single midpoint (tau_d = 0.03 Gyr). The stress-test note is right that these are different quantities, but the paper's own simulations and PDE solution use the same fitting procedure as the observations, and the three agree, so the comparison is not as loose as it first appears. More concerning: the selection correction is only inverse transit probability. The sample comes from the Exoplanet Archive, which mixes Kepler, TESS, and ground-based surveys with different completeness functions; those could bias the slope. The Q0 constraint is normalized to an assumed break age of 2 Gyr, and the quoted range includes only bootstrap scatter. If the true break age is 1 or 4 Gyr, Q0 shifts by a factor of 2. That should be stated as a systematic, not folded into the bootstrap.\n\nMy honest read: the qualitative conclusion — steep, negative alpha, efficient dissipation at long periods, weak at short periods — is likely right. The precise ranges are conditional. This deserves a serious referee. It's a novel method with a clear result and an externally checkable prediction (rare detectable orbital decay). I'd send it out, but I'd ask the authors to either propagate the selection and break-age systematics or explicitly reframe the ranges as conditional.\n\nWho should read it: anyone working on hot Jupiter evolution, tidal theory, or survey completeness. Worth citing once the ranges are stated with their caveats.","headline":"Clever new method for constraining stellar tidal dissipation from the decay-time distribution; the qualitative result is likely right, but the quoted ranges are conditional on selection and break-age assumptions.","tokens_in":29724,"tokens_out":4733,"would_cite":true,"duration_ms":43538,"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":"The distribution of tidal decay times in short-period exoplanet populations constrains the stellar tidal quality factor and favors g-mode resonance locking.","keywords":["tidal dissipation","stellar tidal quality factor","orbital decay timescale","hot Jupiters","resonance locking","continuity equation","exoplanet population statistics","transit timing"],"falsifier":"A concrete test: measure the initial semi-major axis distribution of young hot Jupiters from an unbiased survey of newly formed systems and recompute the steady-state slopes of the paper's four source-distribution models; if the true distribution moves the predicted slope at fixed $\\alpha$ by more than the bootstrap uncertainties, or if transit-timing surveys find many planets near $P \\simeq 1$ day decaying as fast as WASP-12 b ($|P/\\dot P| \\sim 3$ Myr), the inferred $\\alpha$ range and the resonance-locking preference would be overturned.","tokens_in":28576,"feed_emoji":"🪐","tokens_out":9645,"duration_ms":91084,"temperature":0.7,"pith_summary":"This paper claims that the population of short-period exoplanets, viewed through the time each planet has left before tidal inspiral, contains a steady-state power-law segment whose shape encodes the efficiency of tidal dissipation inside the host star. Using a continuity equation for the decay-time distribution, the authors show that the slope of this segment is set by the semi-major-axis dependence of the decay rate, and they measure that slope in a sample of 655 transiting planets, including 252 hot Jupiters. The data favor a stellar tidal quality factor $Q'_\\star = Q_0 (P/2\\ \\mathrm{days})^\\alpha$ with $10^{5.5} \\lesssim Q_0 \\lesssim 10^7$ and $-4.33 \\lesssim \\alpha \\lesssim -2$, meaning dissipation must decrease sharply as orbital period shortens. This rules out constant-$Q'_\\star$ equilibrium tides and, among current mechanisms, best matches resonance locking with gravity modes in the stellar interior.","feed_headline":"Stellar tides weaken steeply with planet period, census shows","feed_subtitle":"A steady-state slice of the decay-time distribution points to g-mode resonance locking in hot-Jupiter host stars.","key_machinery":"The load-bearing object is the tidal decay timescale $\\tau_d$, the time remaining before a planet spirals into its star; under stellar tides it scales as $\\tau_d \\propto a^{1/\\gamma}$ with $\\gamma = 2/(13+3\\alpha)$ for the adopted parameterization $Q'_\\star = Q_0 (P/2\\ \\mathrm{days})^\\alpha$. Planets obey a continuity equation $\\partial\\rho/\\partial t + \\partial(\\rho u)/\\partial x = S(x)$ in $x = \\ln\\tau_d$, whose steady-state solution has a log-log slope $p = [1 - (\\tau_d/\\tau_{d,p})^\\gamma]/[1 - (\\gamma+1)^{-1}(\\tau_d/\\tau_{d,p})^\\gamma]$ that is nearly constant below the break at $\\tau_{d,p} \\sim$ system age. Comparing this predicted slope with the power-law slope fitted to the observed decay-time histograms is what constrains $\\alpha$, while the break location constrains $Q_0$. The low-$\\tau_d$ region is steady because planets with $\\tau_d$ shorter than the population age have time to decay and replace one another, making the shape insensitive to most details of the birth distribution.","core_discovery":"On the paper's own terms, the central discovery is that the tidal decay-time distribution of observed short-period planets has a steady-state power-law region below $\\tau_d \\sim 1$ Gyr, exactly as predicted by a continuity-equation model of planets born and then decaying inward. The slope of that region in the observed sample is too shallow ($p \\simeq 0.6$) for a constant stellar tidal quality factor, which would give $p \\simeq 0.86$; matching the observed slope requires a steep period dependence in $Q'_\\star$. The paper consequently infers $Q'_\\star = Q_0 (P/2\\ \\mathrm{days})^\\alpha$ with $10^{5.5} \\lesssim Q_0 \\lesssim 10^7$ and $-4.33 \\lesssim \\alpha \\lesssim -2$, and finds that the inferred $Q'_\\star$ versus tidal period matches the resonance-locking prediction of g-modes more closely than equilibrium tides, inertial waves, or free-running internal gravity waves. The same analysis shows clean tidal sculpting only for planets around cool stars, with hot-star systems showing no steady-state segment, consistent with the g-mode resonance-locking picture.","pith_inferences":["If the steady-state slope is as robust as claimed, the same method could be applied to subsamples split by planet mass, stellar spectral type, or evolutionary stage to map where one tidal mechanism gives way to another.","The inferred $\\alpha$ range implies $\\tau_d$ depends on semi-major axis almost as weakly as $\\tau_d \\propto a^{1/2}$ at the steep end; future measurements of the birth semi-major axis distribution from young hot Jupiters could directly test whether the steady-state slope shifts as predicted.","Resonance locking predicts that individual systems' decay rates should track the stellar mode-evolution timescale rather than a smooth function of period; long-baseline transit-timing campaigns may detect scatter or clustering in $|P/\\dot P|$ around this prediction.","The paper's relation between the steady-state slope and disruption rate could be turned around: a future sample of planetary engulfment transients would provide an independent check on the inferred $Q'_\\star$ and on the supply of planets into decaying orbits."],"forward_implications":["Stellar tidal dissipation in these systems is strongly frequency dependent, so an inspiraling planet spends more time at short orbital periods than a constant-$Q'_\\star$ model would predict.","The observed slope rules out constant-$Q'_\\star$ equilibrium-tide models and makes resonance locking with g-modes the favored mechanism among those compared.","Planets around cool stars below the Kraft break at $T_{\\rm eff} = 6250$ K show tidal sculpting, while hot-star systems do not, implying the dominant dissipation mechanism depends on stellar interior structure.","Detectable orbital decay should be rare: the inferred $Q'_\\star$ near $P \\simeq 1$ day is large, so WASP-12 b is an outlier while Kepler-1658 b matches the population-level prediction.","The disruption rate of short-period planets can be estimated directly from the power-law slope via $N_{\\rm dis}/N_{<\\tau_{d,\\max}} = (\\Delta t/\\tau_{d,\\max})^p$, giving a population-level link to engulfment events."],"supporting_citations":[{"why":"Supplies the standard equilibrium-tide formula for orbital decay and defines the reduced tidal quality factor $Q'_\\star$ used throughout.","marker":"Goldreich & Soter 1966"},{"why":"Provides the $\\dot a/a$ stellar-tide decay expression and the earlier population-level evidence of tidal sculpting this work builds on.","marker":"Jackson et al. 2009"},{"why":"Gives the empirical $Q'_\\star \\propto P^{-3.1}$ spin-up calibration whose parameterization and favored $\\alpha$ range are directly compared.","marker":"Penev et al. 2018"},{"why":"Predicts the $Q'_\\star$ versus tidal period relation from resonance locking that best matches the paper's constraints.","marker":"Ma & Fuller 2021"},{"why":"Provides the weakly nonlinear g-mode dissipation prediction that is compared against and found inconsistent with the data.","marker":"Essick & Weinberg 2016"},{"why":"Supplies the review of dissipation mechanisms and the critical mass for internal gravity wave breaking used to interpret the constraints.","marker":"Barker 2020"},{"why":"Measured Kepler-1658 b's orbital decay and inferred its $Q'_\\star$, matching the population-level prediction.","marker":"Vissapragada et al. 2022"},{"why":"Argues that nonlinear damping does not suppress resonance locking at short periods, supporting the favored mechanism.","marker":"Zanazzi et al. 2024"},{"why":"Is the source of the observed planetary systems used to build the decay-time histograms.","marker":"NASA Exoplanet Archive 2024"}],"fun_headline_variants":["Tidal decay clock reveals g-mode resonance locking in hot Jupiters","Hot Jupiter infall times show stellar tides weaken steeply with period","Decay-time power law points to g-mode resonance in host stars","Steep stellar Q' drop inferred from planet decay-time distribution","Empirical decay-time census favors resonance locking over other tides"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole inference rests on the assumption that the observed low-$\\tau_d$ planets form a steady-state population whose fitted slope is close to the analytic steady-state solution, so that the unknown initial semi-major axis distribution and residual selection effects do not shift the slope by more than the reported uncertainties.","fun_headline_variants_meta":{"raw":{"variants":["Tidal decay clock reveals g-mode resonance locking in hot Jupiters","Hot Jupiter infall times show stellar tides weaken steeply with period","Decay-time power law points to g-mode resonance in host stars","Steep stellar Q' drop inferred from planet decay-time distribution","Empirical decay-time census favors resonance locking over other tides"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000608,"raw_usage":{"total_tokens":2879,"prompt_tokens":1036,"completion_tokens":1843,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":652,"completion_tokens_details":{"reasoning_tokens":1754}},"tokens_in":652,"tokens_out":1843,"duration_ms":17217,"temperature":1.0,"reasoning_tokens":1754,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:13:38.976772+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test: measure the initial semi-major axis distribution of young hot Jupiters from an unbiased survey of newly formed systems and recompute the steady-state slopes of the paper's four source-distribution models; if the true distribution moves the predicted slope at fixed $\\alpha$ by more than the bootstrap uncertainties, or if transit-timing surveys find many planets near $P \\simeq 1$ day decaying as fast as WASP-12 b ($|P/\\dot P| \\sim 3$ Myr), the inferred $\\alpha$ range and the resonance-locking preference would be overturned.","supporting_citations":[],"review_version":1}