Pith. sign in

REVIEW 3 major objections 6 minor 12 references

Secular Resonances in Planet-Hosting Binary Stars. I. General Theory

T0 review · 3 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read 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…

desk verdict 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. read the letter →

arxiv 2507.17092 v1 pith:R2QFVRBX submitted 2025-07-23 astro-ph.EP

Nader Haghighipour , Michael Andrew This is my paper · ORCID
classification astro-ph.EP
keywords secularresonancesbinarystarsterrestrialplanetformationcircumstellarplanetsperturbationtheoryprotoplanetarydisks
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

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.

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 (3)
  1. [2.2, Eq. (27) and 4.2.4] 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.
  2. [4.1, Figures 5–6] 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.
  3. [4.2.2, Figure 8 and footnote] 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.
minor comments (6)
  1. [Abstract] The phrase 'farther way from the primary' should read 'farther away from the primary.'
  2. [2.2, Eq. (24)] There is a double comma in the displayed equation: 'de_E/dt = ... , , dϖ_E/dt'; one comma should be removed.
  3. [2.1, Eq. (8)] The expression for b(0)_1/2(α_jk) contains a stray 'i' in the second hypergeometric term and unbalanced parentheses; please correct the typography.
  4. [Table 1] 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.
  5. [Section 5] In the concluding paragraph, 'mall bodies' should be 'small bodies.'
  6. [Figure 8 caption] The caption spells 'secondary' as 'secodnary'; please correct.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: resonance locations derive from an independent eigenvalue calculation, and self-citations are contextual only.

full rationale

The derivation of resonance locations is self-contained: g_k are eigenvalues of the matrix [A_jk] in Eqs. (15)-(23), and Ω_E is given by Eq. (27) as a sum of secular precession rates from the two planets and the secondary star; both are evaluated directly from masses and semimajor axes with no fitted parameters. The resonance condition Ω_E = g_k is therefore a genuine calculation, not an input. The validating N-body integrations use the same physical parameters, so visual agreement in Figures 5-6 is a test of the analytic model rather than a circular confirmation. The self-citations (Haghighipour & Winter 2016; Quarles et al. 2020, 2024) supply context for disk mass-loss rate comparisons and planetary stability, but the central location prediction does not depend on them. The omission of disk self-gravity from Eq. (27), despite the 'regressing effect' of the disk invoked in Section 4.2.4, is a physical completeness/correctness concern that could affect the agreement, but it is not a circular reduction because the theory neither fits nor derives the disk contribution from the quantity it predicts. No step reduces to its own input by construction.

Assumptions & free parameters 5 free parameters · 6 assumptions · 0 invented entities

The central theory relies on standard secular perturbation theory (Laplace-Lagrange) with Laplace coefficients and a low-order expansion in eccentricity; these are pulled from prior literature (Ellis & Murray 2000). No numbers are fitted to data. The application to specific binaries uses chosen physical parameters (planet masses and semimajor axes, disk profile) that are not derived from the theory. No new entities such as new particles or forces are introduced.

free parameters (5)
  • inner planet semimajor axis = 1.6 au
    Chosen to keep planets inside the binary stable zone and to mimic the Jupiter-Saturn 5:2 period ratio (Section 3).
  • outer planet semimajor axis = 2.94 au
    Chosen together with the inner value to match Jupiter-Saturn period ratio and stability (Section 3).
  • planet masses = 1 M_Jup, 1 M_Sat
    Set equal to Jupiter and Saturn for comparison with solar system secular resonance studies (Section 3).
  • disk surface density exponent = -1.5
    Adopted for the simulated disk; not derived from theory (Section 4.1).
  • disk inner and outer boundaries = 0.5 au, 1.5 au
    Extent of the planetesimal disk in simulations (Section 4.1).
assumptions (6)
  • domain assumption Coplanarity of binary, planets, and disk
    The paper assumes a coplanar system to apply the lowest-order secular theory (Section 2); inclination effects are postponed to Paper II.
  • domain assumption No mean-motion resonances between the two planets
    Secular terms are selected as independent of mean longitudes (Section 2.1); a MMR would invalidate this separation.
  • standard math Second-order-in-eccentricity, first-order-in-mass-ratio expansion
    The disturbing function is expanded to second order in e and first order in mass ratio (Section 2.1), the standard Laplace-Lagrange approximation.
  • domain assumption Planets remain stable under the secondary's perturbation
    The paper assumes both giant planets maintain stable orbits in the binary (Section 2, citing Quarles et al. 2020, 2024).
  • domain assumption Secondary star treated as a distant exterior perturber
    The secondary enters the disturbing function in the same form as an outer planet (Eqs. 11-13), which is valid for small alpha but ignores higher-order binary effects.
  • domain assumption Negligible disk self-gravity in analytical precession
    Equation (27) for Omega_E omits the disk's gravity, while simulations include the disk; the paper acknowledges a 'regressing effect' of the disk (Section 4.2.4) without quantifying it.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Secular Resonances in Planet-Hosting Binary Stars. I. General Theory." pith.science (2026). https://pith.science/paper/R2QFVRBX

@misc{pith2026250717092,
  author       = {Pith},
  title        = {Pith review of: Secular Resonances in Planet-Hosting Binary Stars. I. General Theory},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/R2QFVRBX}},
  note         = {Machine review of arXiv:2507.17092}
}
read the original abstract

Motivated by the diversity of circumstellar planets in binary stars and the strong effects of the secular resonances of Jupiter and Saturn on the formation and architecture of the inner solar system, we have launched an expansive project on studying the effects of secular resonances on the formation of terrestrial planets around a star of a moderately close binary. As the first phase of our project, we present here the general theory of secular resonances in dual-star systems where the primary hosts two giant planets. Using the concept of generalized disturbing function, we derive the formula for the locations of secular resonances and show that in systems where the perturbation of the secondary star is stronger, the locations of secular resonances are farther way from the primary and closer to the giant planets. The latter implies that in such systems, terrestrial planet formation has a larger area to proceed with more of the protoplanetary disk being available to it. To demonstrate the validity of our theoretical results, we simulated the evolution of a protoplanetary disk interior to the inner giant planet. Results, in addition to confirming our theoretical predictions, pointed to an important finding: In binary stars, the perturbation of the secondary suppresses the secular resonances of giant planets. Simulations also show that as the disk loses material, secular resonances move inward, scattering objects out of the disk and/or facilitating their collisional growth. We present results of our study and discuss their implications for the simulations of terrestrial planet formation.

Figures

Figures reproduced from arXiv: 2507.17092 by the authors.

Figure 1
Figure 1. — Schematic presentation of the system showing the binary star and the two planets [PITH_FULL_IMAGE:figures/full_fig_p016_1.png] view at source ↗
Figure 2
Figure 2. — Graphs of the variation of the frequency [PITH_FULL_IMAGE:figures/full_fig_p017_2.png] view at source ↗
Figure 3
Figure 3. — Locations of the secular resonances of the inner and outer planets, [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
Figures from the paper (6 more)
Figure 4
Figure 4. Figure 4: — Locations of the secular resonances of the inner and outer planets, [PITH_FULL_IMAGE:figures/full_fig_p019_4.png]
Figure 5
Figure 5. Figure 5: — Secular resonances in a GM and GK binary with a semimajor axis of 20 au. The [PITH_FULL_IMAGE:figures/full_fig_p020_5.png]
Figure 6
Figure 6. Figure 6: — Secular resonances in a GG and GF binary with a semimajor axis of 20 au. The [PITH_FULL_IMAGE:figures/full_fig_p021_6.png]
Figure 7
Figure 7. Figure 7: — Graph of the evolution of the primary’s protoplanetry disk in a binary without [PITH_FULL_IMAGE:figures/full_fig_p022_7.png]
Figure 8
Figure 8. Figure 8: — Locations of secular resonances in a 20 au GK (top) and GF (bottom) binary [PITH_FULL_IMAGE:figures/full_fig_p023_8.png]
Figure 9
Figure 9. Figure 9: — The evolution of the mass of the protoplanetary disk in the binaries of Figures 5 [PITH_FULL_IMAGE:figures/full_fig_p024_9.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

12 extracted references · 10 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

  3. [3]

    Chambers, J. E. 1999, , 304, 793

    work page 1999

  4. [4]

    E., Quintana, E

    Chambers, J. E., Quintana, E. V., Duncan, M. J. & Lissauer, J. J. 2003, , 123, 2884

    work page 2003

  5. [5]

    Ellis, K. M. & Murray, C. D. 2000, Icarus 147, 129

    work page 2000

  6. [6]

    & Winter, O

    Haghighipour, N. & Winter, O. C. 2016, Cel. Mech. Dyn. Ast., 124, 235

    work page 2016

  7. [7]

    Levison, H. F. & Agnor, c. 2003, , 125, 2692

    work page 2003

  8. [8]

    & Loeb, A

    Li, G., Naoz, S., Kocsis, B. & Loeb, A. 2014, , 785, 116

    work page 2014

Show all 12 references
  1. [9]

    & Kne z evi \' c , Z

    Milani, A. & Kne z evi \' c , Z. 1992, Icarus, 98, 211

    1992

  2. [10]

    & Kne z evi \' c , Z

    Milani, A. & Kne z evi \' c , Z. 1994, Icarus, 107, 219

    1994

  3. [11]

    & Haghighipour, N

    Quarles, B., Li, G., Kostov, V. & Haghighipour, N. 2020, , 159, 80

    2020

  4. [12]

    Quarles, B., Bhaskar, H. G. & Li, G. 2024, in: Encyclopedia of Astrophysics, Elsevier Pub. Co. (arXiv:2407.13901)

    2024 arXiv

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.