{"id":"2ba53eba-fcd6-4b81-afb7-deb7cdf52adb","arxiv_id":"2607.18399","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Charged Ellis–Bronnikov wormholes develop a previously missed l=2 polar quasinormal-mode instability at sufficiently large masses, in addition to the known radial instability.","lead":"This paper computes how charged Ellis–Bronnikov wormholes vibrate and finds a new instability in their nonradial oscillations. The result matters because it tightens the conditions under which such wormholes could act as black-hole mimics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The l=2 polar instability is plausible but rests on a single unvalidated spectral computation in the least-accurate regime; an independent numerical cross-check is needed before this becomes a secure result.","rationale":"The reader's weakest assumption — numerical reliability of the spectral solver for the l=2 b2 polar branch in the large-parameter regime — is exactly the load-bearing concern. I checked whether the boundary-condition choice e^{iωr*} could be an internal flaw: it is not, because the modified tortoise coordinate r* tends to +∞ on both asymptotic ends, so the factorization implements outgoing waves at both infinities. The uncharged benchmarks and EM isospectrality provide some support for the code, but those checks do not cover the new unstable branch. The numerical crossing is not marginal (ω_I jumps from -0.001 to +0.013 in Table 2, and Table 9 gives ω_I≈0.067), so I would not escalate to REJECT; the result is plausible but unverified. Conditional acceptance, as the reader recommends, is the right level of confidence, and the specific test above (time-domain evolution plus a second frequency-domain solver) would settle the issue. No change to the verdict is required.","tokens_in":24602,"tokens_out":15360,"duration_ms":127850,"concrete_test":"Perform an independent time-domain evolution of the linearized polar l=2 system (finite-difference or discontinuous Galerkin) on the uncharged EB background for Λ=0.8 and Λ=0.9, and for subcritical Λ=0.75, γ1=0.9; fit the late-time signal to e^{-iω t} and compare Im(ω). In parallel, recompute the same eigenvalues with a shooting method (or a second spectral code) at Np=20,40,60 and check that Im(ω) converges to the same positive value and that the b2 branch is continuously connected from the M/rT=0 degenerate limit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim — a previously unknown instability in the fundamental l=2 b2 polar branch (§6) — depends entirely on the Chebyshev spectral solver. The authors state in §4 that accuracy 'remains better than 10^{-3}' but degrades for higher modes and larger Λ and γ1; they give no per-mode error estimates, no Np/convergence data, and no independent method. The crossing in Table 2 (EB: ω_I=-0.001 at Λ=0.7, +0.013 at Λ=0.8) is not marginal at the stated 1e-3 level, and Table 9 (subcritical Λ=0.75, γ1=0.80–0.90) lists ω_I≈0.067, so the instability is not a tiny numerical effect. The concern is rather that this is a single spectral implementation in exactly the parameter region the authors flag as least trustworthy, and the b2 branch is identified by continuity across parameter space with no second solver. If the mode is spurious or the branch is mislabeled at the avoided crossing, the headline instability disappears even though the rest of the QNM catalog may be correct.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives the axial, radial, and polar perturbation equations for static charged Ellis–Bronnikov wormholes in Einstein–Maxwell theory with a phantom scalar, and computes quasinormal-mode spectra with a Chebyshev spectral method. The uncharged limit is benchmarked against known Ellis–Bronnikov results and the expected axial–polar electromagnetic isospectrality is recovered. For charged configurations, the paper tracks axial and polar branches across critical, subcritical, and supercritical families, finding that charge generally increases damping times near the extremal Reissner–Nordstrøm limit. The headline new claim is a previously unknown polar instability: the fundamental l=2 b2 branch crosses ωI=0 at M/rT≈0.3 for uncharged EB wormholes (Sec. 5.1.2, Fig. 3, Table 2) and is also unstable for subcritical wormholes at sufficiently large Λ (Sec. 5.3, Table 9). The authors also note that critical wormholes do not show this instability in the computed range, while supercritical solutions appear stable.","tokens_in":24888,"tokens_out":5223,"duration_ms":48156,"significance":"If the polar l=2 b2 instability is real, it constitutes a genuinely new dynamical-instability channel for charged and uncharged EB wormholes, distinct from the known radial l=0 instability, and it would strengthen the constraints on the viability of these wormhole models. The paper's strengths are its explicit perturbation systems, the closed-form backgrounds, the reproduction of the known uncharged spectrum, and the demonstration of EM isospectrality in the uncharged limit. There are no fitted parameters in the reported spectra. The main result, however, depends on a single spectral implementation in the parameter regime that the authors themselves flag as least accurate, with no independent numerical confirmation. The catalog of stable modes may well be reliable, but the central instability claim requires additional verification before it can be considered secure.","major_comments":[{"comment":"The headline instability rests on one spectral method without independent validation. Sec. 4 states that accuracy 'remains better than 10^-3' but degrades for higher excited modes and larger Λ and γ1, and no per-mode error estimates, Np-convergence data, or second numerical method are provided. The crossing in Table 2 occurs between ωI r0 = -0.001 (Λ=0.7) and +0.013 (Λ=0.8), and Table 9 reports ωI ≈ 0.067 for Λ=0.75 in the same high-parameter regime. While these values exceed 10^-3 in absolute terms, the absence of convergence documentation for exactly this branch and regime leaves the central instability claim under-supported. I request a convergence table at the crossing and at the unstable points, plus at least one independent check (e.g., time-domain integration, a different discretization, or a separate spectral implementation).","section":"Sec. 5.1.2 / Table 2, with Sec. 4"},{"comment":"The b1/b2 labels are introduced because the polar modes are mixed gravitational–scalar, and the branches exhibit crossings and repulsions in the complex-frequency plane. The claim that the unstable mode is 'the fundamental l=2 b2 branch' requires unambiguous branch tracking through these features, especially near M/rT ≈ 0.23 where the imaginary parts of the b1 and b2 branches cross and the real parts repel. The paper does not provide eigenfunction overlaps or any explicit branch-tracking criterion. If the labels switch at an avoided crossing, the instability could belong to a different physical mode, or the continuity of the branch used in Tables 2 and 9 could be broken. Please provide evidence of branch continuity, for instance overlap integrals of eigenvectors between adjacent parameter points.","section":"Sec. 5.1.2 / Figs. 3 and 5"},{"comment":"The conclusion that critical wormholes do not exhibit the polar l=2 instability is stated to hold only within the computed range, and the text concedes that 'it is possible that such an instability appears for large values of the parameters, for which our method loses accuracy.' This is appropriately cautious, but it means the claimed dichotomy (stable critical, unstable subcritical) is not a proven boundary but an upper-bound statement. A quantitative statement of the maximum accessible parameter value and the corresponding accuracy would make the limitation precise.","section":"Sec. 6 / critical and subcritical conclusions"}],"minor_comments":[{"comment":"The appendix tables report frequencies with r0=1, while the main figures scale by the throat radius rT. Since rT ≠ r0 for charged wormholes, the reader cannot directly compare, e.g., the Λ=0.8 crossing in Table 2 with the M/rT ≈ 0.3 statement in Sec. 5.1.2. Please state the relation between rT, M, Qe, and r0, and give the key crossing also in rT units.","section":"Tables 1–12 vs. Fig. 3"},{"comment":"The phrase 'breaking of isospectrality between the b1 and b2 branches' is imprecise: isospectrality was introduced for axial and polar EM modes. The b1/b2 splitting is a degeneracy-breaking of the scalar–gravitational polar branches. Please rephrase to avoid confusion.","section":"Sec. 5.4.2"},{"comment":"The expansion coefficient i0 is easily misread as the imaginary unit multiplied by zero. Please use a distinct symbol such as \\iota_0 or i_0.","section":"Eq. (44)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid QNM catalog with a clean uncharged benchmark, and the derivations are explicit. The novelty, however, is concentrated in the polar instability claim, which is exactly where the numerical evidence is thinnest. The stress-test concern about a single unvalidated spectral computation in the least-accurate regime is justified. I would not recommend acceptance until an independent numerical cross-check and convergence data for the unstable branch are provided. This is fixable within the manuscript's scope, so major_revision rather than rejection seems appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: this is a careful, mostly convincing QNM catalog for charged Ellis–Bronnikov wormholes, and the claimed new l=2 polar instability is probably real, but the paper would be more persuasive if the unstable branch had been checked with an independent method.\n\nWhat's actually new: the authors derive the full axial and polar perturbation equations for the Einstein–Maxwell–phantom system, including coupled gravitational, EM, and scalar degrees of freedom, and compute spectra across all three charged wormhole families (subcritical, critical, supercritical). The uncharged limit reproduces the known EB frequencies and shows EM isospectrality, which is a strong consistency check. They also extend the mass range beyond earlier work. The headline result—a previously unseen instability in the fundamental l=2 polar b2 mode beyond M/rT ≈ 0.3—is genuinely new and, if correct, adds a nonradial instability channel to the known radial one.\n\nThe soft spot is numerical verification. The instability is identified with a single Chebyshev spectral solver, and the authors themselves note that accuracy degrades for larger Λ and γ1—exactly the region where the instability appears. The crossing in Table 2 (from -0.001 to +0.013) is not tiny at their stated 10^-3 accuracy, and Table 9 shows ωI ≈ 0.067, so this isn't just noise. But there are no convergence plots, no per-mode error bars, and no time-domain evolution or second spectral implementation. The branch labeling is done by continuity across parameter space. That's enough to make the result plausible, not enough to make it bulletproof.\n\nI also note the paper honestly acknowledges where the method loses accuracy and where they had to stop. The perturbation equations are given explicitly, so an independent check is feasible. The citation pattern is fine—mostly their own prior work, but the background is closed-form and the benchmarks are independent.\n\nWho this is for: anyone working on wormhole stability or black-hole mimickers. It's a solid subfield contribution, not a paradigm shift.\n\nMy recommendation: it deserves a serious referee. The numerical verification concern is substantial but addressable—ask for a second method or detailed convergence study on the unstable branch before publication, or at least a clear statement that the instability is robust across grid sizes. I'd lean conditional accept.","headline":"Plausible new l=2 polar instability in charged EB wormholes, well-derived but needs an independent numerical check before I'd call it secure.","tokens_in":25356,"tokens_out":2773,"would_cite":true,"duration_ms":23782,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Static charged wormholes carry a previously unknown nonradial instability: the fundamental quadrupolar polar mode grows once the mass is large enough, in addition to the known radial instability.","keywords":["charged Ellis-Bronnikov wormholes","quasinormal modes","nonradial instability","polar perturbations","phantom scalar field","extremal Reissner-Nordström limit","Einstein-Maxwell theory","spectral methods"],"falsifier":"Run an independent time-domain integration of the linearized polar perturbation equations for an uncharged Ellis-Bronnikov wormhole with mass-to-throat ratio about 0.4: if the quadrupolar polar perturbation rings down instead of growing exponentially, the reported instability is a numerical artifact.","tokens_in":24551,"feed_emoji":"🕳️","tokens_out":9759,"duration_ms":74950,"temperature":0.7,"pith_summary":"This paper computes the quasinormal-mode spectrum of static, spherically symmetric charged Ellis-Bronnikov wormholes — solutions of Einstein-Maxwell theory with a phantom scalar field — for both axial and polar perturbations, including gravitational, electromagnetic, and scalar degrees of freedom. Its central result is the discovery of a nonradial instability: for uncharged Ellis-Bronnikov wormholes, the fundamental polar branch labeled b2 at multipole l=2 crosses from damped to growing once the mass-to-throat ratio exceeds about 0.3, and a similar instability appears for subcritical charged wormholes with sufficiently large parameter Λ. This instability is distinct from the well-known radial (l=0) instability and, if real, means that the nonradial sector imposes additional constraints on the dynamical viability of these wormholes. The paper also finds that electric charge generally lengthens damping times, with mode frequencies and damping rates approaching the extremal Reissner-Nordström limit in the most charged configurations.","feed_headline":"A second, nonradial instability hits massive wormholes","feed_subtitle":"Beyond the radial mode, a quadrupolar polar branch grows past a critical mass, tightening constraints on wormhole mimickers.","key_machinery":"The load-bearing object is the b2 polar branch: the second family of polar quasinormal modes that appears because gravitational and phantom-scalar perturbations are coupled, so that each multipole carries two polar branches (b1 and b2) rather than one. The paper follows this branch with a Chebyshev spectral method, a numerical expansion of the perturbation functions in Chebyshev polynomials on a compactified radial grid: the perturbation equations are reduced, via spherical-harmonic decomposition, to a quadratic eigenvalue problem in the complex frequency ω, and the sign of Im(ω) decides stability (negative = damped, positive = growing). The same machinery produces the axial, electromagnetic","core_discovery":"The authors derive the full first-order perturbation equations for charged Ellis-Bronnikov wormholes, decomposing the metric, electromagnetic, and phantom-scalar perturbations into axial and polar sectors. In the uncharged limit they recover the known Ellis-Bronnikov spectrum and confirm the expected isospectrality of electromagnetic modes between axial and polar sectors. Their principal claim is that the fundamental polar branch at multipole l=2 — the branch they call b2, in which gravitational and phantom-scalar perturbations are strongly coupled — becomes unstable for sufficiently massive wormholes. For uncharged Ellis-Bronnikov wormholes its imaginary frequency crosses zero at a mass-to-","pith_inferences":["Editorial inference: The most direct confirmation would be an independent time-domain evolution of the linearized polar equations for an uncharged Ellis-Bronnikov wormhole at M/r_T ≈ 0.4; if the quadrupolar channel does not grow, the reported crossing is a numerical artifact.","Editorial inference: If the instability survives nonlinear evolution, it could set the effective lifetime of massive wormhole black-hole mimickers, possibly making them shorter-lived than the radial instability alone would suggest.","Editorial inference: The authors' analogy between charge and rotation suggests a concrete test: compute the l=2 polar modes of rapidly rotating Ellis-Bronnikov wormholes to see whether rotation similarly suppresses the nonradial instability; the paper only conjectures the radial analogue.","Editorial inference: The finding that charge suppresses the l=2 growth rate raises the possibility of a critical charge-to-mass ratio above which the nonradial instability disappears entirely — a threshold the paper does not map, but a natural next calculation."],"forward_implications":["If the l=2 polar instability is real, stability analyses of Ellis-Bronnikov wormholes that focus only on radial modes are incomplete: sufficiently massive configurations are dynamically unstable in the nonradial sector as well.","Electric charge can substantially reduce damping rates: as the extremal Reissner-Nordström limit is approached, modes become long-lived and all branches accumulate near zero imaginary frequency, so charged wormholes can ring almost like the extremal black hole.","Subcritical charged wormholes inherit the l=2 instability for large Λ, but its growth rate is suppressed by charge; the critical and supercritical families show no l=2 instability in the computed range, with the caveat that accuracy degrades at large parameters.","The radial instability relaxes near the extremal limit (ω_I scales roughly as (1 − M/r_T)^3.1), giving arbitrarily long instability timescales, while supercritical wormholes instead show two unstable radial branches that merge and acquire real frequencies with opposite signs."],"fun_headline_variants":["Massive charged wormholes face new nonradial instability","Quadrupolar mode destabilizes heavy wormholes","Charge shrinks damping but spawns polar instability","Wormhole viability tightened by l=2 instability","Nonradial mode adds new constraint on wormholes"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The result stands or falls on whether the numerical solver's tracking of one particular oscillation branch — the second polar branch of the quadrupolar mode — remains trustworthy in the large-mass regime, where the authors report degraded accuracy, so that its crossing from damped to growing is physical rather than a numerical artifact.","fun_headline_variants_meta":{"raw":{"variants":["Massive charged wormholes face new nonradial instability","Quadrupolar mode destabilizes heavy wormholes","Charge shrinks damping but spawns polar instability","Wormhole viability tightened by l=2 instability","Nonradial mode adds new constraint on wormholes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000149,"raw_usage":{"total_tokens":1034,"prompt_tokens":751,"completion_tokens":283,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":495,"completion_tokens_details":{"reasoning_tokens":208}},"tokens_in":495,"tokens_out":283,"duration_ms":4605,"temperature":1.0,"reasoning_tokens":208,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T15:30:21.417882+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run an independent time-domain integration of the linearized polar perturbation equations for an uncharged Ellis-Bronnikov wormhole with mass-to-throat ratio about 0.4: if the quadrupolar polar perturbation rings down instead of growing exponentially, the reported instability is a numerical artifact.","supporting_citations":[],"review_version":1}