{"id":"71f8a68c-1c47-42aa-815c-d53ea380b1aa","arxiv_id":"2507.18399","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A numerically efficient algorithm and a dimensional-regularization renormalization scheme make the modified Lüscher equation practical for lattice systems with long-range forces.","lead":"This paper works out a practical way to compute a special mathematical function used to extract particle scattering information from lattice simulations when a long-range force is present. In a test model, only the simplest partial wave is needed, and the method avoids the large renormalization problems of earlier approaches.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim rests on perturbativity of the long-range force; the companion assumption that natural-size ERE parameters hold in all partial waves is only checked for repulsive, perturbative couplings, and the ℓ=4 subtraction scheme is not validated for attractive OPE-like channels.","rationale":"The reader identified the same load-bearing assumption: the long-range force is perturbative and does not create bound states/low-lying resonances alone, and the paper itself flags this as needing additional scrutiny. My stress-test confirms this is not a manufactured concern: Section 4.2 states 'only repulsive interactions are considered', and the whole finite-piece construction in Eq. (4.1) assumes the Born series (or the resolvent) is sufficiently convergent that a small number of divergent terms can be separated and the remainder satisfies a well-defined LS equation. For attractive couplings, the perturbative convergence seen in Figs. 7 and 8 is not guaranteed, and physical pion exchange includes attractive partial waves. I do not see an internal inconsistency or a mathematical error in the algebra for repulsive couplings, and the internal checks (µ-independence, Coulomb limit, Hamiltonian spectrum) are genuine supporting evidence. The concern is about scope: the central claim about natural-size ERE parameters in all partial waves and the practical claim about truncation are demonstrated only in the repulsive perturbative regime. That is a qualification, not a refutation, so the CONDITIONAL verdict stands. A cleaner additional point would have been the absence of estimated numerical errors on the VEGAS integrations and the lack of released code, but those are secondary to the perturbativity issue and the reader already noted them.","tokens_in":22381,"tokens_out":2594,"duration_ms":22805,"concrete_test":"On the toy potential of Eq. (3.10), set g negative (attractive Yukawa long-range force), choosing |g| such that the long-range potential alone supports a near-threshold bound or virtual state (e.g., scan g from +0.073m down to the value where the infinite-volume LS equation for V_L alone has a pole). For each g in this scan, compute M_0(q0) and M_4(q0) by (i) the paper's decomposition Eq. (4.1) with the same number of subtracted Born terms and same subtraction scheme, and (ii) a direct solution of the full LS equation with Gaussian quadrature plus a separate subtracted dispersion/spectral representation. If the difference grows dramatically as the bound-state pole approaches threshold, or if the Born series fails to converge, then the decomposition in Eq. (4.1) is not valid for attractive long-range forces and the central claim must be qualified accordingly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claim is that the modified Lüscher approach allows truncating to the S-wave (and adding the G-wave changes nothing), and that its dimensional-regularization renormalization yields modified effective range expansion parameters of natural size in all partial waves. The load-bearing premise is the perturbativity of the long-range force and the absence of bound states/low-lying resonances generated by that force alone. This is stated explicitly in Conclusions item (ii): 'we assumed that the long-range force is perturbative and does not create bound states/low-lying resonances alone. ... this one is more restrictive and might require additional scrutiny.' Section 4.2 reinforces this by stating that only repulsive interactions are considered, since for attractive interactions bound states can emerge and the Born series is no longer convergent. This restriction matters because the toy model uses g = 0.073m, for which the Born series converges rapidly (Figs. 7, 8), and the extraction of the finite piece M_fin^ℓ via solving the Lippmann-Schwinger equation for T_fin^L (Eq. 4.1) relies on the same perturbative split of G_L into G_div^L plus G_0 T_fin^L G_0. If the long-range potential is attractive and strong enough to generate a bound state or near-threshold virtual state, the integral equation for T_fin^L has a pole near the physical region, and the decomposition into a finite number of divergent Born terms plus a well-behaved remainder can fail. Yet the paper's broader motivation is pion exchange in systems such as T_cc and baryon-baryon scattering, where OPE is not purely repulsive and attractive channels exist. Thus the central claim about all partial waves and natural-size ERE parameters is not established outside the demonstrated repulsive perturbative regime.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a numerical implementation of the modified L\\\"uscher equation proposed in Ref. [1] for two-body systems with a known long-range potential plus a short-range interaction. It develops a subtraction procedure for the finite-volume Green function that separates exponentially suppressed corrections from a finite-volume remainder, and it introduces a dimensional-regularization scheme with threshold subtractions for the infinite-volume loop function M_\\ell(q_0), giving explicit Feynman-parameter formulas for \\ell=0 and \\ell=4. In a two-Yukawa toy model restricted to the A_1^+ irrep, the authors claim that the modified quantization condition truncated to the S-wave already reproduces the exact finite-volume spectrum and that including the G-wave changes almost nothing, in contrast to the standard L\\\"uscher equation. They further claim that the proposed renormalization scheme yields modified effective range expansion parameters of natural size in all partial waves. The paper checks the M \\to 0 limit against the known Coulomb result and compares the quantization condition with exact plane-wave Hamiltonian eigenvalues.","tokens_in":22833,"tokens_out":7342,"duration_ms":90473,"significance":"If the central claims hold, the paper would provide a practically useful tool for lattice analyses of systems with one-pion exchange and near-threshold left-hand cuts: partial-wave truncation in the modified L\\\"uscher equation would converge much faster than in the standard approach, and dimensional regularization would avoid the large cutoff subtraction terms that make higher-partial-wave effective range parameters unnatural. The paper's strengths are its explicit checks: the result is \\mu-independent in the tested range, the Coulomb limit is reproduced, and the modified quantization condition matches the exact Hamiltonian spectrum in the toy model. The authors also honestly state in Conclusions item (ii) that the long-range force is assumed to be perturbative and not to create bound states or low-lying resonances by itself. However, the regime of validity is narrower than the abstract suggests, and the 'natural size' assertion is not backed by a quantitative extraction of ERE parameters.","major_comments":[{"comment":"The central claim that the modified quantization condition can be truncated to the S-wave and yields natural-size ERE parameters in all partial waves is established only for repulsive, perturbative long-range forces. Section 4.2 explicitly says that only repulsive interactions are considered and that for attractive interactions bound states can emerge and the Born series is no longer convergent; Conclusions item (ii) flags the perturbativity assumption as 'more restrictive and might require additional scrutiny.' This matters because the decomposition in Eq. (4.1), where G_L is split into a finite number of divergent Born terms plus G_0 T_fin^L G_0, and the subsequent numerical solution of the Lippmann-Schwinger equation for T_fin^L, rely on T_fin^L being regular in the physical region. For an attractive OPE-like potential with a near-threshold bound or virtual state, this split can fail. Since the abstract and Section 3.2 advertise the method without this caveat, the paper should either restrict the headline claims to the perturbative repulsive regime or demonstrate, for example with a weakly attractive Yukawa benchmark, that the method remains accurate when the Born series is not convergent.","section":"§4.2 and Conclusions item (ii)"},{"comment":"The abstract claims that the proposed renormalization scheme gives modified effective range expansion parameters of natural size in all partial waves, but the paper never actually extracts or tabulates those parameters. Figure 9 shows Re K_M^\\ell(q_0^2) only for \\ell=0 and in arbitrary units, while Figure 8 shows M_\\ell(q_0) for \\ell=4 but not the resulting ERE coefficients. 'Natural size' cannot be assessed from these plots without specifying the dimensionless ratios in which the coefficients are measured. Please provide a quantitative statement, e.g. the first few coefficients of K_M^\\ell(q_0^2) for \\ell=0 and \\ell=4 in units of the available scales m and M, or revise the abstract so that it does not overstate what is demonstrated.","section":"§4.4 and §4.2"},{"comment":"The headline result that the S-wave-only modified quantization condition 'already reproduces the exact energy level' and that adding the G-wave changes nothing is presented through D(q_0) curves superimposed on vertical lines rather than through quantitative residuals. The reader cannot judge the size of the residual or the actual effect of truncating at \\ell=4. Moreover, the VEGAS integrations used for M_\\ell(q_0) in Figures 7 and 8 are shown without statistical uncertainties, so the claim that higher loops are 'visually indistinguishable' from the full solution is not quantitatively supported. Please quote explicit numbers, such as |E(\\ell_{\\rm max}=0)-E_{\\rm exact}| and |E(\\ell_{\\rm max}=4)-E_{\\rm exact}| in units of q_0^2/M^2 or M, and report integration uncertainties for the numerical results in Section 4.2.","section":"§3.2 and §4.2"}],"minor_comments":[{"comment":"There is a typo in 'algoritm'; it should read 'algorithm'.","section":"Abstract"},{"comment":"The sentence 'The potential of the toy model, which are used to produce the synthetic lattice data' has a subject-verb agreement error; it should be 'which is used'.","section":"§3.2"},{"comment":"In the sentence 'One could use, for example, MS or MS renormalization scheme', the second 'MS' appears to be missing an overline; please clarify whether the \\overline{MS} scheme is intended.","section":"§4.1.1"},{"comment":"The caption refers to the 'left-land threshold'; this should be 'left-hand threshold'.","section":"Fig. 9 caption"},{"comment":"The warning about subthreshold poles of G_\\mu and G_f that cancel in the sum is important, but the statement that 'one can always adjust the free parameter \\mu' is not supported by a concrete criterion. Please state how \\mu should be chosen in practice or show a scan over \\mu for a representative case.","section":"§3.1"},{"comment":"The values plotted for \\ell=4 in Figure 8 are very small (of order 10^{-8}); please specify the normalization and mass dimensions of M_\\ell(q_0) so that the reader can interpret the 'natural size' of the resulting ERE parameters.","section":"§4.2"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the authors are unusually candid about the main limitation, stating in Conclusions item (ii) that the perturbativity of the long-range force is restrictive and needs scrutiny. That makes this a borderline case rather than a clear reject. My view is that the numerical algorithm is sound for the stated regime, but the abstract and Section 3.2 overgeneralize the claims to all long-range forces and to all partial waves. A major revision that restricts the claims or adds a benchmark for an attractive channel, together with a quantitative extraction of ERE parameters, would make the paper acceptable. The paper's citation of the authors' own previous work is appropriate in this context."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is the numerical follow-up to the authors' modified Lüscher formalism, and the new material is genuinely useful: an efficient algorithm for the modified zeta-function and a dimensional-regularization subtraction scheme that avoids the unpleasantly large cutoff-scale polynomials of cutoff regularization. The paper earns its keep. The checks are real: the Coulomb limit is reproduced, the result is μ-independent, and the modified quantization condition matches the exact finite-volume Hamiltonian spectrum in the toy model. The demonstration that S-wave-only truncation already works, while the standard Lüscher approach converges slowly, is convincing within the stated regime. The appendix on partial-wave contributions to energy shifts is a nice touch and explains the uneven convergence pattern quantitatively.\n\nThe soft spots are real but mostly acknowledged by the authors. The load-bearing assumption is that the long-range force is perturbative and does not create bound states or low-lying resonances by itself. They state this explicitly in Conclusions item (ii), and Section 4.2 restricts to repulsive interactions because the Born series can fail for attractive ones. That restriction matters for the stated physics motivation: one-pion exchange in Tcc and baryon-baryon systems has attractive channels, and the method's accuracy there is not established. The claim about natural-size modified effective range parameters \"in all partial waves\" is also extrapolated from ℓ=0 and ℓ=4; the ℓ=4 subtraction scheme is demonstrated for the repulsive Yukawa case, not for an attractive OPE-like channel. Minor issues: the numerical integrations carry no error bars, and no code or data are shipped, which would have made the algorithm reproducible without re-implementation. The self-citation of Ref. [1] is fine—it is the formalism paper, and this is the implementation.\n\nOverall, the paper is honest, technically careful, and useful to lattice practitioners working on systems with long-range forces. The central derivation holds up; the limitations are explicit and addressable. I would send this to a serious referee. It is not ready as is—the attractive-channel question and the missing error bars should be addressed—but it deserves full review rather than desk rejection.","headline":"Solid implementation paper with a clean renormalization scheme and honest checks, but the headline claim about all partial waves is only demonstrated in a repulsive perturbative regime.","tokens_in":23295,"tokens_out":1118,"would_cite":true,"duration_ms":15489,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["11.15.Ha","12.38.Gc","13.75.Cs"],"model":"deepseek-v4-flash","headline":"Treating the long-range force in a plane-wave basis makes the modified Lüscher equation converge so fast that S-wave-only truncation reproduces exact energy levels, while the proposed renormalization keeps parameters natural.","keywords":["modified Lüscher equation","long-range force","effective range expansion","left-hand cut","dimensional regularization","finite-volume quantization","lattice QCD","Yukawa potential"],"falsifier":"A direct falsifier is the attractive counterpart of the toy model: set $g$ to $-g$ in Eq. (3.10) and check whether the modified S-wave-only quantization condition still reproduces the exact Hamiltonian spectrum. The paper itself predicts it will not, because attractive Yukawa forces produce bound states and the Born series stops converging; observing the failure would confirm the stated limitation, while observing success would contradict it. On the repulsive side, rerunning the $M_S=10M$ comparison and checking that the S-wave-only modified solution agrees with the exact level to within the quoted exponentially small corrections would refute the decoupling claim if it departed visibly.","tokens_in":22182,"feed_emoji":"⚛️","tokens_out":9415,"duration_ms":91298,"temperature":0.7,"pith_summary":"This paper gives a practical numerical recipe for the Lüscher quantization condition when the two particles also feel a known long-range force such as one-pion exchange. It argues that treating this long-range part in a plane-wave basis, and expanding only the short-range part in partial waves, makes the partial-wave expansion converge very fast: in the toy model the modified quantization condition with the S-wave alone already reproduces the exact finite-volume energy level, and adding the G-wave changes nothing. The paper further works out the ultraviolet renormalization of the modified zeta-function in dimensional regularization, subtracting a few terms of the Taylor expansion at threshold so that the modified effective-range parameters stay of natural size in all partial waves. If this carries over to real lattice data, scattering analyses on the lattice would need only a small number of partial waves even when a long-range pion exchange is present.","feed_headline":"S-wave only reproduces lattice levels with long-range force","feed_subtitle":"A renormalized zeta-function tames pion exchange, so lattice analyses need few partial waves and survive the left-hand cut.","key_machinery":"The load-bearing object is the modified Lüscher zeta-function, the finite-volume analogue of the loop function $M_\\ell(q_0)$ formed from the Green function of the long-range potential (a sum of one-pion-exchange ladder diagrams). Two decompositions do the work. The free propagator is split into a subtracted piece with no singular denominator in the physical region and a residual piece; this renders the finite-volume shift $\\Delta H$ ultraviolet-finite and reduces the problem to numerically solving a Lippmann-Schwinger equation in a plane-wave basis, with exponentially small corrections dropped. Then, because each term of the Born series of $G_L$ has a lower divergence index than the previous one, $G_L$ is split into a divergent part containing only $2\\ell+2$ terms and a finite part governed by an integral equation for $T_{\\rm fin}$; the divergent part is handled in dimensional regularization and renormalized by subtracting a polynomial at threshold. That subtraction prescription is what keeps the modified effective-range parameters of natural size in all partial waves.","core_discovery":"The central claim is that the modified Lüscher equation—where the known long-range potential is kept in a plane-wave basis and only the unknown short-range part is expanded in partial waves—can be implemented numerically and converges far better than the standard Lüscher equation. In the toy model with a short-range scale $M_S = 10M$, the modified quantization condition with the S-wave only already reproduces the exact energy level, and adding the G-wave does not change anything; even in the borderline case $M_S = 2M$ its convergence is much better than the standard one. The paper also establishes a renormalization scheme: the modified zeta-function's ultraviolet divergences are isolated in a finite number of Born terms, evaluated in dimensional regularization, and removed by subtracting the first terms of the Taylor expansion at threshold. In this scheme the modified effective range expansion parameters are of natural size in all partial waves, and the modified effective range function is smooth and almost linear across the region of the left-hand cut, where the standard K-matrix becomes singular and complex.","pith_inferences":["The successful Coulomb limit suggests the same subtraction scheme could be applied to other massless-exchange forces, where analytic results exist to benchmark the numerics; the paper only verifies the leading logarithm and a few coefficients.","The restriction to repulsive interactions is the main obstacle to pion-exchange channels with attractive partial waves; extending the method would require resumming or analytically continuing past the breakdown of the Born series, a problem the paper leaves open.","The group-theoretic shell analysis implies a practical selection rule: which excited states can be trusted with S-wave-only truncation depends on the lattice-momentum shell, not just on the energy, so users can choose states whose G-wave mixing coefficients are small.","Because real lattice calculations use moving frames and other irreps, the practical payoff for systems like $T_{cc}$ will depend on how well the decoupling survives those generalizations; the paper only demonstrates the center-of-mass $A_1^+$ case."],"forward_implications":["For lattice analyses of systems with one-pion exchange, the number of partial waves one must keep is small; in the toy model S-wave-only truncation already suffices when the short-range force is not too light.","The left-hand-cut region, where the standard Lüscher method fails, becomes accessible: the modified effective range function stays smooth and real there, so energy levels in that region can be analyzed.","Once the modified zeta-function is tabulated for a known long-range potential, extracting scattering information from lattice data proceeds exactly as in the standard Lüscher approach.","Exponentially suppressed finite-volume corrections are much smaller in the modified equation (e.g., a momentum-scale ratio of 11 versus 1.1 for the ground state at $M_S=10M$), improving accuracy near the left-hand cut.","Dimensional regularization with threshold subtraction avoids the unnaturally large cutoff polynomials that arise in higher partial waves, so the effective-range parameters remain of natural size."],"supporting_citations":[{"why":"Provides the modified Lüscher equation and effective-range formalism whose numerical implementation this paper carries out.","marker":"[1]"},{"why":"Defines the original modified effective range function that the whole approach generalizes to finite volume.","marker":"[24]"},{"why":"Establishes the standard finite-volume quantization condition that serves as the comparison baseline throughout the paper.","marker":"[23]"},{"why":"Supply the analytic Coulomb-scattering result used to verify the massless limit of the computed loop function.","marker":"[30, 31]"},{"why":"Give the VEGAS adaptive integration algorithm used for the numerical evaluation of the Feynman-parameter integrals.","marker":"[28, 29]"},{"why":"Documents the earlier observation that cutoff regularization makes modified effective-range parameters unnaturally large, motivating the dimensional-regularization scheme.","marker":"[27]"}],"fun_headline_variants":["Modified Lüscher zeta-function tames long-range forces","S-wave only suffices for lattice levels with long-range force","Renormalized zeta-function gives natural ERE parameters in all partial waves","Modified Lüscher equation needs few partial waves with long-range force"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes the long-range force is weak enough to be handled as a small perturbation and does not by itself create bound states or near-threshold resonances; the authors flag this as their most restrictive assumption, noting it breaks down for attractive forces.","fun_headline_variants_meta":{"raw":{"variants":["Modified Lüscher zeta-function tames long-range forces","S-wave only suffices for lattice levels with long-range force","Renormalized zeta-function gives natural ERE parameters in all partial waves","Modified Lüscher equation needs few partial waves with long-range force"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000637,"raw_usage":{"total_tokens":2913,"prompt_tokens":902,"completion_tokens":2011,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":518,"completion_tokens_details":{"reasoning_tokens":1938}},"tokens_in":518,"tokens_out":2011,"duration_ms":15094,"temperature":1.0,"reasoning_tokens":1938,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:32:13.706462+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct falsifier is the attractive counterpart of the toy model: set $g$ to $-g$ in Eq. (3.10) and check whether the modified S-wave-only quantization condition still reproduces the exact Hamiltonian spectrum. The paper itself predicts it will not, because attractive Yukawa forces produce bound states and the Born series stops converging; observing the failure would confirm the stated limitation, while observing success would contradict it. On the repulsive side, rerunning the $M_S=10M$ comparison and checking that the S-wave-only modified solution agrees with the exact level to within the quoted exponentially small corrections would refute the decoupling claim if it departed visibly.","supporting_citations":[],"review_version":1}