Planetary Desert around Compact Binaries: Dynamical Instability Triggered by Resonance-Induced Eccentricity Excitation
Pith reviewed 2026-05-15 18:39 UTC · model grok-4.3
The pith
Resonance with a decaying binary pumps planet eccentricities until orbits cross and scatter, emptying inner regions around short-period binaries.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
When an eccentric binary decays via tides, an outer planet can be captured into resonance advection in eccentricity, a state in which its apsidal precession locks with that of the binary, driving extreme eccentricity growth. While such growth can occur in a binary-single planet system, the parameter space is limited and may not necessarily induce instability. In a multi-planet system, however, the excited orbit inevitably crosses those of its neighbors, which triggers violent planet-planet scatterings and produces collisions or ejections. These mutual gravitational interactions amplify the localized instability into a system-wide chain reaction that drastically reshapes the orbital架构 and can
What carries the argument
resonance advection in eccentricity, in which a planet's apsidal precession locks with the binary's during tidal decay, pumping the planet's eccentricity to high values while the binary shrinks.
If this is right
- Mutual planet-planet interactions turn limited single-planet eccentricity growth into system-wide clearing.
- The process operates specifically when the binary is eccentric and decaying through tides.
- Inner regions of multi-planet systems around compact binaries are preferentially emptied by the resulting scatterings.
- Single-planet systems may remain stable even after the resonance-driven eccentricity increase.
Where Pith is reading between the lines
- Surviving planets at larger separations could retain eccentricity or inclination signatures from earlier resonance passages.
- The mechanism predicts a sharp inner cutoff in planet occurrence that could be checked against transit survey statistics.
- Higher-order effects or general-relativistic precession might shift the resonance locations and should be tested in refined models.
- Similar resonance-driven clearing might appear around other classes of evolving binaries with tidal decay.
Load-bearing premise
The single-averaged secular equations remain valid and sufficient throughout the eccentricity excitation and scattering phases without requiring full N-body integration or higher-order corrections.
What would settle it
A full N-body integration of the same multi-planet circumbinary systems that shows no orbit crossings or scatterings despite the predicted resonance capture and eccentricity growth would falsify the instability mechanism.
Figures
read the original abstract
Compact binaries with orbital periods shorter than about 7 days show an absence of transiting planets, a feature known as the ``circumbinary planet desert". The physical mechanism behind this desert remains unclear. We investigate its origin by simulating the long-term dynamics of multi-planet circumbinary systems with evolving inner binaries. Our simulations are based on the single-averaged secular equations that average only over the binary orbital period and fully incorporate planet-planet interactions. When an eccentric binary decays via tides, an outer planet can be captured into resonance advection in eccentricity, a state in which its apsidal precession locks with that of the binary, driving extreme eccentricity growth. While such growth can occur in a binary-single planet system, the parameter space is limited and may not necessarily induce instability. In a multi-planet system, however, the excited orbit inevitably crosses those of its neighbors, which triggers violent planet-planet scatterings and produces collisions or ejections. Crucially, these mutual gravitational interactions amplify the ``localized" instability of a single planet into a system-wide chain reaction, drastically reshaping the orbital architecture and potentially clearing out the inner regions of planetary systems. Our results suggest that the resonance-induced instability provides a natural explanation for the observed circumbinary planet desert.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that the observed circumbinary planet desert around compact binaries (P_bin < 7 days) originates from resonance-induced eccentricity excitation. Using single-averaged secular equations that incorporate planet-planet interactions and tidal binary decay, the authors find that apsidal resonance advection drives extreme eccentricity growth; in single-planet systems this effect is limited, but in multi-planet systems it triggers orbit crossings, violent scatterings, collisions, and ejections that clear the inner regions.
Significance. If the mechanism is robust, the work supplies a dynamical explanation for a clear observational feature in circumbinary exoplanets and illustrates how mutual gravitational interactions can amplify localized instabilities into system-wide clearing. The long-term secular approach is computationally efficient for exploring evolutionary timescales, which is a methodological strength.
major comments (1)
- [Methods (single-averaged secular equations) and abstract] The central results rest on the single-averaged secular equations remaining valid through resonance capture, eccentricity excitation to orbit-crossing values, and subsequent scattering. The abstract and methods description state that these equations average only over the binary period while retaining planet-planet terms, yet standard secular theory assumes moderate eccentricities and slow evolution; at high e the neglected short-period terms and close-encounter dynamics become non-perturbative. No explicit comparison to full N-body integrations is reported for the high-e regime, which is load-bearing for the claim that multi-planet amplification clears the inner disk.
minor comments (2)
- [Abstract] The abstract would benefit from a brief statement of the explored ranges in binary period, planet mass, and number of planets to allow readers to assess the generality of the reported instability.
- [Methods] Notation for the resonance advection state and the definition of the secular Hamiltonian could be clarified with an explicit equation reference in the methods.
Simulated Author's Rebuttal
We thank the referee for their constructive comments on our manuscript. We address the major concern regarding the validity of the single-averaged secular equations below and have revised the manuscript to strengthen the presentation.
read point-by-point responses
-
Referee: [Methods (single-averaged secular equations) and abstract] The central results rest on the single-averaged secular equations remaining valid through resonance capture, eccentricity excitation to orbit-crossing values, and subsequent scattering. The abstract and methods description state that these equations average only over the binary period while retaining planet-planet terms, yet standard secular theory assumes moderate eccentricities and slow evolution; at high e the neglected short-period terms and close-encounter dynamics become non-perturbative. No explicit comparison to full N-body integrations is reported for the high-e regime, which is load-bearing for the claim that multi-planet amplification clears the inner disk.
Authors: We agree that the single-averaged secular equations have limitations at high eccentricities, where short-period terms and close-encounter dynamics are no longer perturbative. Our approach uses these equations to follow the slow secular evolution and resonance advection during binary tidal decay, which drives eccentricity growth until orbits cross in multi-planet systems. Once crossings occur, the resulting instabilities and clearing are a generic outcome of orbit crossing. To address the concern, we will add a new subsection discussing the validity range of the secular approximation and include direct comparisons with full N-body integrations for representative cases at the onset of high eccentricity to confirm the transition to instability. revision: yes
Circularity Check
No significant circularity: forward secular simulations yield independent dynamical outcomes
full rationale
The paper's derivation consists of forward numerical integration of single-averaged secular equations that incorporate planet-planet interactions and binary tidal decay. The resonance-capture and eccentricity-excitation behavior, followed by scattering in multi-planet cases, emerges directly from the time evolution of those equations rather than from any parameter fit, self-definition, or self-citation chain that would force the final architecture. No quoted step equates a claimed prediction to an input by construction; the central claim that such instability clears the inner region is therefore an output of the dynamics, not a renaming or tautology.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Single-averaged secular equations accurately capture the long-term dynamics including planet-planet interactions and resonance effects.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Our simulations are based on the single-averaged secular equations that average only over the binary orbital period and fully incorporate planet-planet interactions.
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
When an eccentric binary decays via tides, an outer planet can be captured into resonance advection in eccentricity...
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Armstrong, D. J., Osborn, H. P., Brown, D. J. A., et al. 2014, MNRAS,444, 1873
work page 2014
-
[2]
Bate, M. R. 2009, MNRAS,392, 590
work page 2009
-
[3]
Bate, M. R. 2019, MNRAS,484, 2341
work page 2019
-
[4]
Baycroft, T. A., Santerne, A., Triaud, A. H. M. J., et al. 2025, MNRAS,541, 2801
work page 2025
- [5]
- [6]
-
[7]
Borkovits, T., Hajdu, T., Sztakovics, J., et al. 2016, MNRAS, 455, 4136
work page 2016
-
[8]
Bryan, M. L., Bowler, B. P., Knutson, H. A., et al. 2016, ApJ, 827, 100
work page 2016
-
[9]
Burgasser, A. J., Simcoe, R. A., Bochanski, J. J., et al. 2010, ApJ,725, 1405
work page 2010
- [10]
-
[11]
Chen, C., Martin, R. G., Lubow, S. H., & Nixon, C. J. 2024, ApJL,961, L5
work page 2024
- [12]
-
[13]
Coleman, G. A. L. 2024, MNRAS,530, 630
work page 2024
-
[14]
Cuello, N., & Giuppone, C. A. 2019, A&A,628, A119
work page 2019
-
[15]
Correia, A. C. M., Udry, S., Mayor, M., et al. 2005, A&A,440, 751 de El´ ıa, G. C., Zanardi, M., Dugaro, A., & Naoz, S. 2019, A&A, 627, A17
work page 2005
- [16]
- [17]
-
[18]
Doolin, S., & Blundell, K. M. 2011, MNRAS,418, 2656
work page 2011
-
[19]
Doyle, L. R., Carter, J. A., Fabrycky, D. C., et al. 2011, Science, 333, 1602
work page 2011
- [20]
- [21]
-
[22]
Eggleton, P. P., & Kiseleva-Eggleton, L. 2001, ApJ,562, 1012
work page 2001
- [23]
- [24]
-
[25]
H., Naoz, S., Wei, L., & Farr, W
Faridani, T. H., Naoz, S., Wei, L., & Farr, W. M. 2022, ApJ,932, 78
work page 2022
-
[26]
Fleming, D. P., Barnes, R., Graham, D. E., et al. 2018, ApJ,858, 86
work page 2018
- [27]
-
[28]
K., Carter, B., King, R., & O’Toole, S
Getley, A. K., Carter, B., King, R., & O’Toole, S. 2017, MNRAS, 468, 2932
work page 2017
- [29]
- [30]
-
[31]
Hamers, A. S., Perets, H. B., & Portegies Zwart, S. F. 2016, MNRAS,455, 3180
work page 2016
- [32]
-
[33]
Hansen, B. M. S., & Naoz, S. 2020, MNRAS,499, 1682
work page 2020
- [34]
-
[35]
Ivanova, N., & Taam, R. E. 2003, ApJ,599, 516
work page 2003
-
[36]
Janson, M., Asensio-Torres, R., Andr´ e, D., et al. 2019, A&A, 626, A99
work page 2019
- [37]
- [38]
-
[39]
Kley, W., & Nelson, R. P. 2012, ARA&A,50, 211
work page 2012
- [40]
- [41]
- [42]
- [43]
- [44]
-
[45]
Leung, G. C. K., & Lee, M. H. 2013, ApJ,763, 107
work page 2013
- [46]
- [47]
- [48]
-
[49]
Lidov, M. L. 1962, P&SS,9, 719
work page 1962
- [50]
- [51]
-
[52]
Liu, B., D’Orazio, D. J., Vigna-G´ omez, A., et al. 2022, PhRvD, 106, 123010
work page 2022
- [53]
- [54]
- [55]
- [56]
- [57]
-
[58]
Moe, M., & Kratter, K. M. 2018, ApJ,854, 44
work page 2018
-
[59]
Mogan, S., & Zanazzi, J. J. 2025, ApJ,988, 150 M¨ uller, S., Baron, J., Helled, R., et al. 2024, A&A,686, A296 Mu˜ noz, D. J., & Lai, D. 2015, PNAS,112, 9264
work page 2025
-
[60]
Naoz, S., & Fabrycky, D. C. 2014, ApJ,793, 137
work page 2014
-
[61]
Naoz, S., Farr, W. M., Lithwick, Y., et al. 2013, MNRAS,431, 2155
work page 2013
- [62]
-
[63]
Nelson, R. P. 2003, MNRAS,345, 233
work page 2003
- [64]
-
[65]
Penzlin, A. B., Kley, W., & Nelson, R. P. 2021, A&A,645, A68
work page 2021
- [66]
- [67]
- [68]
- [69]
- [70]
-
[71]
Standing, M. R., Sairam, L., Martin, D. V., et al. 2023, NatAs,7, 702
work page 2023
- [72]
- [73]
- [74]
-
[75]
R., Tremaine, S., & Kazandjian, M
Touma, J. R., Tremaine, S., & Kazandjian, M. V. 2009, MNRAS, 394, 1085
work page 2009
- [76]
-
[77]
Vick, M., Lai, D., & Anderson, K. R. 2019, MNRAS,484, 5645
work page 2019
-
[78]
Vinson, B. R., & Chiang, E. 2018, MNRAS,474, 4855 von Zeipel, H. 1910, AN,183, 345
work page 2018
-
[79]
Wei, L., Naoz, S., Faridani, T., & Farr, W. M. 2021, ApJ,923, 118
work page 2021
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.