{"id":"680edd6f-01d4-4fa7-8e9f-af3255de56d4","arxiv_id":"2502.01736","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A two-planet merger into Kepler-221 d can explain the non-adjacent (b,c,e) three-body resonance, but only if convergent migration reformed the chain and the expansion avoided a specific (c,d,e) resonance.","lead":"This paper proposes that the Kepler-221 system originally contained five planets and that two of them merged to form the present-day middle planet, which would explain the unusual three-body resonance among the other three. It uses N-body simulations to map the conditions under which this rare architecture can form, and derives constraints on planet masses and tidal parameters that can be tested with new observations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The model's reformation and age requirements demand contradictory tidal strengths: Phase III fails for Qphy < 3, while the <650 Myr age requires Qphy < 2, so the full history works only if an unmodeled phase-dependent dissipation mechanism is invoked.","rationale":"The reader's weakest assumption was the unmodeled torque or enhanced damping needed for resonance reformation; I agree that is a genuine gap, but I see the age-Qphy conflict as at least as decisive because it shows the two required dissipative conditions are nearly mutually exclusive: reformation needs weak damping, while expansion within <650 Myr needs strong damping. The paper's own successful all-phase simulation has an effective age of about 2 Gyr, and the obliquity-tide paragraph is a suggestion rather than a modeled mechanism. I nonetheless keep the reader's CONDITIONAL verdict rather than moving to REJECT: the paper explicitly frames the scenario as a feasibility demonstration, states the tension in conclusion item 6, and offers a clean mass-ratio prediction (acde > 1.56, Equation 11) that future RV follow-up can test. The analytical derivation of Equation 10 is a useful contribution, and including planet d in Phase IV is a genuine improvement over Goldberg & Batygin (2021). The proposed test is inexpensive and settles whether the tidal-history contradiction is real: if Phase III already fails at Qphy = 2, then no single tidal-Q history can satisfy both the reformation and age constraints, and the obliquity-tide resolution must be explicitly demonstrated rather than assumed.","tokens_in":33020,"tokens_out":18595,"duration_ms":188261,"concrete_test":"Run the Phase III parameter grid of Figure 4 with Qphy = 2 (the age-required maximum from Equation 14) and the same outward torque on planet b as in Table 3, and record whether any run reforms the zeroth-order (b,c,e) 3BR.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's own Section 6.3 exposes a quantitative tension that is more load-bearing than any single missing torque: Phase III reformation requires gentle convergent migration corresponding to Qphy > 3 (Table 3 and Figure 4), while the system age (<650 Myr, Berger et al. 2018) combined with the expansion-rate relation texp ≈ 380 Qphy Myr demands Qphy < 2 (Equation 14). The fully successful four-phase simulation (Figure 2) corresponds to a physical age of about 2 Gyr, not <650 Myr. The proposed resolution via obliquity tides (Equation 15) is invoked but not implemented; for it to work, the effective damping must be weak during resonance reformation and strong during orbital expansion, yet the instability and reformation phase is exactly when non-zero obliquity should be acquired and would already speed up damping. Without a self-consistent tidal model that produces this time ordering, the model has not shown that Kepler-221 can reach its observed state within its inferred age. The authors are candid about this limitation, but it is the condition most likely to falsify the central historical claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a four-phase dynamical history for the four-planet system Kepler-221, in which five planets originally formed in a first-order resonance chain; after disk dispersal an instability breaks the chain, two planets merge to form the present planet d, and tidal dissipation plus a prescribed outward torque on planet b reform and then expand the non-adjacent (b, c, e) 6:3:1 three-body resonance to the observed period ratios. The authors support the scenario with REBOUND/REBOUNDx N-body simulations, a single tuned full-phase simulation, phase-by-phase parameter studies, and an analytical derivation of the expansion slope acde in the (c, d, e) period-ratio plane. The paper derives mass-ratio constraints from the requirement acde > 1.56 and discusses applications to K2-138 and Kepler-402.","tokens_in":33259,"tokens_out":7107,"duration_ms":71074,"significance":"If the scenario is correct, the paper provides a plausible formation channel for non-adjacent three-body resonances with an intervening non-resonant planet, a role for the intermediate planet d in destabilizing or preserving the resonance during tidal expansion, and falsifiable mass-ratio predictions for a system with weak TTVs. Strengths include the clean analytical derivation of the expansion slope in Appendix C, the explicit inclusion of planet d in the expansion dynamics (which previous work neglected), the use of public community codes, and the candid discussion of the model's limitations in Section 6.3 and the conclusions. The mass-ratio and density-ratio constraints in Figure 8 are concrete and testable with radial-velocity follow-up. However, as discussed below, the age versus tidal-damping inconsistency acknowledged by the authors is load-bearing and prevents the paper, in its current form, from establishing that the full four-phase history actually reproduces the Kepler-221 system within its inferred age.","major_comments":[{"comment":"The paper contains a quantitative internal inconsistency that is load-bearing for the central historical claim. The resonance-reformation simulations require gentle migration with Qphy > 3 (Section 4.4 and Table 3), while the <650 Myr age estimate (Berger et al. 2018) combined with texp ≈ 380 Qphy Myr (Eq. 14) requires Qphy < 2. The fully successful simulation in Figure 2 corresponds to an expansion time of about 2 Gyr, not <650 Myr. The proposed resolution via obliquity tides (Eq. 15) is qualitative and is not implemented in any simulation; moreover, the period of dynamical instability and resonance reformation is exactly when non-zero obliquity would be excited, so the required sequence of weak damping during Phase III and strong damping during Phase IV lacks a demonstrated physical basis. As it stands, the paper has not shown that the four-phase model can reach the observed configuration within the age of the system.","section":"Section 6.3, Eq. (14), Conclusion item 6"},{"comment":"The reformation of the (b, c, e) zeroth-order 3BR is conditional on a convergent migration between b and c, achieved either by an outward exponentially decaying torque on planet b (Eq. 8) or by strongly enhanced damping on planet c (Qc/Qb = 0.1). The paper parameterizes these mechanisms but does not demonstrate that either operates in Kepler-221; the cited mass-loss and dynamical-tide origins are not modeled from first principles. Since the paper's own simulations show that without such a torque the resonance reforms only as a first-order 3BR or not at all (Figure 3a,d), the physical origin of this torque is a load-bearing assumption. I recommend that the authors either implement a concrete physical mechanism or, if the mechanism is to remain an ansatz, state more prominently that the scenario requires this unmodeled condition and quantify the resulting joint probability of the full sequence rather than only individual phases.","section":"Section 4.3, Eq. (8), Section 4.4, Figure 3"},{"comment":"The mass constraints derived from acde > 1.56 are inverse matching conditions rather than independent predictions. The threshold starts from a specific initial condition, namely the intersection of Pe/Pc = 3 with the (3, 5, 8) resonance line (Appendix C, Figure C.1), and the successful Phase IV simulations use mass model M2a, which was selected after the constraint was derived because it satisfies acde > 1.56. The mass ratios of b, c, and e are therefore constrained only under the assumption that the expansion begins at that particular point and proceeds along the zero-order 3BR line. I ask the authors to clarify explicitly that these are necessary conditions within the assumed scenario, not observational predictions independent of the starting geometry, and to discuss how the constraint changes if the initial Pd/Pc is not exactly at the adopted intersection.","section":"Section 5.2.2, Figure 8, Eq. (10)-(11)"}],"minor_comments":[{"comment":"The caption for Figure 2 states that the orbital expansion phase is represented by 'panels d and h,' but panel d is already described earlier in the same caption as a 3BR-angle panel; the body text says panel d shows the (b, c, e) period ratios and panel h shows the (c, d, e) period ratios. Please correct the caption to avoid this inconsistency.","section":"Figure 2 caption and Section 3"},{"comment":"The sentence 'In panel c and g of Figure 2, the torque applied to planet b amounts to a 2% mass loss in 7.5 Myr' appears to refer to the text of Section 4.3 rather than to anything displayed in Figure 2, whose panels c and g show orbital elements and resonance angles; please clarify the reference.","section":"Section 4.3"},{"comment":"The manuscript contains numerous spacing artifacts from LaTeX or PDF extraction, including 'e fficient,' 'di fferent,' 'V ogt,' and 'Go´ zdziewski'; these should be cleaned before publication.","section":"Throughout"},{"comment":"Table 3 is dense and difficult to parse because the parameter column mixes symbols with ranges and notes; consider splitting the table into per-phase blocks or adding a separate parameter-definition table for symbols such as tb_a0,III and tb_Gamma,IV.","section":"Section 2.3, Table 3"},{"comment":"The sign convention in ΔPe/ΔPc could be made explicit: since Pc decreases while Pe increases during expansion, the raw derivative in Eq. (C.8) is negative, and the positive slope acde in Eq. (C.9) follows after multiplying by -Pc^2/Pd^2; adding one sentence to this effect would prevent reader confusion.","section":"Appendix C, Eq. (C.8)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a well-crafted scenario study with a clean analytical core and honest limitations, but the age versus Qphy tension in Section 6.3 is the crux of the central claim. I believe the authors can address it with a self-consistent implementation of obliquity tides or a comparable phase-dependent dissipation model, or by clearly reframing the result as a conditional scenario with an older age or unknown mechanism. The mass-ratio constraints and the role of planet d are valuable and should survive revision. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a real advance on Goldberg & Batygin (2021). It adds the collision origin of planet d, the post-collision reformation problem, and most usefully the role of d in the expansion phase, where the (c,d,e) (3,5,8) resonance can break the (b,c,e) chain. The analytical slope acde (Eq. C.9) is a clean derivation from angular-momentum and resonance invariants, and the paper is unusually honest about its success rates and assumptions.\n\nWhat the paper does well: the modular N-body campaign covers each phase separately, quantifies where failures happen, and shows that the zeroth-order 3BR does not reform without convergent b-c migration. The result that d must start on the correct side of the (3,5,8) resonance, and that acde > 1.56 is needed to reach the observed period ratios, is new and testable through future TTV/RV mass measurements. The byproduct discussion of first-order 3BR formation (Appendix A) is a nice bonus for K2-138.\n\nSoft spots, in order. The central claim is feasibility, not uniqueness. The full run needs a prescribed kick, an unmodeled outward torque on b (or Qc/Qb = 0.1), a chosen mass model, and a starting location for d that avoids the dangerous resonance; the paper says so, but each tuned ingredient multiplies the prior. The unmodeled torque on b is a real concern, though the paper does show order-of-magnitude consistency with mass-loss and dynamical-tide estimates; the age tension in Section 6.3 is the one that actually bites. Phase III requires Qphy > 3, the <650 Myr system age requires Qphy < 2, and the successful full simulation corresponds to about 2 Gyr. Invoking obliquity tides is plausible, but it is not implemented, and the required time ordering—weak damping during reformation, strong damping afterward—is exactly the order a Cassini-state obliquity would have to explain rather than assume. Also, the mass constraints are inverse requirements: M2a is chosen partly because it gives acde > 1.56. That is acceptable for a scenario paper, but these are constraints, not predictions. No code or data are released, which makes independent checking harder than it should be.\n\nWho benefits: resonant-chain dynamicists and people working on post-disk instability and giant impacts. The paper deserves a serious referee. I would send it to review and ask the authors either to implement a self-consistent tidal/obliquity model for the age problem or to soften the historical claim to a conditional feasibility study. The analytical slope and the warning about d's role will be citable regardless.","headline":"A solid, inventive four-phase merger scenario for Kepler-221 that advances the resonant-chain story, with a real age-vs-tidal-strength tension the authors acknowledge but do not resolve.","tokens_in":33869,"tokens_out":4119,"would_cite":true,"duration_ms":41866,"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":"Kepler-221's odd resonance chain is the fossil of a five-planet system in which two planets merged, and a four-phase simulation can reproduce the observed arrangement.","keywords":["Kepler-221","three-body resonance","planet-planet merger","resonance chain instability","tidal dissipation","orbital migration","exoplanet architecture"],"falsifier":"Measure the masses of planets b, c, d, and e, for example by radial velocities or transit-timing variations. If the mass ratios put the system in the gray region of the paper's diagram, i.e., $a_{cde} \\le 1.56$ with planet e as massive as or heavier than c, the proposed tidal expansion cannot reach the observed period ratios without fine-tuned torques; if planet d turns out to be the lightest planet, the merger origin of d would also not be supported.","tokens_in":32699,"feed_emoji":"🪐","tokens_out":8071,"duration_ms":71008,"temperature":0.7,"pith_summary":"This paper tries to establish that the Kepler-221 system did not always look like it does today. It argues that the present architecture, with planets b, c, and e in a 6:3:1 three-body resonance while the intermediate planet d stays out of it, is the end product of a four-phase history: five planets formed in a chain of first-order resonances, the chain broke after disk dispersal, two planets merged to become d, and tidal dissipation later rebuilt and stretched the (b, c, e) resonance to the observed period ratios. The argument matters because Kepler-221 is a case where simple disk migration cannot explain what is seen, and because the proposed history yields concrete, testable predictions about the planets' masses and tidal properties. A full N-body simulation is shown to pass through all four phases and arrive at the observed period ratios, with the key quantitative condition being a slope $a_{cde} > 1.56$.","feed_headline":"A merged pair of planets can explain Kepler-221's odd resonance","feed_subtitle":"The 6:3:1 b–c–e resonance survives while planet d stays out; the scenario predicts d is a merger remnant.","key_machinery":"The load-bearing objects are the zeroth-order three-body resonance, an angle built from the mean longitudes of b, c, and e, $\\phi_{3BR}=2\\lambda_b-5\\lambda_c+3\\lambda_e$, whose libration defines the 6:3:1 chain, and the expansion slope $a_{cde}$, the ratio in which the period ratios $P_e/P_d$ and $P_d/P_c$ grow during tidal expansion. The resonance angle carries the identity of the chain, while the slope carries the geometry of the (c, d, e) period-ratio plane. The paper derives $a_{cde}$ analytically from angular-momentum conservation and the invariance of the resonance angle, obtaining $a_{cde}=(4.90 m_b + 4.94 m_c)/(7.27 m_e - m_b)$, and shows by N-body integration that trajectories with $a_{cde}>1.56$ can avoid the destructive (3,5,8) cde resonance and land on the observed ratios. The same machinery identifies a favorable mass model, roughly equal masses with $m_c \\approx 1.25 m_b$ and $m_c \\approx 1.25 m_e$, and, when masses are unfavorable, requires a time-dependent outward torque on planet b to steepen the effective slope.","core_discovery":"The paper's central claim is that the Kepler-221 system is the relic of a five-planet resonance chain that was destroyed and partially rebuilt. In the proposed scenario, five planets migrate into a first-order chain of 2:1, 3:2, 4:3, and 3:2 resonances; after the gas disk disperses, an instability (delivered in the simulations as a small kick to one of the two planets d1 and d2) breaks the whole chain, and the two d planets merge into the observed planet d. The (b, c, e) 6:3:1 zeroth-order three-body resonance then reforms only if planets b and c undergo convergent migration, realized either as an outward torque on b (parametrized as an exponentially decaying torque, interpretable as mass loss or dynamical tides) or as anomalously strong tidal damping on c. Once the resonance is re-established, tidal dissipation expands the chain along the (2,3,5) 3BR line toward $P_c/P_b \\approx 2.035$ and $P_e/P_c \\approx 3.228$. Planet d is not passive in this last phase: if the trajectory crosses the (3,5,8) three-body resonance of c, d, and e, that resonance destabilizes and breaks the (b, c, e) chain in more than 95% of the simulations. A successful run therefore requires the expansion slope $a_{cde} = \\Delta(P_e/P_d)/\\Delta(P_d/P_c)$ to exceed 1.56, which in turn constrains the planet mass ratios to values close to the peas-in-a-pod model with a slightly more massive planet c.","pith_inferences":["Going beyond the paper's explicit claims, the slope condition $a_{cde}>1.56$ supplies a generally useful diagnostic: in any compact resonant chain, an interloper planet that is not part of the chain controls whether expansion can proceed, and its observed period ratios encode the relative masses of the chain members.","A clean testable extension would target the predicted mass ordering with radial velocities: if future data place the mass ratios in the gray region of the paper's diagram, e.g., with planet e as massive as or heavier than c, the proposed expansion phase fails without fine-tuned torques, and a merger origin for d becomes much less plausible.","One could also search for analogous systems in the Kepler and TESS samples, pairs of planets straddling a non-resonant third body with the outer pair near a three-body resonance, as candidate products of the same chain-breaking process that need not have started from exactly five planets."],"forward_implications":["If the four-phase scenario is correct, planet d is a merger remnant, so its mass should be comparable to or larger than the other planets rather than following a monotonic mass-radius relation.","The requirement $a_{cde}>1.56$ means planets b, c, and e must have nearly equal masses with planet c somewhat heavier; planets c and e would then be low-density super-puff bodies.","The expansion phase carries an age tension: reaching the observed period ratios within roughly 650 Myr needs effective tidal damping $Q_{\\mathrm{phy}}<2$, while the reformation phase needs $Q_{\\mathrm{phy}}>3$, implying either an older system or an extra dissipation channel such as obliquity tides.","If planet d is initially placed on the wrong side of the (3,5,8) cde resonance, the chain is destroyed in over 95% of simulations, so the present architecture records a narrow avoidance of that resonance.","The same merger-and-reformation logic is offered for Kepler-402 and for the first-order (2,3,4) three-body resonance at the outer edge of K2-138."],"supporting_citations":[{"why":"Identifies the non-adjacent (b, c, e) 6:3:1 three-body resonance in Kepler-221 and provides the baseline case for tidal expansion of the chain.","marker":"Goldberg & Batygin 2021"},{"why":"Supplies the theory that zeroth-order three-body resonances require convergent migration to reform, the central requirement of the post-collision phase.","marker":"Petit 2021"},{"why":"Provides the tidal eccentricity damping timescale used in the simulations and the framework for orbital expansion of resonant chains.","marker":"Papaloizou et al. 2018"},{"why":"Gives the mass-radius relation used for mass model M1, the baseline against which successful mass models are compared.","marker":"Chen & Kipping 2017"},{"why":"Defines the peas-in-a-pod equal-mass framework used for mass model M2 and the successful M2a variant.","marker":"Weiss & Petigura 2020"},{"why":"Sets the system age, stellar parameters, and planet radii that impose the 650 Myr timescale and the tidal-quality constraints.","marker":"Berger et al. 2018"},{"why":"Provides the post-disk instability scenario invoked to break the first-order resonance chain after disk dispersal.","marker":"Raymond et al. 2022"}],"fun_headline_variants":["Kepler-221's odd resonance hints at a lost planet merger","How five planets became four in Kepler-221","Merged planets behind Kepler-221's 6:3:1 resonance","Kepler-221's resonance from a merged planetary pair"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that, after the collision that created planet d, planets b and c moved toward each other convergently for a while, which the model supplies either as an outward force on planet b or as much stronger tidal braking on planet c; if no such mechanism operated in the real system, the b/c/e resonance never reforms and the whole history fails.","fun_headline_variants_meta":{"raw":{"variants":["Kepler-221's odd resonance hints at a lost planet merger","How five planets became four in Kepler-221","Merged planets behind Kepler-221's 6:3:1 resonance","Kepler-221's resonance from a merged planetary pair"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001174,"raw_usage":{"total_tokens":4971,"prompt_tokens":1179,"completion_tokens":3792,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":795,"completion_tokens_details":{"reasoning_tokens":3720}},"tokens_in":795,"tokens_out":3792,"duration_ms":22408,"temperature":1.0,"reasoning_tokens":3720,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T14:39:28.418306+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the masses of planets b, c, d, and e, for example by radial velocities or transit-timing variations. If the mass ratios put the system in the gray region of the paper's diagram, i.e., $a_{cde} \\le 1.56$ with planet e as massive as or heavier than c, the proposed tidal expansion cannot reach the observed period ratios without fine-tuned torques; if planet d turns out to be the lightest planet, the merger origin of d would also not be supported.","supporting_citations":[{"cited_title":"& Batygin , K","cited_arxiv_id":null,"evidence_quote":"Identifies the non-adjacent (b, c, e) 6:3:1 three-body resonance in Kepler-221 and provides the baseline case for tidal expansion of the chain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the theory that zeroth-order three-body resonances require convergent migration to reform, the central requirement of the post-collision phase."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the peas-in-a-pod equal-mass framework used for mass model M2 and the successful M2a variant."},{"cited_title":"N., Izidoro , A., Bolmont , E., et al","cited_arxiv_id":null,"evidence_quote":"Provides the post-disk instability scenario invoked to break the first-order resonance chain after disk dispersal."}],"review_version":1}