pith. sign in

arxiv: 2508.14535 · v2 · submitted 2025-08-20 · ✦ hep-ph · astro-ph.HE

Axion Star Bosenova in Axion Miniclusters

Zihang Wang , Yu Gao This is my paper

Pith reviewed 2026-05-18 22:41 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.HE
keywords axion dark matterminiclustersaxion starsbosenovaself-interactionQCD axionaxion-like particles
0
0 comments X p. Extension

The pith

For the QCD axion, bosenova occurs within the age of the Universe for miniclusters with initial overdensity δ≳100.

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

The paper studies axion stars that form at the centers of axion miniclusters and grow by pulling in axions from the surrounding structure. It incorporates self-interaction effects to find when the star mass exceeds a stability limit and triggers a bosenova. A sympathetic reader would care because the process could happen often in dense enough miniclusters and might leave traces in the distribution or signals from axionic dark matter. The analysis shows self-interaction can control the growth rate, and bosenova becomes common across much of the allowed axion parameter space once the initial overdensity is high enough.

Core claim

Axionic dark matter forms miniclusters that host central axion stars. The star accretes axions from the host until its mass exceeds the stability threshold set by self-interactions, at which point a bosenova occurs. For the QCD axion this happens inside the age of the Universe whenever the initial overdensity satisfies δ ≳ 100. Self-interaction dominates the mass growth for many cases, and bosenova takes place in a large fraction of axion parameter space when the initial overdensity is large.

What carries the argument

the self-interaction-limited stability threshold of dilute axion stars, which fixes the mass at which continuous accretion from the minicluster produces a bosenova

If this is right

  • Bosenova occurs for the QCD axion in miniclusters with initial overdensity δ ≳ 100 within the age of the Universe.
  • Self-interaction can dominate the axion star mass growth.
  • Bosenova occurs in a large fraction of axion parameter space for miniclusters with large initial overdensities.
  • Explicit conditions for bosenova are given for both the QCD axion and temperature-independent axion-like particles.

Where Pith is reading between the lines

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

  • This growth and explosion cycle could reduce the final mass of axion stars that survive in miniclusters.
  • Repeated bosenova events might alter the internal density profiles of axion miniclusters over time.
  • The same accretion-until-threshold logic could apply to other self-interacting scalar dark-matter candidates.

Load-bearing premise

The axion star mass grows continuously by accreting axions from the host minicluster until it exceeds the self-interaction-limited stability threshold, without significant dynamical disruption or mass loss prior to the bosenova.

What would settle it

Finding that axion stars in high-overdensity miniclusters lose mass or are disrupted by other processes before reaching the self-interaction stability limit, or observing no bosenova signatures where the model predicts them.

Figures

Figures reproduced from arXiv: 2508.14535 by Yu Gao, Zihang Wang.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison of the MC mass [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Parameter space for the ALP. The blue line divides [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
read the original abstract

Axionic dark matter can form structures known as miniclusters that host an axion star at their center. The axion star feeds on the host, and the axion star mass may grow beyond its stability limit, leading to a potential bosenova. Since a dilute axion star has a stable mass limit only when self-interaction is considered, we include axion self-interaction effects in this paper, and specify the condition for bosenova in the QCD axion and temperature-independent axion-like particle parameter spaces. We find that self-interaction may dominate the mass growth of the axion star. For a minicluster with a large initial overdensity, bosenova occurs in a large fraction of axion parameter space. For the QCD axion, bosenova occurs within the age of the Universe for miniclusters with an initial overdensity $\delta\gtrsim 100$.

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 models the growth of axion stars embedded in axion miniclusters, incorporating attractive self-interactions to determine the stability limit. It calculates the accretion-driven mass increase from the host minicluster and identifies the conditions under which the star exceeds this limit, triggering a bosenova. For the QCD axion, the authors conclude that bosenova occurs within the age of the Universe for miniclusters with initial overdensity δ ≳ 100, and that self-interactions can dominate the growth in large regions of parameter space for both QCD axions and temperature-independent ALPs.

Significance. If the uninterrupted accretion model holds, the result would imply that bosenovae are a generic outcome for sufficiently overdense miniclusters, providing a potential observational signature or constraint on axion dark matter models. The inclusion of self-interaction effects and explicit parameter-space mapping strengthens the work relative to purely gravitational treatments.

major comments (2)
  1. [§3 (growth modeling) and abstract] The central claim (abstract and §4) that bosenova occurs for δ ≳ 100 rests on the assumption that the axion star accretes mass continuously at a rate set by the minicluster density until it reaches the self-interaction-limited M_max without net mass loss. The manuscript provides no quantitative estimate or simulation check of competing processes such as gravitational scattering off density fluctuations, tidal stripping, or resonant ejection within the inhomogeneous minicluster environment; these could shorten the effective growth time below the Hubble time.
  2. [§2 and Eq. (stability limit)] The stability threshold is taken from prior axion-star literature (cited in §2), but the overlap of the present minicluster environment with the isolated-star assumptions underlying that threshold is not re-derived or tested. If the dense host alters the effective self-interaction or binding, the quoted M_max would shift, directly affecting the δ ≳ 100 threshold.
minor comments (2)
  1. [Introduction] Notation for the overdensity δ and the self-interaction coupling λ should be defined explicitly at first use rather than relying on standard conventions.
  2. [Figure 3] Figure captions for the parameter-space plots should state the precise numerical criterion used to delineate the bosenova region (e.g., growth time < t_universe).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment below and have made revisions to the manuscript where necessary to clarify our assumptions and strengthen the presentation.

read point-by-point responses

  1. Referee: [§3 (growth modeling) and abstract] The central claim (abstract and §4) that bosenova occurs for δ ≳ 100 rests on the assumption that the axion star accretes mass continuously at a rate set by the minicluster density until it reaches the self-interaction-limited M_max without net mass loss. The manuscript provides no quantitative estimate or simulation check of competing processes such as gravitational scattering off density fluctuations, tidal stripping, or resonant ejection within the inhomogeneous minicluster environment; these could shorten the effective growth time below the Hubble time.

    Authors: We agree that competing dynamical processes could affect the net accretion rate. Our calculation provides the growth timescale under the assumption of uninterrupted accretion from the minicluster density profile. In highly overdense miniclusters (δ ≳ 100), the central regions are dense enough that the accretion time is significantly shorter than the age of the Universe, allowing for bosenova even if some mass loss occurs due to scattering or stripping. We have expanded the discussion in §3 to explicitly state this assumption and its limitations, emphasizing that this work focuses on the impact of self-interactions rather than a full dynamical simulation of the minicluster. A quantitative evaluation of loss mechanisms would require high-resolution N-body simulations, which are beyond the scope of the present analytic model but represent an important direction for future work. revision: partial

  2. Referee: [§2 and Eq. (stability limit)] The stability threshold is taken from prior axion-star literature (cited in §2), but the overlap of the present minicluster environment with the isolated-star assumptions underlying that threshold is not re-derived or tested. If the dense host alters the effective self-interaction or binding, the quoted M_max would shift, directly affecting the δ ≳ 100 threshold.

    Authors: The stability criterion for the axion star is based on the point where the attractive self-interaction overcomes the repulsive quantum pressure and gravity within the star. Because the axion star is self-gravitating and its central density greatly exceeds that of the surrounding minicluster (by many orders of magnitude), the internal dynamics and stability are determined locally by the star's own properties. The minicluster acts primarily as a mass reservoir for accretion without substantially modifying the star's binding or self-interaction effects. We have added a brief justification in §2 explaining the applicability of the isolated-star stability limit to the embedded case. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation from dynamical accretion modeling

full rationale

The paper models axion star growth via continuous accretion from the minicluster host until the mass exceeds the self-interaction stability threshold, then integrates the growth timescale against the age of the universe to obtain the δ ≳ 100 condition for QCD axions. This threshold and growth rate are taken from established axion-star equations in the literature rather than being defined in terms of the target bosenova outcome. No self-definitional loops, fitted parameters renamed as predictions, or load-bearing self-citation chains appear in the derivation chain. The central result remains a model-based prediction that can be tested against external simulations or observations and does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis rests on standard axion-field equations and minicluster formation assumptions drawn from prior literature; no new free parameters are introduced in the abstract, and no new particles or forces are postulated.

axioms (2)
  • domain assumption Dilute axion stars possess a stable mass limit only when axion self-interaction is included.
    Explicitly stated in the abstract as the reason for incorporating self-interaction effects.
  • domain assumption The axion star accretes mass from the surrounding minicluster on timescales short enough to reach the stability limit within cosmic time.
    Implicit in the claim that bosenova occurs for δ≳100 within the age of the Universe.

pith-pipeline@v0.9.0 · 5682 in / 1424 out tokens · 39810 ms · 2026-05-18T22:41:08.796787+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

35 extracted references · 35 canonical work pages · 12 internal anchors

  1. [1]

    (43) and Eq

    − e− 1 ν2 (κ2 2+κ2 3) i · (κ1κ2)−1 Z d3⃗ξ (1 + ξ2)−4ϕ⃗ κ1(⃗ξ)ϕ⃗ κ2(⃗ξ)ϕ∗ ⃗ κ3(⃗ξ) 2 ≡ λ2ρ2 g 16π2m5ak2g α(ν, ˜N) , (47) where α(ν, ˜N) is a dimensionless numerical factor and κ3 is given by Eq. (43) and Eq. (44). For ⃗ κ1 and ⃗ κ2, we use the approximation ϕ⃗ κ(⃗ξ) = J0 2 r β ξ − ⃗ n· ⃗ξ ! . (48) For ⃗ κ3, the approximation can still be used with a slight...

    work page

  2. [2]

    R. D. Peccei and H. R. Quinn, CP Conservation in the Presence of Instantons, Phys. Rev. Lett. 38, 1440 (1977)

    work page 1977

  3. [3]

    Weinberg, A New Light Boson?, Phys

    S. Weinberg, A New Light Boson?, Phys. Rev. Lett. 40, 223 (1978)

    work page 1978

  4. [4]

    Wilczek, Problem of Strong P and T Invariance in the Presence of Instantons, Phys

    F. Wilczek, Problem of Strong P and T Invariance in the Presence of Instantons, Phys. Rev. Lett. 40, 279 (1978)

    work page 1978

  5. [5]

    Sikivie, Of Axions, Domain Walls and the Early Uni- verse, Phys

    P. Sikivie, Of Axions, Domain Walls and the Early Uni- verse, Phys. Rev. Lett. 48, 1156 (1982)

    work page 1982

  6. [6]

    Preskill, M

    J. Preskill, M. B. Wise, and F. Wilczek, Cosmology of the Invisible Axion, Phys. Lett. B 120, 127 (1983)

    work page 1983

  7. [7]

    L. F. Abbott and P. Sikivie, A Cosmological Bound on the Invisible Axion, Phys. Lett. B 120, 133 (1983)

    work page 1983

  8. [8]

    Dine and W

    M. Dine and W. Fischler, The Not So Harmless Axion, Phys. Lett. B 120, 137 (1983)

    work page 1983

  9. [9]

    Z. G. Berezhiani, A. S. Sakharov, and M. Y. Khlopov, Primordial background of cosmological axions, Sov. J. Nucl. Phys. 55, 1063 (1992)

    work page 1992

  10. [10]

    Axions In String Theory

    P. Svrcek and E. Witten, Axions In String Theory, JHEP 06, 051, arXiv:hep-th/0605206

    work page internal anchor Pith review Pith/arXiv arXiv

  11. [11]

    The landscape of QCD axion models

    L. Di Luzio, M. Giannotti, E. Nardi, and L. Visinelli, The landscape of QCD axion models, Phys. Rept. 870, 1 (2020), arXiv:2003.01100 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2020

  12. [12]

    E. W. Kolb and I. I. Tkachev, Large amplitude isothermal fluctuations and high density dark matter clumps, Phys. Rev. D 50, 769 (1994), arXiv:astro-ph/9403011

    work page internal anchor Pith review Pith/arXiv arXiv 1994

  13. [13]

    Studies of the motion and decay of axion walls bounded by strings

    S. Chang, C. Hagmann, and P. Sikivie, Studies of the mo- tion and decay of axion walls bounded by strings, Phys. Rev. D 59, 023505 (1999), arXiv:hep-ph/9807374

    work page internal anchor Pith review Pith/arXiv arXiv 1999

  14. [14]

    Fairbairn, D

    M. Fairbairn, D. J. E. Marsh, J. Quevillon, and S. Rozier, Structure formation and microlensing with axion miniclusters, Phys. Rev. D 97, 083502 (2018), arXiv:1707.03310 [astro-ph.CO]

    work page arXiv 2018

  15. [15]

    E. R. Siegel, M. P. Hertzberg, and J. N. Fry, Probing Dark Matter Substructure with Pulsar Timing, Mon. Not. Roy. Astron. Soc. 382, 879 (2007), arXiv:astro- ph/0702546

    work page arXiv 2007

  16. [16]

    V. S. H. Lee, A. Mitridate, T. Trickle, and K. M. Zurek, Probing Small-Scale Power Spectra with Pulsar Timing Arrays, JHEP 06, 028, arXiv:2012.09857 [astro-ph.CO]

    work page arXiv 2012

  17. [17]

    Q. Qiu, Y. Gao, H.-j. Tian, K. Wang, Z. Wang, and X.-M. Yang, Wide binary evaporation by dark solitons: implications from the GAIA catalog, JCAP 02, 001, arXiv:2404.18099 [hep-ph]

    work page arXiv

  18. [18]

    Wang and Y

    Z. Wang and Y. Gao, Axion minicluster halo limits from wide binary disruption, Phys. Rev. D111, 043042 (2025), arXiv:2409.02468 [hep-ph]

    work page arXiv 2025

  19. [19]

    E. W. Kolb and I. I. Tkachev, Axion miniclusters and Bose stars, Phys. Rev. Lett. 71, 3051 (1993), arXiv:hep- ph/9303313

    work page arXiv 1993

  20. [20]

    E. D. Schiappacasse and M. P. Hertzberg, Analysis of Dark Matter Axion Clumps with Spherical Symme- try, JCAP 01, 037, [Erratum: JCAP 03, E01 (2018)], arXiv:1710.04729 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  21. [21]

    Eggemeier and J

    B. Eggemeier and J. C. Niemeyer, Formation and mass 10 growth of axion stars in axion miniclusters, Phys. Rev. D 100, 063528 (2019), arXiv:1906.01348 [astro-ph.CO]

    work page arXiv 2019

  22. [22]

    D. G. Levkov, A. G. Panin, and I. I. Tkachev, Gravita- tional Bose-Einstein condensation in the kinetic regime, Phys. Rev. Lett. 121, 151301 (2018), arXiv:1804.05857 [astro-ph.CO]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  23. [23]

    J. H.-H. Chan, S. Sibiryakov, and W. Xue, Conden- sation and evaporation of boson stars, JHEP 01, 071, arXiv:2207.04057 [astro-ph.CO]

    work page arXiv

  24. [24]

    Budker, J

    D. Budker, J. Eby, M. Gorghetto, M. Jiang, and G. Perez, A generic formation mechanism of ultralight dark matter solar halos, JCAP 12, 021, arXiv:2306.12477 [hep-ph]

    work page arXiv

  25. [25]

    D. G. Levkov, A. G. Panin, and I. I. Tkachev, Relativistic axions from collapsing Bose stars, Phys. Rev. Lett. 118, 011301 (2017), arXiv:1609.03611 [astro-ph.CO]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  26. [26]

    H. Di, L. Shao, Z. Yi, and S.-B. Kong, Novel standard candle: Collapsing axion stars, Phys. Rev. D 110, 103031 (2024), arXiv:2409.05120 [hep-ph]

    work page arXiv 2024

  27. [27]

    P. J. Fox, N. Weiner, and H. Xiao, Recurrent axion stars collapse with dark radiation emission and their cosmo- logical constraints, Phys. Rev. D 108, 095043 (2023), arXiv:2302.00685 [hep-ph]

    work page arXiv 2023

  28. [28]

    Escudero, C

    M. Escudero, C. K. Pooni, M. Fairbairn, D. Blas, X. Du, and D. J. E. Marsh, Axion star explosions: A new source for axion indirect detection, Phys. Rev. D 109, 043018 (2024), arXiv:2302.10206 [hep-ph]

    work page arXiv 2024

  29. [29]

    Di, Stimulated decay of collapsing axion stars and fast radio bursts, Eur

    H. Di, Stimulated decay of collapsing axion stars and fast radio bursts, Eur. Phys. J. C 84, 283 (2024), arXiv:2311.06860 [hep-ph]

    work page arXiv 2024

  30. [30]

    The QCD axion, precisely

    G. Grilli di Cortona, E. Hardy, J. Pardo Vega, and G. Villadoro, The QCD axion, precisely, JHEP 01, 034, arXiv:1511.02867 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv

  31. [31]

    L. D. Landau and E. M. Lifshitz, Course of theoretical physics III: Quantum mechanics, Non-relativistic theory (3rd ed.) (Pergamon Press, 1977)

    work page 1977

  32. [32]

    Axion Cosmology Revisited

    O. Wantz and E. P. S. Shellard, Axion Cosmology Revis- ited, Phys. Rev. D 82, 123508 (2010), arXiv:0910.1066 [astro-ph.CO]

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  33. [33]

    Early seeds of axion miniclusters

    A. Vaquero, J. Redondo, and J. Stadler, Early seeds of axion miniclusters, JCAP 04, 012, arXiv:1809.09241 [astro-ph.CO]

    work page internal anchor Pith review Pith/arXiv arXiv

  34. [34]

    Standard Model thermodynamics across the electroweak crossover

    M. Laine and M. Meyer, Standard Model thermody- namics across the electroweak crossover, JCAP 07, 035, arXiv:1503.04935 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv

  35. [35]

    P. J. McMillan, Mass models of the Milky Way, Mon. Not. Roy. Astron. Soc.414, 2446 (2011), arXiv:1102.4340 [astro-ph.GA]

    work page internal anchor Pith review Pith/arXiv arXiv 2011