{"id":"80c32682-4c26-4431-950f-73e89ae33737","arxiv_id":"2509.04001","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For black holes embedded in Hernquist dark matter halos, the shadow radius and quasinormal mode frequencies are redshifted by a factor 1 - C + C^2/6 in the halo compactness C, with EHT observations implying C <= 0.092.","lead":"This paper computes how a Hernquist-type dark matter halo changes the shadow size and gravitational-wave ringdown frequencies of a supermassive black hole. It finds a simple formula for the frequency shift and an upper bound on the halo compactness from Event Horizon Telescope observations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (5.1) is a leading-eikonal formula; finite-l corrections make it inaccurate for l=2 unless tested at ε→0.","rationale":"The reader's weakest assumption identifies the same structural issue: the eikonal correspondence is used to promote Eq. (5.1) to a statement about low-multipole QNMs. My closer reading of Table 1 and Appendix B shows that the finite-l corrections omitted from Eq. (5.1) produce a ~3% discrepancy at l=2, C=0.3, ε=0.1; this is larger than the numerical/WKB error shown in Figure 3. This does not undermine the shadow bound C≤0.092, which is derived directly from bc in Sec. 6. It does, however, mean the QNM redshift formula should be labeled as a leading-eikonal approximation with a quantified error, or extended to include the c1/κ corrections, before being used as a precise prediction for ringdown. The paper's own limitation that small ε requires very large N and that tables are only given for ε=0.1 strengthens this concern. The verdict CONDITIONAL is therefore appropriate; I would not change it, but I would request the finite-l error quantification as a condition.","tokens_in":12256,"tokens_out":9968,"duration_ms":94032,"concrete_test":"Using the coefficients in Appendix B, evaluate Eq. (4.31) at ε=0 (or ε=10^-3), l=2 (κ=2.5), n=0, and form the ratio Re ω(C,ε)/Re ω(0,0). Plot against 1−C+C^2/6 for C=0.1, 0.2, 0.3. If the corrected ratio deviates from Eq. (5.1) by more than ~2%, Eq. (5.1) is not an accurate description of the l=2 ringdown spectrum; if it agrees to <1%, the eikonal extrapolation is validated. For a fully independent check, run a Leaver continued-fraction solver at the same parameters to confirm the pseudospectral/matrix results.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central QNM claim—Eq. (5.1), ω(C,ε)/ω(0,0)=3√3 M_BH/b_c=1−C+C^2/6—is obtained from the eikonal limit (4.31) by keeping only bcω=κ and neglecting the c1/κ, c3/κ^3, c5/κ^5 terms (Sec. 4.3, Appendix B). For the gravitational-wave-relevant l=2 mode, κ=2.5, so these terms are not negligible. Table 1 (ε=0.1, C=0.3, l=2, n=0) gives the 6th-order eikonal value Reω=0.278094 against the pseudospectral/matrix value 0.276139—a 0.7% discrepancy—whereas Eq. (5.1) predicts a ratio 0.715, while the numerical ratio is 0.739, a ~3% difference. The paper validates Eq. (5.1) visually for ε=10^-1 to 10^-3 and C≤0.5, but the actual numerical tables use only ε=0.1; the ε→0 limit is taken analytically, not checked. If the finite-l correction does not vanish smoothly as ε→0 at large C, then Eq. (5.1) is not a statement about the l=2 ringdown spectrum but an eikonal approximation, and the advertised 'redshift relative to Schwarzschild QNMs' needs a stated error budget. The shadow bound C≤0.092 is independent and unaffected.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a spherically symmetric black hole surrounded by a Hernquist-type dark matter halo, the Schwarzschild-Hernquist solution of Cardoso et al. (2021). It computes the photon sphere and shadow radius as functions of halo compactness C and mass ratio ε, giving small-C analytic expansions. It then computes axial gravitational quasi-normal modes (QNMs) with pseudospectral, matrix, and 6th-order WKB methods, and proposes a redshift formula ω(C,ε)/ω(0,0)=1−C+C²/6+O(C³) in the eikonal and ε→0 limits. Using the EHT shadow-size uncertainty of about 10%, it derives an upper bound C≤0.092.","tokens_in":12642,"tokens_out":10017,"duration_ms":88670,"significance":"If the main claims hold, the paper provides a simple, potentially observable mapping from halo compactness to shadow radius and QNM frequency shift, including a concrete astrophysical bound from EHT observations. The numerical work is carefully cross-validated: pseudospectral and matrix methods agree to six digits, and WKB matches numerical results for large ℓ. The analytic expansions for the photon sphere and impact parameter are useful and reduce correctly to known Schwarzschild limits. The shadow bound is a direct, falsifiable prediction. However, the QNM redshift formula is derived in the eikonal limit, and its accuracy for low multipoles is not quantified; this limits the strength of the claims about observable ringdown frequencies.","major_comments":[{"comment":"The redshift formula is obtained from the eikonal limit (4.31) together with the small-C expansion of the critical impact parameter, dropping all finite-ℓ terms (c1/κ, c3/κ³, ...). At the observationally relevant ℓ=2 mode, these terms are not negligible: for C=0.3, Table 1 gives the 6th-order eikonal Reω=0.278094 versus the numerical value 0.276139 (~0.7% error), while Eq. (5.1) predicts a ratio ω(C)/ω(0)=0.715 compared with the numerical ratio 0.739 (~3% error). The paper states that Eq. (5.1) is \"sufficiently accurate up to C≤0.3\" without providing an error budget or clearly stating that it is an eikonal-limit prediction. Since the abstract and Section 6 present this as the redshift of \"the QNMs\" relative to Schwarzschild, the statement can be misread as applying to the ℓ=2 ringdown spectrum. Please either explicitly identify Eq. (5.1) as an eikonal-limit result with a quantitative err","section":"Section 5, Eq. (5.1)"},{"comment":"The paper claims to calculate axial gravitational QNMs \"up to C∼O(1)\", but Tables 1 and 2 only present results for C≤0.3, and Figure 3 stops at C=0.5. No numerical data are shown in the claimed O(1) range. Given the authors' own statement that the eikonal approximation \"works bad when C>0.3,\" the claim of coverage up to O(1) is not supported by the evidence shown. Please either provide the missing results (e.g., a table for C=0.5, 0.7, 1.0) or amend the claim to reflect the actually presented range.","section":"Abstract and Section 6"}],"minor_comments":[{"comment":"The word \"fit\" is misleading: Eq. (5.1) is derived analytically from the eikonal limit and the impact-parameter expansion, not fitted to numerical data. Please rephrase, e.g., \"we derive\" or \"we obtain\".","section":"Abstract and Section 5"},{"comment":"The formula for Υ is typeset ambiguously: \"arctan r + a0 − MDM√MDMξ\" should have parentheses around the arctangent argument. This is a readability issue.","section":"Eq. (2.8)"},{"comment":"The symbol C is used both for the compactness and as the constant term in the cubic equation (3.2). This can confuse the reader; consider renaming one of them, e.g., using C₀ or a different letter for the cubic coefficient.","section":"Eqs. (2.11) and (3.2)"},{"comment":"Several typos: \"Noth\" should be \"Note\", \"alow\" should be \"allow\", \"paseudospectral\" should be \"pseudospectral\", \"quantity\" should be \"quantify\", \"Schwarschild\" should be \"Schwarzschild\".","section":"Throughout"},{"comment":"The bottom panels show ratios of QNMs to Schwarzschild for ε=10⁻¹, 10⁻², 10⁻³, but the caption does not specify the multipole ℓ and overtone n for these panels. The top caption indicates ℓ=2 and distinguishes n=0 and n=1 by marker shapes; please clarify that the bottom panels follow the same convention.","section":"Figure 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope and the numerical methods are solid. The main concern is that the central QNM redshift claim is an eikonal-limit result but is presented without an error budget for low multipoles, and the abstract overstates the parameter range. These issues are fixable with additional clarification and possibly a short quantitative table. I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper correctly computes the Schwarzschild-Hernquist shadow and QNMs, and the genuinely new item is the second-order compactness correction to the ringdown redshift, omega/omega_Schw = 1 - C + C^2/6, plus an EHT-derived bound C <= 0.092. That bound is robust. The redshift formula, however, is a leading-eikonal statement, and the paper's own l=2 numbers show it drifting to about 3% off at C=0.3. The abstract calls it a fit, which is wrong—it is analytic from b_c—and the claim of validity up to C<=0.3 needs a stated error bar.\n\nNow the details. The spacetime is from Cardoso et al., and the shadow/QNM phenomenology has been covered in several earlier papers. What's new is the explicit second-order term C^2/6, which extends the leading-order result of Refs. [29,30], and the numerical cross-checks between pseudospectral, matrix, and WKB methods that agree to six digits. The analytic expansions for the photon sphere, impact parameter, and the discriminant analysis are internally consistent. The shadow bound from EHT's 10% accuracy is direct and does not rely on the eikonal correspondence.\n\nSoft spots. Equation (5.1) is derived by dropping the subleading eikonal corrections in (4.31), keeping only b_c omega = kappa. For l=2, kappa=2.5, those corrections are large even for Schwarzschild, and they carry C-dependence. The paper's own Table 1 shows the 6th-order eikonal value at C=0.3, epsilon=0.1, l=2, n=0: Re omega=0.278094 vs the numerical 0.276139, a 0.7% discrepancy; the ratio to Schwarzschild is 0.739 vs the formula's 0.715, about 3% off. The bottom panel of Figure 3 makes the mismatch visible, yet the text asserts the formula is \"sufficiently accurate up to C<=0.3\" without an error budget. The epsilon -> 0 limit is taken analytically, not numerically—the tables only use epsilon=0.1, and the small-epsilon curves in Figure 3 come from the same eikonal expansion, so the astrophysically relevant limit is not independently tested. These are fixable: restrict the claim to eikonal l, provide the 1/kappa-corrected version, or at minimum state the C-dependent error at l=2.\n\nMinor: the abstract says \"fit\" where the derivation is analytic, and \"up to C~O(1)\" overstates the tabulated range (0.3). Self-citation is appropriate; the prior leading-order results are credited.\n\nWho this is for: anyone estimating environmental systematics for ringdown tests of GR, or using EHT shadows to bound halo compactness. It deserves a serious referee; with a tightened validity claim and clearer language it would be a useful reference.","headline":"A solid but modest extension of the known environmental QNM redshift, with a clean shadow bound and an over-stated validity range for the new C^2/6 formula.","tokens_in":13137,"tokens_out":4454,"would_cite":true,"duration_ms":39404,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.70.-s","95.35.+d","04.30.-w"],"model":"deepseek-v4-flash","headline":"A black hole embedded in a Hernquist dark-matter halo has a shadow radius that grows with the halo compactness C, and its ringdown frequencies are redshifted by the factor 3√3 M_BH / b_c = 1 - C + C²/6; current shadow-imaging data bound the","keywords":["Schwarzschild-Hernquist black hole","dark matter halo","Einstein cluster","black hole shadow","quasi-normal modes","eikonal correspondence","compactness bound","Event Horizon Telescope"],"falsifier":"Run an independent time-domain evolution code for the axial l=2 fundamental mode of the Schwarzschild-Hernquist black hole at, say, C=0.5 and ε=0.01, and compare Re[ω]/Re[ω_Schwarzschild] with 1-C+C²/6. If the difference exceeds the numerical error (a few tenths of a percent), the claimed C≤0.3 validity range of the redshift formula fails. Observationally, a future shadow measurement of the Galactic center black hole with a few-percent precision would tighten or contradict the C≤0.092 bound: a shadow diameter above 5.716 M_BH would require either a larger compactness or a different halo model.","tokens_in":12167,"feed_emoji":"🕳️","tokens_out":14498,"duration_ms":123239,"temperature":0.7,"pith_summary":"This paper studies a black hole sitting at the center of a Hernquist dark-matter halo, using an exact solution of Einstein's equations built from an Einstein cluster with zero radial pressure. It claims that two observable signals are controlled mainly by the halo compactness C = M_DM/a_0: the shadow radius grows with C, and the gravitational-wave ringdown (quasi-normal mode) frequencies are redshifted relative to vacuum Schwarzschild. The central result is a simple scaling law—the ringdown frequency ratio equals 3√3 M_BH/b_c, which expands as 1 - C + C²/6 + O(C³) for C ≤ 0.3. Taking the Event Horizon Telescope's shadow measurements at their roughly 10% precision, the paper converts the enlarged shadow into an upper bound C ≤ 0.092. This matters because future ringdown observations of galactic black holes must account for this environmental redshift, and because the bound constrains how concentrated dark matter can be around a supermassive black hole.","feed_headline":"Black-hole shadow data cap halo compactness at 0.092","feed_subtitle":"Ringdowns shift as 1-C+C²/6; EHT shadows put the halo compactness below 0.092.","key_machinery":"The load-bearing objects are the compactness parameter C = M_DM/a_0 and the critical impact parameter b_c of the photon sphere. For null geodesics the effective potential is V_L = f(r)/r², and the photon sphere follows from r - 3m(r) = 0; b_c = r_ph/√f(r_ph) then fixes both the shadow radius and, through the eikonal light-ring/QNM correspondence, the ringdown frequency scale. The argument expands b_c in powers of C, inserts it into the eikonal correspondence b_c ω = κ + c_1/κ + ..., and obtains the redshift formula (5.1); the same expansion is checked against direct numerical QNM calculations at l=2 and l=3.","core_discovery":"The paper's central claim is that for the Schwarzschild-Hernquist black hole, the strong-field observables are governed to leading order by a single dimensionless number, the halo compactness C = M_DM/a_0. In the astrophysically relevant limit of small mass ratio ε = M_BH/M_DM, the critical impact parameter of the photon sphere—which determines the shadow radius seen by a distant observer—grows with C, and the quasi-normal mode frequencies are redshifted by the inverse of that same impact parameter: ω(C, ε)/ω(0,0) = 3√3 M_BH / b_c = 1 - C + C²/6 + O(C³). The paper computes axial gravitational QNMs by three independent methods (matrix, pseudospectral, and sixth-order WKB), confirms the redshi","pith_inferences":["Editorial inference: if the same eikonal redshift structure persists for other Einstein-cluster halo profiles, then the ringdown spectrum would encode an integrated compactness rather than the detailed shape of the density profile, making different halo models hard to tell apart by ringdown alone.","Editorial inference: the paper checks the low-multipole validity of the redshift formula for only a few parameter combinations (ε = 0.1, C = 0.1–0.3); a systematic scan over (C, ε) with an independent time-domain code would show whether the formula can serve as a template for actual ringdown observations, which are dominated by l = 2.","Editorial inference: a next-generation very-long-baseline interferometer with shadow precision of a few percent would tighten the bound on C by an order of magnitude, or, if the measured shadow came out larger than 5.716 M_BH, would point to a more compact halo or a breakdown of the Hernquist Einstein-cluster model.","Editorial inference: the same b_c-based redshift factor should also appear in the damping times of the ringdown, so gravitational-wave detectors with good signal-to-noise ratio on the ringdown of a galactic black hole could constrain halo compactness independently of the shadow measurement."],"forward_implications":["If the bound C ≤ 0.092 holds, current shadow images of M87* and the Galactic center black hole cannot distinguish a Schwarzschild-Hernquist black hole from vacuum Schwarzschild, since the shadow shift stays within the observational error.","A future claim of a shadow excess at the roughly 10% level could be interpreted as a halo-compactness measurement rather than a deviation from general relativity.","The ringdown redshift is controlled by b_c, so measuring a galactic black hole's shadow radius and its ringdown frequency together would test the eikonal relation and the halo model simultaneously.","Dark-matter halos with compactness near the allowed range produce enough frequency redshift to act as a systematic error in ringdown-based tests of general relativity, so those tests must marginalize over environmental compactness.","The extension to C ~ O(1) makes the calculation relevant for dense dark-matter spikes or compact halos, where lower-order redshift formulas fail."],"supporting_citations":[{"why":"Supplies the exact Schwarzschild-Hernquist solution: the Einstein-cluster metric, Hernquist density profile, and mass function used throughout.","marker":"[1]"},{"why":"Provides the axial perturbation master equation and effective potential for black holes surrounded by generic dark-matter profiles, on which the QNM calculation is built.","marker":"[16]"},{"why":"Gives a closely related shadow-bound derivation for black holes with dark-matter halos, the result this paper refines toward C ≤ 0.092.","marker":"[25]"},{"why":"Established the leading-order redshift of QNMs proportional to compactness for C ≤ 10^-2, which this paper extends to second order.","marker":"[29]"},{"why":"Independent earlier study of Schwarzschild-Hernquist QNMs and tidal Love numbers, providing the leading-order redshift context and cross-check.","marker":"[30]"},{"why":"Defines the vacuum Schwarzschild Regge-Wheeler equation, the baseline whose QNMs the redshift formula is measured against.","marker":"[38]"},{"why":"Supplies the numerical linearization of the quadratic eigenvalue problem used in the pseudospectral computation of QNMs.","marker":"[39]"},{"why":"Supplies the matrix method used as an independent numerical route to the QNM frequencies.","marker":"[41]"},{"why":"Supplies the sixth-order WKB formula used for the semi-analytic QNM estimates and eikonal expansion.","marker":"[44]"},{"why":"Establishes the eikonal light-ring/QNM correspondence that turns the photon-sphere impact parameter into the ringdown redshift formula.","marker":"[47]"}],"fun_headline_variants":["Shadow radius sets dark-matter halo bound: C≤0.092","BH ringdowns redshift with halo compactness: 1-C+C²/6","EHT caps halo compactness at 0.092 for galactic BH","Schwarzschild-Hernquist BH: compactness controls shadow and QNMs"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The redshift formula (5.1) assumes the eikonal (large-multipole) relation between ringdown frequency and photon-sphere impact parameter stays accurate at the low multipoles l=2 and l=3 for compactness up to C≈0.3, where the paper checks it only for a small set of parameter values.","fun_headline_variants_meta":{"raw":{"variants":["Shadow radius sets dark-matter halo bound: C≤0.092","BH ringdowns redshift with halo compactness: 1-C+C²/6","EHT caps halo compactness at 0.092 for galactic BH","Schwarzschild-Hernquist BH: compactness controls shadow and QNMs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000693,"raw_usage":{"total_tokens":2965,"prompt_tokens":727,"completion_tokens":2238,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":471,"completion_tokens_details":{"reasoning_tokens":2163}},"tokens_in":471,"tokens_out":2238,"duration_ms":18257,"temperature":1.0,"reasoning_tokens":2163,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T10:28:49.412746+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run an independent time-domain evolution code for the axial l=2 fundamental mode of the Schwarzschild-Hernquist black hole at, say, C=0.5 and ε=0.01, and compare Re[ω]/Re[ω_Schwarzschild] with 1-C+C²/6. If the difference exceeds the numerical error (a few tenths of a percent), the claimed C≤0.3 validity range of the redshift formula fails. Observationally, a future shadow measurement of the Galactic center black hole with a few-percent precision would tighten or contradict the C≤0.092 bound: a shadow diameter above 5.716 M_BH would require either a larger compactness or a different halo model.","supporting_citations":[{"cited_title":"Stability of a Schwarzschild singularity,","cited_arxiv_id":null,"evidence_quote":"Defines the vacuum Schwarzschild Regge-Wheeler equation, the baseline whose QNMs the redshift formula is measured against."},{"cited_title":"A Matrix Method for Quasinormal Modes: Schwarzschild Black Holes in Asymptotically Flat and (Anti-) de Sitter Spacetimes","cited_arxiv_id":"1610.08135","evidence_quote":"Supplies the matrix method used as an independent numerical route to the QNM frequencies."},{"cited_title":"Quasinormal behavior of the d-dimensional Schwarzschild black hole and higher order WKB approach,","cited_arxiv_id":null,"evidence_quote":"Supplies the sixth-order WKB formula used for the semi-analytic QNM estimates and eikonal expansion."}],"review_version":1}