The S stars' zone of avoidance in the Galactic center
Pith reviewed 2026-05-23 07:54 UTC · model grok-4.3
The pith
Binary disruption near the supermassive black hole plus orbital relaxation reproduces the S stars' zone of avoidance and thermal eccentricities.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The observed S-star orbital distributions, including the zone of avoidance and their thermal eccentricity distribution, can be largely explained by the continuous disruption of binaries near the central supermassive black hole, followed by orbital relaxation.
What carries the argument
Continuous binary disruption near the supermassive black hole followed by nonresonant and resonant relaxation of the remnant stars.
If this is right
- S stars are the surviving members of binaries that reached the black hole on highly eccentric orbits.
- The zone of avoidance and eccentricity distribution together constrain the initial binary population and the relaxation rates in the central parsec.
- S-star statistics become a diagnostic of the dynamical environment within a few parsecs of the black hole.
- The same process operating continuously can maintain the observed S-star population against orbital decay.
Where Pith is reading between the lines
- Similar orbital gaps may appear around supermassive black holes in other galaxies if binary fractions and relaxation operate on comparable scales.
- Tighter constraints on the present-day S-star sample could back out the binary separation distribution at 5-100 pc.
- Extending the models to include stellar evolution or resonant relaxation with a stellar cusp could predict additional observable signatures such as mass segregation.
Load-bearing premise
Empirical binary-property distributions taken from 5-100 pc scales and the simplified treatments of nonresonant and resonant relaxation are accurate enough that their output matches the observed zone of avoidance.
What would settle it
Discovery of even one S star with pericenter inside the zone of avoidance, or a statistically significant departure from a thermal eccentricity distribution among the full sample, would falsify the explanation.
Figures
read the original abstract
This paper investigates the origin and orbital evolution of S stars in the Galactic center using models of binary disruption and relaxation processes. We focus on explaining the recently discovered ``zone of avoidance'' in S-star orbital parameters, defined as a region where no S stars are observed with pericenters of $\log(r_p / {\rm AU}) \leq 1.57 + 2.6(1 - e)$ pc. We demonstrate that the observed S-star orbital distributions, including this zone of avoidance and their thermal eccentricity distribution, can be largely explained by the continuous disruption of binaries near the central supermassive black hole, followed by orbital relaxation. Our models consider binaries originating from large scales ($5$--$100$ pc) and incorporate empirical distributions of binary properties. We simulate close encounters between binaries and the black hole, tracking the remnant stars' orbits. The initially highly eccentric orbits of disrupted binary remnants evolve due to nonresonant and resonant relaxation in the Galactic center potential. While our results provide insights into the formation mechanism of S stars, there are limitations, such as uncertainties in the initial binary population and mass function and simplifications in our relaxation models. Despite these caveats, our study demonstrates the power of using S-star distributions to probe the dynamical history and environment of the central parsec of our Galaxy.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that the observed S-star orbital distributions, including the recently identified 'zone of avoidance' defined by log(r_p/AU) ≤ 1.57 + 2.6(1-e) and a thermal eccentricity distribution, can be largely reproduced by forward simulations of binary disruptions near the central supermassive black hole (with binaries drawn from empirical properties at 5-100 pc scales) followed by orbital evolution under nonresonant and resonant relaxation.
Significance. If the central claim holds, the work provides a physically motivated explanation for S-star properties via continuous binary disruption plus relaxation, using forward modeling rather than direct fitting to the target data. This is a strength, as it generates testable predictions from dynamical processes with empirical inputs.
major comments (2)
- [Abstract] Abstract: The central claim that the models 'can be largely explained' the zone of avoidance and thermal eccentricity distribution lacks any quantitative goodness-of-fit metrics (e.g., KS statistic, chi-squared, or posterior predictive checks) or error analysis on the simulated distributions; without these, the match cannot be assessed given the listed limitations in binary population, mass function, and relaxation.
- [Abstract] Abstract (relaxation models): The reproduction of the specific avoidance boundary hinges on the simplified nonresonant and resonant relaxation prescriptions; the abstract explicitly flags these simplifications, but no sensitivity tests or comparisons to treatments including resonant relaxation saturation, mass-spectrum dependence, or triaxiality are described, raising the risk that the carved-out region is an artifact of the chosen implementation.
minor comments (1)
- [Abstract] Abstract: The zone-of-avoidance inequality is written as log(r_p/AU) ≤ 1.57 + 2.6(1-e) pc; the trailing 'pc' appears to be a unit error and should be removed for clarity.
Simulated Author's Rebuttal
We thank the referee for their constructive comments on our manuscript. We address each major comment below.
read point-by-point responses
-
Referee: [Abstract] Abstract: The central claim that the models 'can be largely explained' the zone of avoidance and thermal eccentricity distribution lacks any quantitative goodness-of-fit metrics (e.g., KS statistic, chi-squared, or posterior predictive checks) or error analysis on the simulated distributions; without these, the match cannot be assessed given the listed limitations in binary population, mass function, and relaxation.
Authors: We agree that the current version relies on qualitative visual comparison rather than formal statistical metrics. The models reproduce the zone of avoidance boundary and thermal eccentricity distribution, but without KS tests or error bars from parameter variations the strength of the match is difficult to quantify. We will add KS statistics between simulated and observed distributions plus uncertainty estimates in the revised manuscript. revision: yes
-
Referee: [Abstract] Abstract (relaxation models): The reproduction of the specific avoidance boundary hinges on the simplified nonresonant and resonant relaxation prescriptions; the abstract explicitly flags these simplifications, but no sensitivity tests or comparisons to treatments including resonant relaxation saturation, mass-spectrum dependence, or triaxiality are described, raising the risk that the carved-out region is an artifact of the chosen implementation.
Authors: The abstract does note the simplifications in the relaxation treatment. We performed limited internal checks on relaxation timescale variations during model development but did not include a dedicated sensitivity analysis in the manuscript. We will add a short discussion of robustness to changes in the nonresonant and resonant relaxation prescriptions in the revised version; however, a full treatment of resonant relaxation saturation, mass-spectrum effects, and triaxiality lies outside the scope of this work. revision: partial
Circularity Check
No significant circularity; forward modeling with external empirical inputs
full rationale
The paper's central result is obtained via forward Monte Carlo simulations of binary disruption near the SMBH followed by nonresonant and resonant relaxation, using empirical binary-property distributions drawn from 5-100 pc scales as inputs. The zone of avoidance and thermal eccentricity distribution are reported as emergent outcomes rather than quantities to which any model parameter is fitted. The abstract explicitly flags uncertainties in those inputs and simplifications in the relaxation treatment, confirming that the match is not enforced by construction. No self-definitional equations, fitted-input predictions, or load-bearing self-citations appear in the derivation chain.
Axiom & Free-Parameter Ledger
free parameters (1)
- initial binary population and mass function
axioms (2)
- domain assumption Binaries originate from large scales (5-100 pc) with empirical property distributions
- domain assumption Nonresonant and resonant relaxation operate as modeled in the Galactic center potential
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...
-
[2]
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
- [3]
-
[4]
2024, arXiv e-prints, arXiv:2410.19881
Akiba , T., Naoz , S., & Madigan , A.-M. 2024, arXiv e-prints, arXiv:2410.19881
-
[5]
Stellar Dynamics and Stellar Phenomena Near A Massive Black Hole
Alexander , T. 2017, ArXiv e-prints [ [arXiv] 1701.04762 ]
work page internal anchor Pith review Pith/arXiv arXiv 2017
- [6]
- [7]
- [8]
- [9]
-
[10]
Antonini , F., Faber , J., Gualandris , A., & Merritt , D. 2010, , 713, 90
work page 2010
- [11]
-
[12]
Astropy Collaboration , Price-Whelan , A. M., Lim , P. L., et al. 2022, , 935, 167
work page 2022
-
[13]
Astropy Collaboration , Robitaille , T. P., Tollerud , E. J., et al. 2013, , 558, A33
work page 2013
-
[14]
Bahcall , J. N. & Wolf , R. A. 1976, , 209, 214
work page 1976
- [15]
-
[16]
Bar-Or , B. & Alexander , T. 2014, Classical and Quantum Gravity, 31, 244003
work page 2014
- [17]
- [18]
- [19]
-
[20]
Baumgardt , H., Amaro-Seoane , P., & Sch \"o del , R. 2018, , 609, A28
work page 2018
-
[21]
Bonnell , I. A. & Rice , W. K. M. 2008, Science, 321, 1060
work page 2008
- [22]
- [23]
- [24]
-
[25]
Cao , C. Y., Liu , F. K., Li , S., Chen , X., & Wang , K. 2024, arXiv e-prints, arXiv:2411.09278
- [26]
-
[27]
2015, Monthly Notices of the Royal Astronomical Society, 452, 1068
Chen, Y., Bressan, A., Girardi, L., et al. 2015, Monthly Notices of the Royal Astronomical Society, 452, 1068
work page 2015
-
[28]
2014, Monthly Notices of the Royal Astronomical Society, 444, 2525
Chen, Y., Girardi, L., Bressan, A., et al. 2014, Monthly Notices of the Royal Astronomical Society, 444, 2525
work page 2014
- [29]
-
[30]
Collaboration, T. G., Dayem, K. A. E., Abuter, R., et al. 2024, Improving constraints on the extended mass distribution in the Galactic Center with stellar orbits
work page 2024
-
[31]
Dale , J. E., Davies , M. B., Church , R. P., & Freitag , M. 2009, , 393, 1016
work page 2009
-
[32]
B., Blackwell , R., Bailey , V
Davies , M. B., Blackwell , R., Bailey , V. C., & Sigurdsson , S. 1998, , 301, 745
work page 1998
-
[33]
Dessart , L., Ryu , T., Amaro Seoane , P., & Taylor , A. M. 2024, , 682, A58
work page 2024
-
[34]
Do , T., Hees , A., Ghez , A., et al. 2019, Science, 365, 664
work page 2019
- [35]
- [36]
- [37]
-
[38]
Duncan , M. J. & Shapiro , S. L. 1983, , 268, 565
work page 1983
-
[39]
Feldmeier-Krause , A., Kerzendorf , W., Neumayer , N., et al. 2017, , 464, 194
work page 2017
-
[40]
B., Pichon , C., & Magorrian , J
Fouvry , J. B., Pichon , C., & Magorrian , J. 2017, , 598, A71
work page 2017
-
[41]
Fragione , G., Capuzzo-Dolcetta , R., & Kroupa , P. 2017, , 467, 451
work page 2017
- [42]
-
[43]
M., Cheung , P., Portegies Zwart , S
Fregeau , J. M., Cheung , P., Portegies Zwart , S. F., & Rasio , F. A. 2004, , 352, 1
work page 2004
-
[44]
Freitag , M., Amaro-Seoane , P., & Kalogera , V. 2006 a , , 649, 91
work page 2006
-
[45]
Freitag , M., G \"u rkan , M. A., & Rasio , F. A. 2006 b , , 368, 141
work page 2006
-
[46]
K., Gillessen , S., Dodds-Eden , K., et al
Fritz , T. K., Gillessen , S., Dodds-Eden , K., et al. 2011, , 737, 73
work page 2011
-
[47]
Gallego-Cano , E., Fritz , T., Sch \"o del , R., et al. 2024, , 689, A190
work page 2024
-
[48]
Gallego-Cano , E., Sch \"o del , R., Dong , H., et al. 2018, , 609, A26
work page 2018
- [49]
- [50]
-
[51]
Generozov , A., Nayakshin , S., & Madigan , A. M. 2022, , 512, 4100
work page 2022
-
[52]
Generozov , A., Stone , N. C., Metzger , B. D., & Ostriker , J. P. 2018, , 478, 4030
work page 2018
-
[53]
Genzel , R., Eckart , A., Ott , T., & Eisenhauer , F. 1997, , 291, 219
work page 1997
-
[54]
Genzel , R., Thatte , N., Krabbe , A., Kroker , H., & Tacconi-Garman , L. E. 1996, , 472, 153
work page 1996
-
[55]
M., Duch \^e ne , G., Matthews , K., et al
Ghez , A. M., Duch \^e ne , G., Matthews , K., et al. 2003, , 586, L127
work page 2003
-
[56]
Ghez , A. M., Klein , B. L., Morris , M., & Becklin , E. E. 1998, , 509, 678
work page 1998
-
[57]
Gillessen , S., Plewa , P. M., Eisenhauer , F., et al. 2017, , 837, 30
work page 2017
- [58]
- [59]
-
[60]
GRAVITY Collaboration , Abuter , R., Aimar , N., et al. 2022, , 657, L12
work page 2022
-
[61]
GRAVITY Collaboration , Abuter , R., Amorim , A., et al. 2018, , 615, L15
work page 2018
-
[62]
GRAVITY Collaboration , Abuter , R., Amorim , A., et al. 2020, , 636, L5
work page 2020
-
[63]
GRAVITY Collaboration , Abuter , R., Amorim , A., et al. 2021, , 647, A59
work page 2021
-
[64]
Gualandris , A., Mapelli , M., & Perets , H. B. 2012, , 427, 1793
work page 2012
- [65]
-
[66]
G \"u rkan , M. A., Freitag , M., & Rasio , F. A. 2004, , 604, 632
work page 2004
- [67]
-
[68]
Hailey, C. J., Mori, K., Bauer, F. E., et al. 2018, Nature, 556, 70
work page 2018
-
[69]
Hamers , A. S. & Perets , H. B. 2017, , 846, 123
work page 2017
-
[70]
Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357
work page 2020
-
[71]
Heggie , D. C. 1975, , 173, 729
work page 1975
-
[72]
Hills , J. G. 1983, , 88, 1269
work page 1983
-
[73]
Hills , J. G. 1988, , 331, 687
work page 1988
-
[74]
Hills , J. G. 1991, , 102, 704
work page 1991
- [75]
- [76]
-
[77]
Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90
work page 2007
- [78]
-
[79]
2016, in Positioning and Power in Academic Publishing: Players, Agents and Agendas, ed
Kluyver, T., Ragan-Kelley, B., P \'e rez, F., et al. 2016, in Positioning and Power in Academic Publishing: Players, Agents and Agendas, ed. F. Loizides & B. Schmidt, IOS Press, 87 -- 90
work page 2016
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.