{"id":"b0ae8810-cd9f-45e8-be77-b6ed1e183784","arxiv_id":"2607.19755","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Near-horizon analysis yields analytic critical couplings for scalar-tensor instabilities of regular black holes, but the accompanying area-quantization claim is undermined by a factor-of-2 error.","lead":"This paper derives closed-form formulas for the scalar-field coupling at which regular black holes become unstable, using a near-horizon expansion of the perturbation potential. It also claims to recover black-hole area quantization at that critical coupling, but that final step contains an algebraic error.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's central threshold criterion—identifying instability onset with V'(r_h)=0—is supported only by heuristic numerical examples; no proof or precise numerical test establishes that a zero-energy bound state appears exactly at ζ_c, so Eqs. (3.8)–(3.9) may not give the true instability onset.","rationale":"The reader's weakest_assumption correctly identifies the b=0 criterion as the load-bearing bridge between the near-horizon algebra and the instability claim. My stress-test reaches the same point: the paper asserts, but does not prove, that the first unstable mode appears exactly when the near-horizon potential has an extremum on the horizon. The supporting evidence is a single illustrative numerical example (Fig. 1) and a qualitative statement about other cases. Since the central critical-coupling formulas (3.8) and (3.9) rest entirely on this identification, the claim is not yet established. The factor-of-two error in Eq. (3.14) noted by the reader is real and affects the area-quantization section, but it is peripheral to the instability analysis; the threshold criterion is more fundamental. A targeted numerical zero-mode search would settle whether the criterion is correct. Until then, the appropriate verdict remains CONDITIONAL: the paper should either prove the b=0 threshold or validate it with a precise spectral calculation before the critical-coupling formulas are accepted.","tokens_in":13481,"tokens_out":9139,"duration_ms":95945,"concrete_test":"For one background (e.g., NC Schwarzschild with θ=0.2, M=1, ℓ=2, μ=0.5), compute the Regge-Wheeler spectrum for ω²<0 numerically (shooting method on the full radial equation, not the near-horizon expansion), scanning ζ around the value ζ_c=5.673 quoted in Fig. 1. Find the smallest ζ for which a normalizable bound-state solution exists. Also solve the zero-mode problem (ω=0) directly. If the onset ζ differs from ζ_c by more than the numerical tolerance, the b=0 criterion and hence Eqs. (3.8)–(3.9) are not the instability threshold.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. 3, Eq. (3.3) approximates the potential as V(x)≈a x²+b x, and the critical coupling is obtained from b=0. This assumes the onset of tachyonic instability coincides with an extremum of V located exactly on the horizon. The paper provides no analytical argument that the first negative-energy (bound) mode of the Regge-Wheeler operator appears at that parameter value; a negative well outside the horizon can exist for a range of couplings before it is deep enough to support a bound state (or, conversely, a very shallow well may still support a bound state in 1D). The only evidence is Fig. 1 for a single configuration and the statement that other plots behave similarly. For the Einstein-coupled model the near-horizon form (3.3) is an approximation, so Eq. (3.9) inherits the same unverified identification.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyses linear scalar perturbations of spherically symmetric regular black holes (non-commutative Schwarzschild and Bardeen/Hayward/ABG) in two non-minimally coupled models: Ricci coupling and Einstein-tensor coupling. Expanding the metric near the horizon as f ≈ 2κx + ..., the effective potential is approximated by V(x) ≈ a x² + b x. The instability threshold is identified with V'(r_h) = 0, which yields closed-form critical couplings, Eqs. (3.8) and (3.9). The authors further derive near-horizon QNM frequencies, argue that at the critical coupling the modes become purely imaginary, and use the spacing of these modes to reproduce an equidistant area spectrum A = 8πn without the highly-damped approximation. Numerical verification is based on time-domain profiles reported in the authors' earlier work.","tokens_in":13698,"tokens_out":10543,"duration_ms":100264,"significance":"If correct, Eqs. (3.8)-(3.9) would be remarkably simple horizon-data formulas for the onset of instability, with no fitted parameters. The paper performs several sensible consistency checks (Schwarzschild limit, eikonal limits, mass dependence) and the qualitative trends in ℓ, θ, q match the authors' earlier numerics. The derivation of the threshold, however, rests on an unproved equivalence between V'(r_h) = 0 and the appearance of an unstable mode, and the QNM section contains a factor-of-two algebraic error. The manuscript is therefore a promising but incomplete treatment; its main claims need additional spectral justification.","major_comments":[{"comment":"The identification of the instability threshold with b = 0 is the load-bearing step, but it is not proved. b = 0 is only the condition that the near-horizon effective potential has an extremum at the horizon. For a Schrödinger operator with potential a x² + b x (or any generic potential with a negative well), the onset of a zero-energy bound state is controlled by the lowest eigenvalue crossing zero, which generally occurs at a finite negative b, not at b = 0. The evidence offered, Fig. 1 for a single NC-Schwarzschild configuration and a sentence that other plots behave similarly, is heuristic. Please prove the equivalence or provide a controlled numerical test (e.g., compute the fundamental QNM frequency as a function of ζ and show that Im ω → 0 exactly at b = 0), or Eqs. (3.8)-(3.9) remain unjustified.","section":"Sec. 3.1, Eq. (3.3)"},{"comment":"The Gamma-pole algebra is incorrect. Solving the displayed condition gives ω = -b/(2√-a) - iκ(2n+1), not ω = -b/(2√-a) - iκ(n+1/2). Consequently the spacing of near-horizon modes is Δω = 2κ rather than κ, so the area quantization in Sec. 3.3 becomes A = 4πn, not A = 8πn. The qualitative result Re ω = 0 at b = 0 survives, but Eq. (3.17), the area-spectrum claim, and the comparison with the highly-damped literature must be revised.","section":"Sec. 3.2, Eq. (3.14)"},{"comment":"The abstract describes Eqs. (3.8) and (3.9) as 'exact expressions', but the text states that the near-horizon recipe gives the exact result only for the scalar model and that the tensor case requires the near-horizon expansion. Eq. (3.9) is therefore approximate, and no error estimate is given. Please either limit the exactness claim to the scalar-model formula or quantify the accuracy of Eq. (3.9) by comparing with higher-order terms or with direct numerical solution of the full potential.","section":"Abstract, Sec. 3.1, Eqs. (3.8)-(3.9)"}],"minor_comments":[{"comment":"The term 'the so-called RGB' should be 'RBH'.","section":"Sec. 1"},{"comment":"The d-dimensional redefinition A = ψ/r^{(d-2)/2} is stated only for the Ricci-coupled model; clarify whether the tensor-model reduction in Sec. 2.2 has an analogous form.","section":"Sec. 2.1, footnote 1"},{"comment":"The generalized Laguerre function notation is incomplete (missing script L) and parentheses are unbalanced; please rewrite these equations for readability.","section":"Eqs. (3.12)-(3.13)"},{"comment":"'Comparying' should be 'Comparing'.","section":"Fig. 4 caption"},{"comment":"'General area quantization' overstates the result; the derivation applies to the near-horizon sector of spherically symmetric backgrounds and should be phrased accordingly.","section":"Sec. 3.3"}],"recommendation":"major_revision","confidential_remarks":"The central risk is the threshold identification. If the authors can supply a rigorous proof or a high-precision numerical test showing that the onset of instability coincides with b = 0, the paper could be suitable for publication. The QNM algebra error is objective and must be fixed, as it changes the area-spectrum spacing by a factor of two."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dave,\n\nThe useful thing here is the pair of closed-form critical coupling expressions, Eqs. (3.8) and (3.9), for the onset of scalar-field instability on a class of regular black holes. If you take the near-horizon quadratic expansion at face value, the Ricci-coupled formula follows cleanly from V'(r_h)=0, and the tensor formula is at least plausible. That turns a numerical threshold search into a one-line algebraic check. The paper checks the Schwarzschild limit, the eikonal limit, and the dependence on scalar mass; those checks are consistent, and the time-domain profiles in Fig. 1 support the claimed critical values for the examples shown.\n\nThe soft spots are real. First, the central identification - instability starts exactly when the near-horizon potential has an extremum on the horizon, i.e., b=0 - is asserted, not proved. In the pure quadratic model V(x)=a x^2+b x on the half-line x>0, a zero-energy bound state requires b = -2κ sqrt(a) for the lowest mode, not b=0. So the paper's criterion is not self-evident, and one plotted example is weak evidence. The formulas may be good approximations, but calling them exact overstates what is shown. Second, the tensor formula is not derived in the text; the reader is asked to trust the algebra. Third, the QNM section has a clear factor-of-two error: applying the paper's own Gamma-pole condition gives ω = -b/(2 sqrt(-a)) - iκ(2n+1), not -iκ(n+1/2). That changes the frequency spacing from κ to 2κ and undermines the area-quantization claim advertised in the abstract. That section should be corrected or retracted.\n\nThe core instability analysis is a useful heuristic, but it is not a rigorous threshold theorem. The paper would be credible if it presented the criterion as an empirical or approximate condition and supplied a derivation of the tensor formula plus a systematic numerical test across the parameter range. As it stands, the gap between the claim and the evidence is too large for unconditional acceptance.\n\nSend it to a referee who knows black hole perturbation theory; it is not a desk reject. The referee should ask for the tensor derivation, a correction of the QNM factor, and either a proof or an honest downgrading of the threshold criterion. If those are supplied, the formulas are worth publishing.\n\nBest,","headline":"Useful analytic formulas for instability thresholds in regular black holes, but the central threshold criterion is asserted rather than proven, and the QNM/area-quantization section contains a factor-of-two error.","tokens_in":14223,"tokens_out":16396,"would_cite":false,"duration_ms":164362,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper derives closed-form critical coupling constants beyond which regular black holes become unstable under non-minimal scalar-field perturbations, identifying the threshold with an extremum of the near-horizon effective potential loca","keywords":["regular black holes","non-minimal scalar-tensor coupling","tachyonic instability","Regge-Wheeler equation","near-horizon approximation","quasinormal modes","area quantization","critical coupling constant"],"falsifier":"For the NC Schwarzschild tensor model with θ = 0.2, ℓ = 2, μ = 0.5, the paper predicts ζ_c = 5.673. Numerically integrate the exact Regge-Wheeler equation for ζ = 5.6 (just below) and search for a negative well outside r_h or late-time growth in the time-domain profile; finding either would contradict the claim that the threshold sits exactly at b = 0. Alternatively, check directly whether the exact first derivative V'(r_h) vanishes at the predicted ζ_c for the tensor model using the full potential, not the near-horizon approximation.","tokens_in":13318,"feed_emoji":"🕳️","tokens_out":7754,"duration_ms":75444,"temperature":0.7,"pith_summary":"This paper tries to establish that regularity does not protect black holes from instability once a scalar field is allowed to couple to the spacetime curvature. The authors study four families of regular black holes — non-commutative Schwarzschild, Bardeen, Hayward, and ABG — under two non-minimal couplings: the field coupled directly to the Ricci scalar, and its derivatives coupled to the Einstein tensor. They derive closed-form expressions for the critical coupling at which the black hole becomes unstable, essentially by demanding that the effective potential of the Regge-Wheeler equation have an extremum exactly at the horizon. If the derivation holds, the threshold is computable from horizon data alone: horizon radius, surface gravity, and the second derivative of the metric function. At the threshold, the quasinormal ringing gives way to a purely imaginary mode, and the area spectrum A = 8πn follows without the usual highly-damped approximation.","feed_headline":"A horizon condition predicts when regular black holes destabilize","feed_subtitle":"Near-horizon potentials give exact critical couplings in two scalar-tensor models; at the threshold, ringing becomes purely imaginary.","key_machinery":"The machinery is the near-horizon expansion of the effective potential in the Regge-Wheeler equation for a static spherically symmetric metric with f(r_h) = 0. Writing f(r) ≈ 2κ x + (1/2) f''(r_h) x² with x = r − r_h, and expanding the potential as V(x) ≈ a x² + b x, the paper identifies the instability threshold with b = 0, i.e. V'(r_h) = 0. Solving b = 0 for the coupling ζ gives the critical values. In the Ricci-coupled model this is a direct equation; in the Einstein-tensor-coupled model the coefficients a and b themselves come from the near-horizon expansion of the more complicated kinetic coupling, which is an additional approximation the paper notes.","core_discovery":"At the critical value of the non-minimal coupling, the near-horizon effective potential of the Regge-Wheeler equation, V(x) ≈ a x² + b x, has its extremum exactly at the event horizon; the condition is b = 0. For the Ricci-coupled scalar model this condition is exact and yields ζ_c(Ricci) = [r_h² μ² + 2 r_h κ + ℓ(ℓ+1)] / [r_h² f''(r_h) + 8 r_h κ − 2]. For the Einstein-tensor-coupled model, after a near-horizon expansion, it yields ζ_c(Einstein) = 2 r_h [r_h(r_h μ² + 2κ) + ℓ(ℓ+1)] / [(2 r_h κ + ℓ(ℓ+1))(r_h f''(r_h) + 4κ)]. At this threshold the perturbations neither ring nor decay — the late-time tail is a straight line — and the real part of the near-horizon quasi-normal frequency vanishes,","pith_inferences":["Inference: Because the two closed-form formulas depend only on r_h, κ, f''(r_h), μ, and ℓ, they should be directly testable against full numerical integration of the Regge-Wheeler equation for any metric in the same class, including the singular Reissner-Nordström geometry; the paper does not carry out that check.","Inference: The paper's b = 0 criterion, if true, implies a geometric picture of the onset of instability: the negative well is born exactly on the horizon and then migrates outward; this could be checked by tracking the location of the minimum of V(r) as ζ is swept through ζ_c for all four families.","Inference: The area-quantization result at ζ_c suggests that critical coupling could serve as a proxy for highly damped modes in other contexts, such as computing grey-body factors or entropy spectra, though the authors do not explore those.","Inference: A gap in the argument is that absence of instability for ζ < ζ_c is not proven analytically; a full proof would need to show the exact potential has no negative well outside the horizon whenever ζ is below the b = 0 value."],"forward_implications":["For any spherically symmetric black hole whose metric function has a Taylor expansion near the horizon, Eqs. (3.8) and (3.9) give a direct algebraic prediction of the coupling at which scalar perturbations become unstable, without solving the full perturbation equations.","In the Schwarzschild limit both formulas diverge, recovering the known linear stability of Schwarzschild against these scalar perturbations for any finite coupling.","In the tensor model, the critical coupling becomes independent of the multipole ℓ in both the massless (μ = 0) and eikonal (ℓ → ∞) limits, while in the Ricci model ℓ → ∞ removes the instability entirely; these are testable predictions.","At ζ = ζ_c the near-horizon QNM frequencies become purely imaginary, so the onset of instability is accompanied by a mode that neither oscillates nor decays.","The spacing of near-horizon QNM frequencies at criticality is Δω = κ = 2πT_H, reproducing the area spectrum A = 8πn without invoking highly damped modes."],"fun_headline_variants":["Critical coupling makes black hole horizon a potential extremum","When potential peaks at horizon, regular black holes destabilize","Exact critical coupling for black hole instability found analytically","Near-horizon analysis reveals threshold for black hole instability","Purely imaginary ringing marks black hole instability threshold"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise, asserted and supported by one plotted numerical example rather than proved, is that the onset of instability is exactly the condition b = 0, i.e. an extremum of the near-horizon effective potential located on the horizon; if a negative potential well could already exist for smaller coupling without touching the horizon, the derived critical values would be wrong, and the tensor-model formula further inherits the near-horizon approximation the paper i","fun_headline_variants_meta":{"raw":{"variants":["Critical coupling makes black hole horizon a potential extremum","When potential peaks at horizon, regular black holes destabilize","Exact critical coupling for black hole instability found analytically","Near-horizon analysis reveals threshold for black hole instability","Purely imaginary ringing marks black hole instability threshold"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000454,"raw_usage":{"total_tokens":2141,"prompt_tokens":790,"completion_tokens":1351,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":534,"completion_tokens_details":{"reasoning_tokens":1274}},"tokens_in":534,"tokens_out":1351,"duration_ms":11255,"temperature":1.0,"reasoning_tokens":1274,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T11:48:46.367324+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the NC Schwarzschild tensor model with θ = 0.2, ℓ = 2, μ = 0.5, the paper predicts ζ_c = 5.673. Numerically integrate the exact Regge-Wheeler equation for ζ = 5.6 (just below) and search for a negative well outside r_h or late-time growth in the time-domain profile; finding either would contradict the claim that the threshold sits exactly at b = 0. Alternatively, check directly whether the exact first derivative V'(r_h) vanishes at the predicted ζ_c for the tensor model using the full potential, not the near-horizon approximation.","supporting_citations":[],"review_version":1}