{"id":"0549de3f-3d6d-4d7e-b107-7e600cb33565","arxiv_id":"2510.21256","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Superradiant instability of Kerr-Newman black holes is confined to μ > qQ/M and below an analytically determined boundary that disagrees with older numerics.","lead":"Rotating, charged black holes are shown to become superradiantly unstable to a massive charged scalar field in an exactly characterized region of parameter space, with analytic formulas for the instability boundary. The paper also argues that a widely used 2004 numerical boundary is off because it relied on an approximation this exact method avoids.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central quartic (3.11) is the foundation of every exact boundary, but its altered a0 is unexplained and the derivation of μ0 from it is not shown; independent derivation is needed.","rationale":"The reader's weakest_assumption identifies precisely the imported quartic (3.11) with the unexplained a0 change and the unproved λ-independence. My reading confirms that this is the load-bearing point: every claimed exact result—μ0, the quasibound/superradiance overlap regions, and the growth-rate boundary—is derived from this quartic. The additional observation that (3.11) as written is identically zero at ω=μ makes the derivation of (4.8) even less secure and strengthens the reader's concern. I do not see a deeper or more fundamental issue: if the quartic is independently verified, the paper's analytic region mapping and the Furuhashi–Nambu comparison are genuine contributions. Therefore the correct disposition remains CONDITIONAL: the central equation must be derived or matched to [36] before the 'exact' claims can be accepted. No change to the reader's verdict is needed.","tokens_in":14486,"tokens_out":10260,"duration_ms":86616,"concrete_test":"Independently re-derive the characteristic equation by applying the confluent-Heun polynomial (Frobenius truncation) condition directly to the radial equation (2.9) using the transformations (3.5)–(3.8), keeping the angular separation constant λ explicit. Check: (i) whether the resulting a0 equals -(Mμ-Qq)^2D^2μ^2 as in (3.12a) or matches the original coefficient in [36]; (ii) whether λ cancels from the polynomial condition; and (iii) whether substituting ω=μ into the correctly derived polynomial yields (4.8), or whether an additional condition is needed. This can be done with a symbolic algebra system and would settle whether (3.11) is the correct exact spectrum and whether the derived boundaries are valid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—the exact superradiant-instability boundaries (4.14), the quasibound threshold μ0=qQ/M, and the growth-rate boundary (5.8)—all follow from the characteristic quartic (3.11). Two features make this quartic the least secure link.\n\nFirst, the paper states that the coefficients a_n coincide with those of [36] except for a0, where a0=-(Mμ-Qq)^2D^2μ^2 (3.12a). No derivation of this altered a0 is given. Since a0 is the coefficient that singles out the nontrivial quasibound threshold μ0, any error in this coefficient propagates into every region boundary and growth-rate formula.\n\nSecond, (3.11) is written as ∑_{n=0}^4 a_n(ω^n−μ^n)=0. As written, this equation is identically satisfied at ω=μ for every value of μ, so the statement 'ω=μ at the boundary ... in which case the equation reduces to (Mμ−Qq)^2D^2μ^2=0' (around (4.8)) does not follow directly. An additional argument is needed—for example, that μ0 is a double root or that the intended characteristic equation is ∑ a_n ω^n=0—but no such argument is supplied.\n\nFinally, the claimed independence of the resonance equation from the angular separation constant λ is asserted on the authority of [36] but not proved here. If λ does not cancel, the quartic is incomplete and the 'exact' region plots are not exact.\n\nThese are addressable but currently unresolved foundation issues for the paper's headline claims.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the superradiant instability of Kerr-Newman black holes coupled to a massive, charged scalar field. Using the Vieira-Bezerra-Kokkotas (VBK) method, it takes the quartic characteristic equation (3.11) for the quasibound frequency ω as its starting point and derives analytic conditions for the existence of quasibound and superradiant modes. The central claims are: (i) quasibound states exist only for μ > μ0 = qQ/M; (ii) the superradiant instability region is the overlap of the quasibound condition with the superradiance condition μ < ω_c; (iii) the boundary where Im ω changes sign is given by the exact expression (5.8), with perturbative approximations (5.5)-(5.6); and (iv) the numerical boundary of Furuhashi and Nambu [45] is inaccurate because it uses λ = ℓ(ℓ+1). The paper includes region plots in the (a,Q) plane and comparison plots with [45].","tokens_in":14891,"tokens_out":5395,"duration_ms":52613,"significance":"If the central quartic (3.11) is correct, the paper would be a substantial contribution: it provides parameter-free analytic boundaries of the superradiant-instability region without the hydrogenic approximation, makes falsifiable predictions about the parameter space, and offers a concrete explanation for a known numerical discrepancy. The algebraic reductions are mostly transparent, and the figures clearly illustrate the claimed regions and boundaries. However, the entire analysis rests on a characteristic equation that is imported from [36] and modified without derivation, so the significance is conditional on the correctness and provenance of that equation.","major_comments":[{"comment":"The characteristic quartic is the foundation of every subsequent result. The coefficients are said to coincide with those of [36] except for a0, which is changed to a0 = -(Mμ-Qq)^2 D^2 μ^2 without derivation. Since this coefficient is precisely what generates the quasibound threshold μ0 = qQ/M in Eq. (4.9), an unverified or incorrect a0 would invalidate all central claims. Please provide a derivation of (3.11) from the confluent-Heun polynomial condition, or at least a direct verification against [36] explaining why a0 differs from the published expression.","section":"§3, Eqs. (3.11)-(3.12a)"},{"comment":"Equation (3.11) is written as Σ_{n=0}^4 a_n(ω^n - μ^n) = 0. At ω = μ, every term vanishes, so the equation is identically satisfied for all μ. The paper's statement that at the quasibound boundary ω = μ and 'the equation reduces to (Mμ-Qq)^2 D^2 μ^2 = 0' does not follow directly from (3.11). An additional argument is needed, e.g. that the physical mode corresponds to a double root at the boundary or that the original resonance condition is actually Σ a_n ω^n = 0 before rearrangement. Without such an argument, the derivation of μ0 and all subsequent region boundaries is not logically established.","section":"§4.1, Eq. (4.8)"},{"comment":"The instability regions in Fig. 3 are plotted using the inequality μ < ω_c with ω_c = (ma + qQ r_+)/(r_+^2 + a^2), which follows from (4.13) after replacing ω by μ via ω ≲ μ. This is an approximation, not an exact condition, unless Re ω = μ is proven for the relevant modes. Moreover, Section 5 shows that the boundary where Im ω changes sign is μ1, not simply ω_c, so the overlap plotted in Fig. 3b may not be the exact superradiant-instability region claimed in the title and abstract. Please either justify exactness or label Fig. 3 as the approximate/superradiance-region plot with a clear caveat.","section":"§4.2, Eq. (4.13) and Fig. 3"},{"comment":"The exact expression for μ1, the boundary where Im ω = 0, is stated as 'Surprisingly, it is possible to solve for μ1 exactly' but no derivation is given. Since Eq. (5.8) is a headline exact result and is used to benchmark the numerical discrepancy with [45], the reader cannot verify it from (3.11) as presented. Please supply a derivation or a detailed outline, including the branch choices for the square roots.","section":"§5, Eq. (5.8)"},{"comment":"The paper asserts that the characteristic resonance equation is independent of the separation constant λ, but this is not proved here; it is borrowed from [36]. This λ-independence is load-bearing for the argument that the numerics of [45] are inaccurate because they set λ = ℓ(ℓ+1). If λ does not cancel from the resonance condition, the quartic (3.11) is incomplete and the discrepancy explanation collapses. Please include a proof or a more explicit derivation of the λ-independence, rather than relying solely on a statement from [36].","section":"§5.1, last paragraph"}],"minor_comments":[{"comment":"The name 'Viera' should be 'Vieira' (also in the reference to the VBK method).","section":"§2, after Eq. (2.10)"},{"comment":"The near-horizon radial behavior is written as C1 (r - r_+) e^{-i(ω-ω_c)/(2κ_+)}, which appears to be missing the power of (r - r_+) in the exponential. It should presumably be (r - r_+)^{-i(ω-ω_c)/(2κ_+)}.","section":"Eq. (4.1)"},{"comment":"The Furuhashi-Nambu boundary in the figure is described as plotted 'by eye' as approximately linear. Since the comparison with [45] is a central quantitative claim, please provide the actual numerical data points or a reproducible extraction method, rather than an eyeballed line.","section":"Fig. 5"},{"comment":"Reference [48] contains a typo: 'Phys Lett B' should be 'Phys. Lett. B'.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper's main risk is not circularity but reliance on an imported, modified characteristic quartic and a logical gap at Eq. (4.8). These issues are potentially fixable if the authors can supply the missing derivations. If they cannot, the exactness claims should be substantially weakened. I recommend major revision rather than rejection, because the project is valuable and the identified gaps appear addressable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere's the quick take: this paper gives clean analytic maps of the superradiant-instability region for Kerr-Newman black holes with a massive charged scalar, using the Vieira-Bezerra-Kokkotas quartic. The genuinely new things are the threshold μ0=qQ/M for quasibound states, the a-Q overlap diagrams (fig 3b), and the exact boundary (5.8) for the growth rate. The face-off with Furuhashi-Nambu is a real contribution: they argue that the numerics use the hydrogenic approximation λ=ℓ(ℓ+1), which is not the true separation constant, and that many later papers benchmark against that numerical boundary. That point deserves to be heard.\n\nBut there is a load-bearing problem. Everything rests on equation (3.11), imported from Vieira-Bezerra-Muniz with a0 changed and no derivation given. As written, (3.11) is ∑ a_n(ω^n−μ^n)=0, which is identically satisfied at ω=μ for any μ. The step around (4.8) — that at ω=μ the equation reduces to (Mμ−Qq)^2D^2μ^2=0 — simply does not follow from (3.11). You need an extra argument: the intended equation is presumably ∑ a_n ω^n=0, and at the boundary you get ∑ a_n μ^n=0, and then show that reduces to a0=0. That gap is not cosmetic; it supports μ0 and every boundary derived from it. The claim that the quartic is independent of the separation constant λ is also taken on faith from [36]. If λ doesn't cancel, the equation is incomplete and the 'exact' regions are wrong.\n\nTwo more soft spots, in proportion: fig. 3b uses the approximation μ≈ω for the superradiance cut despite the paper's 'exact' framing — they mention it, but it means the regions aren't exactly what they say. And (5.8), the exact boundary, is quoted without derivation. The sign of Imω for qQ<0 isn't discussed either, though the plot shows those cases.\n\nIf the foundations hold, this is a solid, useful contribution. If the quartic derivation is wrong, everything falls. The fix is straightforward: supply the derivation of a0 and μ0 from an unambiguous characteristic equation, and prove or cite a proof for the λ-independence before the 'exact' language is used. I would send this to a competent referee with exactly those questions. I wouldn't cite it in my own work until that's sorted, but I do think it deserves referee time rather than a desk reject.","headline":"A useful analytic map of KN superradiant instability, but the central quartic is borrowed and the derivation of the threshold does not follow as written — referee it, but require the derivation.","tokens_in":15349,"tokens_out":6240,"would_cite":false,"duration_ms":51951,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57"],"pacs":["04.70.Bw"],"model":"deepseek-v4-flash","headline":"This paper claims that superradiant instability of a Kerr-Newman black hole with a massive charged scalar field is exactly the overlap of the quasibound-state condition μ > qQ/M and the superradiant condition, with the instability boundary","keywords":["superradiance","Kerr-Newman black hole","massive charged scalar field","quasibound states","confluent Heun equation","characteristic quartic","growth rate boundary","black hole stability"],"falsifier":"Integrate the full separated Klein-Gordon radial equation (2.9) numerically, using the exact spheroidal eigenvalue λ obtained by solving the angular equation (2.8), and scan across the predicted boundary (5.8); if the sign change of Imω occurs at a value of μM that differs from (5.8), or if the boundary shifts with ℓ, the central claim fails. A cheaper check is to re-derive the quartic from the confluent Heun polynomial condition and verify the coefficient a0 = −(Mμ−Qq)²D²μ².","tokens_in":14402,"feed_emoji":"🕳️","tokens_out":6591,"duration_ms":58650,"temperature":0.7,"pith_summary":"The paper sets out to characterize, exactly and without approximation, where a Kerr-Newman black hole becomes unstable to superradiance when surrounded by a massive charged scalar field. Using the exact polynomial (quartic) resonance condition for the field modes, the authors show that quasibound states exist only when the scalar mass exceeds μ0 = qQ/M, and that the instability region is precisely the overlap of this condition with the superradiant frequency condition. They derive first- and second-order perturbative formulas and an exact closed-form expression for the boundary line where the mode growth rate changes sign. The sharp consequence is that the instability region moves from nearly neutral, Kerr-like black holes for light fields to highly charged, near-extremal Kerr-Newman black holes for heavier fields, with no instability in the Reissner-Nordström limit.","feed_headline":"Exact formula fixes superradiant instability region for Kerr-Newman","feed_subtitle":"A quartic resonance equation replaces the hydrogenic approximation and corrects previous numerical boundaries.","key_machinery":"The carrying object is the characteristic resonance equation (3.11), a quartic polynomial in the mode frequency ω obtained by demanding polynomial (confluent Heun) solutions of the radial Klein-Gordon equation in the Kerr-Newman background. From this quartic the paper derives the quasibound threshold μ0 = qQ/M, the sign of Re√(μ²−ω²), and the closed-form boundary formulas (5.5), (5.6), (5.8) for the growth-rate cutoff. The quartic's noted independence from the angular separation constant λ is what lets the region plots be drawn without first solving the angular equation.","core_discovery":"The central discovery is that the superradiant-instability boundary can be stated exactly: quasibound states require μ > μ0 = qQ/M, superradiance requires μ < ωc, and the line along which the growth rate Im(Mω) turns negative is given by the exact closed-form expression (5.8), together with its perturbative approximations (5.5) and (5.6). The paper obtains these by working directly with the roots of the characteristic quartic resonance equation (3.11), rather than the hydrogenic approximation. It also shows that a widely used numerical boundary of earlier work is shifted, because that work fixed the separation constant at the hydrogenic value rather than letting the exact (λ-independent) res","pith_inferences":["Because the resonance equation is claimed to be independent of λ, the paper implicitly predicts that the instability boundary is identical for all angular quantum numbers ℓ, not just the ℓ=1 mode; an independent numerical solution of the full separated system (including the spheroidal eigenvalue) could test this directly.","The exact boundary (5.8) offers a parameter-free benchmark: a time-domain simulation of a charged scalar field on a Kerr-Newman background near the predicted line should show the onset of exponential growth at the analytic value of μM.","The same quartic machinery could be carried over to dyonic (magnetically charged) black holes, where μ0 would generalize to include a magnetic charge product; if the quartic structure survives, the instability region would again be an overlap of two inequalities.","For ultralight scalar fields (μM ≪ 1) the threshold μ0 = qQ/M confines instability to black holes with tiny charge Q < μ, so axion-like models remain effectively Kerr-like; the paper's exact formulas make this separation quantitative."],"forward_implications":["The instability region in black-hole parameter space becomes a pair of inequalities: μ > qQ/M for quasibound states and μ < ωc for superradiance, with the growth-rate cutoff in closed form (5.8).","Region plots such as fig. 3b give the exact (a,Q) unstable domain for any scalar mass μ, replacing numerical scans.","Existing numerical results used as benchmarks for the instability boundary should be re-evaluated; the paper's exact boundary shifts the point where the growth rate turns negative.","In the Reissner-Nordström limit the instability region collapses to points, reaffirming the known stability of charged, non-rotating black holes against this mechanism.","The method extends to other asymptotically flat black hole families and higher-spin fields without the hydrogenic limit."],"fun_headline_variants":["Exact quartic roots pin down Kerr-Newman superradiance","Superradiant instability boundary now exact, not hydrogenic","Exact formula corrects Kerr-Newman instability map","Instability region shifts when separation constant is exact","Exact expression fixes superradiant boundary for Kerr-Newman"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The whole derivation rests on the characteristic quartic (3.11), imported with a modified leading coefficient a0 and with the asserted λ-independence; if that quartic is not the exact resonance condition, every boundary and region plot in the paper inherits the error.","fun_headline_variants_meta":{"raw":{"variants":["Exact quartic roots pin down Kerr-Newman superradiance","Superradiant instability boundary now exact, not hydrogenic","Exact formula corrects Kerr-Newman instability map","Instability region shifts when separation constant is exact","Exact expression fixes superradiant boundary for Kerr-Newman"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000196,"raw_usage":{"total_tokens":1175,"prompt_tokens":701,"completion_tokens":474,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":445,"completion_tokens_details":{"reasoning_tokens":390}},"tokens_in":445,"tokens_out":474,"duration_ms":4895,"temperature":1.0,"reasoning_tokens":390,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T08:18:40.299766+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Integrate the full separated Klein-Gordon radial equation (2.9) numerically, using the exact spheroidal eigenvalue λ obtained by solving the angular equation (2.8), and scan across the predicted boundary (5.8); if the sign change of Imω occurs at a value of μM that differs from (5.8), or if the boundary shifts with ℓ, the central claim fails. A cheaper check is to re-derive the quartic from the confluent Heun polynomial condition and verify the coefficient a0 = −(Mμ−Qq)²D²μ².","supporting_citations":[],"review_version":1}