{"id":"12536fda-d66d-49cc-9992-8cb472d09118","arxiv_id":"2507.17092","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In a binary star system, the companion star's gravity moves the giant planets' secular resonances farther out and suppresses them, leaving a larger region where terrestrial planets can form.","lead":"Two researchers derive equations for where secular resonances, the slow gravitational alignments that sculpt planetary systems, sit in binary stars that host two giant planets. Their simulations suggest the companion star pushes these resonances outward and weakens them, potentially widening the zone where rocky planets could form.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (27) omits disk self-gravity even though §4.2.4 attributes resonance migration to the disk's 'regressing effect'; the claimed theory-simulation agreement is therefore not a clean test.","rationale":"I read the paper as attempting to derive the secular resonance locations for a two-planet binary system using standard Laplace-Lagrange theory and to establish, with N-body simulations, that a stronger secondary moves these locations outward and suppresses the resonances. The analytic machinery in Section 2 is conventional, and the trend that g_k increases with secondary mass and decreasing binary separation is visible in Tables 1–4. The weakest point is the test-particle precession rate: Eq. (27) contains only the planets and the secondary star, while the simulations include a massive disk, and the paper itself states in §4.2.4 that the disk has a 'regressing effect' responsible for inward resonance migration. This makes the Section 4.1 agreement a comparison between a disk-free prediction and a disk-bearing simulation; without a quantitative estimate of the disk's contribution to Ω_E, the agreement could be coincidental. The suppression claim is similarly vulnerable because Figure 8 compares full-binary and no-secondary runs that used different integrators and computers, and the disk is present in both panels, so the damping could be a disk effect. These weaknesses are fixable with targeted reruns, so rejection is not warranted; the paper needs the additional control to support its general claims.","tokens_in":13253,"tokens_out":8871,"duration_ms":103015,"concrete_test":"Rerun the Section 4.1 validation simulations with the same initial conditions and integrator, but with the mutual gravity of disk particles set to zero, treating embryos and planetesimals as test particles while keeping the planets and secondary active. Overlay the resulting eccentricity excitation on the g1 and g2 locations from Eq. (27). If resonances appear at the predicted locations with disk self-gravity disabled, the omission is not load-bearing; if they shift substantially or disappear, the claimed theory-simulation agreement is an artifact of disk self-gravity rather than a confirmation of the binary secular theory.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The test-particle precession rate in Eq. (27) sums only the two giant planets and the secondary star. The validating simulations, however, include a massive disk of about 100 Moon-to-Mars embryos plus over 400 planetesimals at 0.5–1.5 au, and §4.2.4 explicitly states that the disk exerts a 'regressing effect' whose weakening as the disk loses mass is what moves the secular resonances inward. That is an admission that disk self-gravity contributes to Ω_E at the epochs where the theory is compared with simulation. If the disk contribution is comparable to the g1 and g2 frequencies (roughly 0.1–0.3 deg/yr in Tables 1–4), the resonance locations predicted by Eq. (27) are shifted, so the visual agreement in Figures 5–6 could be coincidental or dominated by disk-driven precession rather than by the binary secular theory. The same omission contaminates the suppression claim in §4.2.2: the disk is present in both panels of Figure 8, and disk self-gravity, not the secondary, might be what damps the eccentricity response. The assumption is unflagged in the analytic derivation, and no run without disk self-gravity is presented.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops a secular perturbation theory for a coplanar binary star system in which the primary hosts two giant planets. The authors derive the eigenvalues g_k of the secular matrix for the four-body system (inner planet, outer planet, secondary star) and the test-particle precession rate Ω_E of a small body (Eq. 27), defining secular resonances by the condition Ω_E = g_k. They apply the theory to binaries with secondary masses from 0.4 to 1.3 solar masses and binary semimajor axes from 20 to 50 au, and validate it with N-body integrations of a planetesimal/embryo disk interior to the inner planet. The paper claims that stronger secondary perturbation moves the secular resonances outward (closer to the giant planets), that the secondary star suppresses the secular resonances of the giant planets, and that disk mass loss causes the resonances to migrate inward.","tokens_in":13501,"tokens_out":5334,"duration_ms":53884,"significance":"If the theory is correct, it provides a parameter-free analytic tool for predicting secular resonance locations in planet-hosting binaries, using only masses and semimajor axes as input. The derivation follows the standard Laplace-Lagrange secular framework and extends it to include the secondary star as a third massive perturber; the eigenvalues are computed in closed form from a 3x3 matrix. The claimed suppression of secular resonances by the secondary, and the inward migration of resonances as the disk loses mass, are potentially important for terrestrial planet formation in binaries. However, the validation is purely qualitative and the analytic precession rate omits disk self-gravity, which the paper itself invokes to explain the resonance migration. These issues must be addressed before the central claims can be considered established.","major_comments":[{"comment":"The test-particle precession rate Ω_E in Eq. (27) includes only the two giant planets and the secondary star, omitting the self-gravity of the planetesimal/embryo disk. Yet Section 4.2.4 states that the disk has a 'regressing effect' on the bodies and that loss of disk mass weakens this effect, causing the secular resonances to migrate inward. This is an explicit admission that disk self-gravity contributes to the secular precession of disk bodies at the epochs where the theory is compared with simulation in Section 4.1. If the disk contribution is comparable to g1 and g2 (order 0.1–0.3 deg/yr in Tables 1–4), the predicted resonance locations are shifted and the visual agreement in Figures 5–6 could be coincidental or dominated by disk-driven precession. Please quantify the disk contribution to Ω_E at the initial and comparison epochs, or include the disk term in Eq. (27), and present a validation run with disk self-gravity disabled.","section":"2.2, Eq. (27) and 4.2.4"},{"comment":"The claimed agreement between theory and simulation is qualitative. The text says secular resonances appear 'either precisely on their predicted locations or in their slight vicinity,' but no quantitative criterion is given for identifying a secular resonance in the simulation, no epoch is stated for the snapshots shown in the bottom panels, and no error bars or scatter measure is provided. Because the resonances migrate inward as the disk loses mass (Section 4.2.4), a comparison between initial theoretical locations and snapshots at unspecified times cannot establish agreement. Please specify the times of the snapshots, define an objective measure of resonance location (e.g., temporary libration of ϖ−ϖ_planet, or eccentricity enhancement above a defined threshold), and compare predicted versus measured locations with uncertainties.","section":"4.1, Figures 5–6"},{"comment":"The suppression claim compares simulations with the secondary star (integrated with the Chambers et al. 2003 code) to simulations without the secondary (integrated with Mercury). Using two different integrators introduces a possible systematic difference, and the footnote attributes differences between the no-secondary panels to running on different computers. More importantly, the disk is present in both panels, so the reduced eccentricity excitation in the with-secondary panels could be caused by disk self-gravity rather than by the secondary's perturbation. Please verify the suppression using a single integrator for both cases and, if possible, with disk self-gravity removed, so that the secondary's effect is isolated.","section":"4.2.2, Figure 8 and footnote"}],"minor_comments":[{"comment":"The phrase 'farther way from the primary' should read 'farther away from the primary.'","section":"Abstract"},{"comment":"There is a double comma in the displayed equation: 'de_E/dt = ... , , dϖ_E/dt'; one comma should be removed.","section":"2.2, Eq. (24)"},{"comment":"The expression for b(0)_1/2(α_jk) contains a stray 'i' in the second hypergeometric term and unbalanced parentheses; please correct the typography.","section":"2.1, Eq. (8)"},{"comment":"The A33 entries for the M-star binary are negative at 20 au but positive at 30, 40, and 50 au. Since A33 is a self-term of the secular matrix it should be positive, and this sign change may indicate a typographical error in the table.","section":"Table 1"},{"comment":"In the concluding paragraph, 'mall bodies' should be 'small bodies.'","section":"Section 5"},{"comment":"The caption spells 'secondary' as 'secodnary'; please correct.","section":"Figure 8 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the scope of the journal and the analytical derivation is standard and internally coherent. The main risk is that the validation is not yet clean: disk self-gravity is omitted from the analytic precession rate despite being invoked to explain resonance migration, and the theory-simulation comparison is qualitative with no epochs or error bars. The suppression claim also relies on a cross-integrator comparison. These issues are addressable with additional analysis or targeted simulations, so major revision rather than rejection seems appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe useful kernel here is a tidy analytical extension: a 3x3 Laplace-Lagrange secular matrix for a binary with two giant planets, yielding resonance locations as intersections of the test-particle precession rate with the g_k eigenvalues. That part is standard but competently assembled, and the specific application to planet-hosting binaries is genuinely new. The qualitative prediction that stronger secondary perturbation moves resonances outward is a direct consequence of the eigenvalues increasing, and it holds up in the algebra.\n\nThe paper does not, however, close the deal on its two more ambitious claims: that the theory is confirmed by N-body simulations, and that the secondary suppresses giant-planet secular resonances. The comparison in Section 4.1 is purely visual—'precisely or in slight vicinity'—with no error bars or quantitative measure of agreement. More troubling, Eq. (27) for the test-particle precession rate includes only the two planets and the secondary star, while the simulations include a massive disk of embryos and planetesimals interior to the inner planet. The paper itself states in §4.2.4 that the disk's 'regressing effect' weakens as mass is lost and that this is what moves the resonances inward. That is an admission that disk self-gravity contributes to Ω_E at the comparison epochs. If that contribution is comparable to g1 or g2 (roughly 0.1–0.3 deg/yr in the tables), the predicted resonance locations shift, and the apparent agreement in Figures 5–6 could be partly coincidental. The suppression claim in §4.2.2 is also confounded: the no-secondary runs used a different integrator (Mercury vs. the special-purpose Chambers et al. code) and different computers, as the footnote admits, so differences in eccentricity excitation could be numeric artifacts. Table typos (e.g., the A13 and A33 signs for G/F binaries) further undermine trust in the numerical results.\n\nThe central analytical mechanism is sound and worth having. The validation is the soft spot, and the disk self-gravity omission is load-bearing for the comparison. Researchers modeling terrestrial planet formation in binaries will get a first-order tool, but they should treat the suppression claim cautiously. This deserves peer review, but I would ask the authors for simulations without disk self-gravity or a quantitative estimate of its contribution to Ω_E, and for a clean control run with the same integrator and machine for the suppression test. With those additions, the theory could stand as a useful tool.","headline":"A standard secular-theory extension with a plausible qualitative prediction, but the simulation validation omits disk self-gravity and the suppression claim rests on uncontrolled comparisons.","tokens_in":14031,"tokens_out":4513,"would_cite":true,"duration_ms":41947,"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":"This paper develops a general theory that predicts where secular resonances from two giant planets fall in a binary star system, showing that a stronger secondary star pushes those resonances outward—and its simulations indicate the…","keywords":["secular resonances","binary stars","terrestrial planet formation","circumstellar planets","secular perturbation theory","protoplanetary disks"],"falsifier":"For a specific binary with known stellar masses, separation, and giant-planet orbits, measure the eccentricity distribution of small bodies in the disk after ~0.1 Myr and compare the location of the eccentricity-excitation band with the predicted $\\Omega_E = g_k$ crossings; alternatively, rerun the same disk simulation with the planetesimal disk's self-gravity artificially turned off and check whether the resonance locations and suppression amplitude change—if they move substantially, the disk's neglected self-gravity is contaminating the agreement.","tokens_in":13032,"feed_emoji":"🪐","tokens_out":3513,"duration_ms":30720,"temperature":0.7,"pith_summary":"The paper asks how the secular resonances of two giant planets—the same kind of Jupiter–Saturn resonances that sculpted the inner solar system—behave when the host star has a close stellar companion. It develops a full four-body secular theory, derives the eigenfrequencies $g_k$ of the planets plus secondary star and the test-particle precession rate $\\Omega_E$, and locates resonances where $\\Omega_E = g_k$. The central predictions are that a more massive or closer secondary star pushes the resonance locations outward, closer to the giant planets and farther from the primary, and that the secondary's perturbation suppresses the dynamical excitation these resonances would otherwise produce. If correct, the theory gives a predictive map of where terrestrial planet formation can proceed in planet-hosting binaries, and it changes how such systems should be simulated.","feed_headline":"A theory predicts where secular resonances lie in binary stars","feed_subtitle":"Stronger secondary stars push the resonances outward and suppress their orbital excitation, shaping where rocky planets can form.","key_machinery":"The central object is the secular perturbation matrix $[A_{jk}]$ of the four-body system, whose eigenvalues $g_k$ are the precession frequencies of the pericenter longitudes of the two giant planets and the secondary star. Its off-diagonal entries mediate the transfer of the secondary's perturbation to the planets, and through the resonance condition $\\Omega_E = g_k$ it sets the radial locations where a small body's pericenter precession matches a planetary mode.","core_discovery":"Working to second order in eccentricity and first order in mass ratio, the paper writes the secular disturbing function of the four-body system (primary, two giant planets, secondary star) in the compact matrix form $R_j = \\tfrac{1}{2} n_j a_j^2 \\sum_k A_{jk} e_j e_k \\cos(\\varpi_j - \\varpi_k)$, with the 3x3 matrix $[A_{jk}]$ given explicitly in terms of Laplace coefficients. The secular precession rates $g_k$ are the eigenvalues of this matrix, and for a test particle the precession rate is $\\Omega_E = \\sum_j \\tfrac{1}{4} n_E (m_j/m_P) \\alpha_{Ej}^2 b_{3/2}^{(1)}(\\alpha_{Ej})$. Secular resonances occur where $\\Omega_E = g_k$. The theory yields two main claims: (1) the magnitude of $g_k$ increases with the secondary's mass and with decreasing binary separation, so stronger secondary perturbation moves the resonance locations to larger semimajor axes—closer to the giant planets and farther from the primary; and (2) in N-body simulations of a protoplanetary disk interior to the inner giant planet, resonances appear at or near their predicted locations, but the secondary star suppresses the resonances' effect, producing weaker eccentricity excitation than the same system without the secondary. The paper also reports that as the disk loses mass, the resonances migrate inward, scattering objects and enhancing collisional growth; this inward migration is attributed to the disk's weakened regressing effect plus slight inward migration of the planets.","pith_inferences":["The theory's resonance-location formula could be tested observationally in systems where the giant planets' eccentricities are growing: the secular resonances' outward shift should manifest as a depletion or excited-eccentricity band whose location tracks the binary's periastron distance, not just its semimajor axis.","The suppression claim, if it holds, implies that terrestrial-planet formation in close binaries may be less impeded by Jupiter-like companions than the solar system's asteroid belt would suggest; one testable extension is to run the same disk simulations with varying secondary eccentricity to quantify how suppression scales with the binary's eccentricity.","The inward migration of resonances during disk mass loss suggests a feedback loop: resonance-driven scattering removes disk mass, which weakens the disk's regressing effect, which moves the resonances inward into denser disk regions, accelerating further mass loss; this loop's strength could be checked by measuring whether the migration stalls when the disk is fully dispersed."],"forward_implications":["Secular resonances in a planet-hosting binary lie at computable semimajor axes, and those locations depend measurably on the binary's mass ratio and separation: stronger secondary perturbation pushes $g_1$ and $g_2$ outward and closer to the giant planets.","More of the protoplanetary disk interior to the inner giant planet remains available for terrestrial planet formation in binaries with a stronger secondary perturbation, because the resonance zones move away from the primary.","The secondary star suppresses the eccentricity excitation of objects captured in the giant planets' secular resonances, so resonance-induced clearing of the disk is weaker than in single-star systems.","As the disk loses mass, secular resonances migrate inward and can scatter planetesimals out of the disk or push them into collisional growth, with the migration rate and timing depending on the disk's initial mass and surface density.","Because the $g_k$ and $\\Omega_E$ are independent of eccentricity, the predicted resonance locations hold for eccentric planets and eccentric binaries as well."],"supporting_citations":[{"why":"supplies the second-order-in-eccentricity secular disturbing function terms with Laplace coefficients used to build the matrix $[A_{jk}]$.","marker":"Ellis & Murray 2000"},{"why":"provides the special-purpose N-body integrator used for the binary-star disk simulations that confirm the predicted resonance locations.","marker":"Chambers et al. 2003"},{"why":"is the single-star benchmark for secular resonance locations and migration; the paper compares its mass-loss and resonance-migration behavior against that work's Figure 8.","marker":"Haghighipour & Winter 2016"},{"why":"supplies the Mercury integrator used for the comparison simulations without the secondary star in Section 4.2.2.","marker":"Chambers 1999"},{"why":"motivates the relevance of secular resonances to the architecture of the asteroid belt and hence to terrestrial-planet formation.","marker":"Milani & Knežević 1992"},{"why":"grounds the stability constraint that keeps the two giant planets inside the binary's stable zone for the chosen orbital parameters.","marker":"Quarles et al. 2020"}],"fun_headline_variants":["Binary companion shifts secular resonances outward","How secondary stars redirect planet-forming resonances","Suppressed resonances: binary stars alter planet formation","Disk mass loss sends resonances inward in binaries","Secular resonance theory for dual-star planet systems"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analytical precession rate $\\Omega_E$ of a small body ignores the self-gravity of the planetesimal disk, even though the simulations include the disk and the paper itself credits the disk with a 'regressing effect' that shifts the resonances inward as mass is lost.","fun_headline_variants_meta":{"raw":{"variants":["Binary companion shifts secular resonances outward","How secondary stars redirect planet-forming resonances","Suppressed resonances: binary stars alter planet formation","Disk mass loss sends resonances inward in binaries","Secular resonance theory for dual-star planet systems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000254,"raw_usage":{"total_tokens":1666,"prompt_tokens":1139,"completion_tokens":527,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":755,"completion_tokens_details":{"reasoning_tokens":458}},"tokens_in":755,"tokens_out":527,"duration_ms":6049,"temperature":1.0,"reasoning_tokens":458,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:57:28.506118+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a specific binary with known stellar masses, separation, and giant-planet orbits, measure the eccentricity distribution of small bodies in the disk after ~0.1 Myr and compare the location of the eccentricity-excitation band with the predicted $\\Omega_E = g_k$ crossings; alternatively, rerun the same disk simulation with the planetesimal disk's self-gravity artificially turned off and check whether the resonance locations and suppression amplitude change—if they move substantially, the disk's neglected self-gravity is contaminating the agreement.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the second-order-in-eccentricity secular disturbing function terms with Laplace coefficients used to build the matrix $[A_{jk}]$."},{"cited_title":"E., Quintana, E","cited_arxiv_id":null,"evidence_quote":"provides the special-purpose N-body integrator used for the binary-star disk simulations that confirm the predicted resonance locations."},{"cited_title":"& Winter, O","cited_arxiv_id":null,"evidence_quote":"is the single-star benchmark for secular resonance locations and migration; the paper compares its mass-loss and resonance-migration behavior against that work's Figure 8."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the Mercury integrator used for the comparison simulations without the secondary star in Section 4.2.2."},{"cited_title":"& Kne z evi \\' c , Z","cited_arxiv_id":null,"evidence_quote":"motivates the relevance of secular resonances to the architecture of the asteroid belt and hence to terrestrial-planet formation."},{"cited_title":"& Haghighipour, N","cited_arxiv_id":null,"evidence_quote":"grounds the stability constraint that keeps the two giant planets inside the binary's stable zone for the chosen orbital parameters."}],"review_version":1}