{"id":"fdbc31d7-2689-46bc-ade5-306eefa1c9da","arxiv_id":"2505.03883","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A new basis-invariant excitation factor for coupled black hole perturbation systems shows resonant amplification at avoided crossings between fundamental modes of different fields in the Einstein-Maxwell-axion theory.","lead":"This paper defines excitation factors for black hole quasinormal modes in systems where several fields are coupled, and applies the definition to an Einstein-Maxwell-axion black hole. It finds that the longest-lived modes from two different degrees of freedom repel each other and their excitation factors are amplified, a resonance-like signature that could appear in ringdown signals.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The basis-invariant excitation factor Bαn in Eq. (21) carries an extra decoupling-limit factor from other fields' reflection coefficients; its growth near the avoided crossing may be a normalization artifact rather than a physical resonance, so a time-domain amplitude check is needed.","rationale":"The paper's central claim is two-part: (i) avoided crossings occur between longest-lived QNMs originating from fundamental modes of different fields, and (ii) the excitation factors defined in Eq. (21) are amplified as a resonance. Part (i) is a well-defined spectral statement, and while the numerical evidence would benefit from convergence data, the paper's frequency-domain results are plausible. The more fragile part is (ii), because the novel basis-invariant definition of Bαn does not coincide with the physical single-field excitation amplitude even in the decoupling limit; it carries an extra factor composed of the other fields' reflection coefficients. As two QNM frequencies approach, that extra factor can grow on its own, so the observed amplification is not by itself evidence of a physical resonance. The paper acknowledges this ambiguity in Sec. II B 2 and footnote *3, where a different normalization (Eq. 23) produces qualitatively different growth. Thus the most load-bearing uncertainty is whether Bαn measures what a gravitational-wave observer would see. A time-domain calculation with a concrete source settles this directly, and the authors already identify such a calculation as necessary future work. Since the concern does not refute the spectral claims but does cast doubt on the physical interpretation of the headline resonance, the verdict should remain CONDITIONAL rather than ACCEPT or REJECT. This is a different specific concern from the reader's focus on master equations and branch labeling, hence partial agreement.","tokens_in":22149,"tokens_out":20660,"duration_ms":225327,"concrete_test":"Perform a time-domain ringdown computation for the axial EMA system: evolve a localized Gaussian pulse (or a delta-source) in the gravitational channel at Q/M = 0.1 for gaγγ = 0, 20, 28 (the avoided crossing), and 40; extract the complex amplitudes of the two longest-lived QNMs from the late-time waveform; compare the amplitude ratios against the predictions from Bg0 and Be0 including source projection. If the extracted amplitude enhancement at gaγγ ≈ 28 is absent or much weaker than |Bg0|/|Bg0| at gaγγ = 0, then Eq. (21) does not measure the physically excited amplitude. If the enhancement matches, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Equation (21) defines Bαn = i detA(out) / ((2iω)^N d/dω detA(in)) at ω = ωαn. In the decoupling limit this reduces to the standard single-field excitation factor only up to the extra factor ∏_{β≠α} A(out,β)_β / (2iω A(in,β)_β), as the paper itself notes in Sec. II B 2. That factor is built from the in/out reflection coefficients of the other fields at the complex QNM frequency. Near the avoided crossing of the gravity- and EM-led fundamental modes, ωg0 approaches ωe0, so A(in,EM)(ωg0), whose zero is the EM QNM frequency, becomes small; this alone inflates Bg0 even in a decoupled theory with a level crossing. The paper never computes an observable amplitude (e.g., the residue of the physical Green's function component in the gravitational channel), so it is open whether the growth in Fig. 5 reflects genuine resonant excitation of the coupled mode or an artifact of the determinant normalization. Footnote *3 shows that a different, channel-adapted normalization gives qualitatively different behavior, underscoring that the chosen measure is not neutral.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":22415,"tokens_out":7680,"duration_ms":80776,"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":[{"comment":"","section":"Sec. II B 2, Eq. (21)"},{"comment":"","section":"Sec. III C and Figs. 7–8"},{"comment":"","section":"Sec. III C 2 and Fig. 1"},{"comment":"","section":"Sec. III D 2, Eq. (75)"}],"minor_comments":[{"comment":"","section":"Eq. (19)"},{"comment":"","section":"Eqs. (39)–(41)"},{"comment":"","section":"Sec. III C 1"},{"comment":"","section":"Sec. III D 1"},{"comment":"","section":"Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well written and the formal development is careful, but the central physical interpretation rests on a normalization-dependent quantity. The authors appear aware of the issue (footnote *3 and the closing paragraph), yet they do not resolve it within the paper. I would ask for either a derivation showing that B_αn controls the physical waveform amplitude or a time-domain/Green's-function computation of an observable amplitude before recommending acceptance. No concerns about attribution or scholarly conduct."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new piece here is the basis-invariant excitation factor B_αn for coupled perturbation systems, Eq. (21), and the demonstration that in the Einstein-Maxwell-axion system the gravity-led and EM-led fundamental modes undergo an avoided crossing. That is a real extension of the Kerr overtone resonance story: it happens for longest-lived modes from different degrees of freedom, not just high overtones of one field. The Wronskian/Green's-matrix derivation is careful, and the authors are appropriately explicit about the tortoise-coordinate phase ambiguity and about why the alternative normalization in Eq. (23) is basis-dependent. Credit is due for that level of care.\n\nThe numerical work also looks broadly plausible: continued fraction for the coupled master system, direct integration for the excitation factors, and consistent behavior across Q/M and coupling values. The scaling relations in Figs. 6–8 are believable as a description of what the code computes.\n\nWhere I part company from their interpretation is the amplification claim. The stress-test concern is on target: B_αn contains (d/dω det A(in))^{-1}, so when two QNM frequencies approach each other, the determinant derivative naturally becomes small and B grows even if no extra physical excitation occurs. The paper itself points out that in the decoupling limit B_αn carries an extra factor built from the other fields' A(in) and A(out), and footnote 3 notes that a different normalization gives qualitatively different behavior. That means the observed growth in Fig. 5 is not yet established as a physical resonance; it could be a normalization artifact. The right check—computing the residue of a physical Green's function component in the gravitational channel, or doing a time-domain waveform with a source—is not done. The authors admit the frequency-domain limitation, but the missing step is exactly what would make the headline claim convincing.\n\nOther soft spots are minor in comparison: no code or convergence tables, and the acknowledged degradation of convergence at large coupling. The branch labeling ('gravity-led', 'EM-led') is reasonable but could become ambiguous right at the crossing, which they partially acknowledge. These do not undermine the formal framework or the existence of the avoided crossing itself.\n\nWho is this for? People working on black hole perturbation theory, modified gravity, and ringdown signatures of extra fields. It deserves a serious referee and, ultimately, publication—but the resonant amplification claim needs to be either reframed as a feature of the chosen measure or tested with a channel-resolved amplitude/time-domain calculation. If I were the editor, I would send it to review and ask for that check.","headline":"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.","tokens_in":22910,"tokens_out":1146,"would_cite":false,"duration_ms":14569,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.70.-s","04.30.-w"],"model":"deepseek-v4-flash","headline":"Coupled black hole quasinormal modes from different fields repel and amplify at avoided crossings.","keywords":["quasinormal modes","black hole perturbation theory","excitation factors","avoided crossing","resonance","Einstein-Maxwell-axion theory","ringdown","continued fraction method"],"falsifier":"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.","tokens_in":21911,"feed_emoji":"🕳️","tokens_out":8928,"duration_ms":78104,"temperature":0.7,"pith_summary":"This paper claims that when black hole perturbation fields are coupled, the longest-lived quasinormal modes of different fields can repel each other in the complex-frequency plane, and that the ease with which these modes are excited is sharply amplified at the repulsion point. It proposes a basis-invariant excitation factor $B_{\\alpha n}$ for coupled systems, Eq. (21), and tests it on an Einstein-Maxwell-axion black hole, where the axion-photon coupling mixes gravitational, electromagnetic, and axion perturbations. The observed amplification is a resonance signature analogous to the one previously seen between Kerr overtones, but here it occurs between fundamental modes of different degrees of freedom, which is a new regime. If the claim stands, ringdown signals from such theories would show a distinctive, parameter-dependent resonant enhancement, giving a concrete way to search for extra fields around black holes.","feed_headline":"Coupled black hole modes resonate at avoided crossings","feed_subtitle":"Gravity and electromagnetic ringdown channels amplify each other exactly where their frequencies repel.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the coupled master equations (48)--(50) for the EMA axial sector and the instability analysis that the paper builds on.","marker":"[42]"},{"why":"Introduces avoided-crossing QNM resonance, the two-level model, and the relation $|B_+-B_-|\\propto|\\omega_+-\\omega_-|^{-1}$ used to interpret the results.","marker":"[21]"},{"why":"Defines the single-field excitation factor whose coupled-system generalization is Eq. (21).","marker":"[22]"},{"why":"Provides the direct integration and matrix-valued continued fraction techniques used for the numerical QNM and excitation-factor computations.","marker":"[44]"},{"why":"Leaver's continued fraction method, extended to matrix form, is used to compute overtone frequencies.","marker":"[55]"},{"why":"The Green's matrix construction used in Appendix A to derive the coupled-system waveform and excitation factor.","marker":"[47]"}],"fun_headline_variants":["Avoided crossings amplify black hole ringdown","Coupled modes resonate at black hole frequency repulsion","Gravity and EM waves resonate on black holes","Excitation factors amplify at black hole avoided crossings","Black hole modes resonate when frequencies repel"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Avoided crossings amplify black hole ringdown","Coupled modes resonate at black hole frequency repulsion","Gravity and EM waves resonate on black holes","Excitation factors amplify at black hole avoided crossings","Black hole modes resonate when frequencies repel"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000991,"raw_usage":{"total_tokens":4170,"prompt_tokens":884,"completion_tokens":3286,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":500,"completion_tokens_details":{"reasoning_tokens":3217}},"tokens_in":500,"tokens_out":3286,"duration_ms":22085,"temperature":1.0,"reasoning_tokens":3217,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:42:52.237610+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Moncrief, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the direct integration and matrix-valued continued fraction techniques used for the numerical QNM and excitation-factor computations."},{"cited_title":"Lee and E","cited_arxiv_id":null,"evidence_quote":"The Green's matrix construction used in Appendix A to derive the coupled-system waveform and excitation factor."}],"review_version":1}