{"id":"e520a9ab-32a6-4329-ab27-0ad28710ea4e","arxiv_id":"2505.02878","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A pole expansion in the uniformization variable expresses two-hadron imaginary-time correlation functions in terms of resonance pole positions and residues, offering a direct lattice QCD analysis method.","lead":"This paper proposes extracting the masses and widths of unstable hadrons, like the rho meson and Lambda(1405), directly from imaginary-time correlation functions in lattice QCD by expressing those correlation functions as a sum of pole terms. The authors demonstrate the expansion on two phenomenological models and argue it can replace multi-step analyses with a direct fit to pole parameters.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mittag-Leffler expansion in Eq. (4) omits the entire function term and assumes away left-hand cuts; the finite-volume gap also blocks direct lattice use, so the central parametrization claim is not yet established.","rationale":"The reader's verdict is CONDITIONAL, and my concern is consistent with that. The paper is honest about its assumptions: it explicitly neglects left-hand cuts, three-body channels, and finite-volume effects. The model demonstrations are internally consistent and show the pole expansion works for the chosen toy interactions, which is real but limited evidence. The most load-bearing gap is not a flaw in the algebra but an unproven assertion that QCD correlators have the required meromorphy and that the entire function in the Mittag-Leffler expansion vanishes. The finite-volume issue is a separate but equally concrete obstruction to the proposed lattice application. Neither gap is shown fatal; both are addressable. Hence the reader's CONDITIONAL verdict stands, and I would not raise it to ACCEPT without either a proof of the asymptotic behavior of D(u) in a realistic model or a synthetic finite-volume test.","tokens_in":13138,"tokens_out":8191,"duration_ms":103753,"concrete_test":"Take the vector-dominance model of Sec. IV, add a one-pion-exchange t-channel diagram so that D(p0) has a left-hand cut, and compute the exact D(k) after uniformization. Then subtract the rho and pi-pi pole contributions of Table I. If the remainder is nonvanishing and does not decrease as additional pole pairs are included, Eq. (4) is missing cut/entire-function contributions, and the pole-only C(τ) parametrization fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formula, Eq. (4), is presented as a Mittag-Leffler expansion, but the theorem states D(u) = sum of principal parts + an entire function E(u). The authors drop E(u) without justification. This is not harmless: even if E(u) is entire in the uniformization variable u, the map u(p0) has square-root branch points, so E(u(p0)) has a nonzero discontinuity across the physical threshold cut and contributes to C(τ) via Eq. (1). Therefore C(τ) is not parametrized only by pole positions and residues unless one proves D(u) -> 0 (or specifies the subtraction constants) as |u|→∞ in the relevant directions. The model demonstrations do not test this: the vector-dominance and chiral-unitary Lagrangians are selected to have no left-hand cuts and only a few poles, so the pole sum reproduces the correlator by construction. In QCD, crossed-channel exchanges produce left-hand cuts in p0, which become branch points in u rather than poles; the paper declares these negligible without a quantitative estimate. In addition, lattice QCD supplies finite-volume correlators with discrete spectra, for which the infinite-volume pole-cut representation of Eq. (5) is not directly valid; the authors concede this at the end but still propose the method. These gaps leave the central claim that C(τ) is exactly a sum of pole terms unverified for the intended application.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes representing the two-hadron imaginary-time correlation function C(τ) in lattice QCD as a sum of pole terms in a uniformization variable. The authors start from the spectral representation of C(τ), introduce a uniformization variable u that makes the two-hadron propagator single-valued, and invoke the Mittag-Leffler theorem to write D(u) as a sum over pole terms. They derive explicit expressions for the single-channel rho-meson case and the two-channel Lambda(1405) case, and they demonstrate the pole expansion in the vector-dominance model and the chiral unitary model, respectively. On this basis, they propose the pole expansion as a new fitting method to extract masses and widths of unstable states from lattice QCD correlation functions. The main claims are that C(τ) is parametrized only by pole positions and residues and that the pole expansion holds for both single-channel and coupled-channel scatterings.","tokens_in":13450,"tokens_out":4715,"duration_ms":66889,"significance":"If the central claim is correct, the paper offers an attractive new route to resonance parameters from Euclidean correlation functions, potentially complementing the Luscher and HAL QCD methods. The uniformization-variable framework is well motivated, and the model comparisons are a useful sanity check. The paper is clearly written and the phenomenological demonstrations are instructive. The significance is presently limited, however, because the derivation relies on an unproven meromorphy assumption, the entire-function term of the Mittag-Leffler expansion is dropped without justification, and the model tests do not exercise the physics that could invalidate the expansion, such as left-hand cuts or finite-volume effects. The paper is honest about the finite-volume issue at the end, but the title and abstract present the method as ready for lattice QCD use.","major_comments":[{"comment":"The Mittag-Leffler expansion of a meromorphic function D(u) is, in general, the sum of principal parts plus an entire function E(u). Eq. (4) omits E(u) entirely, but the paper gives no argument that E(u)=0 or that its contribution to the imaginary-time correlator is negligible. Because u(p0) has square-root branch points at thresholds, even an entire E(u) produces a nonzero discontinuity across the physical cut when composed with u(p0), so E contributes to C(τ) through Eq. (1). Thus the statement that C(τ) is 'parametrized only by the pole positions and residues' is not established. The model demonstrations in Sections 3 and 4 cannot settle this point, since the model amplitudes are pole-dominated by construction. I ask the authors to justify the omission of the entire function, or to show explicitly how it is absorbed or controlled.","section":"Section 2, Eq. (4)"},{"comment":"The demonstrations in the vector-dominance and chiral-unitary models do not test the key assumption that left-hand cuts and three-or-more-body channels can be neglected. In these toy models, the propagators are built from s-channel pole terms plus free two-particle propagation, so the agreement between the direct calculation and the pole sum is essentially a consistency check. To make the central claim convincing, the paper would need a model or a general argument in which left-hand cuts from t- and u-channel exchanges are present and are shown to be representable or negligible in the uniformization variable. A quantitative estimate of the neglected left-hand-cut contribution in a physically motivated model would be a natural addition.","section":"Sections 3 and 4, Eqs. (18) and (23)"},{"comment":"The proposed application is to lattice QCD, where correlation functions are computed in a finite volume and have discrete spectra. The spectral representation in Eq. (1) and the pole expansion in Eq. (4) are infinite-volume statements. The authors acknowledge this at the end and say the finite-volume question should be clarified, but the abstract and title already present the method as a tool for lattice QCD analysis. This gap is load-bearing for the claimed application. The paper would be acceptable either with a finite-volume formulation of the pole expansion or with an explicit statement that the lattice application is conditional on future work on finite-volume effects.","section":"Final paragraph, 'There is one thing which should be clarified'"}],"minor_comments":[{"comment":"Equation (5) as written places the symbol 'Disc' outside the sum and then defines C(τ,u_n) by an integral over p0, which makes the notation ambiguous. Presumably the discontinuity is to be taken on D(u(p0)) before the integration; the formula should be rewritten to make this clear.","section":"Equation (5)"},{"comment":"The two-channel demonstration is shown only for one matrix element, D_{Kbar N, πΣ}, and the agreement is judged visually. A quantitative measure, such as a pointwise relative difference between the direct calculation and the pole-sum result, would strengthen the claim that the expansion 'holds' in the coupled-channel case.","section":"Section 4, Figs. 7 and 8"},{"comment":"The statement that the scattering volume v is 'related to the residues of the poles together with the pole position' is correct, but the step from Eq. (16) to Eq. (17) would benefit from a brief derivation, since the reader must otherwise infer that Im[k^2 D_{ππρ}]/Re[k^2 D_{ππρ}] is being evaluated in the k→0 limit.","section":"Section 3, text near Eq. (17)"},{"comment":"The red lines indicating the 'physical region' are not defined in the captions. Adding a sentence explaining that this line corresponds to the physical p0 values mapped into the k or z plane would improve readability.","section":"Figures 1 and 6"}],"recommendation":"major_revision","confidential_remarks":"The paper is a reasonable contribution to the hep-lat literature, but the central claim is stronger than the evidence. The mathematical omission of the entire function in the Mittag-Leffler expansion and the lack of a quantitative test of the left-hand-cut assumption are the main issues. The authors' own prior work is cited heavily, which is understandable given the continuity, but the novelty relative to Refs. [20-23] should be made more explicit. I would support publication after the major comments are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper does something new: it writes the two-hadron Euclidean correlator as a sum of pole terms in the uniformization variable and gives explicit single- and two-channel residue formulas, tying the residue to the scattering length and reproducing the Maiani-Testa result. The model demonstrations for the rho and Lambda(1405) are clean, and the single-channel one is quantitatively convincing. The authors are also unusually honest about what they assume: no left-hand cuts, no three-body channels, and meromorphy of D(u).\n\nThe soft spots are real. The central formula, Eq. (4), is presented as the Mittag-Leffler expansion, but the theorem gives pole sum plus an entire function. The entire function is dropped without justification. It is not harmless: even if E(u) is entire in u, the map u(p0) has branch points, so E(u(p0)) has a discontinuity across the physical threshold and contributes to C(tau). The model tests cannot see this because vector dominance and the chiral unitary model are chosen so that the correlator is already a sum of poles. In QCD, crossed-channel exchanges generate left-hand cuts, which are branch points in u, not poles; calling them negligible without an estimate is a gap. On top of that, lattice QCD gives finite-volume correlators with discrete spectra; the infinite-volume pole representation is not directly valid there. The authors acknowledge this at the end, but it blocks the proposed fitting method until resolved. The two-channel case is also asserted more than demonstrated: the figures look good, but the text says 'we just point out' and defers details, so there is no quantitative residual to check.\n\nNone of this is fatal if the paper is understood as a proposal. The derivation is standard up to the stated assumptions, and the Maiani-Testa connection is a nice piece of physics. The paper deserves a serious referee, but the referee should ask for a treatment or justification of the omitted entire function, a quantitative estimate of left-hand cuts in the uniformization variable, and a finite-volume analysis or at least a clear statement of when the pole sum is a good approximation in a box. I would want to see those before using the method myself.\n\nMy verdict: send to peer review, conditional. The paper is for lattice hadron spectroscopists and resonance phenomenologists; they will get a clear, useful method proposal with known limitations. I would not cite it as a validated tool, but I might cite it as a proposed alternative.","headline":"A clearly written method proposal that connects pole expansion in the uniformization variable to Euclidean correlators; the main claim is undercut by an unjustified drop of the Mittag-Leffler entire function and an unaddressed finite-volume gap.","tokens_in":13971,"tokens_out":3258,"would_cite":true,"duration_ms":40004,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Gc","13.75.Lb"],"model":"deepseek-v4-flash","headline":"Two-hadron correlation functions are sums of pole terms; resonance masses and widths become direct fit parameters","keywords":["lattice QCD","imaginary-time correlation function","Mittag-Leffler expansion","uniformization variable","hadron resonances","coupled-channel scattering","rho meson","Lambda(1405)"],"falsifier":"Compute the full two-hadron correlation function on a given lattice ensemble, fit it with the pole expansion using a few poles, and compare residuals at short imaginary time: a systematic short-distance deviation, or a fitted pole position that moves with lattice volume or external momentum q, would show that the pole-only parametrization is incomplete.","tokens_in":12948,"feed_emoji":"⚛️","tokens_out":8423,"duration_ms":91543,"temperature":0.7,"pith_summary":"The paper claims that a two-hadron imaginary-time correlation function $C(\\tau)$, of the kind lattice QCD computes, can be represented as a Mittag-Leffler sum over pole terms once the underlying propagator is written in a uniformization variable that makes it single-valued. On this representation, the correlator is parametrized only by the positions and residues of the poles, so fitting lattice data would directly give the complex energies (masses and widths) of unstable states such as the $\\rho$ meson and $\\Lambda(1405)$, along with couplings. The authors verify the expansion numerically in two phenomenological models, the vector-dominance model for $\\rho$ and the chiral unitary model for $\\Lambda(1405)$, and find that the pole sum reproduces the full correlator essentially exactly. If right, this offers a simpler route to resonance parameters from Euclidean correlation functions than multi-step finite-volume analyses.","feed_headline":"Masses and widths become fit parameters in lattice correlators","feed_subtitle":"A uniformized pole expansion makes lattice C(τ) depend only on pole positions and residues, tested on ρ and Λ(1405).","key_machinery":"The load-bearing device is uniformization: for one channel $u = k = \\tfrac{1}{2}\\sqrt{p_0^2 - p^2 - \\varepsilon^2}$, and for two channels $u = z = (k_1 + k_2)/\\Delta$, where $k_i$ are the channel momenta. This variable is chosen so that the threshold branch cuts open and the propagator $D(p_0)$ becomes a single-valued meromorphic function of $u$, allowing the Mittag-Leffler theorem to represent $D(u)$ as a sum of simple pole pairs of the form $(u-u_n)^{-1}$ and $(u+u_n^*)^{-1}$. Substituting this sum into the integral defining $C(\\tau)$ turns the correlator into a sum of known functions $C(\\tau,u_n)$ weighted by residues, so no continuum subtraction or intermediate $K$-matrix parametrization is needed.","core_discovery":"The central discovery is that the imaginary-time correlation function $C_{ij}(\\tau)$, defined by a Laplace transform of the discontinuity of $D_{ij}(p_0)$, can be written as $C_{ij}(\\tau) = \\operatorname{Disc}\\sum_n [r_{ij}^n C(\\tau,u_n) - r_{ji}^{n*} C(\\tau,-u_n^*)]$, where $u$ is the uniformization variable and $C(\\tau,u_n)$ is a known integral. The authors derive this form from the Mittag-Leffler theorem under a meromorphy assumption, explicitly separate dynamical poles ($\\rho$, $\\Lambda(1405)$) from kinematical poles (noninteracting $\\pi\\pi$, $\\bar{K}N$, $\\pi\\Sigma$ pairs), and demonstrate numerically that the sum over just a few poles reproduces both $D$ and $C$ over the plotted range. They conclude that the pole expansion holds for single-channel and coupled-channel (two-channel) scatterings, and they propose it as a method to extract masses and widths of unstable states from lattice QCD correlation functions.","pith_inferences":["Editorial inference: if the meromorphy assumption survives finite-volume checks, the pole-sum fit could replace multi-step analyses for unstable states, reducing model dependence in lattice determinations of widths.","Editorial inference: the same uniformized pole expansion might extend to three-hadron thresholds if a suitable uniformization variable can be constructed, though the paper explicitly leaves three-body channels out.","Editorial inference: a direct numerical test on synthetic lattice correlators generated from a known scattering amplitude could settle whether a few-pole fit recovers the input mass and width at finite volume; the authors themselves flag this as the next step."],"forward_implications":["Lattice QCD analyses of resonances could fit $C(\\tau)$ directly with a handful of pole terms, extracting $m - i\\Gamma/2$ from pole positions without an intermediate finite-volume energy spectrum.","The same parametrization applies to coupled-channel systems, so the several poles of $\\Lambda(1405)$ could be separated from each other and from noninteracting two-hadron contributions.","Large relative momentum between the two hadrons lowers the dynamical pole energy relative to the kinematical pair pole, so such kinematics sharpen the resonance signal in $C(\\tau)$.","Pole residues carry coupling information: in the $\\rho$ case they encode the scattering volume and reproduce the known long-time behavior of the three-point correlation function.","The expansion gives a built-in consistency check: fitted pole positions should not depend on the external momentum $q$ or on which correlation-function component is fitted."],"supporting_citations":[{"why":"It supplies the uniformization procedure that makes the multi-valued propagator D(p0) single-valued.","marker":"[24]"},{"why":"It provides the two-channel uniformization variable z used in the coupled-channel expansion.","marker":"[25]"},{"why":"It establishes the analytical properties of the two-channel S-matrix that justify the single-valued form.","marker":"[26]"},{"why":"It gives the Mittag-Leffler theorem used to expand D as a sum of pole terms.","marker":"[27]"},{"why":"It contains the prior uniformized Mittag-Leffler expansion of the Green's function and T matrix that this paper adapts to correlation functions.","marker":"[21]"},{"why":"It supplies the theorem on the asymptotic behavior of the three-point Euclidean correlation function that the pole expansion reproduces in the q=0 limit.","marker":"[29]"},{"why":"It provides the vector-dominance model used to define the rho interaction Lagrangian for the single-channel demonstration.","marker":"[30]"},{"why":"It gives the model parameters and calculation details for the rho-pole demonstration.","marker":"[31]"},{"why":"It provides the chiral unitary coupled-channel approach for the pi-Sigma and K-bar-N interaction used in the Lambda(1405) demonstration.","marker":"[33]"},{"why":"It supplies the renormalized unitarized interaction without on-shell factorization used for the Lambda(1405) coupled-channel calculation.","marker":"[35]"}],"fun_headline_variants":["Pole expansion extracts rho and Lambda(1405) parameters","Lattice correlators fit pole residues for unstable hadrons","Uniformized pole expansion yields masses and widths","Few poles capture full lattice correlator accurately","New analysis: pole positions from imaginary-time correlators"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation rests on the correlation function being meromorphic in the uniformization variable after left-hand cuts and three-or-more-body channels are neglected; if real correlation functions have those cuts, or if finite-volume lattice data do not obey an infinite-volume pole sum, the proposed fit would fail.","fun_headline_variants_meta":{"raw":{"variants":["Pole expansion extracts rho and Lambda(1405) parameters","Lattice correlators fit pole residues for unstable hadrons","Uniformized pole expansion yields masses and widths","Few poles capture full lattice correlator accurately","New analysis: pole positions from imaginary-time correlators"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000254,"raw_usage":{"total_tokens":1564,"prompt_tokens":935,"completion_tokens":629,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":551,"completion_tokens_details":{"reasoning_tokens":553}},"tokens_in":551,"tokens_out":629,"duration_ms":7848,"temperature":1.0,"reasoning_tokens":553,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:53:16.703604+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full two-hadron correlation function on a given lattice ensemble, fit it with the pole expansion using a few poles, and compare residuals at short imaginary time: a systematic short-distance deviation, or a fitted pole position that moves with lattice volume or external momentum q, would show that the pole-only parametrization is incomplete.","supporting_citations":[{"cited_title":"A New Method to Extract Information of Near-Threshold Resonances: Uniformized Pole-Sum Representation of Green's Function and T-matrix","cited_arxiv_id":"2005.07022","evidence_quote":"It supplies the uniformization procedure that makes the multi-valued propagator D(p0) single-valued."},{"cited_title":"Near-threshold Spectrum from Uniformized Mittag-Leffler Expansion -Pole Structure of $Z(3900)$-","cited_arxiv_id":"2108.11605","evidence_quote":"It provides the two-channel uniformization variable z used in the coupled-channel expansion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives the Mittag-Leffler theorem used to expand D as a sum of pole terms."},{"cited_title":"Kato, Analytical properties of two-channel S-matrix, Annals Phys","cited_arxiv_id":null,"evidence_quote":"It supplies the theorem on the asymptotic behavior of the three-point Euclidean correlation function that the pole expansion reproduces in the q=0 limit."},{"cited_title":"Arfken, H","cited_arxiv_id":null,"evidence_quote":"It provides the vector-dominance model used to define the rho interaction Lagrangian for the single-channel demonstration."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives the model parameters and calculation details for the rho-pole demonstration."}],"review_version":1}