{"id":"9a043a3f-b7fe-49e8-92be-45a9e65a0f09","arxiv_id":"2610.06370","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Linear scalaron quasinormal and quasi-bound spectra of Schwarzschild black holes in metric f(R) gravity depend only on the combination μM, and are recovered independently by continued-fraction and characteristic time-domain methods.","lead":"The paper derives the equation ruling the scalaron (the extra scalar mode of metric f(R) gravity) on a Schwarzschild black hole directly in the Bondi-Sachs null-cone formulation, and shows the ringdown frequencies depend on the theory only through the product of black hole mass and scalaron mass. It then recovers those frequencies two independent ways, and demonstrates pulling both the fast-decaying and the long-lived modes out of a single time-domain waveform.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek/deepseek-v4-flash","headline":"QBS time-domain recovery: cluster MADs ~4e-14 vs TD–CF offsets ~4e-8 (six orders of magnitude), a systematic bias beyond the reported stability, untested for the headline ℓ=2, μM=0.575 simultaneous-extraction case.","rationale":"Central claim analysis: the scalaron obeys the massive scalar master equation with μ² = f'(0)/(3f''(0)); the derivation (trace of eq. (2) linearized, giving (12)–(13)) is textbook and internally consistent. Eq. (25) follows from (20e) by multiplication by r², verified term by term. The tensor reduction via (20b), (20d), ∂r(20d) and the compatibility identity (22) reduces to the GR form (24), with the β,R terms collapsing to the combination (23), which vanishes by (20a) — checked. The μM-degeneracy follows by the scaling symmetry of (25) and is physically correct. The frequency-domain results are validated against independent literature at the 10^-6 level (Table II), the 10^-8 level for QBS (Table III), and 10^-15 for the tensor Nollert check (Table IX). The paper's own limitation statements (b84, b99, b227, b438–b442, b451–b452, b473) are accurate and already cover the reader's weakest_assumption: the analytic f(0)=0, f'(0)≠0, f''(0)>0 class restriction is acknowledged and does not constitute an unresolved flaw, since the resulting equation is benchmarked. The only live concern I found, which the reader did not flag, is the systematic offset of the time-domain quasi-bound extraction relative to its own extremely tight cluster MAD (Table VII ℓ=1: 6.31e-8 offset vs 2.5e-13 MAD; Table VIII ℓ=2: 4.46e-8 offset vs 4.1e-14 MAD), combined with the absence of a robustness scan for the headline ℓ=2, μM=0.575 simultaneous case (Table VI covers only ℓ=1, μM=0.15 QNM). Because the reported TD–CF differences are honest (0.63% QNM damping, 1.2e-3% QBS damping relative errors) and small, the demonstrated recovery claim stands; the concern bears on the interpretation of MADs as stability rather than accuracy, for which the paper already provides the relevant caveat (b441). The concrete test above would settle whether the bias is a finite-domain artifact (case b) or an effective-pole effect requiring restatement (case c). On balance the verdict ACCEPT is unchanged.","tokens_in":36683,"tokens_out":19373,"duration_ms":135892,"concrete_test":"Run a robustness scan for the headline simultaneous-extraction case (ℓ=2, μM=0.575, robs≈8M) analogous to Table VI, varying at minimum the outer domain rmax (2000M, 4000M, 8000M), the inner worldtube xmin (−30, −35, −40), and the initial-data pulse (r0/σ), with Nx and Δu chosen so the reflection time urefl of eq. (55) lies beyond each QBS fit window. Track the extracted QBS pole. (a) If it shifts by much more than its cluster MAD while remaining tight, the observed 4.46e-8 TD–CF offset is a systematic finite-domain bias. (b) If the pole converges toward the continued-fraction value 2Mω_QBS = 1.12493282 − 5.41692409e-5i as rmax and resolution increase, the recovery claim is quantitatively confirmed and the bias is attributed to domain truncation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In the reader's strongest_claim, the characteristic evolution 'independently recovers both mode families from a single waveform.' This is well established for QNMs: the ℓ=1, μM=0.15 pole shifts by ≤1.72e-4 under variations of domain, dissipation, extraction radius and initial data (Table VI), and TD–CF offsets such as 1.4e-5 for (1,0.15) sit inside the cluster MADs (~1.7e-4, Table V). The quasi-bound recovery shows a different pattern that the paper does not explain and does not scan. In Table VII (ℓ=1, μM=0.25) the accepted QBS cluster has MADs (6.22e-13, 2.53e-13) yet the TD median differs from the continued-fraction eigenfrequency by 6.31e-8. In the headline simultaneous case (Table VIII, ℓ=2, μM=0.575) the QBS MADs are (3.7e-14, 4.1e-14) while the TD–CF difference is 4.46e-8 — about six orders of magnitude above the cluster's own spread. A matrix-pencil fit to a window dominated by one damped exponential can lock onto that exponential almost exactly, so an extremely tight cluster certifies only that the late-time waveform is one clean damped exponential; it does not certify accuracy. The systematic offset (relative ~1.2e-5 in the damping) is consistent with a finite-domain 'effective-pole' effect of the initial–worldtube problem. Table VI's robustness scan covers only the ℓ=1, μM=0.15 QNM, not the ℓ=2, μM=0.575 QBS that carries the headline claim of simultaneous recovery, so the magnitude and origin of this bias — and whether it converges away under refinement — is untested. This does not overturn the physics: the massive-scalar identification is independently benchmarked against the literature (Table II, III, IX), and eq. (25) is confirmed to be the standard massive-scalar Schwarzschild equation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek/deepseek-v4-flash","summary":"This paper studies scalaron quasinormal modes and quasi-bound states of Schwarzschild black holes in metric f(R) gravity using the Bondi-Sachs (characteristic) formulation. The authors linearise the f(R) characteristic equations about a Ricci-flat Schwarzschild background and obtain (i) a closed massive-scalar master equation (25) for the Ricci-scalar perturbation, with effective mass mu^2 = f'(0)/(3 f''(0)), and (ii) a tensor master equation (24) that reduces to the standard GR form. They argue that for analytic models with f(0)=0 and f'(0) != 0, the linear scalar spectra depend on the theory only through the dimensionless product mu M. In the frequency domain they implement Leaver-type continued-fraction pole conditions (QNM and quasi-bound branches, distinguished by kappa-sheet) and validate against literature benchmarks (Tables II-III). In the time domain they evolve the scalaron equations with a fourth-order characteristic scheme, verify fourth-order convergence (Table IV), extract complex frequencies with SVD-truncated matrix-pencil stability scans, and demonstrate recovery of QNM and quasi-bound frequencies, including a simultaneous extraction of both mode families from a single waveform at one extraction radius for the case (ell=2, mu M=0.575).","tokens_in":37085,"tokens_out":7309,"duration_ms":50827,"significance":"If the results hold, the paper makes several solid contributions. The explicit reduction of the linearised Bondi-Sachs f(R) system to the massive-scalar master equation, with the parameter-free degeneracy statement that two analytic models sharing f'(0)/f''(0) produce identical linear scalar spectra at fixed M, is clean and testable. The tensor-sector reduction to GR is a useful formalism-level result. The strength of the paper is its validation stack: the continued-fraction implementation reproduces literature QNM benchmarks at the ~1e-6 level or better across massless, low-mass, and moderate-mass regimes (Table II), the time-domain scheme shows verified fourth-order convergence (p = 4.01), the matrix-pencil selection is explicitly insulated from the continued-fraction roots (CF values introduced only after cluster selection), and the paper is notably transparent about the limits of its own extraction procedure: MADs measure fit sensitivity, not total uncertainty; pseudo-spectral stability is outside scope; and the QBS damping is not always resolvable. These honestly stated limitations are part of why the central claims are credible.","major_comments":[{"comment":"The quasi-bound time-domain medians exhibit systematic TD-CF offsets (6.31e-8 for (ell=1, muM=0.25) in Table VII; 4.46e-8 for the headline (ell=2, muM=0.575) case in Table VIII) that exceed the corresponding cluster MADs ((2.53e-13, 6.22e-13) and (3.7e-14, 4.1e-14), respectively) by roughly six orders of magnitude. As the skeptic's note observes, an extremely tight matrix-pencil cluster certifies only that the analysed window is one clean damped exponential; it does not by itself certify accuracy, and the offset is consistent with a finite-domain effective-pole shift of the initial-worldtube problem. The robustness scan of Table VI covers only the (ell=1, muM=0.15) QNM and does not test the (ell=2, muM=0.575) QBS that carries the simultaneous-recovery claim. I request either (i) a domain/refinement scan for the QBS extraction (vary rmax, robs, xmin, and the fitting-window placement) showing whether the offset converges away under refinement, or (ii) an explicit interpretation of the offset as a finite-domain bias with a quantitative estimate of its size. I regard this as a verification request rather than a challenge to the central derivation: the absolute offset (4.46e-8 in 2M omega) is far below physically relevant scales, and the paper already reports these numbers honestly.","section":"Sec. VIII C 4, Tables VII-VIII","authors_acknowledge":"More generally, the reported cluster MADs measure sensitivity to the extraction choices, not total numer- ical uncertainty, and should be interpreted accordingly.","authors_acknowledge_p":0.61}],"minor_comments":[{"comment":"The caption states 'These are the ell = 0 branches shown in Figure 1', but the figure displays ell = 1 branches; it should read 'ell = 1 branches shown in the right panel of Figure 1'.","section":"Fig. 3 caption","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The tensor reduction is presented by asserting that 'standard elimination' yields the solved U, U,r expressions and that all beta- and R-terms collapse into the combination of Eq. (23), which vanishes by Eq. (20a). Since this is the one non-obvious algebraic step in the paper, I suggest including the solved expressions for U and U,r (even in an appendix) so that the GR reduction to Eq. (24) can be verified directly.","section":"Sec. IV A, Eq. (22)-(23)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Two distinct monopole quasi-bound entries at muM = 0.25 appear (Table III: 0.4954575282 - 0.0012730497i; Table VII: 0.48235028 - 0.0090706316i). The footnote clarifies they are different modes, but the overtone / principal-quantum-number labelling (n = 0 vs. N = ell+n+1) should be stated explicitly at first use to prevent reader confusion.","section":"Tables III and VII","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The 1.56e-5 difference for (ell=0, muM=0.05) is explained in Appendix C as a finite-order Hill-determinant value in Ref. [5], with the archived Leaver implementation converging to the present value. This accounting is adequate, and the residual is properly not treated as a discrepancy.","section":"Sec. VIII A, Table II, row (0, 0.05)","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is within scope for this journal and the central claims — the massive-scalar master equation, the tensor GR reduction, and the muM degeneracy — are sound and well validated in the frequency domain. The reader's ACCEPT assessment is, in my view, correct; the skeptic's QBS concern is the one point that deserves author attention, but the reported absolute accuracy (4.46e-8) is far below any physical relevance and the paper is transparent about its extraction statistics. The requested QBS scan is a local addition in Sec. VIII C. The Fig. 3 caption error suggests a light proofread before resubmission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek/deepseek-v4-flash","letter":"What to know: the paper does what it says. The scalaron master equation (25) drops out cleanly from the linearised Bondi–Sachs f(R) system, the tensor sector reduces to the GR form, and the frequency-domain spectra are the known massive-scalar-on-Schwarzschild values parameterised by μM. The literature agreement is genuine — roughly 1e-6 or better in Table II, with the QBS checks in Table III as good as the printed reference precision allows. The characteristic scheme shows true fourth-order convergence (p = 4.01), the continued-fraction work is careful with Nollert termination, and the μM degeneracy is correctly framed as a linear, homogeneous, Ricci-flat result that the authors themselves note higher-order terms and curved backgrounds will break (b451–b452). Time-domain extraction is honestly independent: continued-fraction roots enter only after cluster selection (b357–b366).\n\nThe stress-test note lands on the QBS side, and it is the one soft spot worth taking seriously. In Table VII the ℓ=1, μM=0.25 cluster has MADs of order 1e-13 while its median sits 6.3e-8 from the continued-fraction pole; the headline ℓ=2, μM=0.575 simultaneous case shows the same structure — MADs of 4e-14 against a TD–CF separation of 4.5e-8, roughly 1.2e-5 relative in the damping. The authors say the MADs measure extraction sensitivity, not total uncertainty (b441), which is fair as a caveat but does not answer the pattern. Every QBS case shows the same systematic offset, consistent with a finite-domain effective-pole bias from the initial–worldtube problem, and nothing in the paper explains it or tests whether it converges away under refinement. Table VI's robustness scan covers only the ℓ=1, μM=0.15 QNM — not the quasi-bound pole that carries the simultaneous-recovery claim. To be clear: a tight matrix-pencil cluster certifies that the late-time waveform is one clean damped exponential; it does not certify the accuracy of the pole.\n\nProportion: none of this overturns the physics. The massive-scalar identification is independently benchmarked, and the QBS frequencies themselves come from a well-tested continued-fraction solver. The time-domain QBS section should be read as an illustration of feasibility, not a precision result — the authors mostly frame it that way (b438–b440) — but the claim that time-domain medians agree with continued-fraction values \"within the measured extraction-stability scales\" (b366) does not survive scrutiny for the QBS rows, where cluster spread and TD–CF offset differ by six orders of magnitude.\n\nWho this is for: people working in characteristic-formulation perturbation theory and f(R) black-hole spectra. It deserves a serious referee; the main request should be a convergence/domain-dependence study of the extracted QBS pole and an explicit statement that the tight clusters certify waveform cleanliness, not pole accuracy. I would cite the master-equation derivation and the frequency-domain benchmarks, not the QBS extraction numbers.","headline":"A careful, honest paper whose Bondi-Sachs scalaron derivation and QNM benchmarks are solid; the real soft spot is the quasi-bound time-domain recovery, where ultra-tight clusters sit six orders of magnitude below a systematic TD–CF offset that no scan explains.","tokens_in":37954,"tokens_out":3074,"would_cite":true,"duration_ms":44322,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":["83C57","83D05","83C25","83C35"],"pacs":["04.70.-s","04.50.Kd","04.30.-w"],"model":"deepseek/deepseek-v4-flash","headline":"In metric f(R) gravity, the scalaron around a Schwarzschild black hole is a massive scalar with effective mass $\\mu^2 = f'(0)/(3f''(0))$, so the entire linear quasinormal and quasi-bound spectrum is fixed by the dimensionless product…","keywords":["scalaron","metric f(R) gravity","Schwarzschild black hole","quasinormal modes","quasi-bound states","Bondi-Sachs characteristic formulation","continued-fraction method","massive scalar field perturbation"],"falsifier":"Solve the scalar master equation (25) for the $\\ell=1$, $\\mu M=0.15$ fundamental mode with an independent eigenvalue method and compare against the paper's continued-fraction value $2M\\omega \\approx 0.606080 - 0.183041\\,i$; then re-solve with a second analytic model having the same $f'(0)/f''(0)$ but a non-zero third derivative $f'''(0)$, and verify the root is unchanged. A disagreement beyond roughly $10^{-5}$ in $2M\\omega$ would signal a solver error; a root shift depending on $f'''(0)$ at first order would falsify the claimed degeneracy, while a shift appearing only at second order would confirm the authors' stated first-order scope.","tokens_in":38250,"feed_emoji":"🕳️","tokens_out":23047,"duration_ms":248528,"temperature":0.7,"pith_summary":"The paper sets out to show that in metric f(R) gravity the extra scalar degree of freedom — the scalaron — ringing around a Schwarzschild black hole is not a new kind of perturbation: it obeys the well-studied massive-scalar master equation, with effective mass $\\mu^2 = f'(0)/(3f''(0))$ fixed only by the theory's first two Taylor coefficients at zero curvature, while the tensor perturbations obey the same master equation as in general relativity. The central consequence is a degeneracy: two different analytic f(R) theories with the same ratio $f'(0)/f''(0)$ have identical linear quasinormal-mode and quasi-bound-state spectra, so the scalar sector of any such model is the universal massive-scalar spectrum of Schwarzschild parameterised by the dimensionless product $\\mu M$. Working entirely inside the Bondi–Sachs null-hypersurface formulation, the authors compute both spectral branches with the continued-fraction method, validate them against independent published frequencies, and recover them from fourth-order-convergent characteristic time evolutions by scanning waveforms with matrix-pencil stability analysis. If this is right, ringdown calculations for an entire family of f(R) models collapse to a single effective-mass computation, and distinguishing individual theories at linear order by scalar spectra alone becomes impossible.","feed_headline":"Scalaron spectra of Schwarzschild collapse to one effective mass","feed_subtitle":"One ratio of the theory, f'(0)/f''(0), fixes every linear scalaron frequency; higher-order terms drop out.","key_machinery":"The load-bearing identity is the linearised trace equation $3 f''(0)\\,\\Box R - f'(0)\\,R = 0$, which rearranges into the Klein–Gordon equation $(\\Box - \\mu^2)R = 0$ with $\\mu^2 = f'(0)/(3f''(0))$; this single relation converts the scalaron into a massive scalar on the fixed Schwarzschild background, supplies the closed scalar master equation (25), and makes every computed linear frequency depend on the theory only through $\\mu M$. The companion elimination treats the metric-variable subsystem (20b)–(20d) as an algebraic system for $U$ and its derivatives; imposing the compatibility condition $\\partial_r U - U_{,r} = 0$ collapses, via the hypersurface relation (20a), into the closed tensor master equation (24) in its GR form, so the spin-2 spectrum is untouched by f(R) at this order. Numerically, the scalar QNM and quasi-bound conditions are the vanishing of the two Wronskians $W_+(\\rho)$ and $W_-(\\rho)$ on their respective $\\kappa$-sheets, enforced through three-term continued-fraction recurrences (the standard series-expansion eigenvalue technique, with a large-overtone termination); the independent time-domain check is a fourth-order Runge–Kutta evolution of the reformulated first-order system on a logarithmic radial grid, with Euler–Maclaurin-corrected trapezoidal integration outward from an inner worldtube placed exponentially close to the horizon and with no spectral branch imposed at the outer boundary.","core_discovery":"The core discovery, on the paper's own terms, is a reduction: linearising the Bondi–Sachs characteristic equations of metric f(R) gravity about a Ricci-flat Schwarzschild background splits the system into a closed massive-scalar problem for the Ricci-scalar perturbation $R$ and a tensor problem that reduces to general relativity's. The trace equation gives $(\\Box - \\mu^2)R = 0$ with $\\mu^2 = f'(0)/(3f''(0))$, so the scalar master equation (25) is exactly the massive scalar equation on Schwarzschild; eliminating the metric variables and using the hypersurface relation yields the tensor master equation (24) in its standard GR form. Because coefficients of $R^3$ and higher never enter at first order, the scalar quasinormal-mode and quasi-bound frequencies depend on the theory only through $\\mu M$: at fixed black-hole mass $M$, models sharing $f'(0)/f''(0)$ have identical scalar spectra. The numerically realised spectra are therefore presented as universal massive-scalar Schwarzschild spectra parameterised by $\\mu M$; both the outgoing QNM branch (damped oscillations radiating to infinity) and the decaying quasi-bound branch (long-lived, spatially localised states) are computed as Wronskian zero conditions $W_\\pm(\\rho)=0$ enforced through the continued-fraction method, and a characteristic time evolution with matrix-pencil extraction independently recovers both branches from a single waveform, without imposing either spectral branch at the outer boundary.","pith_inferences":["If the degeneracy holds, scalaron ringdown templates can be shared across entire f(R) model families at linear order; telling individual theories apart would then require sourced perturbations, nonlinear scalaron–tensor coupling, or a non-Ricci-flat background — the directions the authors flag in their concluding remarks.","The demonstrated single-waveform, single-radius recovery of both mode families implies a reusable extraction protocol for massive fields in null-characteristic codes: tune early windows to the ringdown, late windows to the long-lived mode, and let persistence across window, rank, and decimation variations certify the poles — a procedure that could extend directly to massive vector or Dirac perturb","The unresolved $\\ell=2$, $\\mu M=0.25$ quasi-bound decay rate (frequency recovered, damping not) points to a general measurement limit: for extremely long-lived modes a finite waveform constrains the oscillation accurately but the lifetime only weakly, which bounds how finely quasi-bound lifetimes can be probed."],"forward_implications":["Any two analytic f(R) models with the same ratio $f'(0)/f''(0)$ predict identical linear scalar quasinormal and quasi-bound spectra on a Schwarzschild background, so linear-order scalaron ringdown data alone cannot distinguish them.","The tensor master equation reduces exactly to its general-relativistic form, so the spin-2 QNM spectrum of Schwarzschild in metric f(R) gravity is the standard GR spectrum; modified-gravity signatures must show up first in the scalar polarisation.","Along each overtone branch, increasing $\\mu M$ drives the damping toward the quasi-resonant limit $-\\operatorname{Im}(2M\\omega)\\to 0^-$ while leaving the oscillation frequency non-zero, so scalaron modes become arbitrarily long-lived at large effective mass.","In the small-coupling regime the quasi-bound damping rates tend to the matched-asymptotic hydrogenic formula, giving an analytic anchor for the decaying-branch roots as $\\mu M\\to 0$.","A single characteristic waveform at one extraction radius can deliver both the fundamental QNM and quasi-bound complex frequencies from separate fitting windows; the demonstrated $\\ell=2$, $\\mu M=0.575$ case matches the continued-fraction values to $7.35\\times10^{-4}$ (QNM) and $4.46\\times10^{-8}$ (QBS) in $2M\\omega$."],"supporting_citations":[{"why":"Provides the full non-linear Bondi–Sachs characteristic system for metric f(R) gravity that the paper linearises about Schwarzschild; the starting point for both master equations.","marker":"[61]"},{"why":"Supplies the elimination procedure and the auxiliary variable that yield the tensor master equation (24) in the Bondi–Sachs GR framework, to which the f(R) tensor sector is shown to reduce.","marker":"[62]"},{"why":"The continued-fraction, three-term-recurrence eigenvalue method used to impose the QNM and quasi-bound endpoint conditions in the frequency domain.","marker":"[52]"},{"why":"The improved large-overtone termination that stabilises the backward-evaluated continued fractions for highly damped modes.","marker":"[63]"},{"why":"The matched-asymptotic hydrogenic quasi-bound spectrum used as the small-coupling analytic limit the decaying-branch roots must approach.","marker":"[35]"},{"why":"Supplies the low-mass (μM = 0.05) scalar QNM benchmarks for ℓ = 0,...,4 against which the continued-fraction frequencies are validated.","marker":"[5]"},{"why":"Supplies the massless scalar Schwarzschild QNM benchmarks used to fix the frequency convention and anchor the mode sequences.","marker":"[26]"},{"why":"Supplies the ℓ = 1 finite-mass QNM benchmarks at μM = 0.10 and 0.20 in the matched time-domain study.","marker":"[36]"},{"why":"Supplies the ℓ = 2 finite-mass QNM benchmarks at μM = 0.10, 0.15 and 0.20.","marker":"[51]"},{"why":"Supplies the quasi-bound-state reference frequencies at μM = 0.25 used to validate the decaying-branch continued fraction.","marker":"[31]"}],"fun_headline_variants":["Scalaron spectra hinge on one ratio μ M","f(R) scalaron modes collapse to μ M alone","Schwarzschild scalaron frequencies scale with μ M","Bondi–Sachs scalaron spectra depend only on μ M","One mass term sets all scalaron QNM and bound states"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The whole spectral equivalence rests on the linearised trace equation $3f''(0)\\,\\Box R - f'(0)\\,R = 0$, which presumes $f$ analytic near $R = 0$ with $f(0) = 0$, $f'(0) \\ne 0$ and $f''(0) > 0$, and which discards all $O(\\epsilon^2)$ terms — relax any of these conditions and the massive-scalar identification, the $\\mu M$ degeneracy, and the GR tensor reduction no longer follow.","fun_headline_variants_meta":{"raw":{"variants":["Scalaron spectra hinge on one ratio μ M","f(R) scalaron modes collapse to μ M alone","Schwarzschild scalaron frequencies scale with μ M","Bondi–Sachs scalaron spectra depend only on μ M","One mass term sets all scalaron QNM and bound states"]},"model":"deepseek/deepseek-v4-flash","effort":"low","cost_usd":0.000191,"raw_usage":{"cost":0.000190568,"is_byok":false,"cost_details":{"upstream_inference_cost":0.000190568,"upstream_inference_prompt_cost":0.000124768,"upstream_inference_completions_cost":6.58e-05},"total_tokens":1331,"prompt_tokens":1096,"completion_tokens":235,"prompt_tokens_details":{"audio_tokens":0,"video_tokens":0,"cached_tokens":256,"cache_write_tokens":0},"completion_tokens_details":{"audio_tokens":0,"image_tokens":0,"reasoning_tokens":152}},"tokens_in":1096,"tokens_out":235,"cost_source":"openrouter_usage","duration_ms":6234,"image_pages":[],"temperature":1.0,"billed_attempts":1,"reasoning_tokens":152,"transport_provider":"openrouter","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-06T13:54:38.626425+00:00","model_set":{"reader":"deepseek/deepseek-v4-flash"},"falsifier":"Solve the scalar master equation (25) for the $\\ell=1$, $\\mu M=0.15$ fundamental mode with an independent eigenvalue method and compare against the paper's continued-fraction value $2M\\omega \\approx 0.606080 - 0.183041\\,i$; then re-solve with a second analytic model having the same $f'(0)/f''(0)$ but a non-zero third derivative $f'''(0)$, and verify the root is unchanged. A disagreement beyond roughly $10^{-5}$ in $2M\\omega$ would signal a solver error; a root shift depending on $f'''(0)$ at first order would falsify the claimed degeneracy, while a shift appearing only at second order would confirm the authors' stated first-order scope.","supporting_citations":[{"cited_title":"Characteristic Formulation for Metric $f(R)$ Gravity","cited_arxiv_id":"1707.05757","evidence_quote":"Provides the full non-linear Bondi–Sachs characteristic system for metric f(R) gravity that the paper linearises about Schwarzschild; the starting point for both master equations."},{"cited_title":"Quasinormal modes of a Schwarzschild black hole within the Bondi-Sachs framework","cited_arxiv_id":"2407.06636","evidence_quote":"Supplies the elimination procedure and the auxiliary variable that yield the tensor master equation (24) in the Bondi–Sachs GR framework, to which the f(R) tensor sector is shown to reduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The continued-fraction, three-term-recurrence eigenvalue method used to impose the QNM and quasi-bound endpoint conditions in the frequency domain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The improved large-overtone termination that stabilises the backward-evaluated continued fractions for highly damped modes."},{"cited_title":"Detweiler","cited_arxiv_id":null,"evidence_quote":"The matched-asymptotic hydrogenic quasi-bound spectrum used as the small-coupling analytic limit the decaying-branch roots must approach."},{"cited_title":"Critical masses and numerical computation of massive scalar quasinormal modes in Schwarzschild black holes,","cited_arxiv_id":null,"evidence_quote":"Supplies the low-mass (μM = 0.05) scalar QNM benchmarks for ℓ = 0,...,4 against which the continued-fraction frequencies are validated."},{"cited_title":"Scalar/Dirac perturbations of the Reissner-Nordstr\\\"om black hole and Painlev\\'e transcendents","cited_arxiv_id":"2109.06929","evidence_quote":"Supplies the massless scalar Schwarzschild QNM benchmarks used to fix the frequency convention and anchor the mode sequences."},{"cited_title":"Instability of the massive Klein-Gordon field on the Kerr spacetime","cited_arxiv_id":"0705.2880","evidence_quote":"Supplies the ℓ = 1 finite-mass QNM benchmarks at μM = 0.10 and 0.20 in the matched time-domain study."},{"cited_title":"Quasinormal modes of massive fermions in Kerr spacetime: Long-lived modes and the fine structure","cited_arxiv_id":"1712.06667","evidence_quote":"Supplies the ℓ = 2 finite-mass QNM benchmarks at μM = 0.10, 0.15 and 0.20."},{"cited_title":"Quasibound state reminiscent in de-Sitter black holes: Quasinormal modes and the decay of massive fields","cited_arxiv_id":"2401.15156","evidence_quote":"Supplies the quasi-bound-state reference frequencies at μM = 0.25 used to validate the decaying-branch continued fraction."}],"review_version":1}