pith. sign in

arxiv: 2407.11169 · v2 · submitted 2024-07-15 · ✦ hep-ph · astro-ph.GA

Supermassive black hole formation in the initial collapse of axion dark matter

Pith reviewed 2026-05-23 22:37 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.GA
keywords axion dark matterBose-Einstein condensatesupermassive black holesangular momentum transportcosmic dawngravitational collapseQCD axions
0
0 comments X

The pith

Axion dark matter forms supermassive black holes during the collapse of early overdensities by transporting angular momentum outward through rethermalization.

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

Axion dark matter thermalizes by gravitational self-interactions and forms a Bose-Einstein condensate. The paper shows that rethermalization of the axion fluid during the initial collapse of large scale overdensities near cosmic dawn moves angular momentum outward fast enough to allow black holes to form with masses from approximately 10^5 to a few times 10^10 solar masses. This holds for QCD axions and for axion-like particles heavier than 10^{-16} eV/c^2. A sympathetic reader would care because the mechanism supplies a direct route to the supermassive black holes observed at high redshift without intermediate stellar evolution steps.

Core claim

Axion dark matter thermalizes by gravitational self-interactions and forms a Bose-Einstein condensate. The rethermalization of the axion fluid during the initial collapse of large scale overdensities near cosmic dawn transports angular momentum outward sufficiently fast that black holes form with masses ranging from approximately 10^5 to a few times 10^{10} M_⊙. This conclusion holds for QCD axions and for axion-like particles of mass larger than 10^{-16} eV/c^2.

What carries the argument

Rethermalization of the axion Bose-Einstein condensate during gravitational collapse of overdensities, which carries angular momentum outward on the collapse timescale.

Load-bearing premise

Gravitational self-interactions cause axion dark matter to thermalize into a Bose-Einstein condensate whose rethermalization during collapse removes angular momentum on the timescale of the initial overdensity collapse.

What would settle it

A simulation or calculation showing that rethermalization fails to transport angular momentum outward fast enough to permit black hole formation within the collapse time of the overdensity would falsify the claim.

Figures

Figures reproduced from arXiv: 2407.11169 by Pierre Sikivie, Yuxin Zhao.

Figure 1
Figure 1. Figure 1: FIG. 1: Phase space distribution of cold collisionless part [PITH_FULL_IMAGE:figures/full_fig_p014_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Density profile at time [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Black hole mass as a function of [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
read the original abstract

Axion dark matter thermalizes by gravitational self-interactions and forms a Bose-Einstein condensate. We show that the rethermalization of the axion fluid during the initial collapse of large scale overdensities near cosmic dawn transports angular momentum outward sufficientlly fast that black holes form with masses ranging from approximately $10^5$ to a few times $10^{10}~M_\odot$. This conclusion holds for QCD axions and for axion-like particles of mass larger than $10^{-16}$ eV/$c^2$.

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 manuscript claims that axion dark matter thermalizes via gravitational self-interactions to form a Bose-Einstein condensate, and that rethermalization during the initial collapse of large-scale overdensities near cosmic dawn transports angular momentum outward sufficiently rapidly to produce black holes with masses from ~10^5 to a few ×10^{10} M_⊙. The result is asserted to apply to QCD axions and to axion-like particles with m_a ≳ 10^{-16} eV/c².

Significance. If the central claim is substantiated by explicit timescale comparisons, the work would supply a direct-collapse channel for supermassive black holes seeded by axion overdensities, offering a potential explanation for high-redshift quasars that does not rely on stellar-mass seeds or sustained super-Eddington accretion. It would also furnish a new observable link between axion mass and the high-z black-hole mass function.

major comments (2)
  1. [Abstract] Abstract: the central claim that rethermalization transports angular momentum outward on the collapse timescale is stated without any derivation, numerical estimate, or comparison of τ_rel to t_collapse; the manuscript must supply this comparison (including Hubble expansion and the linear-to-nonlinear transition) to establish that the inequality holds across the quoted mass window.
  2. [Derivation of relaxation rate] Section deriving the relaxation rate: the gravitational scattering rate for the condensate must be shown to satisfy τ_rel ≪ t_collapse at the relevant densities, velocities, and redshifts; without an explicit evaluation that accounts for the inhomogeneous density profile during collapse, the transport efficiency remains unverified and the mass range 10^5–10^{10} M_⊙ cannot be justified.
minor comments (1)
  1. [Abstract] Abstract, line 3: 'sufficientlly' is a typographical error and should read 'sufficiently'.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments, which help clarify the presentation of our results on axion dark matter collapse and black hole formation. We address each major comment below.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the central claim that rethermalization transports angular momentum outward on the collapse timescale is stated without any derivation, numerical estimate, or comparison of τ_rel to t_collapse; the manuscript must supply this comparison (including Hubble expansion and the linear-to-nonlinear transition) to establish that the inequality holds across the quoted mass window.

    Authors: We agree that the abstract and main text would benefit from an explicit, consolidated comparison of τ_rel and t_collapse. The manuscript derives the relaxation rate in Section 3 and estimates collapse times near cosmic dawn, but does not tabulate the ratio including Hubble expansion and the linear-to-nonlinear transition for the full mass window. In the revision we will add this comparison (with a new figure) to Section 3 and update the abstract to reference it, confirming the inequality for m_a ≳ 10^{-16} eV/c². revision: yes

  2. Referee: [Derivation of relaxation rate] Section deriving the relaxation rate: the gravitational scattering rate for the condensate must be shown to satisfy τ_rel ≪ t_collapse at the relevant densities, velocities, and redshifts; without an explicit evaluation that accounts for the inhomogeneous density profile during collapse, the transport efficiency remains unverified and the mass range 10^5–10^{10} M_⊙ cannot be justified.

    Authors: We concur that an evaluation using the inhomogeneous density profile is required to rigorously justify the mass range. The present derivation employs a mean-density approximation for the gravitational scattering rate. The revised manuscript will add an appendix containing explicit calculations of τ_rel at multiple points along the collapse trajectory (using the spherical-collapse density profile), evaluated at the relevant redshifts, velocities, and axion masses. These will demonstrate that τ_rel ≪ t_collapse holds across the quoted window. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected; derivation treated as self-contained.

full rationale

No load-bearing steps could be identified that reduce by the paper's own equations or self-citation to inputs by construction, as the full manuscript equations and any cited prior results on rethermalization timescales are not exhibited in the provided source. The abstract states the thermalization premise and the resulting black-hole mass range as a conclusion without showing a fitted parameter renamed as prediction or an ansatz smuggled via self-citation. This is the expected honest non-finding when no explicit reduction is quotable.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on two domain assumptions about axion behavior and one paper-specific assumption about the efficiency of rethermalization; no free parameters or new entities are explicitly introduced in the abstract.

axioms (2)
  • domain assumption Axion dark matter thermalizes by gravitational self-interactions and forms a Bose-Einstein condensate
    Opening sentence of the abstract; treated as established background.
  • ad hoc to paper Rethermalization during initial collapse transports angular momentum outward sufficiently fast for black-hole formation
    The load-bearing step that converts the BEC property into the stated black-hole mass range.

pith-pipeline@v0.9.0 · 5613 in / 1572 out tokens · 25448 ms · 2026-05-23T22:37:19.413194+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Axions explain the formation of supermassive black holes at cosmic dawn

    astro-ph.HE 2026-01 unverdicted novelty 2.0

    Supermassive black holes form naturally at cosmic dawn when dark matter is axions or axion-like particles with mass above 10^{-16} eV/c².

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages · cited by 1 Pith paper

  1. [1]

    Kormendy and L.C

    J. Kormendy and L.C. Ho, Ann. Rev. of Astron. and Astroph. , 51 (2013) 511

  2. [2]

    Kazunori et al

    A. Kazunori et al. (the EHT Collaboration), Ap. J. Lett. 8 75 (2019) L5, and Ap. J. Lett. 930 (2022) L12

  3. [3]

    Frank, A

    J. Frank, A. King and D.J. Raine, Accretion Power in Astrophysics, Third Edition, Cambridge University Press, 2002

  4. [4]

    Agazie et al

    G. Agazie et al. (the NANOGrav. Collaboration), Ap. J. Le tt. 951 (2023) L8; D.J. Reardon et al., Ap. J. Lett. 951 (2023) L6; H. Xu et al., Res. Astron. As troph. 23 (2023) 075024; J. Antoniadis et al. (the EPTA Collaboration), arXiv: 2306.16 227

  5. [5]

    Inayoshi, E

    K. Inayoshi, E. Visbal and Z. Haiman, Ann. Rev. of Astron. and Astroph., 58 (2020) 27

  6. [6]

    Balberg and S.L

    S. Balberg and S.L. Shapiro, Phys. Rev. Lett. 88 (2002) 10 1301; J. Pollack, D.N. Spergel and P. Steinhardt, Ap. J. 804 (2015) 2, 131; W.-X. Feng, H.-B. Yu a nd Y.-M. Zhong, Ap. J. Lett. 914 (2021) 2, L26

  7. [7]

    Rindler-Daller, K

    T. Rindler-Daller, K. Freese, M.H. Montgomery, D. Winge t and B. Paxton, Ap. J. 799 (2015) 210

  8. [8]

    Larson et al., Ap

    R. Larson et al., Ap. J. Lett. 953 (2023) L29; A. Bogdan et a l., Nature Atron. 8 (2024) 126; R. Maiolino et al., arXiv:2308.01230; L.J. Furtak et al., ar Xiv: 2308.05735; R. Miaolino et al., Nature 627 (2024) 59

  9. [9]

    Hooper, A

    D. Hooper, A. Ireland, G. Krujaic and A Stebbins, JCAP 04 ( 2024) 021

  10. [10]

    Sikivie and Q

    P. Sikivie and Q. Yang, Phys. Rev. Lett. 103 (2009) 11130 1

  11. [11]

    Erken, P

    O. Erken, P. Sikivie, H. Tam and Q. Yang, Phys. Rev. D 85 (2 012) 063520

  12. [12]

    Khlopov, B.A

    M.Y. Khlopov, B.A. Malomed and Y.B. Zel’dovich, MNRAS 2 15 (1985) 575

  13. [13]

    Preskill, F

    J. Preskill, F. Wilczek and M. Wise, Phys. Lett. B120 (19 83) 127; L. Abbott and P. Sikivie, Phys. Lett. B120 (1983) 133; M. Dine and W. Fischler, Phys. Le tt. B120 (1983) 137

  14. [14]

    Banik and P

    N. Banik and P. Sikivie, Phys. Rev. D 88 (2013) 123517

  15. [15]

    Bardeen, J.R

    J.M. Bardeen, J.R. Bond, N. Kaiser and A.S. Szalay, Ap. J . 304 (1986) 15

  16. [16]

    Peebles, Ap

    P.J.E. Peebles, Ap. J. 155 (1969) 393

  17. [17]

    Efstatathiou and B.J.T

    G. Efstatathiou and B.J.T. Jones, MNRAS 186 (1979) 133; J. Barnes and G. Efstathiou, Ap. J. 319 (1987) 575

  18. [18]

    Hernanadez, C

    X. Hernanadez, C. Park, B. Cervantes-Sodi and Y.-Y. Cho i, MNRAS 375 (2007) 163

  19. [19]

    C.C. Lin, L. Meistel and F.H. Shu, Ap. J. 142 (1965) 1431; Y.B. Zel’dovich, Astron. and Astroph. 5 (1970) 84; J. Binney, Ap. J. 215 (1977) 492

  20. [20]

    Gebhardt et al., Astron

    K. Gebhardt et al., Astron. J. 122 (2001) 2469

  21. [21]

    Ferrarese and D

    F. Ferrarese and D. Merritt, Ap. J. 538 (2000) L9; K. Gebh ardt et al., Ap. J. (2000) L13

  22. [22]

    Pathria and P.D

    R.K. Pathria and P.D. Beale, Statistical Mechanics, 3rd edition, Elsevier 2011, and references therein

  23. [23]

    Chakrabarty et al., Phys

    S. Chakrabarty et al., Phys. Rev. D 97 (2018) 043531

  24. [24]

    Weinberg, Gravitation and Cosmology , J

    S. Weinberg, Gravitation and Cosmology , J. Wiley and Sons, 1972. 13 x ˙x a) t ≪ tin b) t = tin c) t ≲ tcoll d) t = tcoll FIG. 1: Phase space distribution of cold collisionless part icles during the collpase of a large smooth overdensity near cosmic dawn, at four different times: a) just after the Big Bang, b) when the central part of the overdeisty is at ...

  25. [25]

    0 0 . 5 1 . 0 1 . 5 2 . 0 2 . 5 3 . 0 r

  26. [26]

    2: Density profile at time tin, when the central part of the overdensity is at turnaround, in units where Mf = 1 and R = 1

    25 ρ(r, tin) r|M =1 ¯ρ(tin) FIG. 2: Density profile at time tin, when the central part of the overdensity is at turnaround, in units where Mf = 1 and R = 1. In these units the central density ρ(⃗0, t in) = 3 / 4π and the contemporary average cosmological enegy density ¯ρ(tin) = 4 / 3π 3. The average cosmological enegy density is indicated by the horizontal...

  27. [27]

    0 0 . 2 0 . 4 0 . 6 0 . 8 1 . 0 j 100 102 104 106 108 1010 MBH/M⊙ a) Mf = 3 × 1010 M⊙ b) Mf = 1 × 109 M⊙ c) Mf = 3 × 107 M⊙ FIG. 3: Black hole mass as a function of j ≡ ω intin for three values of Mf . For given Mf there is only a very slight dependence of the black hole mass on tin. The values shown were computed for zcoll = 10 and hence tin = 240 Myr. 16