pith. sign in

arxiv: 1907.01555 · v1 · pith:VVOYBZRYnew · submitted 2019-07-02 · 🌌 astro-ph.GA

An order-disorder phase transition in black-hole star clusters

Jihad Touma , Scott Tremaine , Mher Kazandjian This is my paper

Pith reviewed 2026-05-25 10:50 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords black holesstar clustersphase transitiongalactic nucleitidal disruption eventsgravitational waveseccentric orbitslopsided equilibria
0
0 comments X

The pith

Star clusters around black holes can transition from spherical to lopsided equilibrium under generic evolution.

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

The paper establishes that generic evolutionary processes can drive a star cluster around a massive black hole from a spherical thermal equilibrium into a lopsided equilibrium. In the new state most stars occupy high-eccentricity orbits whose orientations are aligned with one another. The change matters because the spatial structure of the cluster sets the rates of observable transients such as tidal disruptions and gravitational-wave inspirals. Estimates that enforce spherical symmetry will therefore understate those rates if the ordered phase is reached. The authors note that improved models of cluster formation and evolution are needed to determine which phase actually occurs.

Core claim

In the course of generic evolutionary processes, a star cluster surrounding a black hole can undergo a robust phase transition from a spherical thermal equilibrium to a lopsided equilibrium, in which most stars are on high-eccentricity orbits with aligned orientations. The rate of transient events is expected to be much higher in the ordered phase.

What carries the argument

The order-disorder phase transition that converts a spherical thermal equilibrium into a lopsided equilibrium of aligned high-eccentricity orbits.

If this is right

  • Transient event rates rise substantially once the cluster enters the lopsided phase.
  • Rate calculations that assume spherical symmetry become inaccurate for clusters that reach the ordered state.
  • Whether a given cluster is spherical or lopsided depends on the details of its formation and subsequent evolution.

Where Pith is reading between the lines

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

  • The transition may produce detectable anisotropies in the stellar orbits observed near galactic centers.
  • The same mechanism could operate in other self-gravitating systems that possess a dominant central mass.
  • Direct N-body or Monte-Carlo simulations with varied initial conditions could map the range of parameters that trigger the transition.

Load-bearing premise

The evolutionary processes acting on the cluster are generic and sufficient to drive the transition without dependence on special initial conditions or fine-tuned parameters.

What would settle it

A census of stellar orbital eccentricities and apsidal orientations within the inner parsec of a galactic nucleus that shows whether the distribution is isotropic or strongly aligned and eccentric.

Figures

Figures reproduced from arXiv: 1907.01555 by Jihad Touma, Mher Kazandjian, Scott Tremaine.

Figure 1
Figure 1. Figure 1: FIG. 1. Lopsided phase transition in a maximum-entropy [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Secular evolution of a cluster of stars. The initial [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
read the original abstract

The centers of most galaxies contain massive black holes surrounded by dense star clusters. The structure of these clusters determines the rate and properties of observable transient events, such as flares from tidally disrupted stars and gravitational-wave signals from stars spiraling into the black hole. Most estimates of these rates enforce spherical symmetry on the cluster. Here we show that, in the course of generic evolutionary processes, a star cluster surrounding a black hole can undergo a robust phase transition from a spherical thermal equilibrium to a lopsided equilibrium, in which most stars are on high-eccentricity orbits with aligned orientations. The rate of transient events is expected to be much higher in the ordered phase. Better models of cluster formation and evolution are needed to determine whether clusters should be found in the ordered or disordered phase.

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

Summary. The manuscript claims that star clusters surrounding massive black holes undergo a robust order-disorder phase transition, driven by generic evolutionary processes, from a spherical thermal equilibrium to a lopsided equilibrium in which most stars occupy high-eccentricity orbits with aligned orientations; this would substantially increase rates of transient events such as tidal disruptions and extreme-mass-ratio inspirals.

Significance. If substantiated with explicit dynamical evolution, the result would challenge the spherical-symmetry assumption used in most rate calculations for galactic nuclei and could alter predicted event rates by large factors. The identification of a phase transition via statistical mechanics or orbit-averaged methods would constitute a notable conceptual advance if the derivation is parameter-free and the ordered state is shown to be reached from generic initial conditions.

major comments (2)
  1. [Abstract] Abstract: the central claim that the transition occurs 'in the course of generic evolutionary processes' from spherical thermal equilibrium is load-bearing for the paper's astrophysical relevance, yet the provided text supplies no derivation, simulation details, or quantitative evidence that standard two-body and resonant relaxation suffice to drive an ensemble of spherical initial conditions to the lopsided state.
  2. [Abstract] The manuscript appears to establish the existence of a lower-free-energy lopsided equilibrium (via mean-field or orbit-averaged statistical mechanics) but does not demonstrate consistent arrival at that state under the same operators starting from spherical configurations; this leaves the 'generic evolutionary processes' step unverified and undermines the assertion that the ordered phase is reached without special initial conditions or fine-tuned parameters.

Simulated Author's Rebuttal

2 responses · 0 unresolved

Thank you for the opportunity to respond to the referee report. We address the major comments below and propose revisions to clarify the manuscript's results on the phase transition in black-hole star clusters.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that the transition occurs 'in the course of generic evolutionary processes' from spherical thermal equilibrium is load-bearing for the paper's astrophysical relevance, yet the provided text supplies no derivation, simulation details, or quantitative evidence that standard two-body and resonant relaxation suffice to drive an ensemble of spherical initial conditions to the lopsided state.

    Authors: The full manuscript derives the phase transition using orbit-averaged statistical mechanics, showing that the lopsided state is the minimum free-energy configuration for a range of parameters. The generic evolutionary processes are the standard two-body relaxation and resonant relaxation, which are expected to drive the system toward equilibrium. While explicit simulations of the full evolution are beyond the scope of this work, we provide analytical arguments based on relaxation timescales. We will revise the abstract to better reflect the scope of the derivation and add quantitative estimates of the timescales involved. revision: yes

  2. Referee: [Abstract] The manuscript appears to establish the existence of a lower-free-energy lopsided equilibrium (via mean-field or orbit-averaged statistical mechanics) but does not demonstrate consistent arrival at that state under the same operators starting from spherical configurations; this leaves the 'generic evolutionary processes' step unverified and undermines the assertion that the ordered phase is reached without special initial conditions or fine-tuned parameters.

    Authors: We agree that the primary result is the identification of the ordered equilibrium as having lower free energy. The arrival at this state is argued to occur via the same relaxation operators that establish thermal equilibrium in the spherical case. To address the concern, we will include in the revised manuscript a more detailed discussion of the basin of attraction for the ordered state and why generic initial conditions lead to it, based on the mean-field theory. revision: yes

Circularity Check

0 steps flagged

No circularity: equilibria derived independently of evolutionary claim

full rationale

The paper derives the existence of a lopsided equilibrium via orbit-averaged statistical mechanics or mean-field methods applied to the Hamiltonian, showing it has lower free energy than the spherical state. This step is independent of the subsequent assertion that generic relaxation processes will drive the system across the transition; the latter is a statement about dynamics, not a redefinition or fit of the equilibrium itself. No self-citation is used to establish uniqueness of the ordered state, no parameters are fitted to data and then relabeled as predictions, and no ansatz is smuggled through prior work. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The claim rests on the existence and sufficiency of unspecified generic evolutionary processes; no free parameters, new entities, or additional axioms are stated in the abstract.

axioms (1)
  • domain assumption Star clusters around black holes evolve through generic processes that can drive phase transitions
    Invoked directly in the abstract as the driver of the transition.

pith-pipeline@v0.9.0 · 5661 in / 1066 out tokens · 22628 ms · 2026-05-25T10:50:15.846057+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

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

  1. [1]

    An order-disorder phase transition in black-hole star clusters

    M• ≃ 3× 108M⊙(σ/200 km s−1)4.4 over the range 106 ≲ M•/M ⊙ ≲ 1010. The orbits of stars are domi- nated by the gravitational field of the black hole within its sphere of influence, of radius rinfl = GM•/σ2 ≃ 20 pc (M•/108M⊙)0.55. For comparison, the event hori- zon of the black hole is at least five orders of magnitude smaller, r• = 2GM•/c2 = 9.6× 10−6 pc (M•/...

    work page internal anchor Pith review Pith/arXiv arXiv 1907

  2. [2]

    with a softening length b = 0.1a0. At each step in the MCMC simulation we chose a random orbit pair and changed the orbital elements of the pair randomly sub- ject to the constraint that the total angular momentum is conserved. We employed Creutz’s demon algorithm [14] to constrain E at a fixed value. Typically ∼ 106 steps were needed to reach an equilibri...

    work page 2048

  3. [3]

    Kormendy and L

    J. Kormendy and L. C. Ho, Coevolution (or not) of su- permassive black holes and host galaxies. Annu. Rev. As- tron. Astrophys. 51, 511–653 (2013)

    work page 2013

  4. [4]

    K. P. Rauch and S. Tremaine, Resonant relaxation in stellar systems. New Astron. 1, 149–170 (1996)

    work page 1996

  5. [5]

    Sridhar and J

    S. Sridhar and J. R. Touma, Stellar dynamics around a massive black hole - II. Resonant relaxation. Mon. Not. R. Astron. Soc. 458, 4143–4161 (2016)

    work page 2016

  6. [6]

    Bar-Or and J.-B

    B. Bar-Or and J.-B. Fouvry, Scalar resonant relaxation of stars around a massive black hole. Astrophys. J. Lett. 860, L23, 9 pp. (2018)

    work page 2018

  7. [7]

    Alexander, Stellar dynamics and stellar phenomena near a massive black hole

    T. Alexander, Stellar dynamics and stellar phenomena near a massive black hole. Annu. Rev. Astron. Astrophys. 55, 17–57 (2017)

    work page 2017

  8. [8]

    Kocsis and S

    B. Kocsis and S. Tremaine, Resonant relaxation and the warp of the stellar disc in the Galactic Center. Mon. Not. R. Astron. Soc. 412, 187–207 (2011)

    work page 2011

  9. [9]

    Genzel, F

    R. Genzel, F. Eisenhauer, and S. Gillessen, The Galactic Center massive black hole and nuclear star cluster. Rev. Mod. Phys. 82, 3121–3195 (2010)

    work page 2010

  10. [10]

    J. N. Bahcall and R. A. Wolf, Star distribution around a massive black hole in a globular cluster. Astrophys. J. 209, 214–232 (1976)

    work page 1976

  11. [11]

    Cohn and R

    H. Cohn and R. M. Kulsrud, The stellar distribution around a black hole – Numerical integration of the Fokker-Planck equation. Astrophys. J. 226, 1087–1108 (1978)

    work page 1978

  12. [12]

    This behavior is not possible on the resonant-relaxation timescale because semimajor axes are conserved

    The existence of well-defined maximum-entropy states is exceptional in self-gravitating systems, which can in- crease their entropy without limit by placing a small amount of mass on weakly bound or unbound orbits. This behavior is not possible on the resonant-relaxation timescale because semimajor axes are conserved

    work page

  13. [13]

    The following argument may provide some intuition for the physics of the phase transition. Lopsided states are possible because eccentric orbits in a Kepler potential do not precess, so the self-gravity of the stars can corral the apsidal lines of the stellar orbits into an aligned configu- ration. The orbit-averaged gravitational energy between two orbits...

    work page

  14. [14]

    J. R. Touma, S. Tremaine, and M. V. Kazandjian, Gauss’s method for secular dynamics, softened. Mon. Not. R. Astron. Soc. 394, 1085–1108 (2009)

    work page 2009

  15. [15]

    Touma, J., & Tremaine, S., The statistical mechanics of self-gravitating Keplerian disks. J. Phys. A 47, 292001 (2014)

    work page 2014

  16. [16]

    Creutz, Microcanonical Monte Carlo simulation

    M. Creutz, Microcanonical Monte Carlo simulation. Phys. Rev. Lett. 50, 1411–1414 (1983)

    work page 1983

  17. [17]

    Tremaine, Secular stability and instability in stellar systems surrounding massive objects

    S. Tremaine, Secular stability and instability in stellar systems surrounding massive objects. Astrophys. J. 625, 143–155 (2005)

    work page 2005

  18. [18]

    In such a case, lopsidedness may arise on the secular timescale, long before resonant relaxation leads to a maximum- entropy state

    Of course, a cluster could be formed with a non- monotonic phase-space distribution in angular momen- tum that is cold enough for dynamical instability. In such a case, lopsidedness may arise on the secular timescale, long before resonant relaxation leads to a maximum- entropy state

    work page

  19. [19]

    Kormendy and R

    J. Kormendy and R. Bender, The double nucleus and central black hole of M31. Astrophys. J. 522, 772–792 (1999)

    work page 1999

  20. [20]

    M. V. Kazandjian and J. R. Touma, The doubling of stellar black hole nuclei. Mon. Not. R. Astron. Soc. 430, 2732–2738 (2013)

    work page 2013

  21. [21]

    C. K. Brown and J. Magorrian, Three-dimensional Kep- lerian orbit-superposition models of the nucleus of M31. Mon. Not. R. Astron. Soc. 431, 80–91 (2013)

    work page 2013