{"id":"aecf78e9-0d7b-404e-871c-e70f625948c1","arxiv_id":"2505.09866","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A generalized WKB/Bohr-Sommerfeld quantization rule is derived for holographic QCD potentials with an infinite boundary barrier and applied to scalar and vector fields in the soft-wall model, matching shooting-method quasinormal frequencies.","lead":"The authors adapt WKB and Bohr-Sommerfeld quantization rules to compute quasinormal mode frequencies in holographic QCD models, and test them against direct numerical shooting. The method reproduces the shooting results to sub-percent accuracy in most regimes and provides analytic low-temperature formulas for the real part of the frequency.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central quantization rule Eq. (21) as printed is inconsistent with the Gamow formula Eq. (23): the tunneling exponent lacks the factor 2 required to reproduce the tabulated decay rates.","rationale":"I read the paper as making a concrete, falsifiable claim: the generalized Bohr-Sommerfeld rules reproduce the shooting-method quasinormal frequencies across three regimes, with no fitted parameters. The numerical agreement in the tables is internally consistent and genuinely supports that claim; the above-barrier deviations below 1% are especially persuasive. My stress-test did not find a demonstrated physical error in the WKB approach itself. What I did find is that the central formula, as printed, is not self-consistent: Eq. (21) and Eq. (23) disagree by a factor of 2 in the barrier exponent, and Appendix D's displayed Langer potential disagrees with the quoted integral and spectrum. Because the manuscript ships no code or data, these typos are not resolvable from the text, and a reader cannot reproduce the tables from the equations. This reinforces the reader's CONDITIONAL verdict rather than overturning it. The empirical evidence is strong enough that the paper should not be rejected, but the derivation must be corrected and ideally accompanied by the numerical implementation before the central claim can be fully verified.","tokens_in":24069,"tokens_out":25643,"duration_ms":267376,"concrete_test":"Independently re-derive the imaginary part from Eq. (21) exactly as printed: expand the left side to first order in epsilon_I and match the imaginary part with -1/4 exp(i integral sqrt(epsilon_R - V_L) dr). Then evaluate both the resulting exp(-gamma) formula and the paper's Eq. (23) exp(-2 gamma) formula at T = 20 MeV, n = 0, and compare both predictions with Table I. If the printed Eq. (21) is used, the imaginary part will be tens of orders of magnitude too large, confirming the missing factor 2. Separately, recompute the Appendix D integral with the displayed V_L = 2 + 1/(4 z^2) + z^2; if it does not yield -3 pi/2 + m^2 pi/4, then the displayed Langer potential is also incorrect. One of these checks settles whether the formulas as written support the central claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing issue is an internal inconsistency in the central derivation, not the physical plausibility of the Langer matching. Equation (21) states that the Bohr-Sommerfeld integral equals pi(n+1/2) minus (i/4) exp(i integral p_L dr), with no factor 2 in the exponent. Expanding this equation for small imaginary part gives a decay amplitude proportional to exp(-gamma), where gamma is the single-pass barrier action. However, the paper's Gamow formula, Eq. (23), and the numerical values in Tables I-IV require exp(-2 gamma), the standard two-pass tunneling exponent. At T = 20 MeV, the tabulated imaginary parts are around 10^-88; with only exp(-gamma) they would be roughly 10^-44, tens of orders of magnitude larger. Thus the printed generalized Bohr-Sommerfeld rule cannot reproduce the paper's own tables unless a missing factor 2 is silently inserted. Since no code or data accompany the manuscript, a reader cannot determine whether the tables were generated from Eq. (21), from Eq. (23), or from some unstated corrected formula. A second, related typo appears in Appendix D: the displayed Langer potential V_L = 2 + 1/(4 z^2) + z^2 gives the WKB action pi/4 (m^2 - 3), whereas the quoted integral -3 pi/2 + m^2 pi/4 and the claimed spectrum m^2 = 4 c^2 (n+2) correspond to V_L = 2 + 4/z^2 + z^2. These inconsistencies are exactly the places where a reader trying to apply the claimed quantization rules would fail.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives generalized Bohr–Sommerfeld (WKB) quantization rules for quasi-stationary states in holographic QCD, treating the effective Schrödinger potential with an infinite boundary barrier through a Langer-type correction. The rules are applied to scalar and vector fluctuations in the soft-wall AdS black hole model and to a charmonium-like tangent model, and the resulting quasinormal frequencies are compared with a shooting method in three regimes: below the barrier, near the top, and above the barrier. The authors report sub-percent agreement in the real parts and a few-percent agreement in the imaginary parts in most regimes, and they provide analytic low-temperature expansions for the real parts as well as a discussion of dissociation signaled by the above-barrier transition.","tokens_in":24431,"tokens_out":10419,"duration_ms":94165,"significance":"If the claims survive revision, the paper provides a fast semiclassical route to quasinormal-mode spectra in holographic QCD, with analytic control at low temperature and a physical criterion for dissociation. The numerical comparisons are genuine benchmarks: no quantity is fitted to the quasinormal frequencies, and the first-order corrected rules agree with the shooting method to better than about 1% in the real part in nearly all tables (e.g., Tables III, IV, and VII), with imaginary parts agreeing to a few percent in the under-barrier and above-barrier regimes. The extension to vector fields and to a charmonium-like potential demonstrates that the method is not limited to the scalar soft-wall case. However, the printed central quantization rule and several appendix formulas contain internal inconsistencies that currently prevent an independent reader from reproducing the claimed results.","major_comments":[{"comment":"The generalized Bohr–Sommerfeld rule is printed as an equality between the well integral and π(n+1/2) − (i/4) exp(i ∫_{r1}^{r2} p_L dr'), with no factor 2 in the barrier phase. Expanding this equation for small ε_I gives an imaginary part proportional to e^{−γ}, where γ = ∫_{r1}^{r2} √(V_L−ε_R) dr' is the single-pass barrier action. The Gamow formula (23), however, contains e^{2i∫...}, i.e. e^{−2γ}, and the tabulated values require the two-pass exponent: at T = 20 MeV, γ ≈ 100, so e^{−γ} would give |ω_I| ∼ 10^{−44} rather than the printed 1.48 × 10^{−88}. Since Eq. (23) is stated to follow from Eq. (21), the central derivation is internally inconsistent. The manuscript provides no code or data to indicate whether Tables I–IV were generated from Eq. (21), from Eq. (23), or from an unstated corrected rule; this inconsistency must be resolved before the quantization rule can be used as stated.","section":"Eqs. (21) and (23)"},{"comment":"Equation (D3) defines VL(z) = 2c^2 + 1/(4z^2) + c^4 z^2, but adding the Langer term 1/(4r_*^2) to the zero-temperature potential (D1), V(z) = 2c^2 + 15/(4z^2) + c^4 z^2, gives VL(z) = 2c^2 + 4/z^2 + c^4 z^2. The quoted integral −3π/2 + m^2π/4 and the claimed spectrum m^2 = 4c^2(n+2) correspond to the corrected potential, not to the displayed one. With the displayed VL and c = 1, the WKB action is (π/4)(m^2 − 3), which yields m^2 = 4n + 5 rather than the soft-wall spectrum. This error occurs in the place where the paper demonstrates the necessity of the Langer correction and must be corrected.","section":"Appendix D, scalar zero-temperature potential"},{"comment":"The turning point z1 in Eq. (E6) contains √((ε_R − 6c^2)(ε_R − 2c^2)) in the numerator; consistency with the defining equation ε_R − VL = 0 and with the z0 expression in Eq. (E5) requires √((ε_R − 6c^2)(ε_R + 2c^2)). As printed, the two turning points are roots of different quadratic equations, so the subsequent integrals and the analytic formula (E8) cannot be reproduced from the displayed expressions. The corresponding vector-field turning points in Eqs. (E11) and (E12) should be checked for the same type of error.","section":"Appendix E, turning-point expressions"}],"minor_comments":[{"comment":"For T = 55, 56, and 57 MeV, the WKB imaginary parts are printed as 7.96662×10^−3, 1.25368×10^−2, and 1.85427×10^−2, respectively, while the shooting-method values are 8.47488×10^−4, 1.32819×10^−3, and 1.95701×10^−3. The reported deviations of 6.00%, 5.61%, and 5.25% are consistent only if the WKB entries are 7.96662×10^−4, 1.25368×10^−3, and 1.85427×10^−3; the exponents in the table appear to be off by one order of magnitude.","section":"Table VI, n = 1 block"},{"comment":"The text refers to Table X as the 'second mode' and then as the 'ground state' for n = 1, and Table XI as the 'third mode' and then as the 'ground state' for n = 2; these labels should be made consistent.","section":"Sec. IV.A, Tables X and XI"},{"comment":"The exponential factors in Eq. (19) contain ℏ inconsistently: some terms have i/ℏ, while others have i without ℏ. Since the paper sets ℏ = 1 throughout, either remove all ℏ factors or include them uniformly.","section":"Eq. (19)"},{"comment":"The caption assigns the same color (purple) to the n = 6 and n = 7 modes, which makes the figure harder to read; a distinct color would improve clarity.","section":"Fig. 9"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the paper addresses a useful problem, and the benchmark comparisons with the shooting method are valuable. However, the number of typographical and internal inconsistencies in central formulas is unusually high for a journal submission. I recommend inviting a revision that supplies the corrected quantization rule, the corrected appendix expressions, and ideally the code or data used to generate the tables. The absence of code and data is not by itself a reason to reject, but it makes the current inconsistencies unresolvable by the reader."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, honestly, the paper is better than its central typo. The core result—that a generalized Bohr-Sommerfeld rule with a Langer correction reproduces shooting-method QNM frequencies in soft-wall holographic QCD—is demonstrated with real benchmarks. No fitted parameters, three regimes, and the first-order corrections push deviations below 1% in the real part and often below 0.1%. The analytic low-temperature expansions for the real part (E8, E14) are a nice addition. The dissociation-temperature picture is a bonus, not the main point.\n\nThe new content is the adaptation of the atomic-physics BS machinery to potentials with an infinite boundary barrier, which prior black-hole applications did not have. That adaptation is nontrivial because the WKB wave function near the boundary needs the Langer term 1/(4r_*^2), and the paper checks this by matching to the normalizable z^{5/2} falloff. The checks are plausible, and the numerical agreement supports them.\n\nNow the soft spots. Eq (21), the main quantization rule, is missing a factor 2 in the tunneling exponent. As printed, a direct expansion gives an imaginary part proportional to e^{-G}, while the paper's own Gamow formula Eq (23) and all the tables require e^{-2G}. For the T=20 MeV ground state that is the difference between 10^-44 and 10^-88. So the printed rule cannot reproduce the tables. This looks like a typo—Eq (23) is the standard formula and the tables follow from it—but it has to be fixed. Relatedly, Appendix D writes the zero-temperature Langer potential as V_L = 2 + 1/(4z^2) + z^2; the quoted integral and spectrum correspond to V_L = 2 + 4/z^2 + z^2. Table VI also has a factor-10 error in several WKB imaginary parts (the deviations line up with 10^-4, not 10^-3). None of these undermine the method, but each would trip up a reader trying to apply it.\n\nThe Appendix B matching is the part I'd push on: it is done via low-temperature expansion and asymptotic analysis, and the claim that the non-normalizable term is absent for complex omega at finite T is stated more than proven. Given the numerical agreement, I'm willing to believe it, but a referee should ask for a more careful argument or a numerical check at moderate T.\n\nWho is this for? Anyone computing QNM widths in soft-wall or similar holographic models who wants a fast semi-analytic alternative to shooting. It deserves a serious referee; the typos are easy to fix and the core is solid. I'd send it to review, with a request for code/data, not desk-reject.","headline":"Solid WKB/BS benchmark for holographic QNM widths, but Eq (21) has a missing factor 2 and needs a typo-fixing revision before the rules can be used as printed.","tokens_in":24921,"tokens_out":20602,"would_cite":true,"duration_ms":174172,"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":"Generalized Bohr-Sommerfeld rules reproduce quasinormal-mode frequencies in holographic QCD models.","keywords":["quasinormal modes","holographic QCD","soft-wall model","WKB approximation","Bohr-Sommerfeld quantization","glueballs","dissociation","finite temperature"],"falsifier":"A direct check is to compute the ground-state quasinormal frequency at low temperature with a high-precision method that does not rely on WKB matching near the boundary and compare the imaginary part: the zeroth-order Gamow formula deviates from the shooting method by roughly ten percent in $\\omega_I$ while the real part agrees to $10^{-3}$–$10^{-2}$ percent, so a more accurate independent calculation that confirms the shooting values would expose a systematic error in the Gamow exponential. A second check is the analytic low-temperature formula: at $T = 35$ MeV it gives $2.82466$ GeV for the scalar ground state while the shooting and numerical Bohr-Sommerfeld values are $2.82194$ GeV; if this gap does not close with higher-order terms, the claimed range of the analytic expansion fails.","tokens_in":23890,"feed_emoji":"⚛️","tokens_out":7669,"duration_ms":79333,"temperature":0.7,"pith_summary":"This paper seeks to establish that quasinormal modes in holographic QCD—the complex frequencies that describe how thermal fluctuations decay in the dual gauge theory—can be computed from a generalized Bohr-Sommerfeld quantization rule, the same WKB phase-integral condition used for quasi-stationary states in quantum mechanics. The payoff would be a fast semiclassical route to the full complex spectrum: in the soft-wall model, the real part of scalar and vector quasinormal frequencies is reproduced to sub-percent accuracy after first-order corrections, and at low temperature an analytic formula gives that real part directly. The same framework also locates the temperature at which each mode moves above the potential barrier, which the authors read as a dissociation signal consistent with the spectral function's decreasing peak.","feed_headline":"WKB rules match holographic QCD oscillation frequencies","feed_subtitle":"Semiclassical quantization reproduces shooting-method frequencies, with fast analytic thermal mass predictions.","key_machinery":"The machinery is the generalized Bohr-Sommerfeld quantization condition for quasi-stationary states, applied to the Schrödinger-like equation obtained from the holographic field equation. The central object is the Langer-transformed potential $V_L(r_\\ast) = V(r_\\ast) + \\frac{1}{4 r_\\ast^2}$, whose extra term corrects the infinite barrier at the boundary and makes the WKB wave function vanish at $r_\\ast = 0$. For modes near the barrier top, a parabolic approximation produces a gamma-function correction $\\Xi(\\lambda)$; for above-barrier modes, the condition is continued by $\\lambda \\to \\lambda e^{-2\\pi i}$, giving the contour form used in the tables. The imaginary part of the frequency comes from the Gamow formula, an exponential of the barrier integral, and first-order WKB corrections are included through an $I_1$ term. Together these pieces turn the quasinormal-mode boundary-value problem into phase integrals over turning points.","core_discovery":"The paper's central claim is that quasinormal modes in holographic QCD models are reproduced by generalized Bohr-Sommerfeld quantization rules for quasi-stationary states, with an infinite-barrier boundary handled by a Langer-type correction to the effective potential. Three temperature regimes are treated: modes trapped near the potential minimum, modes near the top of the finite potential barrier, and above-barrier modes obtained by analytic continuation. In each regime, the quantization condition yields frequencies that agree with the shooting method; with the first-order WKB correction, deviations in the real part remain below about one percent in most cases. In the soft-wall model at very low temperature, the real part of the scalar and vector quasinormal frequencies is given analytically as a temperature expansion. The shift of modes into the above-barrier region is interpreted as evidence of dissociation of the dual states, in line with the decreasing peak of the spectral function.","pith_inferences":["A natural step the authors do not take is to use the quantization rules as a fitting engine, extracting the soft-wall model parameters such as the dilaton scale directly from measured or lattice thermal masses and widths without repeated full numerical integration.","The predicted low-temperature shift of the real part scales as $T^4$ then $T^8$, which is specific enough that an independent calculation of the thermal pole of the retarded correlator, or a lattice computation of thermal glueball masses, would test the soft-wall model's temperature dependence.","The dissociation temperatures read off from the above-barrier transition could be compared with the inflection point of the spectral-function peak; if the peak persists far beyond those temperatures, the dissociation criterion would need refinement.","The imaginary part agrees only to about ten percent at very low temperature while the real part agrees to parts in $10^{-4}$, suggesting the Gamow exponential factor is the most sensitive ingredient and is the most promising target for a next-order correction.",""],"forward_implications":["The WKB formulas provide a fast semiclassical route to quasinormal frequencies in any holographic model whose effective potential has a well plus a finite barrier, with deviations below one percent for the real part in most tested cases.","At very low temperatures, the analytic expressions for the real part of scalar and vector frequencies allow direct reading of thermal mass shifts without solving the full differential problem.","The temperature at which a mode enters the above-barrier region can be read as a dissociation temperature; the paper lists values such as roughly 77 MeV for the scalar ground state and 105 MeV for the vector ground state.","In the high-temperature above-barrier regime, quasinormal frequencies approach the linear-in-$T$ scaling of conformal AdS$_5$, connecting the low-temperature non-conformal regime to the conformal limit.","Because the quantization rules work for scalar and vector fields in the soft-wall model and for a tangent-model charmonium potential, they should extend to other dilaton profiles and flavors.",""],"supporting_citations":[{"why":"Supplies the generalized Bohr-Sommerfeld quantization rules for quasi-stationary states that the paper adapts to holographic potentials.","marker":"[33–37]"},{"why":"Provides the WKB matching and the parabolic approximation near the potential-barrier top used in the derivation.","marker":"[38]"},{"why":"Gives the Langer transformation whose $1/(4r_\\ast^2)$ correction fixes the WKB wave function near the boundary.","marker":"[40]"},{"why":"Defines the shooting method used as the comparison baseline for the quasinormal frequencies.","marker":"[42]"},{"why":"Supplies the first-order WKB correction term incorporated into the improved quantization conditions.","marker":"[43]"},{"why":"Introduces the soft-wall holographic QCD model whose scalar and vector spectra are computed.","marker":"[8]"},{"why":"Provides the holographic prescription relating quasinormal modes to retarded correlators and the spectral function.","marker":"[5]"},{"why":"Provides the vector-field holographic potential used for the vector meson calculation in the soft-wall model.","marker":"[47]"}],"fun_headline_variants":["Generalized Bohr-Sommerfeld rules reproduce holographic QCD modes","WKB quantization yields analytic holographic QCD frequencies at low T","Holographic QCD quasinormal modes from semiclassical quantization","Bohr-Sommerfeld condition matches shooting for holographic modes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the boundary treatment—the Langer-corrected potential plus the WKB matching that keeps only the normalizable $z^{5/2}$ solution near the boundary—stays valid for complex frequencies at finite temperature; if that matching is wrong, the quantization conditions would select the wrong modes even though the current shooting comparisons agree.","fun_headline_variants_meta":{"raw":{"variants":["Generalized Bohr-Sommerfeld rules reproduce holographic QCD modes","WKB quantization yields analytic holographic QCD frequencies at low T","Holographic QCD quasinormal modes from semiclassical quantization","Bohr-Sommerfeld condition matches shooting for holographic modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000485,"raw_usage":{"total_tokens":2319,"prompt_tokens":799,"completion_tokens":1520,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":415,"completion_tokens_details":{"reasoning_tokens":1442}},"tokens_in":415,"tokens_out":1520,"duration_ms":12580,"temperature":1.0,"reasoning_tokens":1442,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:22:39.594349+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct check is to compute the ground-state quasinormal frequency at low temperature with a high-precision method that does not rely on WKB matching near the boundary and compare the imaginary part: the zeroth-order Gamow formula deviates from the shooting method by roughly ten percent in $\\omega_I$ while the real part agrees to $10^{-3}$–$10^{-2}$ percent, so a more accurate independent calculation that confirms the shooting values would expose a systematic error in the Gamow exponential. A second check is the analytic low-temperature formula: at $T = 35$ MeV it gives $2.82466$ GeV for the scalar ground state while the shooting and numerical Bohr-Sommerfeld values are $2.82194$ GeV; if this gap does not close with higher-order terms, the claimed range of the analytic expansion fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the WKB matching and the parabolic approximation near the potential-barrier top used in the derivation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Langer transformation whose $1/(4r_\\ast^2)$ correction fixes the WKB wave function near the boundary."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the shooting method used as the comparison baseline for the quasinormal frequencies."},{"cited_title":"Scrucca, Quantum physics iii (2025), master program in Physics","cited_arxiv_id":null,"evidence_quote":"Supplies the first-order WKB correction term incorporated into the improved quantization conditions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the vector-field holographic potential used for the vector meson calculation in the soft-wall model."}],"review_version":1}