{"id":"8a32d710-b446-40d1-9a4d-2d436f6aacc4","arxiv_id":"2506.09829","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Purely imaginary gravitational quasinormal modes of Reissner-Nordström-de Sitter black holes are computed and shown to obey a universal small-hole formula and to generate late-time exponential tails.","lead":"A new calculation finds the non-oscillatory vibration modes of the gravitational field around charged black holes in an expanding universe. These modes, inherited from empty de Sitter space, decay slowly and are responsible for the late-time exponential tails seen in simulations of such black holes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper asserts, but never demonstrates, that the computed purely imaginary modes produce the exponential asymptotic tails; Eq. (23), as written, also cannot confirm the universal small-black-hole law.","rationale":"The reader's verdict is CONDITIONAL, and this stress test agrees that the paper should not be fully accepted without additional support. The reader's weakest assumption is numerical reliability; that concern is valid, but the more directly load-bearing gap is the missing demonstration of the tail-attribution claim, which is part of the abstract's central claim. A strong paper could settle this with a single time-domain run and a comparison to the least-damped frequency, so the concern is resolvable and does not warrant rejection. I also note a concrete internal inconsistency in Eq. (23): as typeset it omits the r_h/(2 r_c) term and uses a vague proportionality, so it cannot by itself confirm Eq. (22). This is likely a typographical or presentation error rather than a fatal flaw, because the following sentence states the correct form and Ref. [91] already provides the analytic law. Overall, the physical conclusion is plausible and consistent with the existing scalar-field and Schwarzschild-de Sitter literature, which is why the verdict should remain CONDITIONAL rather than being raised to REJECT.","tokens_in":11737,"tokens_out":9635,"duration_ms":124537,"concrete_test":"Perform one time-domain evolution for a small charged black hole, e.g., r_i = r_h/2, r_c = 10 r_h, l = 2, '-' perturbation: integrate the master equation with the Gundlach-Price-Pullin scheme on two different grid spacings, fit the late-time signal to A exp(-kappa t), and compare kappa with the continued-fraction value of -Im(omega) for the least-damped purely imaginary mode. If the two values do not agree within the estimated discretization error, the tail-attribution claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim has three parts: existence of a purely imaginary de Sitter branch in the coupled gravitational-electromagnetic spectrum, validity of the universal formula (22), and responsibility of these modes for the late-time exponential tails. The first part is credible because it extends known test-field and Schwarzschild-de Sitter results, and the second is supported by an independent analytic result cited as Ref. [91]. The load-bearing weak point is the third: Section III B describes the Gundlach-Price-Pullin time-domain integrator, but the paper presents no time-domain profile, no Prony extraction, and no comparison of an extracted decay rate with Im(omega) of the least-damped purely imaginary mode. The abstract's statement that these modes 'are responsible for the exponential asymptotic tails' is therefore an assertion without displayed evidence. In addition, the numerical confirmation of the universal law is stated through Eq. (23) as omega_n proportional to 1 - r_i/(2 r_c), which omits the r_h/(2 r_c) term appearing in the very next equality and gives no shift for uncharged black holes (r_i = 0), contradicting Fig. 3 and Eq. (22). The absence of truncation orders and convergence checks compounds the problem, since the continued-fraction calculation for purely imaginary modes is exactly the regime where Nollert acceleration is needed and spurious roots are possible. None of this proves the results are wrong, but it means the central claim as written is not substantiated by the material in the paper.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the quasinormal spectrum of gravitational and electromagnetic perturbations of the Reissner-Nordström-de Sitter black hole, focusing on the purely imaginary (non-oscillatory) de Sitter branch. Using the Leaver continued-fraction method with Nollert acceleration and a time-domain Gundlach-Price-Pullin integrator, the authors report that this branch exists for coupled gravitational-electromagnetic perturbations, reduces to empty de Sitter space modes in the small-black-hole limit, and obeys the approximate formula ω_n ∝ 1 - r_i/(2 r_c) for r_h ≪ r_c, which they interpret as confirming the universal law ω_n = ω_n^{(dS)}(1 - M/r_c + O((M/r_c)^2)). They further claim that these modes are responsible for the late-time exponential tails of the perturbation. The paper also describes the interplay between this de Sitter branch and the complex black-hole branch, including the overtaking of the least-damped mode as the black hole grows.","tokens_in":12017,"tokens_out":3736,"duration_ms":46675,"significance":"If the central claims are correct, the paper fills a gap in the quasinormal-mode literature by extending the purely imaginary de Sitter branch, previously known for test fields and for Schwarzschild-de Sitter, to the coupled gravitational-electromagnetic sector of Reissner-Nordström-de Sitter. The result that these modes control late-time exponential tails and that they satisfy the universal small-black-hole formula would be physically relevant for late-time gravitational-wave behavior and for strong cosmic censorship considerations. The authors make appropriate use of standard, well-established numerical tools (Leaver method, Nollert convergence acceleration, and time-domain integration), and the qualitative picture they present is plausible. However, the manuscript as submitted lacks the numerical tables, convergence checks, and time-domain extraction needed to substantiate these quantitative claims, and the approximate formula Eq. (23) contains an internal inconsistency. For these reasons the paper cannot yet be accepted in its present form.","major_comments":[{"comment":"The manuscript provides no numerical tables, no truncation orders, no convergence checks, and no error estimates for the Leaver/Nollert continued-fraction calculations or for the time-domain integration. This is load-bearing because the central existence claim and the quantitative agreement with Eq. (22) rest entirely on these numerical results, and purely imaginary modes are precisely the regime where continued-fraction solutions can be contaminated by spurious roots unless the recurrence reduction and truncation are carefully controlled. Please include representative numerical values of ω_n for a few parameter sets, the number of continued-fraction terms used, and a convergence check (e.g., stability of the root against increasing truncation order).","section":"Sec. III A and III B; Figs. 2-4"},{"comment":"Equation (23) as written is internally inconsistent and cannot serve as the evidence for the universal law. It states ω_n ∝ 1 - r_i/(2 r_c), but the immediately following equality reads 1 - (r_h + r_i)/(2 r_c) + O((r_h/r_c)^2), which contains an additional -r_h/(2 r_c) term. For uncharged black holes (r_i = 0), Eq. (23) predicts no shift from the empty de Sitter value, contradicting Fig. 3 and Eq. (22), which give ω_n/ω_n^{(dS)} = 1 - M/r_c + ... = 1 - r_h/(2 r_c) + ... . The fitted formula must be corrected and supported by a table of numerical data and residuals; without that, the claim that the universal law holds for gravitational perturbations is not substantiated.","section":"Sec. IV B, Eq. (23)"},{"comment":"The abstract and conclusions state that the purely imaginary modes are responsible for the exponential asymptotic tails, but the paper presents no time-domain profile, no Prony extraction, and no comparison of an extracted late-time decay rate with Im(ω) of the least-damped purely imaginary mode. The time-domain method is only described; its results are never displayed or quantified. Since the exponential-tail connection is a central claim, please include an actual time-domain evolution for at least one representative case, extract the late-time decay rate, and compare it with the frequency-domain value.","section":"Sec. III B and Conclusions"},{"comment":"The extension of the universal law (22) from test fields to the coupled gravitational-electromagnetic spectrum is asserted on the basis of the fitted formula (23), but no independent analytic or semi-analytic confirmation is given. Given the problems with Eq. (23) noted above, the reader cannot distinguish a genuine confirmation of the universal law from a fit that merely reproduces the numerical data. Please provide either a derivation of the leading correction for the coupled system or at least a clear table comparing the numerical ω_n with the prediction ω_n^{(dS)}(1 - M/r_c) for several values of r_h/r_c and r_i/r_h.","section":"Sec. IV B"}],"minor_comments":[{"comment":"The title contains a typographical spacing error: \"Reiss ner-Nordström\" should read \"Reissner-Nordström.\"","section":"Title"},{"comment":"The phrase \"multiple number\" in the definition of λ should be \"multipole number.\"","section":"Sec. II"},{"comment":"In the sentence beginning \"As shown for scalar field perturbations in [42] and for gravitational perturbations in [52], in additional to the complex branch,\" the phrase \"in additional to\" should be \"in addition to.\"","section":"Sec. I"},{"comment":"The caption and axes of Fig. 4 are unclear: the top panel shows r_h Re(ω) as a function of r_h/r_c, and the bottom panel shows r_h Im(ω), but the text discusses transitions in the least-damped mode. Please clarify which curves correspond to the de Sitter branch and which to the black-hole branch in both panels, and explain the meaning of the dashed lines (the caption says they represent \"the corresponding modes in the parametric region, when they are not the least damped\").","section":"Fig. 4"},{"comment":"The paper states that ℓ = 2, 3, 4, ... are considered, but only ℓ = 2 results are displayed. If modes with higher ℓ were computed, showing at least a representative table or a statement of their behavior would strengthen the universality claim.","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the qualitative physics is likely correct, but the lack of numerical tables and the inconsistency in Eq. (23) are serious presentation issues. The authors should be encouraged to provide the missing numerical details and to correct the approximate formula, as this would likely make the paper acceptable. I also note that the universal law is cited predominantly from the authors' own previous work; while this is not improper, the new confirmation should be made fully transparent with data tables."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the result is probably right, but the paper as written doesn't prove all of what it claims. The genuinely new piece is the de Sitter branch of purely imaginary quasinormal modes for the coupled gravitational-electromagnetic perturbations of Reissner-Nordström-de Sitter. That is a real gap, and the authors close it using the standard Leaver continued-fraction method with Nollert acceleration, cross-checked with time-domain integration. The physical picture—two branches, least-damped mode crossing from de Sitter to black-hole branch as the hole grows—is clean and consistent with prior scalar-field results.\n\nWhat the paper does well: it sets up the axial/polar isospectrality via the Darboux relation, so the \"±\" treatment is sound, and the small-black-hole limit correctly reproduces the empty de Sitter modes. The claim that the universal law ω = ω^(dS)(1 − M/r_c + ...) holds for gravitational perturbations is plausible because that law was already derived for test fields and is cited from the authors' own prior work.\n\nSoft spots, in order of severity. Most serious: the abstract and conclusions say these modes are responsible for the exponential asymptotic tails, but no time-domain profiles, Prony extractions, or decay-rate comparisons appear anywhere. The tail claim is asserted, not demonstrated. Second, Eq. (23) as printed is wrong: ω ∝ 1 − r_i/(2r_c) gives no shift for uncharged black holes, contradicting Fig. 3 and the very next equality. It should presumably be 1 − (r_h + r_i)/(2r_c) or similar; this looks like a typo, but it sits in the key confirmation step and needs fixing. Third, there are no numerical tables, truncation orders, or convergence checks. For purely imaginary modes the continued fraction is known to need care, and Nollert acceleration is mentioned but not quantified. Without numbers, the work is hard to build on.\n\nNone of this suggests the central computation is wrong; it is likely correct. But the paper is less auditable than it should be, and the tail claim is unsubstantiated as written. Recommendation: send to peer review, with the expectation that the authors add tables, fix Eq. (23), and provide at least one time-domain extraction. The referee will have a clear, doable list.","headline":"First computation of the de Sitter branch for gravitational QNMs of Reissner-Nordström-de Sitter is plausible and likely correct, but the paper overclaims the tail connection and has a misprinted key formula.","tokens_in":12558,"tokens_out":3029,"would_cite":true,"duration_ms":32527,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.30.Nk","04.50.+h"],"model":"deepseek-v4-flash","headline":"Black holes with charge and a cosmological constant host a purely imaginary branch of quasinormal modes that sets the late-time exponential decay.","keywords":["quasinormal modes","Reissner-Nordström-de Sitter","purely imaginary modes","de Sitter branch","gravitational perturbations","exponential tails","continued fractions","cosmological constant"],"falsifier":"Recompute the lowest purely imaginary mode for, say, $\\ell=2$, $r_i=r_h/2$, and $r_c=10 r_h$ with an independent high-precision method, such as direct numerical integration of the radial equation or a spectral collocation scheme, and compare with the continued-fraction value; any mismatch beyond the stated precision would refute the claimed spectrum. Alternatively, measure the late-time decay in a time-domain evolution: if the logarithmic decay rate does not equal the imaginary part of the least-damped purely imaginary mode, the attribution of the exponential tail to this branch fails.","tokens_in":11539,"feed_emoji":"🕳️","tokens_out":15565,"duration_ms":131451,"temperature":0.7,"pith_summary":"This paper asks whether the gravitational and electromagnetic perturbations of a Reissner-Nordström-de Sitter black hole contain the non-oscillatory 'de Sitter branch' of quasinormal modes previously known only for test fields. It answers yes, and shows that these purely imaginary modes reduce to the modes of empty de Sitter space in the small-mass limit, matching the universal formula $\\omega_n = \\omega_n^{(dS)}(1 - M/r_c + O((M/r_c)^2))$. The paper computes the frequencies with a continued-fraction method and demonstrates by time-domain integration that these modes generate the exponential late-time tails of the perturbation. If correct, this completes the quasinormal spectrum of this classical black-hole solution and gives a parameter-free prediction for the late-time decay of small charged black holes in a universe with a cosmological constant.","feed_headline":"Charged black holes gain purely decaying quasinormal modes","feed_subtitle":"These modes follow a universal formula and explain exponential late-time tails of black-hole signals.","key_machinery":"The central object is the purely imaginary quasinormal mode, a solution of the master wave equation (3) whose frequency is purely imaginary, making the mode non-oscillatory and exponentially decaying under the quasinormal boundary conditions of outgoing waves at the cosmological horizon and ingoing waves at the event horizon. The argument is carried by the universal scaling law $\\omega_n = \\omega_n^{(dS)}(1 - M/r_c + O((M/r_c)^2))$, which is tested numerically for the effective potentials (5) and (6). The computational machinery is a frequency-domain continued-fraction construction: the radial function is expanded as a Frobenius series, producing seven- and nine-term recurrence relations that are reduced to three terms by Gaussian elimination, and the quasinormal frequencies are the solutions of the resulting infinite continued-fraction equation, with a convergence-acceleration step applied for the purely imaginary modes. The results are cross-checked by a light-cone time-domain integration scheme that measures the exponential decay rate of the perturbation.","core_discovery":"The paper's discovery is that the coupled gravitational and electromagnetic perturbations of the Reissner-Nordström-de Sitter black hole carry a de Sitter branch of purely imaginary quasinormal frequencies, in addition to the well-known complex black-hole branch. These modes are deformations of the modes of empty de Sitter space and obey the universal law $\\omega_n = \\omega_n^{(dS)}(1 - M/r_c + O((M/r_c)^2))$ in the small-black-hole regime, a formula this paper shows holds for both the '+' and '−' effective potentials. The numerical computation indicates the approximation stays accurate even when the black-hole radius is comparable to the cosmological horizon. Time-domain integration confirms that the exponential asymptotic tails of the perturbation are produced by these purely imaginary modes, establishing them as the least-damped part of the spectrum for small black holes even though they are only weakly excited by typical initial data.","pith_inferences":["The universal law likely holds for any spherically symmetric black hole with a cosmological horizon, independent of the underlying gravitational theory, since the paper's derivation is metric-independent; this extension is testable by computing the de Sitter branch for other backgrounds.","Observational searches for the imprint of a cosmological constant on gravitational-wave signals could target an anomalous exponential-decay phase after the ringdown of a charged black hole rather than fitting a single complex mode, because the de Sitter branch is weakly excited.","The near-linear dependence of the modes on $r_i/r_h$ seen in the paper's Fig. 2 may reflect a hidden symmetry of the perturbation equations, which if identified could yield a closed-form expression for the whole de Sitter branch.","A dedicated high-precision computation in the near-extremal charge limit $r_i \\to r_h$ would test whether the fitted formula (23) remains valid or whether the recurrence reduction changes character there."],"forward_implications":["The quasinormal spectrum of Reissner-Nordström-de Sitter is now known to have two branches, and for small black holes the purely imaginary de Sitter branch is the least-damped set of modes.","The universal law (22) applies to the actual gravitational and electromagnetic perturbations, so the late-time decay rate of a small charged black hole is fixed without free parameters.","Late-time gravitational-wave tails from these black holes are exponential rather than power-law, with the decay rate set by the imaginary part of the least-damped purely imaginary mode.","As the black-hole radius grows, overtones of the complex black-hole branch successively take over as the least-damped modes, making the dominant damping rate a non-monotonic function of the black-hole radius.","In the near-extremal limit the purely imaginary mode persists with a nonzero decay rate, because the black-hole mass always stays below the horizon radius, so the exponential-tail phase never disappears."],"supporting_citations":[{"why":"Shows that purely imaginary modes of Schwarzschild-de Sitter satisfy the ingoing-horizon boundary condition and govern the exponential tails, the phenomenon this paper extends to charged holes.","marker":"[52]"},{"why":"Derives the universal law $\\omega_n = \\omega_n^{(dS)}(1 - M/r_c + O((M/r_c)^2))$ for test fields, which this paper promotes to the gravitational sector.","marker":"[91]"},{"why":"Identifies the purely imaginary modes for scalar-field perturbations of Reissner-Nordström-de Sitter, the starting point for the gravitational analysis.","marker":"[42]"},{"why":"Provide the frequencies of empty de Sitter space that define the zero-mass limit of the de Sitter branch.","marker":"[54, 55]"},{"why":"Supplies the frequency-domain continued-fraction method used to compute the quasinormal frequencies.","marker":"[62]"},{"why":"Supplies the light-cone time-domain integration scheme used to extract the decay rates and verify the exponential tails.","marker":"[63]"},{"why":"Provide the convergence-acceleration step needed to obtain accurate continued-fraction solutions for purely imaginary modes.","marker":"[64, 65]"},{"why":"Gives the coupled perturbation equations and the argument that the axial and polar spectra coincide, reducing the problem to two effective potentials.","marker":"[51]"}],"fun_headline_variants":["Purely imaginary modes found for gravitationally perturbed charged black holes","New branch of purely decaying modes for charged black holes","Exponential tails explained by purely imaginary quasinormal modes","Universal law for purely imaginary black-hole modes","De Sitter branch of modes found for Reissner-Nordstrom black holes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical results assume that the continued-fraction recurrence reduction and the time-domain fitting converge to the true quasinormal frequencies, but the paper reports no convergence checks, truncation orders, or error estimates; if either numerical scheme fails to converge, the claimed frequencies and the fitted universal formula do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Purely imaginary modes found for gravitationally perturbed charged black holes","New branch of purely decaying modes for charged black holes","Exponential tails explained by purely imaginary quasinormal modes","Universal law for purely imaginary black-hole modes","De Sitter branch of modes found for Reissner-Nordstrom black holes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000656,"raw_usage":{"total_tokens":2937,"prompt_tokens":812,"completion_tokens":2125,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":428,"completion_tokens_details":{"reasoning_tokens":2040}},"tokens_in":428,"tokens_out":2125,"duration_ms":17911,"temperature":1.0,"reasoning_tokens":2040,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:38:51.564648+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the lowest purely imaginary mode for, say, $\\ell=2$, $r_i=r_h/2$, and $r_c=10 r_h$ with an independent high-precision method, such as direct numerical integration of the radial equation or a spectral collocation scheme, and compare with the continued-fraction value; any mismatch beyond the stated precision would refute the claimed spectrum. Alternatively, measure the late-time decay in a time-domain evolution: if the logarithmic decay rate does not equal the imaginary part of the least-damped purely imaginary mode, the attribution of the exponential tail to this branch fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows that purely imaginary modes of Schwarzschild-de Sitter satisfy the ingoing-horizon boundary condition and govern the exponential tails, the phenomenon this paper extends to charged holes."},{"cited_title":"Properties of the Reissner-Nordstr\\\"om Spacetimes with a Nonzero Cosmological Constant","cited_arxiv_id":"0803.2685","evidence_quote":"Supplies the light-cone time-domain integration scheme used to extract the decay rates and verify the exponential tails."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the coupled perturbation equations and the argument that the axial and polar spectra coincide, reducing the problem to two effective potentials."}],"review_version":1}