{"id":"c2cbcee8-8d04-4cbf-8fe1-1e51220fdf2f","arxiv_id":"2601.13521","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Every computed fundamental quasinormal frequency for polar metric-dilaton perturbations of dilaton-Euler-Heisenberg black holes has a negative imaginary part, meaning the modes are damped, with ε=+1 and ε=−1 differing in damping near extremality.","lead":"This paper computes the ringdown vibration frequencies of a charged, dilaton-carrying black hole from a string-inspired gravity theory, using two independent numerical methods that mostly agree. All computed modes decay rather than grow, indicating stability, and the damping differs between the two signs of the coupling parameter near the extremal limit.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Master-equation reduction (Eqs. 20-26) is the load-bearing step; the omitted elimination from (15)-(19) and the paper's own over-constraining warning leave the QNM frequencies, and hence the stability claim, vulnerable to a hidden constraint.","rationale":"The reader's weakest assumption is exactly the reduction of the polar sector to (20)-(26); I agree this is load-bearing. The two-method agreement is the paper's main evidence, but it only verifies the reduced system, not the reduction itself. The paper's own over-constraining warning and the absence of a detailed derivation make this the least-secure link. The background metric issue in Eq. (11) is real and should be fixed, but it is likely a typesetting problem; the reduction correctness is the deeper risk to the central stability claim. A symbolic re-derivation or a full-system solve would settle the issue. Since this concern is already reflected in the CONDITIONAL verdict, no change is needed.","tokens_in":20224,"tokens_out":10176,"duration_ms":93646,"concrete_test":"Independently derive the polar perturbation equations from action (1) about the solution (8)-(9) using a symbolic perturbation package (e.g., xAct/xPert), and compare term-by-term with (20)-(26). As a numerical cross-check, solve the original first-order system (15)-(18) with (19) treated as a constraint at (Q_m=0.6, l=2, epsilon=-1), the point where DI and CFM differ by 3.4% in Table III, and compare the fundamental frequency with the table. If the full-system frequency disagrees with the master-equation result, the reduction is over-constrained.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on the two coupled master equations (20) and (23) with potentials (21)-(26). Their derivation from the first-order system (15)-(19) is not shown: Appendix A only cites the Brito-Pacilio notebook [68] and asserts the elimination through (A2)-(A10), while the key algebraic identity (19) is presented without proof. The paper itself warns in Sec. IV.A that the first-order system must not be over-constrained. Both numerical methods solve the same reduced equations, so their mutual agreement does not test the reduction itself. If identity (19) or the elimination imposes an extra constraint—or drops a coupling—every computed frequency, and therefore the all-negative-Im stability conclusion, is in question. Separately, the typeset background (11) is not asymptotically flat as written (A does not tend to 1), so the printed potentials do not uniquely define the claimed spacetime; this is a real reproducibility flaw that should be corrected, but the reduction correctness is the deeper epistemic risk.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies linear polar (even-parity) metric-dilaton perturbations of the magnetically charged dilaton–Euler–Heisenberg (dEH) black hole with dilaton hair. Starting from four radial equations (15)–(18) and an algebraic identity (19), it reduces the polar sector to two coupled second-order master equations for K and S with potentials (21)–(26). The authors compute fundamental quasinormal frequencies using direct integration and matrix-valued continued-fraction/AIM methods for l=0,...,3 and ϵ=±1, compare against Schwarzschild QNMs for Q_m=0.01, and report that all computed modes have negative imaginary parts, from which they conclude stability. They also identify different near-extremal damping behavior for ϵ=+1 and ϵ=−1.","tokens_in":20426,"tokens_out":6016,"duration_ms":57977,"significance":"If correct, the paper would be one of the first stability/ringdown studies of the dEH family and would strengthen the case for using ringdown measurements to constrain dilaton-EH couplings. The numerical work has clear strengths: no output parameters are fitted; the Q_m=0.01 limit reproduces the standard Schwarzschild values (Tables III–IV; e.g. 0.373681−0.0889637i vs 0.3737−0.08896i for l=2 gravitational); and two independent solvers agree to sub-percent in most cases. The l=0,1 potential plots are positive definite. However, the central reduction is asserted rather than demonstrated, the background metric as printed is not asymptotically flat, and the stability conclusion exceeds what a finite set of n=0 modes can show.","major_comments":[{"comment":"The load-bearing step is the reduction of the four-equation first-order system to the two master equations (20) and (23). The paper states identity (19) without proof and Appendix A only lists the ansatz coefficients, citing the Mathematica notebook of Ref. [68]. The numerical agreement between DI and CFM does not test this reduction because both methods solve the same reduced equations. Please provide the full derivation or the notebook, including how H0 and R1 are eliminated and how over-constraining of the first-order system is avoided. Without this, the frequencies—and hence the stability conclusion—cannot be independently checked.","section":"§III and Appendix A, Eqs. (19)–(26)"},{"comment":"As printed, the background metric function A(r) is not asymptotically flat and does not reduce to Schwarzschild: for r→∞, A(r) ≈ 1 − 4M^2/Q_m^2 + 2Mr + O(r^{-5}), and the Q_m→0 limit is singular. The original metric in Eqs. (8)–(9) is asymptotically flat, so this appears to be a coordinate/transcription error, likely a missing denominator involving Q_m^2 + sqrt(Q_m^4+4M^2r^2). Since the potentials in Eqs. (21)–(26) are constructed from these A and B, the tables cannot be reproduced from the printed equations. Please correct Eq. (11) and state the coordinate relation to Eqs. (8)–(9).","section":"§II, Eq. (11)"},{"comment":"The inference that 'all negative imaginary quasinormal frequencies imply stability' is too strong. The work computes only n=0 fundamental frequencies for l≤3; for l=0,1 potential positivity is argued, but for the coupled l=2,3 system no such argument is given. A finite sample of damped modes is evidence, not proof, of linear stability. Please either restrict the claim to 'no unstable modes found in the computed set' or add a stability proof (e.g., S-deformation, potential-positivity, or absence of unstable eigenvalues) for the coupled l≥2 system.","section":"§V and Conclusion"}],"minor_comments":[{"comment":"The table captions use 'bEH' instead of 'dEH'; also the discrepancy is called ΔDA for l=0,1 and ΔDC for l≥2. Standardize the notation.","section":"Tables I–II"},{"comment":"Tables I–II use the asymptotic iteration method (AIM), but Section IV only describes direct integration and matrix-valued continued fraction. Clarify that the l=0,1 results use AIM (Appendix B) while l≥2 uses DI/CFM.","section":"Sec. IV and Appendix B"},{"comment":"The legends and line styles for gravitational versus dilaton modes are hard to follow, especially in Figs. 5(b) and 6(b). Please use explicit labels such as 'Im ω' and distinguish modes more clearly.","section":"Figs. 5–6"},{"comment":"The phrase 'All negative imaginary quasinormal frequencies' should be qualified as 'all computed fundamental quasinormal frequencies' to avoid overstating the result, consistent with the major comment above.","section":"Abstract and Sec. I"},{"comment":"Since the reduction is central, please provide a stable reference or link to the exact version of the Mathematica notebook used, and state which simplifications were applied to eliminate exponential factors and derivatives of A and B.","section":"Appendix A, Ref. [68]"}],"recommendation":"major_revision","confidential_remarks":"The numerical computations appear competent and the Schwarzschild benchmarks are good, but the omitted master-equation reduction and the incorrect printed metric are blockers for reproducibility. If the authors supply the full reduction/code and correct Eq. (11), and temper the stability claim, the paper could be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: if you need quasinormal frequencies for the 2024 dilaton-Euler-Heisenberg black hole family, this is the first place to look — and the paper also shows why you should double-check the printed background before using any of its numbers.\n\nWhat is genuinely new: the first QNM computation for the dEH solution family, with a coupled polar metric-dilaton system that, if the reduction is right, is a reusable tool. The paper does a number of things well. It checks against Schwarzschild at Q_m=0.01 and matches to better than 0.01%; the DI and CFM methods agree to under 1% for almost all data points; and the ε=±1 scan is honest parameter variation, not a fit. The l=0 and l=1 potentials are shown to be positive definite, which is a solid stability argument for those sectors. The authors also flag the over-constraining risk in Sec. IV.A and admit in Sec. VI that axial and electromagnetic perturbations are still open.\n\nThe soft spots are real. The background metric in Eq. (11) is not asymptotically flat: A(r) grows like 2Mr rather than tending to 1, so the equations as printed do not define the spacetime whose modes the tables report. This looks like a transcription error rather than a conceptual one, but it is load-bearing: no one can reproduce the computation from the paper. Second, the reduction from the first-order system to the two master equations (20) and (23) is never actually shown. Appendix A leans on a Mathematica notebook and an unproved algebraic identity (19). Since both numerical methods solve the same reduced system, their cross-agreement does not test the reduction. The authors' own warning about over-constraining is exactly the risk that needs to be closed. Third, the stability claim outruns the evidence: a finite set of n=0, l≤3 polar modes is consistent with stability, not a proof of it. The positive potentials for l=0,1 are stronger evidence than the negative imaginary parts in Tables III–IV. Finally, the one place the two methods disagree most (3.4% for ε=−1, l=2, Q_m=0.6) is also the headline near-extremal behavior, so that qualitative feature is not as solid as the rest.\n\nWho gets value: people doing QNM tables for modified-gravity black holes, or anyone testing the dEH solution's viability. It deserves a serious referee, but the referee should send it back for correction of the metric, a full derivation or a repository, and a moderated stability claim.","headline":"First QNM dataset for dEH black holes, but a wrong printed metric and an under-documented master-equation reduction mean the numbers should not be trusted until corrected.","tokens_in":21030,"tokens_out":6085,"would_cite":false,"duration_ms":65398,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C35"],"pacs":["04.70.Bw","04.30.-w"],"model":"deepseek-v4-flash","headline":"This paper claims the dilaton-Euler-Heisenberg black hole is stable against polar metric-dilaton perturbations: every fundamental quasinormal frequency computed for dilaton modes l=0–3 and gravitational modes l=2–3 has a negative imaginary","keywords":["quasinormal modes","polar perturbations","dilaton-Euler-Heisenberg black holes","stability","dilaton hair","Euler-Heisenberg electrodynamics","direct integration method","matrix-valued continued fraction"],"falsifier":"Solve the original first-order system (15)-(18) with the same horizon and infinity boundary conditions but without using the algebraic identity (19); if any eigenvalue with positive imaginary part appears, the stability claim fails. Separately, verify that the background metric printed in Eq. (11) is asymptotically flat — as r grows it reads A(r)≈2Mr/Q_m^2, which cannot reproduce the Schwarzschild frequencies listed for Q_m=0.01.","tokens_in":20004,"feed_emoji":"🕳️","tokens_out":7569,"duration_ms":73185,"temperature":0.7,"pith_summary":"This paper tries to establish that the dilaton-Euler-Heisenberg black hole — a magnetically charged, spherically symmetric solution of Einstein–Maxwell–dilaton theory with a dilaton-coupled Euler–Heisenberg nonlinear term — is dynamically stable against polar (even-parity) metric and dilaton perturbations. The authors compute the fundamental quasinormal-mode frequencies for dilaton modes with angular momentum l=0,1,2,3 and gravitational modes with l=2,3, using two independent numerical methods whose results agree to better than a few percent. Every computed frequency has a negative imaginary part, meaning every mode decays in time rather than growing. They also find that the damping rates behave differently depending on the sign of the coupling parameter ε=α−β, especially near extremality.","feed_headline":"Dilaton-magnetized black holes: every computed mode is damped","feed_subtitle":"All computed polar modes decay, and the sign of the dilaton coupling changes how fast near extremality.","key_machinery":"The load-bearing object is the algebraic identity (19), obtained from the (r,r) component of the linearized Einstein equations, which eliminates one of the polar metric variables and allows the four coupled first-order equations to be reduced to two coupled second-order master equations (20) and (23) in the tortoise coordinate. These master equations, with potentials (21)–(26), turn the stability question into a complex eigenvalue problem whose eigenvalues are the quasinormal frequencies: a positive imaginary part would signal growth, a negative imaginary part means damping. Two independent numerical schemes — direct integration and matrix-valued continued fraction — are used to solve that e","core_discovery":"The paper studies the linearized polar (even-parity) perturbations of the spherically symmetric dilaton-Euler-Heisenberg black hole with dilaton hair, a magnetically charged solution of Einstein–Maxwell–dilaton theory with a dilaton-coupled Euler–Heisenberg term. In the Regge–Wheeler gauge the polar sector forms four coupled equations for the metric variables and the dilaton; using an algebraic identity and a reduction procedure, the authors obtain two coupled second-order master equations, one for the gravitational channel and one for the dilaton channel. They compute the fundamental (n=0) quasinormal frequencies by direct integration and by a matrix-valued continued fraction method, report","pith_inferences":["If the reduction is exact, the same two-channel scheme should reproduce the ε=0 limit of this family smoothly; the paper's tables approach Schwarzschild at small Q_m but do not tabulate ε=0, so an independent check at ε=0 would test continuity.","The paper computes only n=0 modes, but the l=0,1 effective potentials are positive definite; a spectral argument would likely extend dilaton-mode stability to higher overtones, a step the paper does not take.","The sign asymmetry near extremality tracks the different horizon structures of ε=+1 (single horizon) and ε=−1 (two horizons up to an extremal limit); connecting the damping reversal near Q_m≈0.826 to the disappearance of the second horizon is a natural next calculation.","A direct waveform prediction is within reach: the two coupled master equations give not only frequencies but also channel mixing, so ringdown templates could be generated to assess how much ε and Q_m are actually measurable in a single event."],"forward_implications":["The dEH black hole remains stable against monopole and dipole dilaton perturbations (l=0,1) and against quadrupole and octupole gravitational and dilaton perturbations (l=2,3), at least at the fundamental-mode level and for the parameter range studied.","The real part of the dilaton-mode frequency grows with magnetic charge Q_m, so stronger magnetic hair makes the black hole ring at a higher frequency; the gravitational l=2 mode behaves differently for ε=1 (peaking and declining) than for ε=−1 (monotonic rise).","Near extremality, the damping rate of the ε=−1 gravitational l=2,3 modes reverses and grows, while the ε=1 modes keep declining — the sign of the dilaton coupling leaves an observable fingerprint in the ringdown.","Because all computed imaginary parts are negative, the spacetimes are dynamically stable in this sector, and deviations of their quasinormal frequencies from Schwarzschild values could be used to constrain Q_m and ε with future gravitational-wave ringdown observations."],"fun_headline_variants":["All polar metric-dilaton modes damped for dEH black holes","dEH black holes: every computed polar mode decays","Dilaton-EH black holes stable against all studied polar modes","Coupling sign flips damping behavior near extremal dEH holes","Polar quasinormal spectra: dEH black holes are stable"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The stability conclusion rests on the reduction of the four linearized polar equations to the two master equations being complete and not over-constrained — a step the paper itself flags as risky.","fun_headline_variants_meta":{"raw":{"variants":["All polar metric-dilaton modes damped for dEH black holes","dEH black holes: every computed polar mode decays","Dilaton-EH black holes stable against all studied polar modes","Coupling sign flips damping behavior near extremal dEH holes","Polar quasinormal spectra: dEH black holes are stable"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000303,"raw_usage":{"total_tokens":1578,"prompt_tokens":742,"completion_tokens":836,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":486,"completion_tokens_details":{"reasoning_tokens":747}},"tokens_in":486,"tokens_out":836,"duration_ms":9404,"temperature":1.0,"reasoning_tokens":747,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T09:31:42.235341+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the original first-order system (15)-(18) with the same horizon and infinity boundary conditions but without using the algebraic identity (19); if any eigenvalue with positive imaginary part appears, the stability claim fails. Separately, verify that the background metric printed in Eq. (11) is asymptotically flat — as r grows it reads A(r)≈2Mr/Q_m^2, which cannot reproduce the Schwarzschild frequencies listed for Q_m=0.01.","supporting_citations":[],"review_version":1}