arxiv: 2604.21239 · v1 · submitted 2026-04-23 · 🌀 gr-qc · astro-ph.CO· hep-ph

Recognition: unknown

Relativistic frequency shifts in gravitational waves from axion clouds

Takuya Takahashi

Pith reviewed 2026-05-09 21:40 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.COhep-ph

keywords axion cloudssuperradiancegravitational wavesfrequency shiftsblack holesperturbation theoryself-interactions

0 comments

The pith

A bilinear-form relativistic perturbation theory unifies calculations of frequency shifts in gravitational waves from axion clouds around black holes, including cases with multiple excited modes.

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

The paper seeks accurate predictions for the frequencies and their time evolution in continuous gravitational waves emitted by axion clouds that form around rotating black holes through superradiant instability. Axion self-interactions can populate several superradiant modes at once, opening multiple emission channels whose frequency shifts must be tracked precisely if detectors are to identify the signals. The authors apply relativistic perturbation theory built on a bilinear form to derive these shifts in a single framework that covers both self-interaction effects and the cloud's own gravity. This matters because next-generation gravitational-wave observatories will need reliable templates to search for ultralight bosons in the data.

Core claim

Relativistic perturbation theory based on a bilinear form calculates the frequency shifts of gravitational waves emitted by axion clouds, supplying a simple unified treatment that works for self-interaction-driven multiple modes and that revisits the self-gravity contribution.

What carries the argument

Relativistic perturbation theory based on a bilinear form, which computes frequency shifts by treating the interactions among superradiant modes and the cloud's gravitational field.

If this is right

Frequency shifts can be computed uniformly for gravitational-wave emission channels that arise when several superradiant modes are populated.
The same framework incorporates the cloud's self-gravity without separate ad-hoc adjustments.
Predictions become available for the long-term frequency evolution of signals that future detectors must match.
The method extends directly to any situation in which multiple bosonic modes interact inside a black-hole cloud.

Where Pith is reading between the lines

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

The framework could be reused for other ultralight scalar or vector fields that form clouds around black holes.
Accurate frequency templates derived this way would allow continuous-wave searches to place tighter bounds on boson masses and couplings.
If the bilinear approach holds, it offers a computationally lighter alternative to full nonlinear simulations for initial signal modeling.

Load-bearing premise

The bilinear-form perturbation theory remains valid and accurate when self-interactions excite multiple modes and the cloud's self-gravity is included, without higher-order nonlinear effects becoming dominant.

What would settle it

A high-resolution numerical evolution of an axion cloud with self-interactions that tracks the actual gravitational-wave frequencies and checks whether they match the shifts predicted by the bilinear theory or deviate once multiple modes and self-gravity are fully nonlinear.

Figures

Figures reproduced from arXiv: 2604.21239 by Takuya Takahashi.

**Figure 3.** Figure 3: FIG. 3. GW frequencies including frequency shifts, together [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗

read the original abstract

Superradiant instability of ultralight bosons can produce clouds around rotating black holes, whose continuous gravitational wave (GW) emission is a promising observational target. Precise predictions of the signal frequency and its evolution are essential for detecting such continuous GWs. For axions, self-interactions can populate multiple superradiant modes via nonlinear couplings, and GW emission can occur through various channels. To calculate the frequency shifts of GWs emitted through these channels, we employ relativistic perturbation theory based on a bilinear form. We apply this framework to self-interaction effects for the first time, and also revisit the treatment of the self-gravity contribution. Our results provide a simple and unified framework for calculating frequency shifts, including cases in which multiple modes are excited, and are relevant for next-generation GW observations.

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.

Desk Editor's Note private letter to a colleague

The paper gives a unified bilinear perturbation framework for GW frequency shifts from multi-mode axion clouds but leaves the handling of nonlinear self-interaction terms without clear error bounds or checks.

read the letter

The core advance here is applying the bilinear form of relativistic perturbation theory to self-interaction effects in axion clouds for the first time, while folding self-gravity into the same treatment for cases with several superradiant modes excited. That produces a single set of expressions for the frequency shifts across different emission channels, which the author positions as simpler than prior piecemeal calculations. If the derivations hold, this could feed directly into template banks for continuous-wave searches targeting ultralight axion signals. The paper earns credit for identifying the observational need for precise frequency evolution and for trying to keep the framework compact enough to be usable. The approach is presented as a derivation rather than a fit, which is the right direction. The main limitation is that the abstract and the stress-test description give no indication of explicit convergence tests, comparisons against known single-mode limits, or estimates of the size of neglected quadratic and higher terms once multiple modes are populated. The claim that the bilinear truncation remains accurate under simultaneous self-interactions and self-gravity therefore sits on an unverified assumption. If the full text contains those checks and they are clean, the concern shrinks; if they are absent or only sketched, the central result stays provisional. This is the kind of technical methods paper that belongs in the gravitational-wave and ultralight dark-matter literature. Readers who already work with superradiant clouds and need updated frequency predictions would get the most out of it, provided the math checks out. It is worth sending to peer review so that specialists can examine the perturbation construction and the truncation justification in detail.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a relativistic perturbation theory framework based on a bilinear form to compute frequency shifts in gravitational waves emitted by axion clouds around Kerr black holes. It applies the method to self-interaction effects (claimed for the first time) that can excite multiple superradiant modes via nonlinear couplings, revisits the self-gravity contribution, and asserts that the approach supplies a simple unified framework for frequency-shift calculations even in multi-mode cases, with relevance to next-generation GW observations.

Significance. If the bilinear-form truncation is shown to be consistent and higher-order nonlinear back-reaction remains negligible, the unified treatment would supply a practical tool for precise frequency predictions needed to search for continuous GW signals from superradiant axion clouds. The absence of free parameters in the core derivation and the explicit inclusion of multi-mode self-interactions would strengthen its utility for observational modeling.

major comments (2)

[Abstract and multi-mode results] Abstract and the multi-mode application section: the central claim that the bilinear-form relativistic perturbation theory supplies an accurate unified framework when self-interactions populate multiple superradiant modes rests on the assertion that all relevant frequency-shift contributions are captured at linear order. No explicit error bounds, convergence tests, or demonstration that quadratic and higher nonlinear terms remain negligible at the retained order are provided for the multi-mode case, leaving the truncation unverified.
[Self-gravity treatment] Self-gravity revisit section: the paper states that self-gravity is revisited within the same bilinear structure, yet no direct comparison is shown between the new bilinear treatment and prior calculations to confirm consistency of the frequency-shift correction at the order relevant for GW emission.

minor comments (2)

[Method introduction] Notation for the bilinear form is introduced without an explicit definition of the inner product or the precise perturbation variable in the opening paragraphs of the method section; adding this would improve readability.
[Derivation of frequency shift] The abstract claims the framework is 'simple,' but several intermediate steps in the derivation of the frequency-shift formula are only sketched; expanding one or two key algebraic steps would aid verification.

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 revised the manuscript to incorporate additional discussion and comparisons where the points identify gaps in the original presentation.

read point-by-point responses

Referee: [Abstract and multi-mode results] Abstract and the multi-mode application section: the central claim that the bilinear-form relativistic perturbation theory supplies an accurate unified framework when self-interactions populate multiple superradiant modes rests on the assertion that all relevant frequency-shift contributions are captured at linear order. No explicit error bounds, convergence tests, or demonstration that quadratic and higher nonlinear terms remain negligible at the retained order are provided for the multi-mode case, leaving the truncation unverified.

Authors: We agree that the original manuscript lacks explicit error bounds or convergence tests for the multi-mode truncation. The bilinear framework is constructed so that frequency shifts arise at linear order in the metric and axion perturbations, with self-interactions entering through the quadratic source terms in the mode equations. Higher-order contributions are parametrically suppressed by the smallness of the cloud amplitude relative to the black-hole mass and by the slow secular evolution. In the revised manuscript we have added a dedicated paragraph in the multi-mode section that supplies order-of-magnitude estimates for the size of the neglected quadratic and cubic terms under typical superradiant-cloud parameters, together with a brief consistency argument showing that these terms remain below the precision required for next-generation continuous-wave searches. revision: yes
Referee: [Self-gravity treatment] Self-gravity revisit section: the paper states that self-gravity is revisited within the same bilinear structure, yet no direct comparison is shown between the new bilinear treatment and prior calculations to confirm consistency of the frequency-shift correction at the order relevant for GW emission.

Authors: The referee is correct that no side-by-side comparison was presented. The bilinear derivation recovers the standard self-gravity frequency shift in the single-mode limit by construction. To make the agreement explicit, the revised manuscript now includes a short subsection that analytically reduces the bilinear result to the known expression from earlier literature and tabulates numerical values of the frequency correction for representative cloud parameters, confirming agreement at the order relevant for gravitational-wave emission. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation from established perturbation theory is self-contained

full rationale

The paper derives GW frequency shifts via relativistic perturbation theory based on a bilinear form, applying the framework to self-interactions for the first time while revisiting self-gravity. No quoted equations or steps reduce by construction to fitted inputs, self-definitions, or load-bearing self-citations whose validity depends on the present work. The central claim of a unified framework for multi-mode cases rests on the bilinear perturbation structure rather than renaming or reparameterizing known results. The approach is presented as an application of prior theory with novel extensions, consistent with an independent derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the applicability of relativistic perturbation theory and the bilinear form to the axion-cloud system; no new free parameters or invented entities are introduced in the abstract.

axioms (1)

domain assumption Relativistic perturbation theory based on a bilinear form accurately captures frequency shifts from self-interactions and self-gravity in axion clouds.
Invoked to derive the unified framework for GW frequency shifts.

pith-pipeline@v0.9.0 · 5426 in / 1150 out tokens · 33667 ms · 2026-05-09T21:40:12.346280+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

66 extracted references · 52 canonical work pages · 3 internal anchors

[1]

(4), which we write asOϕ= 0

Bilinear product As the nonperturbative background, we begin with the Lagrangian density for the scalar field, L=− 1 2 √−g gµν∇µϕ∇νϕ+µ 2ϕ2 .(19) The equation of motion is given by Eq. (4), which we write asOϕ= 0. In what follows, we work with the constant-thypersurface in Boyer–Lindquist coordinates. Here, following Ref. [56], we define a bilinear product...
[2]

Perturbation theory For the development of perturbation theory, it is con- venient to move to the Hamiltonian formulation. Here, we introduce the canonical momentum conjugate to the scalar fieldϕ, π= ∂L ∂ ˙ϕ = √−g n µ∇µϕ ,(24) the Hamiltonian density is given by H=π ˙ϕ− L =− 1 2√−gg tt π2 − gti gtt π∂iϕ + 1 2 √−g Aij∂iϕ∂jϕ+µ 2ϕ2 .(25) whereA ij =g ij −g t...

2025
[3]

R. D. Peccei and H. R. Quinn,Phys. Rev. Lett.38, 1440 (1977), [,328(1977)]

1977
[4]

Weinberg,Phys

S. Weinberg,Phys. Rev. Lett.40, 223 (1978)

1978
[5]

Wilczek,Phys

F. Wilczek,Phys. Rev. Lett.40, 279 (1978)

1978
[6]

J. E. Kim and G. Carosi,Rev. Mod. Phys.82, 557 (2010), [erratum: Rev. Mod. Phys.91,no.4,049902(2019)], arXiv:0807.3125 [hep-ph]

work page arXiv 2010
[7]

Dine and W

M. Dine and W. Fischler,Phys. Lett.120B, 137 (1983)

1983
[8]

Preskill, M

J. Preskill, M. B. Wise, and F. Wilczek,Phys. Lett. 120B, 127 (1983)

1983
[9]

L. F. Abbott and P. Sikivie,Phys. Lett.120B, 133 (1983)

1983
[10]

L. Hui, J. P. Ostriker, S. Tremaine, and E. Witten, Phys. Rev. D95, 043541 (2017), arXiv:1610.08297 [astro- ph.CO]

work page Pith review arXiv 2017
[11]

String Axiverse

A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper, and J. March-Russell,Phys. Rev.D81, 123530 (2010), arXiv:0905.4720 [hep-th]

work page internal anchor Pith review arXiv 2010
[12]

Axions In String Theory

P. Svrcek and E. Witten,JHEP06, 051 (2006), arXiv:hep-th/0605206 [hep-th]

work page Pith review arXiv 2006
[13]

V. M. Mehta, M. Demirtas, C. Long, D. J. E. Marsh, L. McAllister, and M. J. Stott,JCAP07, 033 (2021), arXiv:2103.06812 [hep-th]

work page arXiv 2021
[14]

W. H. Press and S. A. Teukolsky,Nature238, 211 (1972)

1972
[15]

T. J. M. Zouros and D. M. Eardley,Annals Phys.118, 139 (1979)

1979
[16]

S. L. Detweiler,Phys. Rev.D22, 2323 (1980)

1980
[17]

S. R. Dolan,Phys. Rev.D76, 084001 (2007), arXiv:0705.2880 [gr-qc]

work page Pith review arXiv 2007
[18]

Brito, V

R. Brito, V. Cardoso, and P. Pani,Lect. Notes Phys.971, pp.1 (2020), arXiv:1501.06570 [gr-qc]

work page arXiv 2020
[19]

Exploring the String Axiverse with Precision Black Hole Physics

A. Arvanitaki and S. Dubovsky,Phys. Rev.D83, 044026 (2011), arXiv:1004.3558 [hep-th]

work page Pith review arXiv 2011
[20]

Brito, V

R. Brito, V. Cardoso, and P. Pani,Class. Quant. Grav. 32, 134001 (2015), arXiv:1411.0686 [gr-qc]

work page arXiv 2015
[21]

M. J. Stott, arXiv:2009.07206 [hep-ph]

work page arXiv 2009
[22]

Baumann, H.S

D. Baumann, H. S. Chia, and R. A. Porto,Phys. Rev. D 99, 044001 (2019), arXiv:1804.03208 [gr-qc]

work page arXiv 2019
[23]

Baumann, H.S

D. Baumann, H. S. Chia, R. A. Porto, and J. Stout,Phys. Rev. D101, 083019 (2020), arXiv:1912.04932 [gr-qc]

work page arXiv 2020
[24]

Baryakhtaret al., in2022 Snowmass Summer Study (2022) arXiv:2203.07984 [hep-ph]

M. Baryakhtaret al., in2022 Snowmass Summer Study (2022) arXiv:2203.07984 [hep-ph]

work page arXiv 2022
[25]

Dyson, T

C. Dyson, T. F. M. Spieksma, R. Brito, M. van de Meent, and S. Dolan,Phys. Rev. Lett.134, 211403 (2025), arXiv:2501.09806 [gr-qc]

work page arXiv 2025
[26]

R. Roy, S. Vagnozzi, and L. Visinelli,Phys. Rev. D105, 083002 (2022), arXiv:2112.06932 [astro-ph.HE]

work page arXiv 2022
[27]

Y. Chen, J. Shu, X. Xue, Q. Yuan, and Y. Zhao,Phys. Rev. Lett.124, 061102 (2020), arXiv:1905.02213 [hep- ph]

work page arXiv 2020
[28]

Y. Chen, C. Li, Y. Mizuno, J. Shu, X. Xue, Q. Yuan, Y. Zhao, and Z. Zhou,JCAP09, 073 (2022), arXiv:2208.05724 [hep-ph]

work page arXiv 2022
[29]

Brito, S

R. Brito, S. Ghosh, E. Barausse, E. Berti, V. Cardoso, I. Dvorkin, A. Klein, and P. Pani,Phys. Rev. Lett.119, 131101 (2017), arXiv:1706.05097 [gr-qc]

work page arXiv 2017
[30]

Tsukada, T

L. Tsukada, T. Callister, A. Matas, and P. Meyers, Phys. Rev. D99, 103015 (2019), arXiv:1812.09622 [astro- ph.HE]

work page arXiv 2019
[31]

Gavilan-Martin, O

D. Gavilan-Martin, O. Simon, D. Krishna, D. F. J. Kim- ball, D. Budker, and A. Wickenbrock, arXiv:2602.23415 [hep-ph]

work page arXiv
[32]

Abbottet al.(LIGO Scientific, Virgo, KAGRA),Phys

R. Abbottet al.(LIGO Scientific, Virgo, KAGRA),Phys. Rev. D105, 102001 (2022), arXiv:2111.15507 [astro- ph.HE]

work page arXiv 2022
[33]

Abbottet al.(KAGRA, LIGO Scientific, VIRGO), Phys

R. Abbottet al.(KAGRA, LIGO Scientific, VIRGO), Phys. Rev. D106, 042003 (2022), arXiv:2204.04523 [astro-ph.HE]

work page arXiv 2022
[34]

A. G. Abacet al.(LIGO Scientific, VIRGO, KAGRA), arXiv:2509.07352 [gr-qc]

work page arXiv
[35]

Yoshino and H

H. Yoshino and H. Kodama,PTEP2014, 043E02 (2014), arXiv:1312.2326 [gr-qc]

work page arXiv 2014
[36]

Yoshino and H

H. Yoshino and H. Kodama,PTEP2015, 061E01 (2015), arXiv:1407.2030 [gr-qc]

work page arXiv 2015
[37]

Brito, S

R. Brito, S. Ghosh, E. Barausse, E. Berti, V. Cardoso, I. Dvorkin, A. Klein, and P. Pani,Phys. Rev.D96, 064050 (2017), arXiv:1706.06311 [gr-qc]

work page arXiv 2017
[38]

Palombaet al.,Phys

C. Palombaet al.,Phys. Rev. Lett.123, 171101 (2019), arXiv:1909.08854 [astro-ph.HE]

work page arXiv 2019
[39]

S. J. Zhu, M. Baryakhtar, M. A. Papa, D. Tsuna, N. Kawanaka, and H.-B. Eggenstein,Phys. Rev. D102, 063020 (2020), arXiv:2003.03359 [gr-qc]

work page arXiv 2020
[40]

Yoshino and H

H. Yoshino and H. Kodama,Prog. Theor. Phys.128, 153 (2012), arXiv:1203.5070 [gr-qc]

work page arXiv 2012
[41]

Yoshino and H

H. Yoshino and H. Kodama,Class. Quant. Grav.32, 214001 (2015), arXiv:1505.00714 [gr-qc]

work page arXiv 2015
[42]

Takahashi, H

T. Takahashi, H. Omiya, and T. Tanaka,Phys. Rev. D 110, 104038 (2024), arXiv:2408.08349 [gr-qc]

work page arXiv 2024
[43]

Collaviti, L

S. Collaviti, L. Sun, M. Galanis, and M. Baryakhtar, Class. Quant. Grav.42, 025006 (2025), arXiv:2407.04304 11 [gr-qc]

work page arXiv 2025
[44]

Omiya, T

H. Omiya, T. Takahashi, T. Tanaka, and H. Yoshino, Phys. Rev. D110, 044002 (2024), arXiv:2404.16265 [gr- qc]

work page arXiv 2024
[45]

Gruzinov, arXiv:1604.06422 [astro-ph.HE]

A. Gruzinov, arXiv:1604.06422 [astro-ph.HE]

work page arXiv
[46]

Baryakhtar, M

M. Baryakhtar, M. Galanis, R. Lasenby, and O. Simon, Phys. Rev. D103, 095019 (2021), arXiv:2011.11646 [hep- ph]

work page arXiv 2021
[47]

Omiya, T

H. Omiya, T. Takahashi, T. Tanaka, and H. Yoshino, JCAP06, 016 (2023), arXiv:2211.01949 [gr-qc]

work page arXiv 2023
[48]

S. J. Witte and A. Mummery,Phys. Rev. D111, 083044 (2025), arXiv:2412.03655 [hep-ph]

work page arXiv 2025
[49]

Punturoet al.,Class

M. Punturoet al.,Class. Quant. Grav.27, 194002 (2010)

2010
[50]

Cosmic Explorer: The U.S. Contribution to Gravitational-Wave Astronomy beyond LIGO

D. Reitzeet al.,Bull. Am. Astron. Soc.51, 035 (2019), arXiv:1907.04833 [astro-ph.IM]

work page internal anchor Pith review arXiv 2019
[51]

Kawamuraet al.,Class

S. Kawamuraet al.,Class. Quant. Grav.28, 094011 (2011)

2011
[52]

Takahashi, H

T. Takahashi, H. Omiya, and T. Tanaka,PTEP2022, 043E01 (2022), arXiv:2112.05774 [gr-qc]

work page arXiv 2022
[53]

Siemonsen, T

N. Siemonsen, T. May, and W. E. East, arXiv:2211.03845 [gr-qc]

work page arXiv
[54]

T. May, W. E. East, and N. Siemonsen,Phys. Rev. D 111, 044062 (2025), arXiv:2410.21442 [gr-qc]

work page arXiv 2025
[55]

Omiya, T

H. Omiya, T. Takahashi, and T. Tanaka,PTEP2021, 043E02 (2021), arXiv:2012.03473 [gr-qc]

work page arXiv 2021
[56]

Omiya, T

H. Omiya, T. Takahashi, and T. Tanaka,PTEP2022, 043E03 (2022), arXiv:2201.04382 [gr-qc]

work page arXiv 2022
[57]

S. R. Green, S. Hollands, L. Sberna, V. Toomani, and P. Zimmerman,Phys. Rev. D107, 064030 (2023), arXiv:2210.15935 [gr-qc]

work page arXiv 2023
[58]

Cannizzaro, L

E. Cannizzaro, L. Sberna, S. R. Green, and S. Hollands, Phys. Rev. Lett.132, 051401 (2024), arXiv:2309.10021 [gr-qc]

work page arXiv 2024
[59]

Cannizzaro, M

E. Cannizzaro, M. Palleschi, L. Sberna, R. Brito, and S. Green, arXiv:2512.15878 [gr-qc]

work page arXiv
[60]

Bilinear products and the orthogonality of quasinormal modes on hyperboloidal foliations

M. Minucci, R. Panosso Macedo, C. Pantelidou, and L. Sberna, arXiv:2604.13182 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv
[61]

D. R. Brill, P. L. Chrzanowski, C. M. Pereira, E. D. Fack- erell, and J. R. Ipser,Phys. Rev.D5, 1913 (1972)

1913
[62]

Baumann, H.S

D. Baumann, H. S. Chia, J. Stout, and L. ter Haar,JCAP 12, 006 (2019), arXiv:1908.10370 [gr-qc]

work page arXiv 2019
[63]

L. Fu, H. Omiya, T. Tanaka, X. Tong, Y. Wang, and H.-Y. Zhu, arXiv:2512.06790 [gr-qc]

work page arXiv
[64]

E. W. Leaver,Journal of Mathematical Physics27, 1238 (1986), https://doi.org/10.1063/1.527130

work page doi:10.1063/1.527130 1986
[65]

Hollands and V

S. Hollands and V. Toomani,Class. Quant. Grav.43, 055001 (2026), arXiv:2405.18604 [gr-qc]

work page arXiv 2026
[66]

Wardell, C

B. Wardell, C. Kavanagh, and S. R. Dolan,Class. Quant. Grav.42, 205007 (2025), arXiv:2406.12510 [gr-qc]

work page arXiv 2025