pith. sign in

arxiv: 1907.04330 · v1 · pith:LXRZWC3Bnew · submitted 2019-07-09 · 🌌 astro-ph.GA

The evolution of kicked stellar-mass black holes in star cluster environments II. Rotating star clusters

Jeremy J. Webb , Nathan W. C. Leigh , Roberto Serrano , Jillian Bellovary , K. E. Saavik Ford , Barry McKernan , Mario Spera , Alessandro A. Trani This is my paper

Pith reviewed 2026-05-25 00:02 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords black holesstar clustersdynamical frictionorbital decayrotating clustersN-body simulationstidal captureangular momentum
0
0 comments X

The pith

Kicked black holes in rotating star clusters have longer orbital decay times due to orbit circularization that halts dynamical friction.

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

This paper studies the trajectories of stellar-mass black holes that receive velocity kicks in rotating star clusters using N-body simulations and an analytic framework. It establishes that as a kicked black hole's orbit decays, it gains angular momentum from interactions with high-rotation-frequency stars, causing the orbit to circularize once the black hole's orbital frequency matches the local stars'. At circularization, dynamical friction becomes ineffective due to low relative velocities, resulting in longer orbital decay times compared to non-rotating clusters. The circularization timescale varies with the cluster's rotation rate and kick velocity, allowing decay into the core in slowly rotating clusters before circularization. This increases the likelihood of tidal capture events, potentially helping form binaries and high-mass black holes.

Core claim

For a black hole kicked outside the cluster's core in a rotating cluster, as its orbit decays the black hole will quickly gain angular momentum as it interacts with stars with high rotational frequencies. Once the black hole decays to the point where its orbital frequency equals that of local stars, its orbit will be circular and dynamical friction becomes ineffective since local stars will have low relative velocities. After circularization, the black hole's orbit decays on a longer timescale than if the host cluster was not rotating.

What carries the argument

Angular momentum gain from rotating stars leading to orbit circularization at frequency matching, which suppresses dynamical friction.

If this is right

  • BHs in rotating clusters have longer orbital decay times than in non-rotating ones.
  • The timescale for orbit circularization depends strongly on the cluster's rotation rate and the initial kick velocity.
  • Kicked BHs in slowly rotating clusters can decay into the core before circularization occurs.
  • The probability of a BH undergoing a tidal capture event increases in rotating clusters, possibly aiding in the formation of binaries and high-mass BHs.

Where Pith is reading between the lines

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

  • This mechanism may increase the retention of black holes in rotating globular clusters by delaying their sinking to the center where they might be ejected.
  • Observations of black hole distributions in clusters with measured rotation could test for this circularization effect.
  • Including gas or stellar evolution in future models might alter the angular momentum exchange and thus the decay times.

Load-bearing premise

The N-body simulations and analytic framework correctly capture the angular-momentum exchange between the kicked black hole and the rotating stellar population without significant numerical artifacts or missing physics such as stellar evolution or gas.

What would settle it

Performing identical N-body simulations of a kicked black hole in a rotating cluster and in an otherwise identical non-rotating cluster and finding that the decay times are the same would falsify the claim of longer decay times due to circularization.

Figures

Figures reproduced from arXiv: 1907.04330 by Alessandro A. Trani, Barry McKernan, Jeremy J. Webb, Jillian Bellovary, K. E. Saavik Ford, Mario Spera, Nathan W. C. Leigh, Roberto Serrano.

Figure 1
Figure 1. Figure 1: The initial mean angular frequency profile of star clus￾ters with q=0.0 (blue), q=0.3 (orange), 0.6 (green), and 1.0 (red) as a function of three dimensional clustercentric distance. time t = 0, we assume a Plummer model for the underlying gravitational potential with scale radius a. We further as￾sume the cluster rotates with rotational frequency Ω = Ωˆ ~ z, which is equal to the mean Ω of stars within a … view at source ↗
Figure 2
Figure 2. Figure 2: The evolution of each kicked BH through its host clus￾ter as functions of R (top panels) and z (bottom panels) for all three kick scenarios described in Section 3. In all three cases, the BH’s orbit eventually circularizes in the xy-plane due to the cluster rotating about the z-axis. rotation rates, BH masses, and kick velocities influence our findings. 4.1 The trajectory of the BH through its host cluster… view at source ↗
Figure 3
Figure 3. Figure 3: The orbital frequency of each kicked BH (blue) and of nearby stars (orange) over a period > 100 Myr for the non￾rotating cluster case (top panel) and all three kick scenarios de￾scribed in Section 3. In the rotating cases, the BH gains angular momentum until its orbit circularizes and it reaches the same orbital frequency as local stars. 4.2 Time evolution of each component of the kicked BH’s oscillatory a… view at source ↗
Figure 5
Figure 5. Figure 5: The same as [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: illustrates that for a range of BH masses, ro￾tation is still able to circularize the BH’s orbit and slow its decay. For the 50 M⊙ and 100 M⊙ BHs, the departure from the non-rotating cases is evident, with the decay of the BH’s amplitude stalling for several tens of Myr before slowly decaying into the cluster’s core. For the 10 M⊙ BH, while the BH’s orbit circularizes around 90 Myr it still reaches the cor… view at source ↗
Figure 7
Figure 7. Figure 7: The amplitudes for simulated kicked BHs with kick velocities ranging from 0.75 to 2 × σc. Solid lines mark clusters rotating with q=1. and dash lines mark non-rotating clusters with q=0. BHs will decay quickly to the clusters centre and do not have time to gain angular momentum from rotating stars. Hence their orbital decay is similar to the non-rotating case described in Webb et al. (2018). However for in… view at source ↗
Figure 8
Figure 8. Figure 8: The radial and tangential kinetic energies for simulated kicked BHs with kick velocities ranging from 0.75 to 2 × σc. Solid lines mark clusters rotating with q=1. and dash lines mark non-rotating clusters with q=0. orbital decay of kicked BHs [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The amplitudes for simulated kicked BHs in clusters rotating with q=0 (blue), 0.3 (orange), 0.6 (green) and 1.0 (red). decrease the rate at which the BH gains angular momentum Gualandris & Merritt (2008). While the decay of the BH in the q=0.3 model cluster is unaffected by the fact that the cluster is rotating, it is still gaining angular momentum as its orbit decays. There is simply an interplay between … view at source ↗
Figure 10
Figure 10. Figure 10: The radial (blue) and tangential (orange) kinetic en￾ergies over a period of 100 Myr for kicked BHs in the non-rotating case (left panel) as well as all three kick scenarios described in Section 3, namely initial kicks aligned with the x-axis (left center panel), the z-axis (right center panel) and initial kicks inclined at a 45◦ angle relative to the x- and z-axes (right panel). of simulations to GC-like… view at source ↗
read the original abstract

In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic framework to study their motion as a function of the host cluster and the kick itself. Our simulations indicate that for a BH that is kicked outside of the cluster's core, as its orbit decays in a rotating cluster the BH will quickly gain angular momentum as it interacts with stars with high rotational frequencies. Once the BH decays to the point where its orbital frequency equals that of local stars, its orbit will be circular and dynamical friction becomes ineffective since local stars will have low relative velocities. After circularization, the BH's orbit decays on a longer timescale than if the host cluster was not rotating. Hence BHs in rotating clusters will have longer orbital decay times. The timescale for orbit circularization depends strongly on the cluster's rotation rate and the initial kick velocity, with kicked BHs in slowly rotating clusters being able to decay into the core before circularization occurs. The implication of the circularization phase is that the probability of a BH undergoing a tidal capture event increases, possibly aiding in the formation of binaries and high-mass BHs.

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 / 2 minor

Summary. The paper continues a study of kicked stellar-mass black holes in star clusters, now focusing on rotating hosts. Using N-body simulations and an analytic framework, it claims that BHs kicked outside the core gain angular momentum through interactions with the rotating stellar population; once the BH orbital frequency matches the local stellar rotation, the orbit circularizes, dynamical friction becomes ineffective due to low relative velocities, and the subsequent decay timescale lengthens relative to non-rotating clusters. The circularization timescale depends on rotation rate and kick velocity; slowly rotating clusters allow decay into the core before circularization. This increases the probability of tidal capture events and binary formation.

Significance. If the central mechanism holds, the result implies that cluster rotation can substantially alter BH retention, orbital decay, and the likelihood of forming BH binaries or high-mass BHs via tidal captures. The combination of direct N-body trajectories with an analytic model that derives the frequency-matching condition from angular-momentum exchange (without additional free parameters) is a strength, as is the falsifiable prediction that decay times lengthen once circularization occurs.

major comments (2)
  1. [Results section (trajectories and decay times)] The central claim that decay times are longer in rotating clusters than in non-rotating ones is load-bearing, yet the manuscript provides no direct side-by-side comparison (e.g., identical initial conditions with rotation turned off) or quantitative ratio of decay timescales. Without this, the lengthening cannot be isolated from other simulation differences.
  2. [Analytic framework] The analytic framework asserts that dynamical friction becomes ineffective after frequency matching, but the text does not show an explicit derivation or test that the torque or friction coefficient drops to zero (or near-zero) at that point; a plot of relative velocity or friction force versus radius would be required to confirm the mechanism is not an artifact of the N-body integrator.
minor comments (2)
  1. [Abstract] The abstract states that 'the timescale for orbit circularization depends strongly on the cluster's rotation rate,' but no scaling relation or example numbers are given; a brief quantitative illustration would improve clarity.
  2. [Figure captions] Figure captions should list the exact number of particles, softening length, and rotation parameter values used in each run to allow reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive report. We address the two major comments below and will revise the manuscript to strengthen the presentation of our results.

read point-by-point responses

  1. Referee: [Results section (trajectories and decay times)] The central claim that decay times are longer in rotating clusters than in non-rotating ones is load-bearing, yet the manuscript provides no direct side-by-side comparison (e.g., identical initial conditions with rotation turned off) or quantitative ratio of decay timescales. Without this, the lengthening cannot be isolated from other simulation differences.

    Authors: We agree that a controlled side-by-side comparison with identical initial conditions is required to isolate the effect of rotation. The present manuscript compares to the non-rotating results of Paper I, but those runs do not share the exact same initial conditions. In the revision we will add new N-body simulations with rotation disabled using the same initial conditions as the rotating cases, together with a figure and table that report the decay-time ratios. revision: yes

  2. Referee: [Analytic framework] The analytic framework asserts that dynamical friction becomes ineffective after frequency matching, but the text does not show an explicit derivation or test that the torque or friction coefficient drops to zero (or near-zero) at that point; a plot of relative velocity or friction force versus radius would be required to confirm the mechanism is not an artifact of the N-body integrator.

    Authors: The analytic model derives the circularization condition from angular-momentum exchange, but we acknowledge that an explicit demonstration of the friction coefficient vanishing at frequency matching is not provided. In the revised manuscript we will expand the analytic section with the derivation of the torque and friction force, and we will add a figure showing relative velocity and estimated friction force versus radius along representative trajectories. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper derives its central claims about BH orbital circularization and extended decay times in rotating clusters directly from N-body simulation trajectories and an analytic model of angular-momentum exchange with the rotating stellar background. The frequency-matching condition and resulting dynamical-friction suppression follow from the modeled physical interactions rather than from any fitted parameter renamed as a prediction, self-citation chain, or definitional equivalence. No load-bearing step reduces by construction to the paper's own inputs; the results are self-contained against the external benchmark of the performed simulations.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review prevents exhaustive enumeration; the central claim rests on standard assumptions of Newtonian N-body dynamics and the validity of the Chandrasekhar dynamical-friction formula in a rotating background.

pith-pipeline@v0.9.0 · 5822 in / 1074 out tokens · 13682 ms · 2026-05-25T00:02:41.199196+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

22 extracted references · 22 canonical work pages · 2 internal anchors

  1. [1]

    2003, Gravitational N -body Simulations: Tools and Algorithms (Cambridge Monographs on Mathe- matical Physics)

    Aarseth, S.J. 2003, Gravitational N -body Simulations: Tools and Algorithms (Cambridge Monographs on Mathe- matical Physics). Cambridge University Press, Cambridge

    work page 2003

  2. [2]

    2016, ApJ, 833, 252 Antonini F., Capuzzo-Dolcetta R., Mastrobuono-Battisti A., Merritt D

    Vesperini, E. 2016, ApJ, 833, 252 Antonini F., Capuzzo-Dolcetta R., Mastrobuono-Battisti A., Merritt D. 2012, ApJ, 750, 111

    work page 2016

  3. [3]

    2016, MNRAS, 455, 25 Bahcall J

    Arca-Sedda, M. 2016, MNRAS, 455, 25 Bahcall J. N., Wolf R.A., 1976, ApJ, 209, 214 Bahcall J. N., Wolf R.A., 1977, ApJ, 216, 883

    work page 2016

  4. [4]

    L., Anderson, J., Piot to, G., van der Marel, R

    Bellini, A., Bianchini, P., Varri, A. L., Anderson, J., Piot to, G., van der Marel, R. P., Vesperini, E., Watkins, L. L. 2017, ApJ, 844, 167

    work page 2017

  5. [5]

    2016, MNRAS, 458, 3644

    Varri, A.L. 2016, MNRAS, 458, 3644

    work page 2016

  6. [6]

    The internal rotation of globular clusters revealed by Gaia DR2

    Bianchini, P., van der Marel, R. P., del Pino, A., Watkins, L. L., Bellini, A., Fardal, M. A., Libralato, M., Sills, A. 2018, MNRAS, Submitted, arXiv:1806.02580 Binney J. 1977, MNRAS, 181, 735 Binney J., Tremaine S., 1987, Galactic Dynamics (Prince- ton: Princeton University Press) Blecha L., Cox T. J., Loeb A., Hernquist L. 2011, MNRAS, 412, 2154 Chandra...

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  7. [7]

    1999, ApJ, 525, 720

    Colpi, M., Mayer, L., Governato, F. 1999, ApJ, 525, 720

    work page 1999

  8. [8]

    2007, MNRAS, 377, 465 Favata M., Hughes S

    Ernst, A., Glaschke, P., Fiestas, J., Just, A., Spurzem, R. 2007, MNRAS, 377, 465 Favata M., Hughes S. A., Holz D. E. 2004, ApJ, 607, L5 Feldmeier A., et al., 2014, A&A, 570, A2 Frank F., Rees F. M., 1976, MNRAS, 176, 633 Generozov A., Stone N. C., Metzger B. D., Ostriker J. P. 2018, MNRAS, 478, 4030 Gualandris A., Merritt D. 2008, ApJ, 678, 780 H´ enault...

    work page 2007

  9. [9]

    & Lagoute, C

    Longaretti, P.-Y. & Lagoute, C. 1996, A&A, 308, 453

    work page 1996

  10. [10]

    2013, A&A, 558, A117

    Lutzgendorf, N., Baumgardt, H., Kruijssen, J.M.D. 2013, A&A, 558, A117

    work page 2013

  11. [11]

    1960, MNRAS, 120, 204 Mackey A

    Lynden-Bell, D. 1960, MNRAS, 120, 204 Mackey A. D., Da Costa G. S., Ferguson A. M. N., Yong D., 2013, ApJ, 762, 65 Merritt D., Milosavljevic M., Favata M., Hughes S. A., Holz D. E. 2004, ApJL, 607, L9 Mezza A. 2002, A&A, 395, 25M Nguyen D. D., et al., 2018, ApJ, 858, 118 Peebles P. J. E., 1972, ApJ, 178, 371

    work page 1960

  12. [12]

    1992, MNRAS, 256, 368

    Pesce, E., Capuzzo-Dolcetta R., Vietri, M. 1992, MNRAS, 256, 368

    work page 1992

  13. [13]

    2016, MNRAS, 462, 2333 Repetto S., Davies, M.B., Sigurdsson, S

    Peuten, M., Zocchi, A., Gieles, M., Gualandris, A., H´ enault-Brunet, V. 2016, MNRAS, 462, 2333 Repetto S., Davies, M.B., Sigurdsson, S. 2012, MNRAS, 425, 2799

    work page 2016

  14. [14]

    L., Morscher, M., Pattabiraman, B., Chat- terjee, S., Haster, C.-J., and Rasio, F

    Rodriguez, C. L., Morscher, M., Pattabiraman, B., Chat- terjee, S., Haster, C.-J., and Rasio, F. A. 2015, Physical Review Letters, 115, 051101

    work page 2015

  15. [15]

    L., Chatterjee, S

    Rodriguez, C. L., Chatterjee, S. and Rasio, F. A. 2016, Physical Review Letters, 93, 084029

    work page 2016

  16. [16]

    Mapping the stability of stellar rotating spheres via linear response theory

    Rozier, S., Fouvry, J.B, Breen, P.G., Varri, A.L., Pi- chon, C., Heggie, D. C., 2019, MNRAS, Submitted, arXiv:1902.09299

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  17. [17]

    2015, MNRAS, 451, 4086

    Spera, M., Mapelli, M., Bressan, A. 2015, MNRAS, 451, 4086

    work page 2015

  18. [18]

    2016, MNRAS,460, 1

    Spera, M., Mapelli, M., Jeffries, R.D. 2016, MNRAS,460, 1

    work page 2016

  19. [19]

    & Mapelli, M

    Spera, M. & Mapelli, M. 2017, MNRAS,470.4739 Tiongco M. A., Vesperini E., Varri A. L., 2017, MNRAS, 469, 683 Tremaine S. D. & Weinberg, M.D. 1984, MNRAS, 209, 729

    work page 2017

  20. [20]

    2010, ApJ, 708, 1598 Trenti M., van der Marel, R

    Trenti, M., Vesperini, E., Pasquato, M. 2010, ApJ, 708, 1598 Trenti M., van der Marel, R. 2013, MNRAS, 435, 3272 Valtonen M., Karttunen H. 2006, The Three-Body Problem (Cambridge: Cambridge University Press) Vicari A., Capuzzo-Dolcetta R., Merritt, D. 2007, ApJ, 662 797

    work page 2010

  21. [21]

    & Vesperini, E

    Webb, J.J. & Vesperini, E. 2016, MNRAS, 463, 2383 Webb J. J., Leigh N. W. C., Singh A., Ford K. E. S., McK- ernan B., Bellovary J. 2018, MNRAS, 474, 3835

    work page 2016

  22. [22]

    Weinberg, M. D. 1986, ApJ, 300, 93 This paper has been typeset from a T EX/ LATEX file prepared by the author. c⃝ 2018 RAS, MNRAS 000, 1–13

    work page 1986