{"id":"3d1e1731-b130-490d-a4e9-598a704c2358","arxiv_id":"1908.06117","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Charged scalar perturbations of higher-dimensional Reissner-Nordström-de Sitter black holes are superradiantly unstable, with the growth rate increasing with spacetime dimension.","lead":"This paper studies how charged scalar fields can grow unstably around higher-dimensional Reissner-Nordström-de Sitter black holes, extracting energy through superradiance. It shows the instability becomes stronger as spacetime dimension increases and traces it to a family of modes tied to the cosmological horizon.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim rests on an unpublished spectral solver with no shown convergence data; the reported unstable l=0 branch could be a numerical artifact, so the instability claim is not yet quantitatively verified.","rationale":"The reader's weakest assumption was exactly the reliability of the numerical frequency-domain solver, and my stress-test pass finds this to be the most load-bearing point. The paper's central claim is a numerical one: positive imaginary parts of specific quasinormal frequencies. All other elements—the master equation, the superradiant scattering relation, the existence of dS modes, and the qualitative potential-well argument—are standard or analytically supported. The single point where the conclusion can break is the numerical solver: if the unstable l=0 branch is misidentified or contaminated, the entire claim of superradiant instability in higher dimensions collapses. The paper does cite cross-checks with WKB, a matrix method, and time-evolution data, which are real mitigations and deserve credit; however, none of these checks are shown, so the reader cannot verify that the reported positive Im(omega) values are physical rather than numerical noise. The small magnitude of the growth rates makes this concern concrete rather than hypothetical. I therefore agree with the reader's assessment and see no reason to change the CONDITIONAL verdict: the paper is likely correct but should be accepted only after the numerical uncertainties are quantified and the verification data are made available.","tokens_in":12966,"tokens_out":10141,"duration_ms":110702,"concrete_test":"For the headline unstable configuration in each dimension (e.g., d=4 with Lambda=0.005, qQ=0.5; d=5 with Lambda=0.25, qQ=3; d=6 with Lambda=1, qQ=7), recompute the l=0 QNM frequency with the package of Ref. [54] at Chebyshev resolutions N=30, 60, 120 and with the independent matrix method of Ref. [60], and also run a time-domain evolution in the same background. Require that Re(omega) and Im(omega) converge to at least 1% and that Im(omega) remains positive across all three methods. Additionally, recompute the two mu=0.4 entries in Table II to resolve the duplicate-row ambiguity. If any method gives Im(omega) <= 0 or the values do not converge, the central instability claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is the existence and dimension-dependence of a charged-scalar superradiant instability in d=4,5,6 RNdS. The evidence is a frequency-domain calculation using the Mathematica package of Ref. [54]. The text says the results were checked with WKB and a matrix method, and acknowledgments mention time-evolution data, but no convergence study, error bars, comparison plots, code, or data are presented. This matters because the unstable modes are obtained by continuously tracking the l=0 pure-dS 'zero-mode' as qQ is increased: exactly the regime where a spectral solver can jump branches or converge to a spurious mode if the grid resolution or initial guess is inadequate. The reported growth rates are small (Im(omega) of order 10^-4 to 10^-3 in Tables I and II), so even modest numerical error could change the sign of Im(omega) and invert the stability conclusion. An apparent internal inconsistency in Table II, where the same row mu=0.4 is listed twice with different frequencies, further suggests the numerical output was not fully vetted. Without quantitative error control, the paper's strongest quantitative statements—that the modes are unstable and that higher dimensions amplify the growth rate—are not established at the claimed precision.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies charged scalar-field perturbations on higher-dimensional Reissner-Nordström–de Sitter (RNdS) black holes, specializing to spacetime dimensions d=4,5,6. It derives a radial master equation, obtains a Wronskian-based criterion for superradiance, and computes quasinormal mode (QNM) frequencies with a frequency-domain spectral package. The central claims are that l=0 charged scalar perturbations become superradiantly unstable in a finite region of the charge-coupling parameter space, that this instability is tied to the pure de Sitter 'zero-mode' family of QNMs, that adding a scalar-field mass stabilizes the system, and that increasing the spacetime dimension amplifies the instability and enlarges the unstable region. The paper reports unstable and stable modes in Tables I and II and illustrates the behavior in Figures 1–5.","tokens_in":13197,"tokens_out":4708,"duration_ms":49428,"significance":"If the numerical results are reliable, the paper is a useful extension of the known four-dimensional superradiant instability to higher-dimensional RNdS spacetimes, and the connection to the pure de Sitter QNM family provides a physically interesting interpretation. The derivation of the master equation and the Wronskian superradiance condition is clean and appears correct, and the explicit tables give concrete frequencies that can be checked by other methods. The paper also makes the falsifiable prediction that all unstable modes satisfy the superradiant condition, while some stable modes satisfy it as well. The main weakness is that the central quantitative claims rest on a spectral solver with no reported convergence study, error bounds, or shown cross-checks, which is especially concerning because the unstable imaginary parts are small and could be affected by numerical error.","major_comments":[{"comment":"The central instability claim is not yet quantitatively verified because no convergence study or error estimates are reported for the spectral solver. The unstable modes have imaginary parts of order 10^-4 to 10^-3, so even modest numerical error could change the sign of Im(ω) and alter the stability conclusion. The manuscript states that the results were checked with WKB and a matrix method, and the acknowledgments mention time-evolution data, but none of these checks are shown quantitatively for the parameter points used. Since the unstable branch is obtained by tracking the l=0 pure de Sitter zero-mode as qQ increases, the possibility of branch jumping or a spurious numerical mode is not excluded. Please provide a convergence test (e.g., frequencies versus grid size/order), show direct comparisons with at least one independent method for representative rows of Tables I and II, and report error bars or an equivalent accuracy estimate.","section":"Sec. V, Tables I–II, Figs. 3–5"},{"comment":"The same mass value μ=4×10^-1 is listed twice in the d=6 block with two different frequencies, ω=0.04496−0.02141i and ω=0.03005−0.10116i. This internal inconsistency suggests that the numerical output was not fully vetted and makes the table unreliable as a reference. Please correct the typo or explain which mass value was actually used, and re-check the corresponding entries.","section":"Table II, d=6 block"},{"comment":"The claim that 'the increment of dimensions amplifies the growth rate and enlarges the region of the parameter space' is not established by the data as presented. The d=4, 5, and 6 panels in Figures 3 and 5 use different cosmological constants (Λ=0.005, 0.25, 1) and different qQ ranges, so the larger qQ range at higher d may simply reflect the larger available parameter space rather than a dynamical amplification. The peak imaginary parts in Table I are not monotonic in d under the chosen parameter sets (d=4: ~0.0016 at qQ=0.5; d=5: ~0.0007 at qQ=1; d=6: ~0.0025 at qQ=10). To support the dimension-dependence claim, the comparison should be made at a controlled physical scale, for example fixed dimensionless quantities such as ΛM^2 or fixed horizon-radius ratios, with equivalent resolution in all dimensions.","section":"Abstract, Sec. V, Figs. 3 and 5"}],"minor_comments":[{"comment":"The notation qQ is used as a dimensionless coupling, but the relation between the scalar-field charge q, the black-hole charge parameter q0, and the physical charge Q is only implicit; a short definition of the dimensionless combination would improve clarity.","section":"Sec. II, Eq. (8)"},{"comment":"The quantity Qmax appearing in the caption is never defined; please define it explicitly, for example via the extremal charge for the given M and Λ.","section":"Fig. 3 caption and Sec. V"},{"comment":"The potential-well argument is used as evidence for instability, but the existence of a potential well is only heuristic and does not by itself prove the QNM spectral property; the text should state more clearly that the instability is established by the frequency-domain calculation rather than by the shape of V(r).","section":"Sec. IV, Fig. 2"},{"comment":"Several typos should be corrected, including 'in expense of' (should be 'at the expense of'), 'Quiet strikingly' (should be 'Quite strikingly'), and 'a d-dimensional' (should be 'a d-dimensional' with the article adjusted).","section":"Introduction and Acknowledgments"},{"comment":"The table footnote says all modes satisfy the superradiant condition, which is a useful check, but it would help to also list the values of ω−Φ(rc) and Φ(r+)−ω explicitly so the reader can verify the inequalities without recomputation.","section":"Sec. V, Table I"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the analytical setup is sound. The main risk is numerical reliability: the central sign-changing result relies on a spectral solver with no convergence data, and the duplicate row in Table II indicates that the numerical output was not fully checked. I would request a revised version with a convergence study and explicit cross-checks before publication; if those checks fail, the instability claim would need to be reconsidered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a sound, incremental extension of the known 4D superradiant instability of charged scalars in RNdS to d=5 and d=6, with a clean derivation of the superradiant condition and a sensible interpretation linking the unstable branch to the pure de Sitter mode family. The main quantitative claim—that higher dimensions amplify the instability—is physically credible but rests on a numerical solver whose convergence is not demonstrated.\n\nWhat the paper does well: the master equation, Wronskian argument, and boundary conditions are all handled correctly. The connection to the zero-mode of the BH-dS family is a nice conceptual step, and the paper honestly notes that superradiance is necessary but not sufficient. The tables show the mode frequencies satisfy the threshold condition, which is a genuine check.\n\nThe soft spots are mostly about reproducibility and numerical hygiene. The results come from an unpublished Mathematica package; the WKB and matrix-method checks are only mentioned, not shown; no error bars or convergence data are provided. The stress-test note is right that the unstable branch is exactly where a spectral solver can jump or return spurious modes. The growth rates are small (Im omega ~ 1e-4 to 1e-3), so modest numerical error could change the sign of the imaginary part. The acknowledgment that time-evolution data confirmed the spectral results is helpful, but those data are not included, so the paper as written does not deliver the evidence one would want for a quantitative claim.\n\nThere is also an apparent typo in Table II: for d=6, Lambda=1, qQ=7, mu=0.4 appears twice with different frequencies. That is minor and doesn't affect the central argument, but it suggests the tables weren't fully proofread.\n\nI think the central claim is probably right—it is a natural extension of an established result, and the effective-potential argument supports it—but the paper as written doesn't yet vindicate the claimed quantitative amplification.\n\nWho is this for: people working on black hole superradiance or dS mode families. It deserves peer review, but a serious referee should ask for convergence studies, the actual WKB/matrix comparisons, and a corrected Table II.","headline":"A plausible and cleanly derived extension of the 4D charged-scalar superradiant instability to d=5,6, but the quantitative claim rests on a spectral solver with no shown convergence data.","tokens_in":13716,"tokens_out":2089,"would_cite":false,"duration_ms":20543,"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.70.Bw"],"model":"deepseek-v4-flash","headline":"Charged scalar fields superradiantly destabilize Reissner-Nordström-de Sitter black holes in four to six dimensions, with higher dimensions amplifying the growth and widening the unstable parameter region.","keywords":["superradiance","charged scalar fields","quasinormal modes","Reissner-Nordström-de Sitter black holes","higher-dimensional black holes","de Sitter spacetime","black hole instability","cosmological horizon"],"falsifier":"Run an independent time-domain integration of the charged Klein-Gordon equation on a fixed $d=5$ RNdS background with $M=1$, $Q=0.5$, $\\Lambda=0.25$, $qQ=1$, $\\mu=0$: the reported dominant mode has $\\omega \\approx 0.01656 + 0.00073 i$, so the field must grow. If it instead decays, or if a second independent frequency-domain method does not reproduce the unstable branch under grid refinement, the central claim fails.","tokens_in":12772,"feed_emoji":"🕳️","tokens_out":10718,"duration_ms":89715,"temperature":0.7,"pith_summary":"The paper claims that charged scalar fields on Reissner-Nordström–de Sitter (RNdS) black holes in four, five, and six spacetime dimensions are superradiantly unstable: the $\\ell=0$ mode grows exponentially in time, extracting energy from the black hole's electromagnetic field. The growing mode belongs to the 'BH de Sitter' family of quasinormal modes, seeded by the purely imaginary zero-mode of pure de Sitter space, rather than to the usual photon-sphere family. Higher dimensions deepen the potential well that traps the field, which both amplifies the growth rate and enlarges the region of the $qQ$ and $\\Lambda$ parameter space where instability occurs. Adding a mass $\\mu$ to the scalar field shifts the modes downward and eventually restores stability, although unstable massive modes still satisfy the superradiant condition. A fair reading of the evidence is that superradiance is necessary but not sufficient for the instability.","feed_headline":"Charged fields grow faster around higher-dimensional black holes","feed_subtitle":"Deeper trapping wells make the superradiant instability set in faster and over a wider parameter range.","key_machinery":"The load-bearing object is the 'BH de Sitter' family of quasinormal modes: modes of asymptotically de Sitter black holes whose frequencies are well approximated by the pure de Sitter spectrum $\\omega/\\kappa_c^{\\mathrm{dS}} = -i(\\ell+2n)$ and $\\omega/\\kappa_c^{\\mathrm{dS}} = -i(\\ell+2n+d-1)$. For $\\ell=0$ this family contains a zero-mode, and turning on the scalar field charge $qQ$ deforms it into a growing mode. The mechanism that allows the growth is the effective potential $V(r)$ in the radial master equation: for $\\ell=0$ it develops a potential well between the photon sphere and the cosmological horizon, trapping waves where the superradiant amplification condition $\\Phi(r_c)<\\omega<\\Phi(r_+)$ can act; higher dimensions make this well deeper.","core_discovery":"The central discovery is that in $d=4,5,6$ non-extremal RNdS spacetimes, the $\\ell=0$ charged scalar mode — the one that becomes the zero-mode $\\omega=0$ of pure de Sitter space when the charge coupling vanishes — turns into a complex quasinormal frequency with $\\mathrm{Im}(\\omega)>0$ as soon as $qQ$ is switched on. Its real part grows monotonically with $qQ$; its imaginary part first rises, peaks, then falls below zero, so stability is restored for sufficiently large charge coupling or scalar mass. Every unstable mode found lies inside the superradiant window $\\Phi(r_c)<\\omega<\\Phi(r_+)$, but so do many stable modes, which the paper reads as proof that superradiance is necessary but not sufficient. As the spacetime dimension increases, the effective potential well deepens, making the instability stronger and extending it across a wider region of parameter space.","pith_inferences":["If the de Sitter zero-mode is the seed, the same mechanism should operate for charged fermion or vector perturbations on RNdS, and possibly for rotating de Sitter black holes whenever a de Sitter mode lies inside the superradiant frequency window; this is a testable extension, not a claim of the paper.","The dimension dependence visible in the tables suggests a scaling law for the peak growth rate and the critical couplings in terms of horizon radii or the cosmological surface gravity; extracting it would give an analytic estimate beyond the three dimensions computed here.","Because the instability strengthens with dimension, nonlinear evolutions in $d=5$ or $d=6$, where the growth timescale is shorter, may be the practical route to determine whether the endpoint is a scalarized black hole or evacuation of the field — a question the paper leaves open."],"forward_implications":["For $d=4,5,6$ Reissner-Nordström-de Sitter black holes, $\\ell=0$ charged scalar perturbations grow exponentially for intermediate values of the charge coupling $qQ$, extracting electromagnetic energy from the black hole.","The growth rate peaks at a finite $qQ$ and then falls; increasing the scalar mass $\\mu$ lowers $\\mathrm{Im}(\\omega)$, so a finite critical mass exists beyond which the instability disappears.","Because the instability rides on the de Sitter zero-mode rather than the photon-sphere family, any asymptotically de Sitter black hole supporting a zero-mode is a candidate for the same superradiant instability.","Higher spacetime dimensions amplify the effect: as $d$ goes from 4 to 6, the unstable region of the $(qQ,\\Lambda)$ parameter space widens and the peak growth rate grows.","The superradiant condition $\\Phi(r_c)<\\omega<\\Phi(r_+)$ is necessary but not sufficient, so detecting a mode inside this window does not by itself prove instability."],"supporting_citations":[{"why":"grounds the original 4D discovery that l=0 charged scalar perturbations of RNdS are superradiantly unstable, the result this paper extends to higher dimensions.","marker":"[39]"},{"why":"establishes that the superradiant condition is necessary but not sufficient in 4D RNdS, the interpretive scheme the paper carries into d=5,6.","marker":"[40]"},{"why":"supplies the frequency-domain numerical method used to compute the quasinormal frequencies reported in the paper.","marker":"[54]"},{"why":"identifies the BH dS family of quasinormal modes and its approximation by pure de Sitter modes, which the paper links to the instability.","marker":"[55]"},{"why":"extends the BH dS family to higher-dimensional scalar perturbations, providing the higher-d mode structure the paper builds on.","marker":"[57]"},{"why":"gives the pure dS quasinormal spectrum used to approximate the BH dS modes and locate the zero-mode.","marker":"[70, 71]"},{"why":"provides the analytic result for 4D RNdS that the paper's Lambda-dependence of the instability matches.","marker":"[74]"}],"fun_headline_variants":["Higher dimensions intensify black-hole superradiance","Charge-driven instability scales with dimension","Superradiant growth gets a boost in higher dimensions","Higher-D black holes speed up charged-field blowup","Dimension hikes destabilize charged scalar fields"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the numerical frequency-domain solver returns accurate quasinormal frequencies for the chosen parameter ranges; the paper reports no convergence study or error bounds, so a spurious mode misidentified as the dominant $\\ell=0$ branch would invalidate the instability claim.","fun_headline_variants_meta":{"raw":{"variants":["Higher dimensions intensify black-hole superradiance","Charge-driven instability scales with dimension","Superradiant growth gets a boost in higher dimensions","Higher-D black holes speed up charged-field blowup","Dimension hikes destabilize charged scalar fields"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000219,"raw_usage":{"total_tokens":1426,"prompt_tokens":909,"completion_tokens":517,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":525,"completion_tokens_details":{"reasoning_tokens":450}},"tokens_in":525,"tokens_out":517,"duration_ms":5787,"temperature":1.0,"reasoning_tokens":450,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:54:39.221189+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run an independent time-domain integration of the charged Klein-Gordon equation on a fixed $d=5$ RNdS background with $M=1$, $Q=0.5$, $\\Lambda=0.25$, $qQ=1$, $\\mu=0$: the reported dominant mode has $\\omega \\approx 0.01656 + 0.00073 i$, so the field must grow. If it instead decays, or if a second independent frequency-domain method does not reproduce the unstable branch under grid refinement, the central claim fails.","supporting_citations":[],"review_version":1}