REVIEW 4 major objections 5 minor 1 cited by
Resonance of black hole quasinormal modes in coupled systems
T0 review · 4 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Coupled black hole quasinormal modes from different fields repel and amplify at avoided crossings.
desk verdict Careful, honest formalism for coupled-field QNM excitation, with a real avoided-crossing example in EMA theory—but the claimed amplification may be partly a normalization artifact and needs a time-domain check. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Three ingredients carry the argument. The basis-invariant excitation factor $B_{\alpha n}=i\det A^{(out)}(2i\omega_{\alpha n})^{-N}\left[d\det A^{(in)}/d\omega\right]^{-1}$ in Eq. (21), built from the coefficient matrices of in-going and outgoing waves, generalizes the single-field excitation factor while removing the field-basis ambiguity. The coupled master equations (48)--(50) for $\psi$, $Z_+$, and $Z_-$ describe the axial sector of the Einstein-Maxwell-axion system, with coupling terms proportional to $g_{a\gamma\gamma}Q$. A matrix-valued Leaver continued fraction computes the QNM frequencies, and a two-level non-Hermitian model $\omega_{\pm}^2=E_c\pm\sqrt{E_d^2+\Delta^2}$ explains why the excitation-factor difference follows a square-root Lorentzian near the crossing.
What would settle it
Time-evolve the axial perturbation system (48)--(50) with a localized broadband source for $Q/M=0.1$ and $g_{a\gamma\gamma}$ on both sides of 28, then decompose the late-time ringdown; the claimed resonance is ruled out if the gravity and EM channel amplitudes do not peak where $|\omega_{g0}-\omega_{e0}|$ is smallest. Alternatively, fix $g_{a\gamma\gamma}Q/M\simeq 2.83$ and reduce $Q/M$; the peak of $|B_{g0}-B_{e0}|$ should grow roughly as $(Q/M)^{-1}$ until the large-coupling instability intervenes.
Extended reading notes
Core claim
On a Reissner-Nordström black hole with $l=2$ and $Q/M=0.1$, the gravity-led and EM-led fundamental quasinormal frequencies approach each other as the axion-photon coupling $g_{a\gamma\gamma}$ grows, then repel in both real and imaginary parts near $g_{a\gamma\gamma}\simeq 28$; the same repulsion appears when $Q/M$ is varied at fixed $g_{a\gamma\gamma}=10$, near $Q/M\simeq 0.26$. The paper shows that the excitation factors $B_{g0}$ and $B_{e0}$ of the participating modes are amplified at these parameter values, with $|B_{g0}-B_{e0}|\propto |\omega_{g0}-\omega_{e0}|^{-1}$, and reproduces the peak shape with a square-root Lorentzian model. This establishes, in the authors' formulation, a new resonance phenomenon specific to coupled perturbation systems, distinct from the overtone resonance found in Kerr black holes.
Load-bearing premise
The load-bearing premise is that the coupled master equations (48)--(50) faithfully represent the axial perturbations of the Einstein-Maxwell-axion black hole and that each QNM branch can be unambiguously labeled by its decoupling-limit origin; if either part fails, the resonance between fundamental modes of different fields is not established.
Editorial extensions
If this is right
- Resonant amplification can occur between the longest-lived quasinormal modes of different fields, not only between highly damped overtones of a single field, widening the observational window for black hole spectroscopy.
- The definition of $B_{\alpha n}$ in Eq. (21) applies to any system of $N$ coupled perturbation variables with a common horizon radius, so it provides a general diagnostic for resonance in modified-gravity and dark-matter-inspired models.
- In the Einstein-Maxwell-axion system, the resonance location is controlled by the combination $g_{a\gamma\gamma}Q/M$ for small $Q/M$, with the spectral-repulsion width narrowing as $Q$ decreases and the peak excitation growing as $Q^{-1}$.
- The same system becomes linearly unstable for large coupling, with a purely imaginary mode turning unstable, so the resonant behavior is confined to the stable, moderate-coupling regime.
- Because the avoided crossing occurs at nearly the same parameter values across overtones, the resonance is not a fine-tuned single-mode effect but a robust spectral feature.
Reading between the lines
- If the basis-invariant excitation factor is adopted as the standard measure, the resonance could be used to constrain axion-photon couplings from ringdown amplitude ratios, not just frequencies, once gravitational and electromagnetic channels are separately resolved.
- The authors' two-level model suggests that the square-root-Lorentzian peak shape is generic for coupled fields, while the lemniscate shape in the Kerr case is specific to overtone coupling; a systematic survey of two-field systems could test this distinction.
- The $Q\to 0$ divergence of the peak excitation hints that extremely small charges with large couplings could produce very narrow resonances, but the simultaneous appearance of instability sets an upper bound on the attainable amplification.
- A direct time-domain simulation of a pulse in the axial sector around an EMA black hole would convert the frequency-domain resonance prediction into a waveform prediction, making the effect testable against template searches.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a basis-independent definition of excitation factors for quasinormal modes of black holes whose perturbation equations form a system of N coupled second-order ODEs. The construction is based on the Green's matrix and Wronskian formalism, and the proposed quantity B_alpha n is designed to reduce to the standard single-field excitation factor in the decoupled limit. The formalism is then applied to the axial sector of a Reissner-Nordström black hole in Einstein-Maxwell-axion theory, where three degrees of freedom (gravitational, electromagnetic, and axion) are coupled. Using a continued-fraction method, the authors report avoided crossings between gravity-led and EM-led fundamental QNMs, accompanied by amplified excitation factors, and interpret this as a resonant excitation. A two-level non-Hermitian model with aligned phases is used to reproduce the spectral trajectories and the scaling of the excitation factors.
Significance. If the claimed effect is genuine, the paper identifies a new type of QNM resonance that is distinct from the overtone avoided crossings in Kerr: it occurs between the longest-lived modes that originate from fundamental modes of different fields. The Green's matrix/Wronskian derivation is careful, the basis-independence argument is clearly presented, and the EMA system provides a concrete, well-motivated example. The explicit recurrence coefficients in Appendix B and the simple two-level model with a falsifiable square-root Lorentzian scaling are useful assets. However, the central claim is not yet established at the level of a physical observable: the excitation factor B_alpha n carries a determinant normalization whose growth near a near-degeneracy may be partly a normalization artifact, and the paper itself notes in footnote *3 that a different channel-adapted normalization gives qualitatively different behavior. No code, convergence tables, error bars, or time-domain waveforms are provided, so the numerical demonstration and its interpretation remain under-supported.
major comments (4)
- [Sec. II B 2, Eq. (21)]
- [Sec. III C and Figs. 7–8]
- [Sec. III C 2 and Fig. 1]
- [Sec. III D 2, Eq. (75)]
minor comments (5)
- [Eq. (19)]
- [Eqs. (39)–(41)]
- [Sec. III C 1]
- [Sec. III D 1]
- [Conclusions]
Circularity Check
The growth of Bαn near the avoided crossing is built into the determinant normalization of Eq. (21), so the 'resonant excitation' confirmation is partly self-definitional; the avoided-crossing spectrum itself is an independent numerical result.
-
self definitional
[Sec. II B 2, Eq. (21); Sec. III D 2, Eq. (75) and Fig. 8; footnote *3]
"Bαn = i detA(out) / ((2iωαn)^N) (d/dω detA(in))^{-1}|_{ω=ωαn} ... Written explicitly, the extra factor takes the form of (2iω)^{-(N-1)} ∏_{β≠α} A^{(out,β)}_β / A^{(in,β)}_β. ... In Ref. [21], it was demonstrated that |B_+−B_-| ∝ |ω_+−ω_-|^{-1}. ... They show excellent agreement ... This result confirms the validity of the model and the resonant excitation."
By construction, Bαn contains the inverse derivative (d/dω detA(in))^{-1} evaluated at the QNM pole. Since detA(in) vanishes at both ωg0 and ωe0 near an avoided crossing, one has detA(in) ≈ C(ω−ωg0)(ω−ωe0), so (d/dω detA(in))^{-1}|_{ωg0} ∝ |ωg0−ωe0|^{-1}. The amplification of B and the Fig. 8 agreement with |ωg0−ωe0|^{-1} therefore follow algebraically from the chosen determinant normalization, not from an independent calculation of an observable amplitude. The paper itself notes that in the decoupling limit Eq. (21) carries the extra factor ∏_{β≠α} A^{(out,β)}_β / A^{(in,β)}_β, which already diverges when another field's A^{(in,β)}_β vanishes at a level crossing, and footnote *3 reports that a different normalization gives qualitatively different behavior.
full rationale
The QNM frequencies and the avoided-crossing structure in the EMA system are obtained by a direct continued-fraction solution of the coupled master equations (48)-(50), which are taken from Ref. [42]; that spectral part is not circular. The circularity concerns the resonance claim: the excitation factor Bαn in Eq. (21) is defined with the inverse derivative of detA(in), so its growth at near-degenerate QNMs is a built-in property of the determinant normalization. The paper's confirmation of resonance by fitting |Bg0−Be0| to |ωg0−ωe0|^{-1} (or a square-root Lorentzian) is therefore a check of the definition rather than a prediction from an independent observable amplitude. The citation of Ref. [21] for Eq. (75) reinforces this issue, but the relation is derivable from the definition itself, so the self-citation is not the sole load-bearing element. Because the central spectral finding remains an independent numerical result, the overall circularity is partial rather than total.
Assumptions & free parameters
free parameters (4)
- p =
not quoted
- q =
not quoted
- ω0 and ϑ =
not quoted
- f0, x0, γ =
not quoted
assumptions (4)
- domain assumption The coupled master equations (48)-(50), taken from Ref. [42], correctly describe linear axial perturbations of the EMA black hole.
- ad hoc to paper Each QNM of the coupled system can be unambiguously labeled by its decoupling-limit parent as gravity-led, EM-led, or axion-led.
- domain assumption The fields are massless, share the same horizon, and the potential matrix vanishes at both boundaries.
- ad hoc to paper The two-level model (71)-(73) with aligned phases captures the relevant dynamics near the avoided crossing.
Cite this review
Pith. "Pith review of Resonance of black hole quasinormal modes in coupled systems." pith.science (2026). https://pith.science/paper/PFRY6H4I
@misc{pith2026250503883,
author = {Pith},
title = {Pith review of: Resonance of black hole quasinormal modes in coupled systems},
year = {2026},
howpublished = {\url{https://pith.science/paper/PFRY6H4I}},
note = {Machine review of arXiv:2505.03883}
}
read the original abstract
Black hole quasinormal modes (QNMs) can exhibit resonant excitations associated with avoided crossings in their complex frequency spectrum. Such resonance phenomena can serve as novel signatures for probing new physics, where additional degrees of freedom are commonly introduced. Motivated by this possibility, we investigate QNMs in systems where multiple degrees of freedom are coupled with each other, and introduce a definition of excitation factors suitable for such systems. To demonstrate our formulation, we apply it to a black hole in the Einstein-Maxwell-axion theory, where we find that avoided crossings can appear even between longest-lived modes originating from the fundamental modes of different degrees of freedom, in contrast to the Kerr case in General Relativity. We show that the excitation factors are indeed amplified as a manifestation of resonance at parameter values corresponding to the avoided crossings.
Figures
Figures from the paper (5 more)
Forward citations
Cited by 1 Pith paper
-
Exceptional line and pseudospectrum in black hole spectroscopy
A continuous line of exceptional points exists in the three-parameter space of a Gaussian-bump-perturbed Regge-Wheeler potential, with pseudospectral contour sizes scaling as ε^{1/2} at second-order EPs.
Reference graph
Works this paper leans on
-
[1]
For this purpose, we consider Eq
Waveform Next, we introduce excitation factors for coupled systems, which characterize how easily each QNM is excited in response to an external source for each QNM. For this purpose, we consider Eq. (3) with an added source term, d2 dr2∗ +ω2 ⃗ˆΨ(ω,r )−V (r)⃗ˆΨ(ω,r ) = ⃗S(r) . (9) *2 Of course, it may in principle be possible that some QNMs (dis)appear at...
-
[2]
Definition With the setup above, let us define the excitation factors in coupled systems. In the case of a single degree of freedom, i.e., N = 1, the definition of the excitation factor is given by [22, 23] Bsingle n = A(out) 2ωn d dωA(in) −1 ω=ωn . (20) In linear perturbation theory, there is an ambiguity in the overall normalization of the wavefunction....
-
[3]
Continued fraction method First, we describe the computation method used in this paper. A representative calculation method is the direct integration method [44], which is widely applicable in many cases. However, it is not suitable for extracting the *5 The differences from the notation in Ref. [42] are as follows: One of our master variables, hodd, corr...
-
[4]
Numerical results Here, we show numerical results of the QNM frequencies of the BH in the EMA system. In particular, we present results only for l = 2, as it is the smallest l where all the degrees of freedom are coupled and is also considered observationally important. This system has three degrees of freedom, which correspond to purely gravitational, EM...
-
[5]
Instability In the regime where gaγγQ/M is relatively large, the real part of the gravity-led QNM frequency tends to decrease beyond the point where the avoided crossing occurs. In this regime, the convergence of the continued fraction method with respect to the truncation order becomes increasingly poor, especially for higher overtones. Moreover, as inve...
-
[6]
First, we briefly describe the numerical calculation method
Numerical results To obtain the excitation factor, it is necessary to compute the coefficient matrix, which in turn requires solving the homogeneous solutions under each boundary condition. First, we briefly describe the numerical calculation method. As discussed in the previous subsection, the resonant excitation is expected among the QNMs originating fr...
-
[7]
Key feature of the resonant excitation Finally, we discuss the resonance behaviors of QNM spectra and the excitation factors. Although QNMs are difficult to treat analytically, their features can be captured by the following simple model. As proposed in Ref. [21], the theory of QNMs can be well understood in terms of quantum mechanics for non-Hermitian sy...
-
[8]
7, it can be inferred that q is proportional to Q
shown in Fig. 7, it can be inferred that q is proportional to Q. This is due to the indirect coupling between gravity-led and EM-led modes through the axion-led mode. We show in Fig. 8 the difference between excitation factors corresponding to the gravity-led and EM-led QNMs along with the inverse of the difference between QNM frequencies |ωg0−ωe0|−1. We ...
Show all 67 references
-
[9]
B. P. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. Lett. 116, 061102 (2016), arXiv:1602.03837 [gr-qc]
2016 arXiv
-
[10]
Abbott et al
R. Abbott et al. (KAGRA, VIRGO, LIGO Scientific), Phys. Rev. X 13, 041039 (2023), arXiv:2111.03606 [gr-qc]
2023 arXiv
-
[11]
C. V. Vishveshwara, Nature 227, 936 (1970)
1970
-
[12]
W. H. Press, Astrophys. J. Lett. 170, L105 (1971)
1971
-
[13]
S. A. Teukolsky, Astrophys. J. 185, 635 (1973)
1973
-
[14]
Chandrasekhar and S
S. Chandrasekhar and S. L. Detweiler, Proc. Roy. Soc. Lond. A 344, 441 (1975)
1975
-
[15]
K. D. Kokkotas and B. G. Schmidt, Living Rev. Rel. 2, 2 (1999), arXiv:gr-qc/9909058
1999 arXiv
-
[16]
Dreyer, B
O. Dreyer, B. J. Kelly, B. Krishnan, L. S. Finn, D. Garrison, and R. Lopez-Aleman, Class. Quant. Grav. 21, 787 (2004), arXiv:gr-qc/0309007
2004 arXiv
-
[17]
Berti, V
E. Berti, V. Cardoso, and C. M. Will, Phys. Rev. D 73, 064030 (2006), arXiv:gr-qc/0512160
2006 arXiv
-
[18]
Berti, V
E. Berti, V. Cardoso, and A. O. Starinets, Class. Quant. Grav. 26, 163001 (2009), arXiv:0905.2975 [gr-qc]
2009 arXiv
-
[19]
R. A. Konoplya and A. Zhidenko, Rev. Mod. Phys. 83, 793 (2011), arXiv:1102.4014 [gr-qc]. 21
2011 arXiv
-
[20]
Baibhav, M
V. Baibhav, M. H.-Y. Cheung, E. Berti, V. Cardoso, G. Carullo, R. Cotesta, W. Del Pozzo, and F. Duque, Phys. Rev. D 108, 104020 (2023), arXiv:2302.03050 [gr-qc]
2023 arXiv
-
[21]
Cardoso, M
V. Cardoso, M. Kimura, A. Maselli, E. Berti, C. F. B. Macedo, and R. McManus, Phys. Rev. D 99, 104077 (2019), arXiv:1901.01265 [gr-qc]
2019 arXiv
-
[22]
McManus, E
R. McManus, E. Berti, C. F. B. Macedo, M. Kimura, A. Maselli, and V. Cardoso, Phys. Rev. D 100, 044061 (2019), arXiv:1906.05155 [gr-qc]
2019 arXiv
-
[23]
Mukohyama, K
S. Mukohyama, K. Takahashi, K. Tomikawa, and V. Yingcharoenrat, JCAP 07, 050 (2023), arXiv:2304.14304 [gr-qc]
2023 arXiv
-
[24]
R. P. Kerr, Phys. Rev. Lett. 11, 237 (1963)
1963
- [25]
-
[26]
Berti, V
E. Berti, V. Cardoso, K. D. Kokkotas, and H. Onozawa, Phys. Rev. D 68, 124018 (2003), arXiv:hep-th/0307013
2003 arXiv
- [27]
-
[28]
G. B. Cook and M. Zalutskiy, Phys. Rev. D 90, 124021 (2014), arXiv:1410.7698 [gr-qc]
2014 arXiv
-
[29]
Motohashi, Phys
H. Motohashi, Phys. Rev. Lett. 134, 141401 (2025), arXiv:2407.15191 [gr-qc]
2025 arXiv
-
[30]
E. W. Leaver, Phys. Rev. D 34, 384 (1986)
1986
- [31]
-
[32]
Giesler, M
M. Giesler, M. Isi, M. A. Scheel, and S. Teukolsky, Phys. Rev. X 9, 041060 (2019), arXiv:1903.08284 [gr-qc]
2019 arXiv
-
[33]
G. B. Cook, Phys. Rev. D 102, 024027 (2020), arXiv:2004.08347 [gr-qc]
2020 arXiv
- [34]
-
[35]
S. Ma, K. Mitman, L. Sun, N. Deppe, F. H´ ebert, L. E. Kidder, J. Moxon, W. Throwe, N. L. Vu, and Y. Chen, Phys. Rev. D 106, 084036 (2022), arXiv:2207.10870 [gr-qc]
2022 arXiv
-
[36]
Takahashi and H
K. Takahashi and H. Motohashi, Class. Quant. Grav. 41, 195023 (2024), arXiv:2311.12762 [gr-qc]
2024 arXiv
-
[37]
M. H.-Y. Cheung, E. Berti, V. Baibhav, and R. Cotesta, Phys. Rev. D 109, 044069 (2024), arXiv:2310.04489 [gr-qc]
2024 arXiv
- [38]
- [39]
- [40]
-
[41]
Y. Yang, E. Berti, and N. Franchini, arXiv:2504.06072 [gr-qc]
-
[42]
R. K. L. Lo, L. Sabani, and V. Cardoso, arXiv:2504.00084 [gr-qc]
-
[43]
Reissner, Annalen Phys
H. Reissner, Annalen Phys. 355, 106 (1916)
1916
-
[44]
Moncrief, Phys
V. Moncrief, Phys. Rev. D 9, 2707 (1974)
1974
-
[45]
Moncrief, Phys
V. Moncrief, Phys. Rev. D 10, 1057 (1974)
1974
-
[46]
Moncrief, Phys
V. Moncrief, Phys. Rev. D 12, 1526 (1975)
1975
-
[47]
Lee and E
K.-M. Lee and E. J. Weinberg, Phys. Rev. D 44, 3159 (1991)
1991
-
[48]
Filippini and G
F. Filippini and G. Tasinato, Class. Quant. Grav. 36, 215015 (2019), arXiv:1903.02950 [gr-qc]
2019 arXiv
-
[49]
De Felice and S
A. De Felice and S. Tsujikawa, Phys. Lett. B 855, 138808 (2024), arXiv:2402.08868 [gr-qc]
2024 arXiv
-
[50]
Boskovic, R
M. Boskovic, R. Brito, V. Cardoso, T. Ikeda, and H. Witek, Phys. Rev. D 99, 035006 (2019), arXiv:1811.04945 [gr-qc]
2019 arXiv
-
[51]
Cardoso, S
V. Cardoso, S. Mukohyama, N. Oshita, and K. Takahashi, Phys. Rev. D 109, 124036 (2024), arXiv:2404.05790 [gr-qc]
2024 arXiv
- [52]
-
[53]
Chaverra, J
E. Chaverra, J. C. Degollado, C. Moreno, and O. Sarbach, Phys. Rev. D 93, 123013 (2016), arXiv:1605.04003 [gr-qc]
2016 arXiv
-
[54]
W. T. Reid, American Journal of Mathematics 53, 443 (1931)
1931
-
[55]
T. C. Sisman and B. Tekin, J. Phys. A 43, 125205 (2010), arXiv:0908.0223 [math-ph]
2010 arXiv
-
[56]
Ould El Hadj, Phys
M. Ould El Hadj, Phys. Rev. D 111, 044055 (2025), arXiv:2411.16538 [gr-qc]
2025 arXiv
-
[57]
P. G. S. Fernandes, C. A. R. Herdeiro, A. M. Pombo, E. Radu, and N. Sanchis-Gual, Phys. Rev. D 100, 084045 (2019), arXiv:1908.00037 [gr-qc]
2019 arXiv
-
[58]
Burrage, P
C. Burrage, P. G. S. Fernandes, R. Brito, and V. Cardoso, Class. Quant. Grav. 40, 205021 (2023), arXiv:2306.03662 [gr-qc]
2023 arXiv
-
[59]
Melis, F
M. Melis, F. Corelli, R. Croft, and P. Pani, Phys. Rev. D 111, 064072 (2025), arXiv:2412.14259 [gr-qc]
2025 arXiv
-
[60]
Regge and J
T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957)
1957
-
[61]
Chandrasekhar, The mathematical theory of black holes (1985)
S. Chandrasekhar, The mathematical theory of black holes (1985)
1985
-
[62]
J. Y. J. Burton and T. Osburn, Phys. Rev. D 102, 104030 (2020), arXiv:2010.12984 [gr-qc]
2020 arXiv
-
[63]
E. W. Leaver, Proc. Roy. Soc. Lond. A 402, 285 (1985)
1985
-
[64]
E. W. Leaver, Phys. Rev. D 41, 2986 (1990)
1990
-
[65]
Nomura and D
K. Nomura and D. Yoshida, Phys. Rev. D 105, 044006 (2022), arXiv:2111.06273 [gr-qc]
2022 arXiv
-
[66]
Molina, P
C. Molina, P. Pani, V. Cardoso, and L. Gualtieri, Phys. Rev. D 81, 124021 (2010), arXiv:1004.4007 [gr-qc]
2010 arXiv
-
[67]
W. D. Heiss, Phys. Rev. E 61, 929 (2000)
2000
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.