{"id":"c3b3be86-dd28-4c96-bed5-a456b02e2ab0","arxiv_id":"2607.08544","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"A matrix-based extension of Papaloizou (2015) gives the tidal separation timescale T for N-planet chains and converts observed offsets into effective Q' bounds, with special sensitivity to the second and outermost planet masses.","lead":"The authors extend Papaloizou's analytical model of tidal separation from three-planet to N-planet resonant chains of adjacent first-order mean-motion resonances, yielding an explicit formula for the time evolution of period-ratio offsets. This lets observers place bounds on the hard-to-measure planetary tidal quality factor Q' even when stellar age or outer masses are poorly known.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The paper's strongest claim is an explicit, invertible timescale for the innermost period-ratio evolution under planetary tides for arbitrary N, provided only adjacent first-order 2P-MMRs are present. That derivation (Appendix B, matrix A/B construction, solution for T, Eq. 12) is internally consistent, recovers the three-planet Papaloizou (2015) limit, and matches the N-body experiments of Fig. 2 across N=3–6. The applications in Sections 4.1–4.3 correctly treat incomplete data and outer-planet omissions as producing upper bounds rather than point estimates, and the sensitivity to m2 and mN is a useful byproduct. The domain restriction flagged by the reader is real and already acknowledged by the authors; it does not undermine the mathematics or the validation inside the stated domain. No additional load-bearing inconsistency appears. Therefore the CONDITIONAL verdict with high confidence remains appropriate; no adjustment is required.","tokens_in":28117,"tokens_out":537,"duration_ms":5592,"concrete_test":"Re-run the six-planet N-body validation of Fig. 2 after inserting a non-adjacent first-order (or second-order) resonance between planets 4–6 while keeping all other parameters fixed; if the analytic n1/n2(t) curve from Eq. 12 still tracks the numerical separation to within a few percent, the collapse-to-Δ12 approximation is more robust than claimed; if it diverges, the upper-bound procedure of Section 4.2 is confirmed as necessary.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (explicit T from the matrix solution of the secular equations, inverted for effective Q' bounds even with incomplete data) holds under the paper's stated assumptions. The derivation in Appendix B correctly collapses all Δij to a single Δ12 only when adjacent first-order 2P-MMRs are present; the authors openly convert applications to real systems (K2-138) into upper bounds by dropping outer planets, and Section 3.2 plus Fig. 3 quantify when initial disk offsets can be neglected. Validation against N-body runs (Fig. 2) and recovery of prior numerical Q' estimates further support the result inside its domain. The reader's weakest_assumption correctly flags the domain restriction, but that restriction is already treated as converting the tool into an upper-bound estimator rather than a hidden flaw in the math.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper extends Papaloizou (2015) to N-planet resonant chains of adjacent first-order 2P-MMRs by casting the secular equations into a matrix form (Appendix B). Solving for the resonant variables yields an explicit timescale T such that the innermost period ratio evolves as n1/n2(t)=(q12+1)/q12 * (1+(1/100)(t/T)^{1/3}) (Eq. 12). Inverting this relation (Eq. 16) produces an effective planetary Q' that can be bounded even when stellar age or an outer mass is unknown. Analytic curves match multi-Gyr N-body runs for N=3–6 (Fig. 2); application to the five inner planets of K2-138 recovers upper bounds consistent with earlier numerical work (Cerioni & Beaugé 2023). The authors also show that a more massive second planet boosts separation while a more massive outermost planet inhibits it.","tokens_in":28290,"tokens_out":894,"duration_ms":8803,"significance":"If the derivation holds inside its stated domain, the work supplies a practical, inexpensive tool for placing quantitative bounds on planetary Q'—a parameter that is otherwise poorly constrained for rocky exoplanets—from the observed architecture of resonant chains. The matrix construction is transparent and immediately generalizable; the explicit T and the inversion formula (Eqs. 12, 16) are falsifiable against N-body integrations and against future systems that satisfy the adjacency/first-order assumptions. Recovery of the prior K2-138 numerical estimate and the clear conversion of incomplete systems into upper bounds are concrete strengths. The paper therefore advances both the analytic toolkit and the observational interpretation of multi-planet resonant systems.","major_comments":[{"comment":"Section 2.3 and Appendix A argue that planetary tides dominate once e1 exceeds a critical value ~10^{-3}, yet the analytic model of Appendix B and Eq. (12) permanently omit stellar tides. For systems that spend many Gyr near the observed offsets (or that have small m2/m0 or large R0/R1), the neglected stellar contribution can become comparable (Fig. A.1, System A). A short quantitative estimate of the fractional error in T (or in the inverted Q') when stellar tides are restored would strengthen the claim that the bounds remain meaningful for the systems of interest.","section":null},{"comment":"Section 3.3 and Eq. (16) treat T as a linear combination of the individual 1/Q'i and then replace all Q'i by a single effective Q'. Matrix M itself depends weakly on the circularization timescales, so the linearity is only approximate. The Monte-Carlo experiments of Sections 4.1–4.2 recover the expected weighted average, but an explicit statement of the residual error introduced by freezing M (or a brief numerical check that the approximation remains <10–20 % across the explored mass/radius range) would make the inversion more robust.","section":null}],"minor_comments":[{"comment":"Equation (12) and the surrounding text use both “1 % separation” and the relative offset f12(t); a single consistent definition early in Section 3 would reduce ambiguity.","section":null},{"comment":"Figure 2 caption should state the precise initial conditions (a1, Q'i vector, stellar mass) so that the comparison can be reproduced without hunting through the text.","section":null},{"comment":"Table 1 mixes fixed and uncertain parameters; a footnote clarifying which quantities were held fixed in the Monte-Carlo draws would help.","section":null},{"comment":"A few typographical slips remain (e.g., “useing” in the Introduction, “off 12(0)” in Section 3.2). A careful proof-read is warranted.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The domain restriction (adjacent first-order MMRs only) is already handled honestly by converting real-system applications into upper bounds; I do not regard it as a hidden flaw. The paper is a solid, incremental advance that fits the journal’s scope."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a practical extension of Papaloizou (2015) that actually delivers an explicit timescale T and a closed-form n1/n2(t) for arbitrary N under adjacent first-order 2P-MMRs. The matrix setup in Appendix B is transparent, the analytic curves track multi-Gyr N-body runs for N=3–6 (Fig. 2), and the K2-138 exercise recovers the earlier numerical Q'~50 result as an upper bound once the outer planet is dropped. That is real utility: you can invert observed offsets for effective Q' even when age or one mass is missing, and the mass-sensitivity results (second planet boosts, last planet inhibits) are clean and new.\n\nWhat the paper does well is stay inside its assumptions and convert the limitations into upper bounds rather than overclaiming. Section 3.2 and Fig. 3 quantify when disk-driven initial offsets can be neglected; the stellar-tide argument is given with a critical-e threshold and an appendix check. Self-citations to their own K2-138 numerics are used only as external benchmarks, not as circular inputs. The free parameters (effective Q', age, outer mass) are exactly the ones the inversion is built to handle.\n\nSoft spots are real but proportionate. The method still requires adjacent first-order resonances only, so most of the longest known chains force you to drop outer planets and report Q'_max. That is already stated clearly; it is a domain restriction, not a hidden flaw in the math. Stellar tides are neglected after the critical-e argument; for the systems they care about the approximation holds, but it is an approximation. No load-bearing circularity or invented entities.\n\nThis is for people who work on resonant-chain dynamics or who need quick Q' bounds from architecture. The math is solid enough that a serious editor should send it to referees. I would cite the T formula and the mass-sensitivity results if I were estimating dissipation in a chain that fits the assumptions. Bring it to reading group if the group cares about analytic tools for tides; otherwise it is a solid methods paper that does what it claims.","headline":"Clean matrix extension of Papaloizou 2015 that turns resonant offsets into usable Q' bounds, with honest domain limits and solid N-body checks.","tokens_in":28923,"tokens_out":529,"would_cite":true,"duration_ms":5992,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"An analytical timescale for tidal opening of long resonant chains lets observers bound planetary Q' even when age or mass is missing.","keywords":["resonant chains","tidal dissipation","modified quality factor Q'","mean-motion resonances","period-ratio offsets","N-planet systems"],"falsifier":"Apply the formula to a well-dated multi-planet chain whose Q' has already been measured independently by tidal circularization or spin evolution; a systematic mismatch larger than the quoted uncertainties would falsify the analytic T.","tokens_in":28961,"feed_emoji":"🌍","tokens_out":556,"duration_ms":5981,"temperature":0.7,"pith_summary":"Resonant chains of three or more planets sit slightly off exact period commensurabilities. That offset is a fossil of long-term tidal energy loss. Direct N-body tidal runs are too slow for broad parameter surveys, so the authors extend a three-planet analytic model to arbitrary length. They obtain a single timescale T that governs how the innermost period ratio grows as the cube root of time. Inverting the formula yields an effective planetary modified quality factor Q' from the observed offset, the masses and radii, and the stellar age. The same formula still returns useful upper or lower bounds when age or an outer-planet mass is unknown, and it shows that the second planet's mass accelerates the opening while the outermost planet's mass resists it. The method therefore turns already-measured orbital architectures into quantitative constraints on poorly known interior dissipation.","feed_headline":"Tidal offsets of resonant chains yield planetary Q' bounds","feed_subtitle":"A cube-root law converts observed period ratios into dissipation limits even without a precise age.","key_machinery":"The matrix M = A^{-1}B that converts the vector of inverse circularization times into the resonant forcing terms Xij; those terms determine the single constant T that appears in the cube-root law for period-ratio growth.","core_discovery":"The secular equations of an N-planet chain of adjacent first-order mean-motion resonances can be cast as a linear matrix system whose solution supplies an explicit tidal-separation timescale T. Once T is known, the innermost period ratio evolves as n1/n2(t)=(q12+1)/q12 * (1+(1/100)(t/T)^{1/3}). Inverting that relation recovers an effective planetary Q' (or a rigorous upper bound) even when the stellar age or one planetary mass is poorly constrained.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Cube-root law converts resonant offsets into planetary Q' bounds","Secular matrix yields tidal timescale for long resonant chains","Resonant chain separations constrain Q' without stellar age","Second-planet mass boosts tidal separation in resonant chains","Analytical tool recovers Q' bounds from poorly known chains"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The derivation requires that only adjacent first-order resonances are present and that any residual offsets left by the disk are negligible compared with the offsets observed today.","fun_headline_variants_meta":{"raw":{"variants":["Cube-root law converts resonant offsets into planetary Q' bounds","Secular matrix yields tidal timescale for long resonant chains","Resonant chain separations constrain Q' without stellar age","Second-planet mass boosts tidal separation in resonant chains","Analytical tool recovers Q' bounds from poorly known chains"]},"model":"grok-4.5","effort":"low","cost_usd":0.007024,"raw_usage":{"total_tokens":1812,"prompt_tokens":866,"num_sources_used":0,"completion_tokens":81,"cost_in_usd_ticks":70240000,"prompt_tokens_details":{"text_tokens":866,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":865,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":866,"tokens_out":81,"duration_ms":8475,"temperature":1.0,"reasoning_tokens":865,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-10T05:38:05.212333+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Apply the formula to a well-dated multi-planet chain whose Q' has already been measured independently by tidal circularization or spin evolution; a systematic mismatch larger than the quoted uncertainties would falsify the analytic T.","supporting_citations":[],"review_version":1}