{"id":"3b5feca0-3d56-4895-b284-6724de3e79d3","arxiv_id":"2412.15037","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Axial gravitational quasinormal modes are computed for Schwarzschild and hairy black holes in Einstein-Weyl gravity, showing stability and new massive spin-2 tones.","lead":"Researchers computed the ringdown frequencies of black holes in a modified theory of gravity with extra curvature terms. They found axial stability and new heavy gravity modes that could leave an observable imprint in gravitational wave signals.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Hairy-branch QNM results rest solely on the 4th-order semi-analytic continued-fraction background, with no independent QNM cross-check or error estimate for that branch; this is a real verification gap, though the qualitative stability and spectrum structure likely survive.","rationale":"The reader's weakest-assumption analysis identifies exactly the same load-bearing point: the hairy-branch QNM computation depends entirely on the accuracy of the semi-analytic continued-fraction background, with no independent cross-check. I agree with this assessment. The paper does several things well: it derives the axial perturbation equations, reproduces known Schwarzschild results, cross-validates the Schwarzschild massive modes with both direct integration and continued fractions, shows convergence of the DI method with increasing nmax, and explicitly restricts itself to axial parity and to modes with omega > mu. These are genuine supporting elements. The remaining soft spot is not an internal inconsistency in the derivation but a missing verification step for the hairy branch. The paper itself acknowledges expected numerical errors and disclaims high precision, which tempers but does not remove the concern. Since the reader's conditional verdict already reflects this gap, no change to the verdict is needed. The qualitative claims — axial stability and the presence of massless-led and massive-led mode families — are likely robust to a few-percent background error, but the quantitative ringdown predictions and the precise statement that all hairy modes differ from GR should be read as contingent on the semi-analytic background being accurate enough. The concrete test proposed would settle this by direct comparison with a fully numerical background and by convergence checks on the continued-fraction order.","tokens_in":24803,"tokens_out":7648,"duration_ms":58303,"concrete_test":"For representative parameter values p = 0.95, 1.0, and 1.1, recompute the fundamental l = 2 vector and tensor QNMs using the fully numerical hairy background obtained by the shooting method described in Sec. III C (integrating Eqs. (13)-(14) from the horizon with the exponentially diverging mode suppressed at infinity), instead of the 4th-order continued-fraction fit. Keep the same DI QNM code and boundary-condition treatment, and compare Re(omega rh) and Im(omega rh) with Figs. 5-8. In parallel, increase the continued-fraction order from 4 to 6 and 8 and monitor convergence of the frequencies. If the frequency shifts exceed roughly 2-3% for any mode, the claimed few-percent deviations from GR and the unqualified statement that all hairy modes differ from GR require explicit qualification; if the shifts are negligible, the concern is settled.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims about the hairy-branch spectrum are computed by solving the three-function system (60)-(61) on the semi-analytic background of Appendix A (Eqs. (A1)-(A12)), which is a 4th-order continued-fraction fit from Ref. [22]. Section IV B states explicitly that a continued-fraction QNM method is not straightforwardly applicable to this background and that numerical errors are expected. The DI method is calibrated only on the Schwarzschild case, where the background is exact; that calibration validates the integrator and boundary-condition treatment, but not the accuracy of the approximate hairy background. Because the perturbation matrices P^h and V^h in Eq. (61) involve A, B and their derivatives, any pointwise error in the fit propagates directly into the QNM frequencies. The paper's quantitative phenomenology — massless-led modes differing from GR by a few percent and hence being potentially measurable by third-generation detectors — operates at the same few-percent level, so a background-induced shift of that size could affect whether the predicted deviations are meaningful. The qualitative stability result (no mode with positive imaginary part) is less sensitive to such errors, but it too is obtained only on this approximate background. This is a verification gap rather than demonstrated error, but it is the most load-bearing assumption behind the hairy-branch part of the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies axial gravitational perturbations of static, spherically symmetric black holes in Einstein-Weyl quadratic gravity, using the auxiliary-field formulation of the theory. It derives the linearized perturbation equations for both the Schwarzschild background and the non-Schwarzschild ('hairy') background, reduces them to first-order-in-derivative systems for three perturbation functions, and computes quasinormal-mode frequencies with direct integration and, for Schwarzschild, continued fractions. For Schwarzschild, the massless axial modes are found to coincide with the Regge-Wheeler spectrum of general relativity, while additional vector and tensor modes associated with the massive spin-2 field appear; for the hairy background, the paper reports massless-led and massive-led mode families that are all different from GR, and it finds no unstable axial modes in either family. The paper closes with a phenomenological discussion of whether these extra modes could be observed in ringdown signals.","tokens_in":24980,"tokens_out":8977,"duration_ms":86889,"significance":"If the numerical results hold, this is a useful first computation of gravitational axial quasinormal modes for the hairy black-hole branch of quadratic gravity, and it provides concrete evidence that massive spin-2 degrees of freedom generate additional ringdown content even when the stationary solution is Schwarzschild. The Schwarzschild-sector computation is strengthened by the use of two independent methods (direct integration and continued fractions) with explicit convergence checks, and the QNM calculation has no free parameters fitted to the target frequencies. The main limitations are that the hairy-branch results rest on a single numerical method applied to an approximate background, and that the stability conclusion covers only the ω>μ QNM sector, which the paper itself notes is not the full massive-field spectrum.","major_comments":[{"comment":"The quantitative claims for the hairy background—in particular the mode frequencies shown in Figs. 5–8 and the statement in Sec. V B that deviations of a few percent are measurable by third-generation detectors—rest entirely on the direct-integration method applied to the approximate continued-fraction background of Appendix A, with no independent verification for this branch. Since the coefficient matrices P^h and V^h in Eq. (61) contain A, B and their derivatives, any pointwise error in the fit propagates directly into the QNM frequencies, and the authors themselves note in Sec. IV B that numerical errors are expected in this computation. The few-percent deviations highlighted in the phenomenology are of the same order as that expected uncertainty, so the present calculation does not demonstrate that the hairy-branch deviations are physical rather than artifacts of the approximate background. The frequencies should be recomputed on the fully numerical background, or at least an error estimate should be obtained by varying the truncation order of the continued-fraction fit, before the quantitative phenomenological claims are made.","section":"Sec. IV B and Appendix A"},{"comment":"The stability claim ('no unstable modes were found') covers only the quasinormal sector with ω > μ: the boundary conditions in Eqs. (50)–(52) and (65)–(67) use k = sqrt(ω^2 − μ^2), which restricts the search to ω > μ, while the paper itself notes in Sec. VI that the massive spin-2 equation also admits quasi-bound states with ω < μ. An instability, if present, could live in that sector, so the conclusion that the solutions are stable under axial perturbations is not established by the present calculation. The text should either extend the analysis to the quasi-bound spectrum or explicitly restrict the stability statement to the ω > μ QNM sector.","section":"Secs. IV A, IV B and VI"},{"comment":"The matrix V^h in Eq. (61) is displayed with the (3,2) entry written as V^h_{23}, which duplicates the (2,3) entry; this should presumably be V^h_{32}. More importantly, the derivation leading to Eq. (61) is only sketched: the text states that the constraint equation and the (θφ) component of Eq. (29) are used to eliminate h0 and F0, but the resulting expressions for P^h and V^h are not given. Given that the entire hairy-branch computation depends on these matrices, at least the general structure of the reduction (or a brief verification) should be provided, or the explicit matrices should be included in an appendix.","section":"Sec. IV B, Eq. (61)"}],"minor_comments":[{"comment":"The sentence 'By substituting (47)-(48) into (40)' appears to refer to the wrong equation: Eq. (40) is the two-function Schrödinger-type system for Q and Z, whereas Eqs. (47)-(49) are the horizon expansions for h1, F1 and F2; the substitution should be into the full three-function system (45).","section":"Sec. IV A"},{"comment":"The paragraph on boundary conditions states that the three forward integrations fix the horizon parameters '(h1^{(0)}, F1^{(0)}, F1^{(0)})'; the third entry should be F2^{(0)}, as given earlier in the same subsection.","section":"Sec. IV B"},{"comment":"At the bifurcation point p ≈ 0.876 the hairy solution approaches the Schwarzschild solution, so the statement that 'all modes are different from those of general relativity' should be qualified, and a numerical check that the hairy modes reduce to the Schwarzschild modes as p → p_min would be a useful consistency test.","section":"Sec. V B"},{"comment":"Fig. 1 compares the numerical and semi-analytical metric functions for p = 1.1, but it does not show the pointwise fractional difference; plotting that difference would directly quantify the accuracy of the approximate background used in the QNM computation.","section":"Fig. 1"},{"comment":"A table with representative QNM frequencies and estimated errors for both backgrounds and for several values of p would make the paper substantially more useful for quantitative comparisons than the figures alone.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for a general relativity journal, and I see no issue with the citation pattern or novelty disclosure. My recommendation is driven by verification gaps in the hairy-branch computation and by the restricted domain of the stability claim, both of which are addressable within the scope of the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The genuinely new content is the axial gravitational QNM spectrum for the hairy (non-Schwarzschild) branch of Einstein-Weyl quadratic gravity, including the coupled massless/massive spin-2 system on that background. The Schwarzschild massive spin-2 modes were already computed by Brito-Cardoso-Pani and Ould El Hadj; the paper says so, recovers them, and cross-validates with two independent methods. That part is solid.\n\nThe paper does a lot right. The perturbation equations are derived cleanly, reduced to a three-function ODE system, and the Schwarzschild branch is checked with both direct integration and continued fractions, with the massless modes matching GR. The DI/CF agreement on Schwarzschild calibrates the integrator. The search for unstable axial modes is a genuine stability test, and finding none is a useful result. The phenomenological section is careful: the excitation of the massive modes is explicitly heuristic, and the axial-only restriction is stated up front. The citation of [59] and [73] for the prior Schwarzschild massive modes is proper.\n\nThe soft spot is exactly where the stress-test puts it. The hairy-branch frequencies are computed only on the 4th-order continued-fraction background from [22], with the DI method and no independent cross-check for that branch. The paper admits this: it says numerical errors are expected and the aim is the general structure rather than high-accuracy predictions. That is an honest caveat, but it matters because the few-percent deviations from GR in Figs. 5-8 are in the same range as the expected background error. A background-induced shift of that size could change whether the predicted deviations are meaningful. The qualitative stability result, no positive imaginary parts, is much less sensitive and probably fine. Also, the claim that all modes are different from GR should be read with the bifurcation endpoint in mind: at p about 0.876 the hairy branch merges with Schwarzschild, so the spectra must coincide there; this is a wording slip, not a substantive error.\n\nNo code, data, or tabulated frequencies are provided, which limits reproducibility but is common for this kind of paper. No circularity: the QNMs are computed from derived equations with the coupling as input, not fitted.\n\nWho is this for? People working on modified-gravity ringdown and black hole spectroscopy. It deserves a serious referee; the referee should push for an error estimate on the hairy background, or a check at one or two p values with the numerical background directly, and tabulated frequencies. I would not desk-reject.","headline":"The genuinely new result is the axial QNM spectrum for the hairy branch in quadratic gravity, but the hairy-branch frequencies rest on an approximate background with no independent cross-check; the Schwarzschild massive modes and the qualitative stability result are on firmer ground.","tokens_in":25581,"tokens_out":2778,"would_cite":true,"duration_ms":24043,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C35","83D05"],"pacs":["04.30.-w","04.70.-s","04.50.Kd"],"model":"deepseek-v4-flash","headline":"Quadratic gravity adds massive spin-2 quasinormal modes to black hole ringdowns and leaves both black hole families stable.","keywords":["quasinormal modes","quadratic gravity","Einstein-Weyl gravity","black hole stability","massive spin-2 modes","ringdown","hairy black holes","Schwarzschild black holes"],"falsifier":"Compute the ℓ = 2 axial quasinormal-mode frequencies of the hairy branch on a fully numerical background, rather than the continued-fraction fit, at several values of p in [0.876, 1.143] and compare with the paper's results; disagreement at the few-percent level would indicate that the reported hairy spectrum is an artifact of the approximate background. Alternatively, any found mode with Im(ω) > 0 in this parameter range would refute the stability claim.","tokens_in":24534,"feed_emoji":"🕳️","tokens_out":5486,"duration_ms":47864,"temperature":0.7,"pith_summary":"The paper studies how black holes ring in Einstein-Weyl quadratic gravity, where the Einstein-Hilbert action gains a Weyl-squared term that adds a massive spin-2 field. It computes the axial quasinormal modes of both known static black hole families: the Schwarzschild solution and a 'hairy' solution whose metric differs from Schwarzschild, with the hair here being purely gravitational. The central result is that the perturbation spectrum decomposes into a massless and a massive spin-2 sector: Schwarzschild keeps the usual general-relativity modes and gains new massive vector and tensor modes, while the hairy black hole has massless-led and massive-led modes that all differ from those of general relativity. No unstable axial modes were found in the radially stable parameter range, providing numerical evidence that both families are stable. The authors argue that the massive modes can be excited during ringdown, making them a possible observational signature of quadratic gravity.","feed_headline":"Massive spin-2 modes join black hole ringdowns in quadratic gravity","feed_subtitle":"Schwarzschild and hairy black holes stay stable, and the new modes hand observers a direct test of modified gravity.","key_machinery":"The argument is carried by the auxiliary tensor field fμν, which rewrites the quadratic-gravity action in a form with second-order field equations while keeping the massive spin-2 degree of freedom explicit. Perturbations are expanded in axial tensor spherical harmonics, giving metric functions h0, h1 and massive-field functions F0, F1, F2; after imposing gauge and constraint conditions the system reduces to a coupled set of Schrödinger-type equations for h1, F1, F2. Quasinormal-mode frequencies are found by imposing ingoing waves at the horizon and outgoing waves at infinity, using direct forward integration for both backgrounds and a continued-fraction method as an independent benchmark for the Schwarzschild case. The hairy background itself is represented by a semi-analytical continued-fraction parametrization, which is what makes the full computation feasible.","core_discovery":"In Einstein-Weyl quadratic gravity, axial gravitational perturbations of static, spherically symmetric black holes are governed by a massless spin-2 field, the usual graviton, coupled to a massive spin-2 field. On a Schwarzschild background the massive field decouples, so the spectrum consists of the standard general-relativity quasinormal modes plus two new series of massive modes, classified as vector and tensor. On the hairy background the fields remain coupled, producing massless-led and massive-led modes, and within the explored parameter range all of these frequencies differ from the corresponding general-relativity values. Searching for modes with positive imaginary frequency in the domain where both backgrounds are radially stable, 0.876 ≲ p ≲ 1.143, the authors found none, which they take as strong numerical evidence of axial stability for both families.","pith_inferences":["The polar sector, not computed here, will likely show similar spectral deformation and would bring in the massive scalar degree of freedom from the βR² term, so a complete ringdown template needs both parities.","The heuristic estimate places the massive-mode excitation at order α in the metric, but a field-redefinition argument suggests it might appear only at order α²; computing the actual excitation by an inspiralling particle would settle which order is physical.","The paper notes the massive spin-2 equation also admits quasi-bound states with frequency below μ; those could source metric perturbations and produce a distinct low-frequency ringdown component that the current analysis does not cover.","If hairy black holes are stable, their larger horizon for a given mass may yield measurable differences in tidal deformability or accretion-disk spectra, extending tests beyond the ringdown signal."],"forward_implications":["Schwarzschild black holes in quadratic gravity, even with the same stationary metric as in general relativity, should emit massive spin-2 modes during ringdown, giving a direct observational handle on the theory.","For hairy black holes the massless-led and massive-led mode frequencies shift by up to a few percent relative to general relativity, in principle measurable by third-generation ground-based detectors if such light black holes exist.","Both black hole families appear stable under axial perturbations in the parameter range where they are radially stable, supporting their classical physical viability.","Observable stationary deviations require black hole masses below Mc ≈ 0.438/μ, so current observations favor small couplings and leave a narrow window for detecting hairy solutions.","The massive modes provide a qualitative signature that does not require the stationary solution to differ from general relativity, so ringdown can reveal quadratic gravity even when the spacetime is Schwarzschild."],"supporting_citations":[{"why":"Establishes that static, asymptotically flat solutions of quadratic gravity are Ricci-scalar flat and that the Einstein-Weyl subclass captures the relevant physics.","marker":"[20]"},{"why":"Provides the semi-analytical continued-fraction hairy black hole background used for all hairy-branch quasinormal-mode computations.","marker":"[22]"},{"why":"Time-domain non-radial evolutions that motivate the stability expectation the paper tests for the gravitational sector.","marker":"[25]"},{"why":"Fixes the radial-stability domain 0.876 ≲ p ≲ 1.143 and the mass bound used to define the parameter range.","marker":"[33]"},{"why":"Supplies the framework for extra non-GR ringdown modes that the paper uses to argue the massive modes can be observed.","marker":"[57]"},{"why":"Derives the massive spin-2 axial perturbation equations and the vector/tensor mode decomposition on Schwarzschild that the paper extends.","marker":"[59]"},{"why":"Provides the axial harmonic decomposition and gauge underlying the perturbation scheme.","marker":"[72]"},{"why":"Continued-fraction method used as the benchmark for the Schwarzschild quasinormal-mode frequencies.","marker":"[75]"}],"fun_headline_variants":["Gravitational ringdown gets a massive spin-2 twist","Quadratic gravity adds massive modes to black hole echoes","Hairy black holes ring with extra graviton modes","New spin-2 modes could test modified gravity","Black hole vibrations reveal massive spin-2 partner"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The hairy-background results assume the semi-analytical continued-fraction metric accurately represents the true solution for every p between 0.876 and 1.143, and this assumption is not cross-checked with an independent method.","fun_headline_variants_meta":{"raw":{"variants":["Gravitational ringdown gets a massive spin-2 twist","Quadratic gravity adds massive modes to black hole echoes","Hairy black holes ring with extra graviton modes","New spin-2 modes could test modified gravity","Black hole vibrations reveal massive spin-2 partner"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000439,"raw_usage":{"total_tokens":2216,"prompt_tokens":917,"completion_tokens":1299,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":533,"completion_tokens_details":{"reasoning_tokens":1223}},"tokens_in":533,"tokens_out":1299,"duration_ms":6586,"temperature":1.0,"reasoning_tokens":1223,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:40:54.549857+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the ℓ = 2 axial quasinormal-mode frequencies of the hairy branch on a fully numerical background, rather than the continued-fraction fit, at several values of p in [0.876, 1.143] and compare with the paper's results; disagreement at the few-percent level would indicate that the reported hairy spectrum is an artifact of the approximate background. Alternatively, any found mode with Im(ω) > 0 in this parameter range would refute the stability claim.","supporting_citations":[],"review_version":1}