pith. sign in

arxiv: 2412.02752 · v2 · submitted 2024-12-03 · 🌌 astro-ph.GA

The S stars' zone of avoidance in the Galactic center

Pith reviewed 2026-05-23 07:54 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords S starsGalactic centersupermassive black holebinary disruptionorbital relaxationzone of avoidanceS-star orbitseccentricity distribution
0
0 comments X

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.

S stars orbit the Milky Way's central black hole but avoid a specific wedge in pericenter-eccentricity space defined by log(r_p / AU) ≤ 1.57 + 2.6(1 - e). Models that start with binaries drawn from observed properties at 5-100 pc, let the black hole tear them apart, and then let the remnant stars undergo nonresonant and resonant relaxation recover both this empty zone and the observed thermal eccentricity distribution. A reader would care because the match supplies a concrete channel that turns distant binaries into the observed close S-star population and turns their orbital statistics into a record of the central parsec's dynamical history.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2412.02752 by Aleksey Generozov, Antonia Drescher, Diogo C. Ribeiro, Felix Mang, Frank Eisenhauer, Guillaume Bourdarot, Hagai B. Perets, Matteo S. Bordoni, Rainer Sch\"odel, Reinhard Genzel, Stefan Gillessen, Thomas Ott.

Figure 1
Figure 1. Figure 1: Flowchart summarizing the key steps in our model for the pro￾duction of S-stars, including references and software used. (b) The pericenter of the encounter is sampled uniformly be￾tween 16 gravitational radii (∼ the tidal radius of a single star) and three times the maximum tidal radius of the bi￾nary population (equation 2). Thus, we implicitly assume binaries are in the full loss cone regime, such that … view at source ↗
Figure 2
Figure 2. Figure 2: Initial (left) and final (right) pericenters and eccentricities for a subsample of model stars (Note the different x-scales in the two panels.) The left panel shows all stars, independently of when they are deposited in the Galactic Center. The right panel shows the 71 model stars surviving to the present day with K ≤ 18 and a ≤ 0.05 pc. Early (late) stars are shown in purple (red). Model stars initially h… view at source ↗
Figure 3
Figure 3. Figure 3: Top panel: Cumulative eccentricity distribution from the model (blue, solid), compared to the observed S-star eccentricity distribution (orange, dashed). For reference, we also show a thermal eccentricity distribution (green, dashed-dotted). Bottom panel: Cumulative pericen￾ter distribution from our model and observations. The numbers in the legend are the p-values from a two-sample KS-test with the model … view at source ↗
Figure 4
Figure 4. Figure 4: Cumulative eccentricity, pericenter, and semi-major axis distributions for bright (K ≤ 16, left panels) and faint (K > 16, right panels) early S-stars. reproducing the observed thermal eccentricity distribution would be more challenging in such scenarios. 5. There is also tension with the measured ages of the bright S￾stars, with the model predicting older ages. Again, this may Article number, page 9 of 15… view at source ↗
Figure 5
Figure 5. Figure 5: Top panel: K-band luminosity function of all model (blue) and observed (orange, dashed) S-stars. Bottom panel: K-band luminosity functions for model (blue) and observed (orange, dashed) early S-stars (with orbits). The red line is the K-band luminosity function of early stars in the central 0.5” from Schödel et al. (2020) (red line; see Ap￾pendix B for details). The luminosity functions are normalized to m… view at source ↗
Figure 7
Figure 7. Figure 7: Top panel: Corner plot showing distributions of semi-major axis (a), eccentricity (e), and K magnitude from theory (orange) and obser￾vations (blue). Bottom panel: Eccentricity versus semi-major axis for the observed S-stars, and for the model stars that are closest in semi-major axis. Article number, page 11 of 15 [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The light green line shows the luminosity function, accounting for the dependence of the binary disruption probability on stellar life￾time. (Effectively, the progenitor binaries are drawn from a present-day mass function rather than an initial mass function). The other lines are the same as in [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: K-band luminosity functions, after discarding stars formed more than 108 yr ago. The blue (light green) line corresponds to binaries drawn from an initial (present-day) mass function. 12 13 14 15 16 17 K 10 1 10 0 10 1 10 2 10 3 C D F × C o n s t Theory Observations Theory PDMF Schodel et al 2020 K-band LF rproj< 0.5'' Early Stars [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

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)
  1. [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.
  2. [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)
  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

2 responses · 0 unresolved

We thank the referee for their constructive comments on our manuscript. We address each major comment below.

read point-by-point responses
  1. 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

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

0 steps flagged

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

1 free parameters · 2 axioms · 0 invented entities

Based on abstract only: the central claim rests on assumed empirical binary distributions from 5-100 pc and on the validity of the relaxation approximations used; no invented entities are introduced.

free parameters (1)
  • initial binary population and mass function
    Abstract notes uncertainties in these inputs that affect the modeled S-star distributions.
axioms (2)
  • domain assumption Binaries originate from large scales (5-100 pc) with empirical property distributions
    Stated directly in the abstract as the source population for the disruption events.
  • domain assumption Nonresonant and resonant relaxation operate as modeled in the Galactic center potential
    Abstract invokes these processes to evolve the initially eccentric orbits.

pith-pipeline@v0.9.0 · 5812 in / 1431 out tokens · 33418 ms · 2026-05-23T07:54:23.309013+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

139 extracted references · 139 canonical work pages · 1 internal anchor

  1. [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. [2]

    write newline

    " 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. [3]

    & Perets , H

    Aharon , D. & Perets , H. B. 2015, , 799, 185

  4. [4]

    2024, arXiv e-prints, arXiv:2410.19881

    Akiba , T., Naoz , S., & Madigan , A.-M. 2024, arXiv e-prints, arXiv:2410.19881

  5. [5]

    Stellar Dynamics and Stellar Phenomena Near A Massive Black Hole

    Alexander , T. 2017, ArXiv e-prints [ [arXiv] 1701.04762 ]

  6. [6]

    & Hopman , C

    Alexander , T. & Hopman , C. 2009, , 697, 1861

  7. [7]

    & Morris , M

    Alexander , T. & Morris , M. 2003, , 590, L25

  8. [8]

    & Pfuhl , O

    Alexander , T. & Pfuhl , O. 2014, , 780, 148

  9. [9]

    2023, , 947, 8

    Amaro Seoane , P. 2023, , 947, 8

  10. [10]

    2010, , 713, 90

    Antonini , F., Faber , J., Gualandris , A., & Merritt , D. 2010, , 713, 90

  11. [11]

    & Merritt , D

    Antonini , F. & Merritt , D. 2013, , 763, L10

  12. [12]

    M., Lim , P

    Astropy Collaboration , Price-Whelan , A. M., Lim , P. L., et al. 2022, , 935, 167

  13. [13]

    P., Tollerud , E

    Astropy Collaboration , Robitaille , T. P., Tollerud , E. J., et al. 2013, , 558, A33

  14. [14]

    Bahcall , J. N. & Wolf , R. A. 1976, , 209, 214

  15. [15]

    2013, , 434, L26

    Balberg , S., Sari , R., & Loeb , A. 2013, , 434, L26

  16. [16]

    & Alexander , T

    Bar-Or , B. & Alexander , T. 2014, Classical and Quantum Gravity, 31, 244003

  17. [17]

    & Alexander , T

    Bar-Or , B. & Alexander , T. 2016, , 820, 129

  18. [18]

    & Fouvry , J.-B

    Bar-Or , B. & Fouvry , J.-B. 2018, , 860, L23

  19. [19]

    K., et al

    Bartko , H., Martins , F., Fritz , T. K., et al. 2009, , 697, 1741

  20. [20]

    2018, , 609, A28

    Baumgardt , H., Amaro-Seoane , P., & Sch \"o del , R. 2018, , 609, A28

  21. [21]

    Bonnell , I. A. & Rice , W. K. M. 2008, Science, 321, 1060

  22. [22]

    2017, , 469, 2042

    Bradnick , B., Mandel , I., & Levin , Y. 2017, , 469, 2042

  23. [23]

    2012, , 427, 127

    Bressan , A., Marigo , P., Girardi , L., et al. 2012, , 427, 127

  24. [24]

    Burkert , A., Gillessen , S., Lin , D. N. C., et al. 2023, arXiv e-prints, arXiv:2306.02076

  25. [25]

    Y., Liu , F

    Cao , C. Y., Liu , F. K., Li , S., Chen , X., & Wang , K. 2024, arXiv e-prints, arXiv:2411.09278

  26. [26]

    & Amaro-Seoane , P

    Chen , X. & Amaro-Seoane , P. 2014, , 786, L14

  27. [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

  28. [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

  29. [29]

    2016, , 823, 102

    Choi , J., Dotter , A., Conroy , C., et al. 2016, , 823, 102

  30. [30]

    G., Dayem, K

    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

  31. [31]

    E., Davies , M

    Dale , J. E., Davies , M. B., Church , R. P., & Freitag , M. 2009, , 393, 1016

  32. [32]

    B., Blackwell , R., Bailey , V

    Davies , M. B., Blackwell , R., Bailey , V. C., & Sigurdsson , S. 1998, , 301, 745

  33. [33]

    Dessart , L., Ryu , T., Amaro Seoane , P., & Taylor , A. M. 2024, , 682, A58

  34. [34]

    2019, Science, 365, 664

    Do , T., Hees , A., Ghez , A., et al. 2019, Science, 365, 664

  35. [35]

    2015, , 809, 143

    Do , T., Kerzendorf , W., Winsor , N., et al. 2015, , 809, 143

  36. [36]

    R., Ghez, A

    Do, T., Lu, J. R., Ghez, A. M., et al. 2013, , 764, 154

  37. [37]

    2016, , 222, 8

    Dotter , A. 2016, , 222, 8

  38. [38]

    Duncan , M. J. & Shapiro , S. L. 1983, , 268, 565

  39. [39]

    2017, , 464, 194

    Feldmeier-Krause , A., Kerzendorf , W., Neumayer , N., et al. 2017, , 464, 194

  40. [40]

    B., Pichon , C., & Magorrian , J

    Fouvry , J. B., Pichon , C., & Magorrian , J. 2017, , 598, A71

  41. [41]

    2017, , 467, 451

    Fragione , G., Capuzzo-Dolcetta , R., & Kroupa , P. 2017, , 467, 451

  42. [42]

    & Sari , R

    Fragione , G. & Sari , R. 2018, , 852, 51

  43. [43]

    M., Cheung , P., Portegies Zwart , S

    Fregeau , J. M., Cheung , P., Portegies Zwart , S. F., & Rasio , F. A. 2004, , 352, 1

  44. [44]

    2006 a , , 649, 91

    Freitag , M., Amaro-Seoane , P., & Kalogera , V. 2006 a , , 649, 91

  45. [45]

    A., & Rasio , F

    Freitag , M., G \"u rkan , M. A., & Rasio , F. A. 2006 b , , 368, 141

  46. [46]

    K., Gillessen , S., Dodds-Eden , K., et al

    Fritz , T. K., Gillessen , S., Dodds-Eden , K., et al. 2011, , 737, 73

  47. [47]

    2024, , 689, A190

    Gallego-Cano , E., Fritz , T., Sch \"o del , R., et al. 2024, , 689, A190

  48. [48]

    2018, , 609, A26

    Gallego-Cano , E., Sch \"o del , R., Dong , H., et al. 2018, , 609, A26

  49. [49]

    2021, , 501, 3088

    Generozov , A. 2021, , 501, 3088

  50. [50]

    & Madigan , A.-M

    Generozov , A. & Madigan , A.-M. 2020, , 896, 137

  51. [51]

    Generozov , A., Nayakshin , S., & Madigan , A. M. 2022, , 512, 4100

  52. [52]

    C., Metzger , B

    Generozov , A., Stone , N. C., Metzger , B. D., & Ostriker , J. P. 2018, , 478, 4030

  53. [53]

    1997, , 291, 219

    Genzel , R., Eckart , A., Ott , T., & Eisenhauer , F. 1997, , 291, 219

  54. [54]

    Genzel , R., Thatte , N., Krabbe , A., Kroker , H., & Tacconi-Garman , L. E. 1996, , 472, 153

  55. [55]

    M., Duch \^e ne , G., Matthews , K., et al

    Ghez , A. M., Duch \^e ne , G., Matthews , K., et al. 2003, , 586, L127

  56. [56]

    M., Klein , B

    Ghez , A. M., Klein , B. L., Morris , M., & Becklin , E. E. 1998, , 509, 678

  57. [57]

    M., Eisenhauer , F., et al

    Gillessen , S., Plewa , P. M., Eisenhauer , F., et al. 2017, , 837, 30

  58. [58]

    & Loeb , A

    Ginsburg , I. & Loeb , A. 2007, , 376, 492

  59. [59]

    & Quillen , A

    Gould , A. & Quillen , A. C. 2003, , 592, 935

  60. [60]

    2022, , 657, L12

    GRAVITY Collaboration , Abuter , R., Aimar , N., et al. 2022, , 657, L12

  61. [61]

    2018, , 615, L15

    GRAVITY Collaboration , Abuter , R., Amorim , A., et al. 2018, , 615, L15

  62. [62]

    2020, , 636, L5

    GRAVITY Collaboration , Abuter , R., Amorim , A., et al. 2020, , 636, L5

  63. [63]

    2021, , 647, A59

    GRAVITY Collaboration , Abuter , R., Amorim , A., et al. 2021, , 647, A59

  64. [64]

    Gualandris , A., Mapelli , M., & Perets , H. B. 2012, , 427, 1793

  65. [65]

    & Ramirez-Ruiz , E

    Guillochon , J. & Ramirez-Ruiz , E. 2013, , 767, 25

  66. [66]

    A., Freitag , M., & Rasio , F

    G \"u rkan , M. A., Freitag , M., & Rasio , F. A. 2004, , 604, 632

  67. [67]

    2017, , 847, 120

    Habibi , M., Gillessen , S., Martins , F., et al. 2017, , 847, 120

  68. [68]

    J., Mori, K., Bauer, F

    Hailey, C. J., Mori, K., Bauer, F. E., et al. 2018, Nature, 556, 70

  69. [69]

    Hamers , A. S. & Perets , H. B. 2017, , 846, 123

  70. [70]

    R., Millman, K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357

  71. [71]

    Heggie , D. C. 1975, , 173, 729

  72. [72]

    Hills , J. G. 1983, , 88, 1269

  73. [73]

    Hills , J. G. 1988, , 331, 687

  74. [74]

    Hills , J. G. 1991, , 102, 704

  75. [75]

    & Alexander , T

    Hopman , C. & Alexander , T. 2006 a , , 645, 1152

  76. [76]

    & Alexander , T

    Hopman , C. & Alexander , T. 2006 b , , 645, L133

  77. [77]

    Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90

  78. [78]

    S., Figer , D

    Kim , S. S., Figer , D. F., & Morris , M. 2004, , 607, L123

  79. [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

  80. [80]

    & Tremaine , S

    Kocsis , B. & Tremaine , S. 2011, , 412, 187

Showing first 80 references.