REVIEW 3 major objections 2 minor 41 references
Results on meson-meson scattering at large $N_\text{c}$
T0 review · 3 major / 2 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read This lattice study of meson-meson scattering at Nc=3-6 finds no tetraquark resonance in the AA channel, but a virtual bound state at Eb/Mpi=1.741(13) at Nc=3, with the expected 1/Nc suppression and visible subleading corrections.
desk verdict Solid preliminary large-Nc scattering study whose headline virtual bound state is a parametrization-dependent hint rather than an established result. 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 argument is carried by a modified effective range expansion with a fixed Adler zero, $k/M_\pi\,\cot\delta_0 = (M_\pi E/(E^2-2z))(B_0+B_1 k^2/M_\pi^2)$, with $z=M_\pi^2$, used to fit the finite-volume energy levels through Lüscher's quantization condition. The energies themselves come from a matrix of correlation functions built from two-pion, two-vector-meson, and local tetraquark operators, solved with a generalized eigenvalue problem; the large-$N_c$ power counting of eqs. (1)-(2), expressing the $SS$/$AA$ amplitudes as an opposite-sign $1/N_c$ term plus subleading $N_f$-dependent corrections and the $AS$ amplitude as $O(1/N_c^2)$, organizes the comparison across $N_c$.
What would settle it
Rerun the $AA$-channel analysis on the same finite-volume energies using a different amplitude parametrization, for example one with a $B_2$ term or with a free Adler zero, and check whether a pole at $E_b/M_\pi\approx 1.741$ remains; if the pole moves by more than the quoted uncertainty or disappears, the virtual-bound-state claim is falsified.
Extended reading notes
Core claim
In the theory with $N_f=4$ degenerate flavors and $M_\pi\approx 590$ MeV, the $AA$-channel pion-pion amplitude has no pole above the two-pion threshold for any $N_c$ studied, and at $N_c=3$ it is consistent with a virtual bound state at $E_b/M_\pi=1.741(13)$ beneath threshold. The $SS$ and $AA$ phase shifts, after removing the leading $1/N_c$ factor, are of similar magnitude across $N_c$, with visible subleading corrections appearing as a linear dependence of the scattering length and effective range on $1/N_c$ for $N_c=4,5,6$. The $AS$ channel interacts very weakly, in line with the expected $O(1/N_c^2)$ suppression. The paper is explicit that the virtual bound state needs confirmation with other amplitude parametrizations.
Load-bearing premise
The analysis assumes that the effective-range formula with the Adler zero fixed at $z=M_\pi^2$ and truncated after the $B_1$ term is the true shape of the pion-pion amplitude over the fitted range, so that the virtual bound state is a real feature of the theory and not an artifact of the fitting function.
Editorial extensions
If this is right
- If the result is right, the $AA$ channel at $N_c=3$ has a subthreshold pole at $E_b/M_\pi=1.741(13)$, i.e. a virtual bound state rather than a resonance above threshold.
- The $SS$ and $AA$ amplitudes follow the leading $1/N_c$ scaling, and the residual $N_c$ dependence can be used to determine the $N_c$ scaling of low-energy constants in one-loop chiral perturbation theory.
- The absence of an above-threshold resonance in the $AA$ channel at this pion mass tells experiment and phenomenology that a tetraquark signal, if present, is a subthreshold or narrower effect in this setup.
- The very weak $AS$ interactions confirm the kinematic $O(1/N_c^2)$ suppression, so that channel is a clean laboratory for $N_c$ counting rather than for resonance hunting.
Reading between the lines
- If the virtual bound state survives alternative parametrizations, following its pole trajectory in $N_c$ would distinguish a subleading-effect state (moving toward threshold and disappearing as $N_c\to\infty$) from an intrinsic tetraquark that stays fixed.
- A robust pole would also allow a compositeness analysis of its residue, separating a meson-meson molecular component from an elementary tetraquark component.
- The apparent $N_c=6$ separation in the $SS$ channel, which the authors flag as 'unnaturally separated', suggests a systematic effect such as the pion-mass mismatch between ensembles; correcting for the $M_\pi$ dependence would sharpen the $N_c\to\infty$ extrapolation.
- The same operator set and analysis pipeline could be applied to the matching doubly charmed channels at the physical pion mass to test whether the LHCb $T_{cs0}(2900)$ states have the same origin as the $N_c=3$ virtual state seen here.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a preliminary lattice study of meson-meson scattering in a theory with N_f=4 degenerate quark flavors and N_c=3--6, at a pion mass of about 590 MeV. Finite-volume energies are extracted from a large operator basis that includes two-pion, two-vector-meson, and local tetraquark operators, and infinite-volume phase shifts are obtained via the L\"uscher quantization condition. The s-wave phase shifts in the SS and AA channels are fitted to a modified effective-range expansion with an Adler zero fixed at z=M_pi^2. The main physics results are: (i) the AA channel at N_c=3 is consistent with a virtual bound state at E_b/M_pi=1.741(13), while no tetraquark resonance is found above threshold; (ii) the amplitudes show the expected 1/N_c suppression with visible subleading corrections; and (iii) the AS channel is very weakly interacting. The paper also studies the N_c dependence of the scattering length and effective range.
Significance. If established, these results would be an important step toward understanding whether the exotic tetraquark candidates seen by LHCb survive in the large-N_c limit, and they would provide a first multi-color lattice determination of subleading 1/N_c corrections to meson-meson scattering beyond threshold. The analysis is methodologically careful in several respects: it uses a large operator set, a GEVP with AIC-based averaging, multiple moving frames and irreps, and a consistent action for valence and sea quarks. The authors are also appropriately cautious in their wording, explicitly flagging the main model-dependence. However, the headline virtual-bound-state claim currently rests on a single two-parameter fit form and is not yet supported by robustness checks.
major comments (3)
- [Sec. 3, Eq. (7) and Fig. 2] The virtual bound state at E_b/M_pi=1.741(13) is inferred from the modified effective-range expansion of Eq. (7) truncated after the B1 term, with the Adler zero fixed at z=M_pi^2. The pole is located at k^2/M_pi^2=-0.24, roughly halfway between the imposed zero at k^2/M_pi^2=-0.5 and the two-pion threshold, a region not directly constrained by any lattice energy level. The manuscript provides no fit-quality information (e.g., chi2/dof or the fitted energies) and no test with alternative parametrizations; the authors themselves state in Sec. 4 that the existence of the state needs to be established with other parametrizations. Because this is the sole evidence for the central result, the paper should either add such robustness tests or present the virtual bound state explicitly as a model-dependent indication rather than as a result.
- [Sec. 3, Fig. 4] The linear extrapolation in 1/N_c for the scattering length and effective range uses the N_c=4--6 points while the text attributes the discrepancy between the SS and AA channels at large N_c to 'mismatches in the pion masses between ensembles'. Since M_pi enters both the ChPT normalization and the two-body kinematics, fitting in 1/N_c without including the pion-mass dependence can bias the extracted subleading coefficients. The authors note this limitation and plan a one-loop ChPT fit, but as it stands the quantitative claim about subleading corrections is not established.
- [Sec. 3, Fig. 3 and Fig. 4] The text states that in the SS channel 'results for N_c=6 seem to be unnaturally separated from the rest, which we are currently investigating', yet this N_c=6 point is included in the linear extrapolation to N_c=infinity in Fig. 4. The sensitivity of the extrapolated values to excluding N_c=6, or to the choice of fit range, should be quantified. Without this, the claimed O(1/N_c^2) corrections are not robust against the anomalous N_c=6 behavior.
minor comments (2)
- [Figures 2 and 3] The panel labels '(a) AA channel' and '(b) SS channel' appear to be swapped relative to the text, which describes figs. 2a and 2b as the SS and AA channels, respectively; the signs of the phase shifts (negative for SS, positive for AA) confirm the swap.
- [Eq. (7)] The term 'pole' in 'we keep the location of the pole fixed at z=M_pi^2' is confusing, since z is an Adler zero of the parametrization, not the pole that would correspond to a bound state. Calling z an Adler zero explicitly would improve clarity.
Circularity Check
No significant circularity: the finite-volume energies are independent lattice inputs, and the virtual bound state is a fitted output whose parametrization dependence is explicitly flagged.
full rationale
The paper's derivation chain is data-driven rather than self-referential. Lattice correlation functions are computed with a large operator set, finite-volume energies are extracted via GEVP fits, and these energies are then converted into infinite-volume scattering amplitudes through the Lüscher quantization condition, Eq. (5), and its partial-wave form, Eq. (6). The modified effective range expansion, Eq. (7), is an assumed parametrization with the Adler zero fixed at z = M_pi^2 and higher-order terms set to zero; the parameters B0 and B1 are fitted to the finite-volume energy data. The resulting phase shifts, scattering lengths, effective ranges, and the Nc = 3 virtual bound state are all outputs of these fits, not inputs used to define the parametrization. The virtual bound state is a pole of the fitted amplitude below threshold, so it depends on the choice of parametrization, but that is a model-robustness issue rather than circularity. The paper explicitly acknowledges this: 'its existence needs to be established also with other parametrizations of the scattering amplitude' (Sec. 4), which is an honest limitation, not evidence that the claim is built into the assumptions. The self-citations to refs. [2-5] provide the simulation setup, ensembles, and a previous near-threshold analysis, but the new finite-volume spectra, phase shifts, and Nc dependence presented here are newly computed quantities. No load-bearing step reduces to a self-cited uniqueness theorem or to a renamed known result. Therefore the paper receives a circularity score of 0.
Assumptions & free parameters
free parameters (1)
- B0, B1 (eq. 7 coefficients) for each channel and Nc =
Not reported in the paper; derived AA pole at Nc=3: E_b/M_pi = 1.741(13)
assumptions (5)
- standard math Lüscher quantization condition, eqs. (5)-(6), maps finite-volume energies to infinite-volume phase shifts.
- domain assumption Single-channel elastic pi-pi scattering is valid for the energy levels used in the amplitude fits.
- domain assumption The N_f=4 degenerate flavor theory with SU(4)_f irreps is a useful proxy for mapping to the LHCb tetraquark channels.
- ad hoc to paper The modified effective range expansion, eq. (7), truncated after B1 with the Adler zero fixed at z = M_pi^2, describes the scattering amplitude over the fitted energy range.
- domain assumption Pion masses are sufficiently matched across Nc ensembles so that Nc scaling is not contaminated by M_pi differences.
invented entities (1)
-
AA-channel virtual bound state near E_b/M_pi = 1.741(13) at Nc=3
Cite this review
Pith. "Pith review of Results on meson-meson scattering at large $N_\text{c}$." pith.science (2026). https://pith.science/paper/BNLVDFN3
@misc{pith2026250119115,
author = {Pith},
title = {Pith review of: Results on meson-meson scattering at large $N_\textc$},
year = {2026},
howpublished = {\url{https://pith.science/paper/BNLVDFN3}},
note = {Machine review of arXiv:2501.19115}
}
abstract
We present results on the large $N_\text{c}$ scaling of meson-meson scattering amplitudes. We work in a theory with $N_\text{f}=4$ degenerate quark flavors and run lattice simulations with $N_\text{c}=3-6$ and pion mass $M_\pi\approx 590$ MeV. We focus on three different scattering channels, two of which have the same quantum numbers as some tetraquark candidates recently found at LHCb. Finite-volume energies are extracted using a large set of operators, containing two-particle operators corresponding to two pions or two vector mesons, and local tetraquark operators. Using L\"uscher's quantization condition, we constrain the infinite-volume scattering amplitudes and investigate subleading $N_\text{c}$ corrections to the large $N_\text{c}$ limit. For one of the channels, we find indications of a virtual bound state at $N_\text{c}=3$, which may be related to one of the aforementioned exotic states.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[5]
A lattice study of $\pi\pi$ scattering at large $N_\text{c}$
J. Baeza-Ballesteros, P. Hernández and F. Romero-López,A lattice study of𝜋𝜋 scattering at large N𝑐, JHEP06(2022) 049 [2202.02291]
work page Pith review arXiv 2022
-
[1]
’t Hooft,A Planar Diagram Theory for Strong Interactions, Nucl
G. ’t Hooft,A Planar Diagram Theory for Strong Interactions, Nucl. Phys. B72(1974) 461
1974
-
[2]
P.HernándezandF.Romero-López, Thelarge 𝑁𝑐 limitofQCDonthelattice ,Eur.Phys.J.A 57 (2021) 52 [2012.03331]
arXiv 2021
-
[3]
Large $N_c$ scaling of meson masses and decay constants
P. Hernández, C. Pena and F. Romero-López,Large𝑁𝑐 scaling of meson masses and decay constants,Eur. Phys. J. C79(2019) 865 [1907.11511]
work page Pith review arXiv 2019
-
[4]
Dissecting the $\Delta I= 1/2$ rule at large $N_c$
A. Donini, P. Hernández, C. Pena and F. Romero-López,Dissecting theΔ𝐼 = 1/2 rule at large𝑁𝑐,Eur. Phys. J. C80 (2020) 638 [2003.10293]
work page Pith review arXiv 2020
-
[6]
Curve collapse for the isospin-2 pion scattering length from QCD with 3, 4, and 5 colors
T. DeGrand,Curve collapse for the isospin-2 pion scattering length from QCD with 3, 4, and 5 colors, 2409.02242
-
[7]
J.R. Pelaez and G. Rios,Nature of the f0(600) from its N(c) dependence at two loops in unitarized Chiral Perturbation Theory, Phys. Rev. Lett.97(2006) 242002 [hep-ph/0610397]
arXiv 2006
-
[8]
Truong,Chiral Perturbation Theory and Final State Theorem, Phys
T.N. Truong,Chiral Perturbation Theory and Final State Theorem, Phys. Rev. Lett.61 (1988) 2526
work page 1988
Show all 41 references
-
[9]
Dobado, M.J
A. Dobado, M.J. Herrero and T.N. Truong,Unitarized Chiral Perturbation Theory for Elastic Pion-Pion Scattering,Phys. Lett. B235 (1990) 134
1990
-
[10]
Witten,Baryons in the 1/n Expansion, Nucl
E. Witten,Baryons in the 1/n Expansion, Nucl. Phys. B160 (1979) 57
1979
-
[11]
Coleman,1/N, in17th International School of Subnuclear Physics: Pointlike Structures Inside and Outside Hadrons, 3, 1980
S.R. Coleman,1/N, in17th International School of Subnuclear Physics: Pointlike Structures Inside and Outside Hadrons, 3, 1980
1980
-
[12]
Weinberg,Tetraquark Mesons in Large𝑁 Quantum Chromodynamics, Phys
S. Weinberg,Tetraquark Mesons in Large𝑁 Quantum Chromodynamics, Phys. Rev. Lett. 110 (2013) 261601 [1303.0342]
2013 arXiv
-
[13]
Knecht and S
M. Knecht and S. Peris,Narrow Tetraquarks at Large N, Phys. Rev. D88(2013) 036016 [1307.1273]
2013 arXiv
-
[14]
Cohen and R.F
T.D. Cohen and R.F. Lebed,Are There Tetraquarks at Large𝑁𝑐 in QCD(F)?,Phys. Rev. D 90 (2014) 016001 [1403.8090]
2014 arXiv
-
[15]
Bijnens and J
J. Bijnens and J. Lu,Meson-meson Scattering in QCD-like Theories, JHEP03(2011) 028 [1102.0172]
2011 arXiv
-
[16]
LHCb collaboration, A model-independent study of resonant structure in𝐵+→𝐷+𝐷−𝐾+ decays,Phys. Rev. Lett.125 (2020) 242001 [2009.00025]. 9 Results on meson-meson scattering at large𝑁c Jorge Baeza-Ballesteros
2020
-
[17]
LHCb collaboration, Amplitude analysis of the𝐵+→𝐷+𝐷−𝐾+ decay, Phys. Rev. D102 (2020) 112003 [2009.00026]
2020
-
[18]
LHCb collaboration, First Observation of a Doubly Charged Tetraquark and Its Neutral Partner, Phys. Rev. Lett.131 (2023) 041902 [2212.02716]
2023
-
[19]
LHCb collaboration, Amplitude analysis of B0→D¯0Ds+𝜋- and B+→D-Ds+𝜋+ decays, Phys. Rev. D108(2023) 012017 [2212.02717]
2023
-
[20]
Molina and E
R. Molina and E. Oset,Tcs¯(2900) as a threshold effect from the interaction of the D*K*, Ds*𝜌 channels, Phys. Rev. D107 (2023) 056015 [2211.01302]
2023 arXiv
-
[21]
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. 2. Scattering States,Commun. Math. Phys.105(1986) 153
1986
-
[22]
Luscher,Two particle states on a torus and their relation to the scattering matrix,Nucl
M. Luscher,Two particle states on a torus and their relation to the scattering matrix,Nucl. Phys. B354 (1991) 531
1991
-
[23]
Del Debbio, A
L. Del Debbio, A. Patella and C. Pica,Higher representations on the lattice: Numerical simulations. SU(2) with adjoint fermions,Phys. Rev. D81 (2010) 094503 [0805.2058]
2010 arXiv
-
[24]
Del Debbio, B
L. Del Debbio, B. Lucini, A. Patella, C. Pica and A. Rago,Conformal versus confining scenario in SU(2) with adjoint fermions, Phys. Rev. D80(2009) 074507 [0907.3896]
2009 arXiv
-
[25]
Hadron Spectrumcollaboration,Energy dependence of the𝜌 resonance in𝜋𝜋 elastic scattering from lattice QCD,Phys. Rev. D87(2013) 034505 [1212.0830]
2013 arXiv
-
[26]
Detmold, D.J
W. Detmold, D.J. Murphy, A.V. Pochinsky, M.J. Savage, P.E. Shanahan and M.L. Wagman, Sparsening algorithm for multihadron lattice QCD correlation functions, Phys. Rev. D104 (2021) 034502 [1908.07050]
2021 arXiv
-
[27]
Dudek, R.G
J.J. Dudek, R.G. Edwards and C.E. Thomas,S and D-wave phase shifts in isospin-2 pi pi scattering from lattice QCD,Phys. Rev. D86(2012) 034031 [1203.6041]
2012 arXiv
-
[28]
Michael and I
C. Michael and I. Teasdale,Extracting Glueball Masses From Lattice QCD,Nucl. Phys. B 215 (1983) 433
1983
-
[29]
M.LuscherandU.Wolff, HowtoCalculatetheElasticScatteringMatrixinTwo-dimensional Quantum Field Theories by Numerical Simulation, Nucl. Phys. B339 (1990) 222
1990
-
[30]
Jay and E.T
W.I. Jay and E.T. Neil,Bayesian model averaging for analysis of lattice field theory results, Phys. Rev. D103 (2021) 114502 [2008.01069]
2021 arXiv
-
[31]
Luscher,Signatures of unstable particles in finite volume,Nucl
M. Luscher,Signatures of unstable particles in finite volume,Nucl. Phys. B364 (1991) 237
1991
-
[32]
Rummukainen and S.A
K. Rummukainen and S.A. Gottlieb,Resonance scattering phase shifts on a nonrest frame lattice,Nucl. Phys. B450(1995) 397 [hep-lat/9503028]
1995 arXiv
-
[33]
Kim, C.T
C.h. Kim, C.T. Sachrajda and S.R. Sharpe,Finite-volume effects for two-hadron states in moving frames, Nucl. Phys. B727 (2005) 218 [hep-lat/0507006]. 10 Results on meson-meson scattering at large𝑁c Jorge Baeza-Ballesteros
2005 arXiv
-
[34]
S. He, X. Feng and C. Liu,Two particle states and the S-matrix elements in multi-channel scattering,JHEP 07(2005) 011 [hep-lat/0504019]
2005 arXiv
-
[35]
Bernard, M
V. Bernard, M. Lage, U.-G. Meissner and A. Rusetsky,Resonance properties from the finite-volume energy spectrum,JHEP08 (2008) 024 [0806.4495]
2008 arXiv
-
[36]
Luu and M.J
T. Luu and M.J. Savage,Extracting Scattering Phase-Shifts in Higher Partial-Waves from Lattice QCD Calculations, Phys. Rev. D83 (2011) 114508 [1101.3347]
2011 arXiv
-
[37]
Briceno and Z
R.A. Briceno and Z. Davoudi,Moving multichannel systems in a finite volume with application to proton-proton fusion, Phys. Rev. D88(2013) 094507 [1204.1110]
2013 arXiv
-
[38]
Briceno,Two-particle multichannel systems in a finite volume with arbitrary spin,Phys
R.A. Briceno,Two-particle multichannel systems in a finite volume with arbitrary spin,Phys. Rev. D89(2014) 074507 [1401.3312]
2014 arXiv
-
[39]
Gockeler, R
M. Gockeler, R. Horsley, M. Lage, U.G. Meissner, P.E.L. Rakow, A. Rusetsky et al., Scattering phases for meson and baryon resonances on general moving-frame lattices, Phys. Rev. D86(2012) 094513 [1206.4141]
2012 arXiv
-
[40]
Yndurain,Low-energy pion physics, hep-ph/0212282
F.J. Yndurain,Low-energy pion physics, hep-ph/0212282
-
[41]
Pelaez, A
J.R. Pelaez, A. Rodas and J. Ruiz De Elvira,Global parameterization of𝜋𝜋 scattering up to 2 GeV, Eur. Phys. J. C79(2019) 1008 [1907.13162]. 11
2019 arXiv
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.