REVIEW 4 major objections 4 minor 236 references
The paper argues that a massive accretion disk torques binary stars into the Hills mechanism, producing bursts of S-stars and hypervelocity stars, and that the Galactic Center's young star disk drove about a hundred such disruptions.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-01 00:29 UTC pith:JQL5TOWB
load-bearing objection Well-validated new mechanism for disk-driven binary disruptions; the central ~10^2 Galactic Center estimate is soft because the evaporation model is a step function. the 4 major comments →
Binary disruptions driven by massive disks around massive black holes
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that accretion-disk-driven disruption (ADDD) is a generic, rich channel for binary destruction around massive black holes. For binaries whose outer orbits are nearly perpendicular to the disk, the disk torque drives their angular momentum down while conserving its component along the disk axis, pushing many of them into the loss cone on a timescale comparable to disk growth. The paper derives a semi-analytical formula for the fraction of binaries disrupted as a function of distance from the black hole, validates it numerically, and estimates that the disk that formed the Galactic Center's young stars disrupted on the order of 100 binaries. The byproducts—implanted S-star
What carries the argument
The key object is the 'librating island' of phase space, a set of orbits in the loss wedge whose argument of pericenter librates rather than circulates under the disk potential. Orbits in this island are the ones the disk torque can drive to the loss cone. The paper uses a fitting formula for the separatrix action of this island, which generalizes previous results to arbitrary disk and cluster density slopes, to write the disrupted-binary distribution as an integral over the loss wedge and island.
Load-bearing premise
The estimate that ~100 binaries were disrupted assumes that binaries survive unchanged until a single evaporation time and then vanish completely; the paper's own variant that lets hard binaries live longer increases the count dramatically, so the number is sensitive to this idealized binary depletion model.
What would settle it
Measuring the ejection-age distribution of hypervelocity stars: the model predicts a burst of ejections about 5 Myr ago, roughly isotropic, whereas steady-state mechanisms predict a flat age distribution. If a complete sample of HVSs shows no burst at that epoch, the ADDD burst estimate would be contradicted.
If this is right
- If a massive disk formed in the Galactic Center ~5 Myr ago, the model predicts a burst of ~100 disrupted binaries, producing an excess of S-stars and hypervelocity stars above the steady background.
- Hypervelocity stars from ADDD are ejected roughly isotropically and in temporal bursts, a combination unique among known mechanisms; this signature can be searched for in future surveys.
- In other galaxies, the S-star analogues implanted by ADDD can later be driven to tidal disruption, producing (possibly repeating) TDEs in the ~0.1–1 Gyr after an AGN phase, which may help explain TDE excesses in post-starburst galaxies.
- The disk torque also drives secular chaos in binaries that are not disrupted, potentially shrinking stellar binaries and merging compact-object binaries via gravitational-wave emission.
Where Pith is reading between the lines
- If ~100 binaries were disrupted in a burst 5 Myr ago, the present-day S-star cluster should contain a population of older, low-mass stars with a distinctive age gap; future faint-object observations could test this by detecting the low-mass tail of implanted stars.
- Because ADDD ejects stars nearly isotropically, surveys that look for a planar or axial signature might miss this channel; joint analysis of ejection times and velocities may be more discriminating than sky position alone.
- The paper's evaporation model may be the main lever on the rate: if hard binaries survive many evaporation times, the number of disruptions in old nuclear clusters could be an order of magnitude higher than the fiducial estimate, making ADDD competitive with two-body relaxation as a steady HVS source.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper extends the Kaur & Stone (2025) disk-driven disruption mechanism from single stars to binaries. It characterizes the outer-orbit conditions for binary disruptions (loss wedge and librating island), proposes a fitting formula for the separatrix action (Eq. 4), constructs a toy binary population (§4.3), and derives semi-analytical plus Monte Carlo estimates of the number of disruptions (Eqs. 7–8). These estimates are compared with N-body integrations in Fig. 6. The application to the Galactic Center's ~5-Myr-old clockwise disk yields a burst of ~10^2 binary disruptions, which the authors connect to S-cluster stars, hypervelocity stars, and possibly S5-HVS1. They also discuss implications for TDEs in post-starburst galaxies and argue that ADDD produces near-isotropic, bursty HVS ejections.
Significance. If the central estimate holds, ADDD would be a new, observationally relevant channel for Hills-mechanism disruptions, with a distinctive combination of bursty and near-isotropic ejection. The paper's strengths include a generic semi-analytical framework that extends KS25 to arbitrary disk/cluster density slopes, a careful Monte Carlo-vs-N-body comparison in Fig. 6, and falsifiable predictions for S-star magnitudes, HVS sky positions, and TDE delay timescales. The S5-HVS1 comparison is post-hoc but not used to fit parameters. The main weakness is that the headline number ~10^2 depends on a step-function binary-evaporation model whose sensitivity the authors themselves demonstrate in Fig. 7. The mechanism is well founded; the quantitative population-level prediction is not yet robust.
major comments (4)
- [§4.3 (items 5–6), Eq. (12), Fig. 7] The step-function evaporation model — binaries are kept unchanged until t_age = t_evap and then removed — is the single most load-bearing population assumption in the paper. Equations (7)–(8) inherit it through r_m in Eq. (10), and Eq. (12) sets t_evap. The authors themselves show in Fig. 7 (purple points) that replacing this by 'long-lived hard binaries' raises N_d by more than an order of magnitude for the Chen+23 SFH, and the text acknowledges that this 'improved physics can cause a drastic increase.' Since the N-body validation in §5 only checks which orbits reach r_t for a prescribed binary population, it does not constrain this survival assumption. The headline '~10^2' is therefore not pinned by the validated orbital dynamics. I request either a calibration of the survival model against the observed binary-fraction profile, or a reframing of the central claim as a conditional estim
- [§4.3, Eq. (10), Eq. (14)] The binary population model takes f_b,h=1 and λ=0 'throughout this work', yet the observational constraints are admittedly sparse. N_d is linear in f_b,h (see the d f_b,h/dain factor in Eq. (14)) and λ enters the radial exponent Λ in Eq. (B11); taking f_b,h~0.3–0.5 and λ>0 would reduce N_d by factors of 2–5 or more. The paper does not explore this in Fig. 7, so the plotted range understates the model uncertainty. The central estimate should be accompanied by a calculation of its sensitivity to these parameters.
- [§5, Fig. 6] The 'excellent agreement' between Monte Carlo and N-body validates the orbit-averaged loss-wedge/librating-island selection, not the full binary-disruption count. The simulations assign a_in and m_b from the population model and record whether the outer orbit reaches r_t; they do not evolve the binary population, evaporate binaries, or model binary-single interactions. Thus the validation is narrower than the abstract's 'validated with numerical simulations' implies. Please state explicitly that the population-level assumptions in §4.3 are not validated by Fig. 6.
- [§6.2, Fig. 7] The text states estimates range 'from a few to a few hundred' and 'our most detailed models suggest the number was 50–100,' while the abstract and conclusions state '~10^2.' Given the factor-of-several spread from β, μ, and IMF, plus the >10 effect of the long-lived-hard-binary variant, calling this 'order-of-magnitude-accurate' is not supported by the model's own sensitivity. I recommend a more cautious summary of the quantitative payload, and a discussion of which future observations (e.g., HVS ejection-time distribution, low-mass S-star population) could discriminate among the variants.
minor comments (4)
- [Eq. (6) and Fig. 4] The piecewise scaling for f_d(a) is stated without derivation in the main text; the Appendix B derivation is clear, but a one-sentence justification in §3.3.1 would help. Also, Fig. 4's dashed lines are arbitrarily normalized; please state the normalization in the caption.
- [§6.2 vs Fig. 7 caption] The text says 'our most detailed models suggest the number was 50–100' while the Fig. 7 caption says 'between ~10 and ~10^2'; make these consistent.
- [References] 'von Ziepel, H. V. 1909' should be 'von Zeipel, H. V. 1909' (standard spelling of the ZLK mechanism).
- [Appendix C] The 'factor 10' in p(shrink|enter) is introduced as the only fine-tuning; since it is set to match the Dodici et al. (2026) simulations, please report its sensitivity (e.g., how N_d changes if this factor is 3 or 30).
Circularity Check
No significant circularity: the ~10^2 estimate is a conditional prediction from an explicit population model, not a fitted or self-referential quantity.
full rationale
The central derivation is self-contained. The orbital-dynamics input—loss wedge, librating island, separatrix action—is taken from Kaur & Stone (2025) and generalized; the separatrix fit (eq. 4) and its J_sat, χ_sat parameters are fitted to numerically evaluated phase-space separatrices, not to the target disruption count. The population integrals (eqs. 7–8) convolve this dynamics with an explicit toy binary population, with the paper stating 'We take f_b,h = 1 and λ = 0 throughout' and 'Observational constraints for our model are sparse'—these are stated assumptions, not fitted values. The '~10^2' estimate is therefore a conditional prediction under those assumptions. The paper itself flags the key sensitivity ('This improved physics can cause a drastic increase in the number of disruptions from older populations') and lists 'Primary sources of uncertainty are the velocity anisotropy, the star-formation history, and the initial mass function slope.' The S5-HVS1 coincidence is post-hoc and is not used to calibrate any model parameter. Self-citations to Dodici et al. (2026) appear in Appendix C and Section 6.2 for binary shrinking, but the paper states this 'only leads to a small increase in the number of disruptions,' so the self-citation is not load-bearing for the central claim. No constructed equality between inputs and predictions was found.
Axiom & Free-Parameter Ledger
free parameters (8)
- mu (disk-to-MBH mass ratio) =
0.1 (fiducial); varied 0.05-1
- gamma_d (disk density slope) =
5/2 (fiducial)
- beta (velocity anisotropy) =
0 (fiducial)
- J_sat, chi_sat (separatrix-action fit parameters) =
Table 2 values per gamma_d (e.g., 0.685±0.005, 3.58±0.14 for gamma_d=2.5)
- alpha_m (IMF slope) =
1.7 (fiducial); varied 1.5, 2.3
- f_b,h (binary fraction at r_h) =
1
- lambda (binary formation efficiency power) =
0
- factor 10 in p(shrink|enter) =
10
axioms (6)
- domain assumption KS25 axisymmetric potential-density pairs represent AGN disk torques
- domain assumption Spherical power-law cluster potential (gamma_c)
- domain assumption Secular orbit-averaged treatment for analytical estimates
- domain assumption Binary disrupts when outer orbit reaches tidal radius r_t
- domain assumption Milky Way's young stellar disk formed in situ from a massive gas disk ~5 Myr ago
- domain assumption NSC star formation histories (three models)
read the original abstract
Accretion disks around massive black holes (MBHs) in galactic nuclei can contain significant mass, in which case their non-spherical potentials exert significant torques on single and binary stars orbiting the MBHs. These torques can drive orbits to very large eccentricities. Previous works have shown that this driving can cause disruptions of single stars by the tidal gravity of the MBH. Here, we characterize the ability of these torques to drive binary stars to the point of disruption via the Hills mechanism. We derive semi-analytical estimates, validated with numerical simulations, of the number of binary disruptions driven by a generic disk in a generic nucleus. Using these results, we estimate that the formation of the ~5-Myr-old disk of stars in the Galactic Center drove a burst of ~10^2 binary disruptions. These disruptions produced an excess of S-cluster stars and hypervelocity stars --- possibly including S5-HVS1, the fastest-known hypervelocity star. In other galaxies, analogues of S-cluster stars produced by this process may evolve into tidal disruption events following a disk phase. Hypervelocity star observations may help clarify the importance of disk-driven disruptions, which eject stars nearly isotropically and in temporal bursts --- a unique combination among processes driving binary disruptions.
Figures
Reference graph
Works this paper leans on
-
[1]
The Loss-cone Problem in Axisymmetric Nuclei. , keywords =. doi:10.1088/0004-637X/774/1/87 , archivePrefix =. 1301.3150 , primaryClass =
-
[2]
Disruptions of stars and binary systems on chaotic orbits in an axisymmetric Milky Way centre. , keywords =. doi:10.1093/mnras/staf1237 , archivePrefix =. 2505.06344 , primaryClass =
-
[3]
doi:10.1093/mnras/stw1220 , eprint =
, keywords =. doi:10.1093/mnras/stw1220 , eprint =
-
[4]
arXiv , author =:1903.00010 , journal =
doi:10.3847/1538-4357/ab1e4d , eid =. arXiv , author =:1903.00010 , journal =
Pith/arXiv arXiv 1903
-
[5]
arXiv , author =:2312.17319 , journal =
doi:10.3847/1538-4357/ad61e1 , eid =. arXiv , author =:2312.17319 , journal =
-
[6]
arXiv , author =:1910.02079 , journal =
doi:10.3847/1538-4357/ab8461 , eid =. arXiv , author =:1910.02079 , journal =
Pith/arXiv arXiv 1910
-
[7]
arXiv , author =:2303.16231 , journal =
doi:10.3847/1538-4357/acf65b , eid =. arXiv , author =:2303.16231 , journal =
-
[8]
arXiv , author =:1509.03633 , journal =
doi:10.1051/0004-6361/201526757 , eid =. arXiv , author =:1509.03633 , journal =
-
[9]
, keywords =. doi:10.1086/177864 , eprint =
-
[10]
doi:10.1038/s41586-019-1556-x , eprint =
, keywords =. doi:10.1038/s41586-019-1556-x , eprint =
-
[11]
arXiv , author =:2002.08967 , journal =
doi:10.1051/0004-6361/202037610 , eid =. arXiv , author =:2002.08967 , journal =
Pith/arXiv arXiv 2002
-
[12]
doi:10.1038/s41586-021-03394-6 , eprint =
, keywords =. doi:10.1038/s41586-021-03394-6 , eprint =
-
[13]
arXiv , author =:2110.10786 , journal =
doi:10.3847/2041-8213/ac313b , eid =. arXiv , author =:2110.10786 , journal =
Pith/arXiv arXiv 2041
-
[14]
arXiv , author =:2207.07511 , journal =
doi:10.1051/0004-6361/202244512 , eid =. arXiv , author =:2207.07511 , journal =
-
[15]
doi:10.1038/s41550-023-02073-y , eprint =
Nature Astronomy , keywords =. doi:10.1038/s41550-023-02073-y , eprint =
-
[16]
arXiv , author =:2009.03321 , journal =
doi:10.3847/1538-4357/abe38d , eid =. arXiv , author =:2009.03321 , journal =
Pith/arXiv arXiv 2009
-
[17]
arXiv , author =:2208.12452 , journal =
doi:10.1051/0004-6361/202244805 , eid =. arXiv , author =:2208.12452 , journal =
-
[18]
arXiv , author =:2209.07538 , journal =
doi:10.3847/2041-8213/ac9f36 , eid =. arXiv , author =:2209.07538 , journal =
Pith/arXiv arXiv 2041
-
[19]
doi:10.1093/mnrasl/slaa020 , eprint =
, keywords =. doi:10.1093/mnrasl/slaa020 , eprint =
-
[20]
doi:10.1093/mnrasl/slad113 , eprint =
, keywords =. doi:10.1093/mnrasl/slad113 , eprint =
-
[21]
arXiv , author =:2001.03504 , journal =
doi:10.3847/1538-4357/abb3cc , eid =. arXiv , author =:2001.03504 , journal =
Pith/arXiv arXiv 2001
-
[22]
doi:10.1093/mnras/176.3.633 , journal =
-
[23]
doi:10.1038/331687a0 , journal =
-
[24]
arXiv , author =:2203.08162 , journal =
doi:10.3847/2041-8213/ac6021 , eid =. arXiv , author =:2203.08162 , journal =
Pith/arXiv arXiv 2041
-
[25]
doi:10.48550/arXiv.2401.12555 , eid =. arXiv , author =:2401.12555 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2401.12555
-
[26]
arXiv , author =:2303.16977 , journal =
doi:10.3847/1538-4357/acc93e , eid =. arXiv , author =:2303.16977 , journal =
-
[27]
doi:10.1086/172607 , journal =
-
[28]
doi:10.1093/mnras/173.3.729 , journal =
-
[29]
doi:10.48550/arXiv.2311.04352 , eid =. arXiv , author =:2311.04352 , journal =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2311.04352
- [30]
-
[31]
doi:10.1093/mnras/stad3470 , eprint =
, keywords =. doi:10.1093/mnras/stad3470 , eprint =
-
[32]
doi:10.1088/0004-637X/700/2/1933 , eprint =
, keywords =. doi:10.1088/0004-637X/700/2/1933 , eprint =
-
[33]
Dynamics of recaptures, ejections, and mergers of stellar mass binaries over multiple encounters with Sgr A*. , keywords =. doi:10.1093/mnras/staf1766 , archivePrefix =. 2505.08499 , primaryClass =
-
[34]
The Two Young Star Disks in the Central Parsec of the Galaxy: Properties, Dynamics, and Formation. , keywords =. doi:10.1086/503273 , archivePrefix =. astro-ph/0601268 , primaryClass =
-
[35]
The Young Stars in the Galactic Center. , keywords =. doi:10.3847/2041-8213/ac68ef , archivePrefix =. 2205.07595 , primaryClass =
Pith/arXiv arXiv 2041
-
[38]
Enhanced activity of massive black holes by stellar capture assisted by a self-gravitating accretion disc. , keywords =. doi:10.1051/0004-6361:20066068 , archivePrefix =. 0704.2781 , primaryClass =
-
[39]
arXiv , author =:2101.05825 , journal =
doi:10.3847/1538-4357/abdc25 , eid =. arXiv , author =:2101.05825 , journal =
-
[40]
doi:10.1093/mnras/stab609 , eprint =
, keywords =. doi:10.1093/mnras/stab609 , eprint =
-
[41]
doi:10.1093/mnras/stad2203 , eprint =
, keywords =. doi:10.1093/mnras/stad2203 , eprint =
-
[42]
arXiv , author =:2210.09945 , journal =
doi:10.3847/1538-4357/acbfb6 , eid =. arXiv , author =:2210.09945 , journal =
-
[43]
arXiv , author =:2304.00775 , journal =
doi:10.1051/0004-6361/202346565 , eid =. arXiv , author =:2304.00775 , journal =
-
[44]
arXiv , author =:2211.09851 , journal =
doi:10.3847/1538-4357/acbd3d , eid =. arXiv , author =:2211.09851 , journal =
-
[45]
arXiv , author =:2102.08135 , journal =
doi:10.3847/1538-4357/ac05c6 , eid =. arXiv , author =:2102.08135 , journal =
-
[46]
arXiv , author =:2110.10855 , journal =
doi:10.3847/2041-8213/ac31aa , eid =. arXiv , author =:2110.10855 , journal =
Pith/arXiv arXiv 2041
-
[47]
arXiv , author =:2404.08138 , journal =
doi:10.3847/1538-4357/ad5cf2 , eid =. arXiv , author =:2404.08138 , journal =
-
[48]
arXiv , author =:2001.03626 , journal =
doi:10.1007/s00159-020-00125-0 , eid =. arXiv , author =:2001.03626 , journal =
Pith/arXiv arXiv 2001
-
[49]
doi:10.1103/RevModPhys.82.3121 , eprint =
Reviews of Modern Physics , keywords =. doi:10.1103/RevModPhys.82.3121 , eprint =
-
[50]
arXiv , author =:1403.6657 , journal =
doi:10.1051/0004-6361/201423481 , eid =. arXiv , author =:1403.6657 , journal =
-
[51]
arXiv , author =:2307.11821 , journal =
doi:10.1051/0004-6361/202347416 , eid =. arXiv , author =:2307.11821 , journal =
-
[52]
doi:10.1146/annurev-astro-091916-055306 , eprint =
, keywords =. doi:10.1146/annurev-astro-091916-055306 , eprint =
-
[53]
doi:10.1103/RevModPhys.15.1 , journal =
-
[54]
doi:10.1016/S1384-1076(96)00012-7 , eprint =
, keywords =. doi:10.1016/S1384-1076(96)00012-7 , eprint =
-
[56]
arXiv , author =:1805.10313 , journal =
doi:10.3847/1538-4357/aadae2 , eid =. arXiv , author =:1805.10313 , journal =
-
[57]
arXiv , author =:1802.08890 , journal =
doi:10.3847/2041-8213/aac88e , eid =. arXiv , author =:1802.08890 , journal =
Pith/arXiv arXiv 2041
-
[58]
doi:10.1086/118405 , journal =
-
[59]
doi:10.1093/mnras/stz1730 , eprint =
, keywords =. doi:10.1093/mnras/stz1730 , eprint =
-
[60]
doi:10.1093/mnras/stz2026 , eprint =
, keywords =. doi:10.1093/mnras/stz2026 , eprint =
-
[61]
arXiv , author =:1812.07053 , journal =
doi:10.3847/1538-4357/ab2f78 , eid =. arXiv , author =:1812.07053 , journal =
-
[62]
doi:10.1093/mnras/stv057 , eprint =
, keywords =. doi:10.1093/mnras/stv057 , eprint =
-
[63]
arXiv , author =:2306.03703 , journal =
doi:10.3847/1538-4357/ad0be2 , eid =. arXiv , author =:2306.03703 , journal =
-
[64]
doi:10.1088/0004-6256/136/6/2552 , eprint =
, keywords =. doi:10.1088/0004-6256/136/6/2552 , eprint =
-
[65]
, month = jul, pages =
-
[66]
doi:10.1093/mnras/staa2720 , eprint =
, keywords =. doi:10.1093/mnras/staa2720 , eprint =
-
[67]
arXiv , author =:1710.02542 , journal =
doi:10.3847/1538-3881/aaa970 , eid =. arXiv , author =:1710.02542 , journal =
-
[68]
arXiv , author =:2203.01947 , journal =
doi:10.3847/1538-4357/ac8aff , eid =. arXiv , author =:2203.01947 , journal =
-
[69]
arXiv , author =:1110.6655 , journal =
doi:10.1088/0004-637X/747/1/4 , eid =. arXiv , author =:1110.6655 , journal =
-
[70]
arXiv , author =:1405.6029 , journal =
doi:10.1088/0004-637X/799/2/118 , eid =. arXiv , author =:1405.6029 , journal =
- [71]
-
[72]
doi:10.1146/annurev-astro-081915-023315 , eprint =
, keywords =. doi:10.1146/annurev-astro-081915-023315 , eprint =
-
[73]
doi:doi/10.1002/asna.19091832202 , journal =
-
[74]
doi:10.1016/0032-0633(62)90129-0 , journal =
-
[75]
doi:10.1086/108790 , journal =
-
[76]
doi:10.1086/111815 , journal =
-
[77]
arXiv , author =:2002.10547 , journal =
doi:10.3847/1538-4357/ab94bc , eid =. arXiv , author =:2002.10547 , journal =
Pith/arXiv arXiv 2002
-
[78]
A Stream of Hypervelocity Stars from the Galactic Center. , keywords =. doi:10.3847/1538-4357/abbf4f , archivePrefix =. 2005.10267 , primaryClass =
Pith/arXiv arXiv 2005
-
[79]
The distribution of stars around the Milky Way's central black hole. II. Diffuse light from sub-giants and dwarfs. , keywords =. doi:10.1051/0004-6361/201730452 , archivePrefix =. 1701.03817 , primaryClass =
-
[80]
The distribution of stars around the Milky Way's central black hole. I. Deep star counts. , keywords =. doi:10.1051/0004-6361/201730451 , archivePrefix =. 1701.03816 , primaryClass =
-
[81]
The old nuclear star cluster in the Milky Way: dynamics, mass, statistical parallax, and black hole mass. , keywords =. doi:10.1093/mnras/stu2452 , archivePrefix =. 1403.5266 , primaryClass =
-
[82]
The nuclear star cluster of the Milky Way: proper motions and mass. , keywords =. doi:10.1051/0004-6361/200810922 , archivePrefix =. 0902.3892 , primaryClass =
-
[83]
Three-dimensional Stellar Kinematics at the Galactic Center: Measuring the Nuclear Star Cluster Spatial Density Profile, Black Hole Mass, and Distance. , keywords =. doi:10.1088/2041-8205/779/1/L6 , archivePrefix =. 1311.0886 , primaryClass =
Pith/arXiv arXiv 2041
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.