REVIEW 3 major objections 4 minor 37 references
Pole-Expansion of Two-Hadron Imaginary-Time Correlation Function -a new method of analysis for unstable states in lattice QCD-
T0 review · 3 major / 4 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read Two-hadron correlation functions are sums of pole terms; resonance masses and widths become direct fit parameters
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Section 2, Eq. (4)] 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.
- [Sections 3 and 4, Eqs. (18) and (23)] 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.
- [Final paragraph, 'There is one thing which should be clarified'] 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.
minor comments (4)
- [Equation (5)] 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 4, Figs. 7 and 8] 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 3, text near Eq. (17)] 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.
- [Figures 1 and 6] 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.
Circularity Check
No significant circularity: Eq. (4)-(5) are a direct Mittag-Leffler consequence under an explicitly stated meromorphy assumption; the model comparisons are independent checks, and the remaining gaps are validity assumptions rather than circular reductions.
full rationale
The derivation chain is: assume meromorphy of D(u) after uniformization (explicitly stated), apply the Mittag-Leffler theorem to write D as a sum of pole terms, then substitute into the spectral representation to obtain C(tau). This is a direct mathematical consequence and does not define an output in terms of a fitted input. The numerical demonstrations compare the pole-sum form against separately computed model correlators (vector-dominance and chiral-unitary), not against pole parameters obtained by fitting C(tau), so there is no fitted-input-called-prediction. The paper does cite the authors' prior works [20-23] for the method and uses Ref. [21] for calculational details, which is self-citation, but it is not load-bearing here because the pole-expansion formulas and the comparisons are presented in this paper. The principal limitations (the dropped entire function in Eq. (4), neglected left-hand cuts and three-body channels, and the conceded finite-volume question) are assumptions or validity gaps for the proposed lattice-QCD application, not circular reductions. No quoted equality reduces to itself by construction, so no significant circularity is found.
Assumptions & free parameters
free parameters (5)
- m_pi = 0.140 GeV =
0.140 GeV
- m_rho = 0.760 GeV, Gamma_rho = 0.160 GeV =
0.760 GeV, 0.160 GeV
- g = 6.05 =
6.05
- f = 0.104 GeV =
0.104 GeV
- C matrix for chiral unitary model =
[[3, -sqrt(3/2)], [-sqrt(3/2), 4]]
assumptions (5)
- domain assumption Neglect of left-hand cuts and three-or-more-body channels
- domain assumption D(p0) is a meromorphic function of the uniformization variable u
- standard math Mittag-Leffler theorem for meromorphic functions
- standard math Spectral representation C(τ) = ∫ e^{-p0 τ} Disc D(p0)
- domain assumption P and CP symmetry impose relations between residues
Cite this review
Pith. "Pith review of Pole-Expansion of Two-Hadron Imaginary-Time Correlation Function -a new method of analysis for unstable states in lattice QCD-." pith.science (2026). https://pith.science/paper/4NXMFXZB
@misc{pith2026250502878,
author = {Pith},
title = {Pith review of: Pole-Expansion of Two-Hadron Imaginary-Time Correlation Function -a new method of analysis for unstable states in lattice QCD-},
year = {2026},
howpublished = {\url{https://pith.science/paper/4NXMFXZB}},
note = {Machine review of arXiv:2505.02878}
}
abstract
We analyze the pole expansion of the two-hadron imaginary-time correlation function. We first explain the general idea that the imaginary-time correlation function is expressed as a sum of the pole terms, the Mittag-Leffler expansion, in terms of the uniformization variable, which makes the S-matrix single-valued. We then derive explicit expressions of the pole expansion for the single-channel ($\rho$ meson) and two-channel ($\Lambda(1405)$) examples and demonstrate that the pole expansion actually holds employing phenomenological models, the vector-dominance model for the $\rho$ meson and the chiral unitary model for $\Lambda(1405)$. From this observation we propose the pole expansion as a method to extract information of unstable states such as masses and widths from the two-hadron imaginary-time correlation functions obtained by lattice QCD simulations.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
R. A. Briceno, J. J. Dudek, and R. D. Young, Scattering processes and resonances from lattice QCD, Rev. Mod. Phys. 90, 025001 (2018), arXiv:1706.06223 [hep-lat]
arXiv 2018
-
[2]
are masses of hadrons in the channel 1 (2): u =z = k1 +k2 ∆ = 1 2 (p2 0− p2−ε2 1)1/2 + (p2 0− p2−ε2 2)1/2 ∆ , (3) where ∆ = 1 2 p ε2 2−ε2 1 [25, 26]. In terms of the uniformization variable, u, Dij(u) can be expressed as a sum of the pole terms, the Mittag-Leffler expansion [27, 28], as Dij(u) = X n rij n u−un − rji∗ n u +u∗ n . (4) Substituting Eq.(4) in...
-
[3]
S. Aoki et al. (CP-PACS), Lattice QCD Calculation of the rho Meson Decay Width, Phys. Rev. D 76, 094506 (2007), arXiv:0708.3705 [hep-lat]
arXiv 2007
-
[4]
X. Feng, K. Jansen, and D. B. Renner, Resonance Parameters of the rho-Meson from Lattice QCD, Phys. Rev. D 83, 094505 (2011), arXiv:1011.5288 [hep-lat]
arXiv 2011
-
[5]
S. Aoki et al. (CS), ρ Meson Decay in 2+1 Flavor Lattice QCD, Phys. Rev. D 84, 094505 (2011), arXiv:1106.5365 [hep-lat]
arXiv 2011
-
[6]
C. B. Lang, D. Mohler, S. Prelovsek, and M. Vidmar, Coupled channel analysis of the rho meson decay in lattice QCD, Phys. Rev. D 84, 054503 (2011), [Erratum: Phys.Rev.D 89, 059903 (2014)], arXiv:1105.5636 [hep-lat]
arXiv 2011
-
[7]
J. J. Dudek, R. G. Edwards, and C. E. Thomas (Hadron Spectrum), Energy dependence of the ρ resonance in ππ elastic scattering from lattice QCD, Phys. Rev. D 87, 034505 (2013), [Erratum: Phys.Rev.D 90, 099902 (2014)], arXiv:1212.0830 [hep-ph]
arXiv 2013
-
[8]
C. Alexandrou, L. Leskovec, S. Meinel, J. Negele, S. Paul, M. Petschlies, A. Pochinsky, G. Ren- don, and S. Syritsyn,P -waveππ scattering and theρ resonance from lattice QCD, Phys. Rev. D 96, 034525 (2017), arXiv:1704.05439 [hep-lat]
arXiv 2017
Show all 37 references
-
[9]
Andersen, J
C. Andersen, J. Bulava, B. H¨ orz, and C. Morningstar, The I = 1 pion-pion scattering am- plitude and timelike pion form factor from Nf = 2 + 1 lattice QCD, Nucl. Phys. B 939, 145 (2019), arXiv:1808.05007 [hep-lat]
2019 arXiv
-
[10]
D. J. Wilson, R. A. Briceno, J. J. Dudek, R. G. Edwards, and C. E. Thomas, Coupled ππ,K ¯K scattering in P -wave and theρ resonance from lattice QCD, Phys. Rev. D 92, 094502 (2015), arXiv:1507.02599 [hep-ph]
2015 arXiv
-
[11]
J. J. Dudek, R. G. Edwards, and D. J. Wilson (Hadron Spectrum), An a0 resonance in strongly coupled πη, KK scattering from lattice QCD, Phys. Rev. D 93, 094506 (2016), arXiv:1602.05122 [hep-ph]
2016 arXiv
-
[12]
J. J. Dudek, R. G. Edwards, C. E. Thomas, and D. J. Wilson (Hadron Spectrum), Resonances in coupledπK−ηK scattering from quantum chromodynamics, Phys. Rev. Lett. 113, 182001 (2014), arXiv:1406.4158 [hep-ph]
2014 arXiv
-
[13]
Bulava et al
J. Bulava et al. (Baryon Scattering (BaSc)), Two-Pole Nature of the Λ(1405) resonance from Lattice QCD, Phys. Rev. Lett. 132, 051901 (2024), arXiv:2307.10413 [hep-lat]
2024
-
[14]
Bulava et al
J. Bulava et al. (Baryon Scattering (BaSc)), Lattice QCD study of πΣ-K¯N scattering and the Λ(1405) resonance, Phys. Rev. D 109, 014511 (2024), arXiv:2307.13471 [hep-lat]
2024
-
[15]
G. Moir, M. Peardon, S. M. Ryan, C. E. Thomas, and D. J. Wilson, Coupled-Channel Dπ, Dη and Ds ¯K Scattering from Lattice QCD, JHEP 10, 011, arXiv:1607.07093 [hep-lat]
-
[17]
Stable Particle States, Commun. Math. Phys. 104, 177 (1986)
1986
-
[18]
Luscher, Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories
M. Luscher, Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories
-
[19]
Scattering States, Commun. Math. Phys. 105, 153 (1986)
1986
-
[20]
Rummukainen and S
K. Rummukainen and S. A. Gottlieb, Resonance scattering phase shifts on a nonrest frame lattice, Nucl. Phys. B 450, 397 (1995), arXiv:hep-lat/9503028
1995 arXiv
-
[21]
Ishii, S
N. Ishii, S. Aoki, and T. Hatsuda, The Nuclear Force from Lattice QCD, Phys. Rev. Lett. 99, 022001 (2007), arXiv:nucl-th/0611096
2007 arXiv
-
[22]
S. Aoki, T. Hatsuda, and N. Ishii, Theoretical Foundation of the Nuclear Force in QCD and its applications to Central and Tensor Forces in Quenched Lattice QCD Simulations, Prog. Theor. Phys. 123, 89 (2010), arXiv:0909.5585 [hep-lat]
2010 arXiv
-
[23]
Humblet and L
J. Humblet and L. Rosenfeld, Theory of nuclear reactions: I. resonant states and collision matrix, Nuclear Physics 26, 529 (1961)
1961
-
[24]
Yamada and O
W. Yamada and O. Morimatsu, New method to extract information of near-threshold reso- nances: Uniformized Mittag-Leffler expansion of Green’s function and T matrix, Phys. Rev. C 102, 055201 (2020), arXiv:2005.07022 [hep-ph]
2020 arXiv
-
[25]
W. A. Yamada, O. Morimatsu, T. Sato, and K. Yazaki, Near-threshold spectrum from a uniformized Mittag-Leffler expansion: Pole structure of the Z(3900), Phys. Rev. D105, 014034 (2022), arXiv:2108.11605 [hep-ph]
2022 arXiv
-
[26]
W. A. Yamada, O. Morimatsu, T. Sato, and K. Yazaki, Survival probability of unstable states in coupled-channels: Nonexponential decay of threshold cusp, Phys. Rev. D 108, L071502 (2023), arXiv:2305.11695 [hep-ph]
2023 arXiv
-
[27]
R. G. Newton, Scattering Theory of Waves and Particles (Springer, 1982)
1982
-
[28]
R. G. Newton, Inelastic scattering, Annals of Physics 4, 29 (1958)
1958
-
[29]
Kato, Analytical properties of two-channel S-matrix, Annals Phys
M. Kato, Analytical properties of two-channel S-matrix, Annals Phys. 31, 130 (1965)
1965
-
[30]
Arfken, H
G. Arfken, H. Weber, and F. Harris,Mathematical Methods for Physicists: A Comprehensive Guide (Elsevier Science, 2013)
2013
-
[31]
E. T. Whittaker and G. N. Watson,A Course of Modern Analysis: an Introduction to the General Theory of Infinite Processes and of Analytic Functions; with an Account of the Principal Transcendental Functions (University Press, 1928)
1928
-
[32]
Maiani and M
L. Maiani and M. Testa, Final state interactions from Euclidean correlation functions, Phys. Lett. B 245, 585 (1990)
1990
-
[33]
J. J. Sakurai, Theory of strong interactions, Annals Phys. 11, 1 (1960)
1960
-
[34]
Klingl, N
F. Klingl, N. Kaiser, and W. Weise, Effective Lagrangian approach to vector mesons, their structure and decays, Z. Phys. A 356, 193 (1996), arXiv:hep-ph/9607431
1996 arXiv
-
[35]
Kaiser, T
N. Kaiser, T. Waas, and W. Weise, SU(3) chiral dynamics with coupled channels: Eta and kaon photoproduction, Nucl. Phys. A 612, 297 (1997), arXiv:hep-ph/9607459
1997 arXiv
-
[36]
Oset and A
E. Oset and A. Ramos, Nonperturbative chiral approach to s wave anti-K N interactions, Nucl. Phys. A 635, 99 (1998), arXiv:nucl-th/9711022
1998 arXiv
-
[37]
Hyodo and D
T. Hyodo and D. Jido, The nature of the Lambda(1405) resonance in chiral dynamics, Prog. Part. Nucl. Phys. 67, 55 (2012), arXiv:1104.4474 [nucl-th]
2012 arXiv
-
[38]
Morimatsu and K
O. Morimatsu and K. Yamada, Renormalization of the unitarized Weinberg-Tomozawa inter- action without on-shell factorization and I=0 K ¯N–πΣ coupled channels, Phys. Rev. C 100, 025201 (2019), arXiv:1903.12380 [hep-ph]
2019 arXiv
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.