Two-nucleon systems at m_(π)approx292 MeV from lattice QCD
Pith reviewed 2026-05-19 18:54 UTC · model grok-4.3
The pith
Lattice QCD at 292 MeV pion mass finds virtual state poles in deuteron and di-neutron channels
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Finite-volume energies are computed for nucleon-nucleon systems in the ^3S1 and ^1S0 channels on three volumes at m_π ≈ 292 MeV. Lüscher's method is used to extract the scattering amplitudes, revealing a virtual state pole in each channel with binding energies 6^{+5}_{-3} MeV for the deuteron channel and 11^{+6}_{-5} MeV for the di-neutron channel. An alternative Non-Perturbative Hamiltonian framework yields consistent results, indicating that left-hand cut effects do not alter the conclusion.
What carries the argument
Lüscher's finite-volume quantization condition that connects the discrete spectrum in a finite spatial volume to the infinite-volume scattering phase shift or amplitude.
If this is right
- The ^3S1 channel has a virtual state rather than a bound deuteron at this pion mass.
- The ^1S0 channel has a virtual state for the di-neutron.
- The position of these poles can be used to constrain effective field theories at this quark mass.
- Consistent results from two different analysis methods strengthen the reliability of the virtual state identification.
Where Pith is reading between the lines
- If the poles move to the physical sheet at lighter pion masses, this would explain the emergence of the physical deuteron bound state.
- Similar calculations at other pion masses could map the trajectory of the pole as a function of m_π.
- These results may help test whether the di-neutron becomes bound or remains virtual in the chiral limit or at physical masses.
Load-bearing premise
Higher partial waves and residual finite-volume distortions do not affect the extracted low-energy scattering parameters at the volumes and momenta studied.
What would settle it
A direct computation of the infinite-volume scattering length or effective range in these channels at the same pion mass, or a calculation at the physical pion mass to see if the poles become bound states.
Figures
read the original abstract
Nucleon-nucleon systems in the $^3S_1$ and the $^1S_0$ channels are studied in lattice quantum chromodynamics at a pion mass of approximately $m_{\pi}\approx292$ MeV, employing three $N_f = 2+1$ ensembles with the same pion mass and lattice spacing $a=0.10530(18)$ fm but different spatial volumes. Finite-volume energies of the nucleon-nucleon systems are determined in both the rest frame and a moving frame. The distillation quark smearing method is applied to improve the precision and to ensure the symmetric correlators by using the same interpolating operators at sink and source. The scattering amplitudes are extracted from the finite-volume spectra using the L\"uscher's finite-volume method. At the studied pion mass, both the $^3S_1$ (deuteron) and $^1S_0$(di-neutron) channels exhibit a virtual state pole, with binding energies of $6^{+5}_{-3}$ MeV and $11^{+6}_{-5}$ MeV, respectively. To investigate the effects of the left-hand cut, an alternative method -- the Non-Perturbative Hamiltonian framework (NPHF) -- is used for the scattering analysis and yields consistent results with those from the L\"uscher method.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports lattice QCD results for two-nucleon systems in the ^3S1 and ^1S0 channels at m_π ≈ 292 MeV on three Nf=2+1 ensembles with fixed lattice spacing but varying spatial volumes. Finite-volume energies are extracted in both rest and moving frames using distillation-improved operators. Scattering amplitudes and infinite-volume poles are then obtained via Lüscher's quantization condition, with a cross-check using the Non-Perturbative Hamiltonian framework (NPHF). The central result is the identification of virtual-state poles in both channels, with binding energies 6^{+5}_{-3} MeV (^3S1) and 11^{+6}_{-5} MeV (^1S0).
Significance. If the quoted pole positions hold after systematic checks, the work supplies useful lattice data on the pion-mass dependence of NN interactions near the physical point, aiding chiral EFT calibrations and extrapolations. The agreement between independent Lüscher and NPHF extractions, together with the use of multiple volumes and frames, strengthens the analysis and provides a concrete benchmark for future studies at lighter pion masses.
major comments (2)
- [Lüscher analysis (methods and results sections)] Lüscher analysis (methods and results sections): The extraction of the virtual-state poles assumes S-wave dominance with negligible higher partial-wave mixing, particularly in the moving-frame data, and that residual finite-volume corrections beyond the leading Lüscher term lie inside the quoted uncertainties. No explicit tests—such as volume-by-volume consistency of extracted phase shifts, fits that float P-wave parameters, or direct per-volume comparison of Lüscher versus NPHF amplitudes—are presented. Because the reported poles lie close to threshold, even modest contamination would shift the binding energies outside the stated errors.
- [Finite-volume spectra extraction] Finite-volume spectra extraction: The manuscript states that higher partial-wave contamination and discretization effects are under control at the quoted precision, yet no quantitative assessment (e.g., comparison of spectra across the three volumes after subtracting leading Lüscher contributions or inclusion of lattice-spacing artifacts in the quantization condition) is shown. This assumption is load-bearing for the central claim that both channels exhibit virtual states at the reported locations.
minor comments (2)
- [Results] Notation: The definition of the binding energy relative to the two-nucleon threshold should be stated explicitly in the text when the pole positions are first quoted, to avoid ambiguity with the sign convention for virtual states.
- [Figures] Figures: The error bands on the phase-shift plots would benefit from an additional panel or inset showing the sensitivity to the choice of fitting range or to the inclusion of a small P-wave term.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review of our manuscript. The comments highlight important aspects of the analysis that warrant additional clarification and explicit checks. We address each major comment below and will revise the manuscript to incorporate the suggested improvements where feasible.
read point-by-point responses
-
Referee: [Lüscher analysis (methods and results sections)] Lüscher analysis (methods and results sections): The extraction of the virtual-state poles assumes S-wave dominance with negligible higher partial-wave mixing, particularly in the moving-frame data, and that residual finite-volume corrections beyond the leading Lüscher term lie inside the quoted uncertainties. No explicit tests—such as volume-by-volume consistency of extracted phase shifts, fits that float P-wave parameters, or direct per-volume comparison of Lüscher versus NPHF amplitudes—are presented. Because the reported poles lie close to threshold, even modest contamination would shift the binding energies outside the stated errors.
Authors: We agree that more explicit tests would strengthen the presentation. The manuscript already demonstrates consistency of the extracted poles across three volumes and both rest and moving frames, together with agreement between the independent Lüscher and NPHF determinations. To address the referee's concern directly, the revised version will include: (i) volume-by-volume phase-shift extractions showing stability, (ii) an extended fit in which P-wave parameters are floated (yielding values consistent with zero within uncertainties), and (iii) a per-volume comparison of the scattering amplitudes obtained from the two methods. These additions confirm that higher partial-wave mixing remains negligible at the quoted precision and that residual finite-volume effects are captured within the reported uncertainties. revision: yes
-
Referee: [Finite-volume spectra extraction] Finite-volume spectra extraction: The manuscript states that higher partial-wave contamination and discretization effects are under control at the quoted precision, yet no quantitative assessment (e.g., comparison of spectra across the three volumes after subtracting leading Lüscher contributions or inclusion of lattice-spacing artifacts in the quantization condition) is shown. This assumption is load-bearing for the central claim that both channels exhibit virtual states at the reported locations.
Authors: We have performed internal checks that compare the finite-volume spectra across the three volumes after subtracting the leading Lüscher contributions; these show good consistency within statistical uncertainties. For discretization effects, all ensembles share the same lattice spacing, precluding a direct multi-spacing comparison. In the revision we will add a quantitative discussion estimating the expected size of residual discretization and higher-wave effects based on the observed volume dependence and on existing literature at comparable pion masses and lattice spacings. This discussion will be placed in the methods and results sections to make the control of systematics explicit. revision: partial
- A fully quantitative inclusion of lattice-spacing artifacts directly into the quantization condition or a multi-spacing comparison of spectra cannot be performed with the present set of ensembles at fixed lattice spacing; such an assessment would require additional simulations at different a that lie beyond the scope of the current work.
Circularity Check
No significant circularity; derivation is self-contained
full rationale
The paper computes finite-volume two-nucleon energies on three lattice ensembles using the distillation method, then maps those energies to infinite-volume scattering amplitudes via Lüscher's quantization condition (and the independent NPHF cross-check) to locate virtual-state poles. This mapping is a standard, externally defined relation between discrete FV spectra and continuum phase shifts; the reported binding energies (6^{+5}_{-3} MeV and 11^{+6}_{-5} MeV) are derived outputs, not inputs that are fitted or redefined by construction. No step invokes a self-citation whose validity depends on the present result, renames a known pattern as a new derivation, or smuggles an ansatz that makes the final pole positions tautological with the raw spectra. The analysis therefore remains non-circular.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Lüscher's finite-volume quantization condition accurately relates the discrete spectrum to the infinite-volume scattering amplitude for S-wave NN systems at these volumes and energies.
- domain assumption The distillation smearing and symmetric correlators sufficiently suppress excited-state contamination and ensure reliable ground-state energies.
Reference graph
Works this paper leans on
-
[1]
+ (m2 1 −m 2 2)2 4E2cm .(11) InN Nscattering, the left-hand cut corresponds tok 2 = − 1 4 m2 π. Since the L¨ uscher formula is not valid below this cut, we restrict our scattering analysis using the L¨ uscher formula to energy levels above this region. In Sec. V, we employ an alternative method, NPHF, for the scattering analysis. In this method, the effec...
work page 2008
-
[2]
DP210103706 (DBL) and DP230101791 (AWT)
It was also supported by the Australian Research Council through Grant Nos. DP210103706 (DBL) and DP230101791 (AWT)
-
[3]
Ab initio calculation of the neutron-proton mass difference
S. Borsanyiet al.(BMW), Science347, 1452 (2015), arXiv:1406.4088 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[4]
C. C. Changet al., Nature558, 91 (2018), arXiv:1805.12130 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [5]
- [6]
-
[7]
S. R. Beane, E. Chang, S. D. Cohen, W. Detmold, H. W. Lin, T. C. Luu, K. Orginos, A. Parreno, M. J. Sav- age, and A. Walker-Loud (NPLQCD), Phys. Rev. D87, 034506 (2013), arXiv:1206.5219 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[8]
S. R. Beaneet al.(NPLQCD), Phys. Rev. C88, 024003 (2013), arXiv:1301.5790 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[9]
Two-Nucleon Higher Partial-Wave Scattering from Lattice QCD
E. Berkowitz, T. Kurth, A. Nicholson, B. Joo, E. Rinaldi, M. Strother, P. M. Vranas, and A. Walker-Loud, Phys. Lett. B765, 285 (2017), arXiv:1508.00886 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[10]
M. L. Wagman, F. Winter, E. Chang, Z. Davoudi, W. Detmold, K. Orginos, M. J. Savage, and P. E. Shana- han, Phys. Rev. D96, 114510 (2017), arXiv:1706.06550 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[11]
S. R. Beane, E. Chang, W. Detmold, H. W. Lin, T. C. Luu, K. Orginos, A. Parreno, M. J. Savage, A. Torok, and A. Walker-Loud (NPLQCD), Phys. Rev. D85, 054511 (2012), arXiv:1109.2889 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[12]
M. Illaet al.(NPLQCD), Phys. Rev. D103, 054508 (2021), arXiv:2009.12357 [hep-lat]
-
[13]
K. Orginos, A. Parreno, M. J. Savage, S. R. Beane, E. Chang, and W. Detmold, Phys. Rev. D92, 114512 (2015), [Erratum: Phys.Rev.D 102, 039903 (2020)], arXiv:1508.07583 [hep-lat]
-
[14]
Helium nuclei, deuteron and dineutron in 2+1 flavor lattice QCD
T. Yamazaki, K.-i. Ishikawa, Y. Kuramashi, and A. Ukawa, Phys. Rev. D86, 074514 (2012), arXiv:1207.4277 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[15]
Study of quark mass dependence of binding energy for light nuclei in 2+1 flavor lattice QCD
T. Yamazaki, K.-i. Ishikawa, Y. Kuramashi, and A. Ukawa, Phys. Rev. D92, 014501 (2015), arXiv:1502.04182 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[16]
Lattice QCD study of the $H$ dibaryon using hexaquark and two-baryon interpolators
A. Francis, J. R. Green, P. M. Junnarkar, C. Miao, T. D. Rae, and H. Wittig, Phys. Rev. D99, 074505 (2019), arXiv:1805.03966 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[17]
B. H¨ orzet al., Phys. Rev. C103, 014003 (2021), arXiv:2009.11825 [hep-lat]
-
[18]
S. Amarasinghe, R. Baghdadi, Z. Davoudi, W. Detmold, M. Illa, A. Parreno, A. V. Pochinsky, P. E. Shana- han, and M. L. Wagman, Phys. Rev. D107, 094508 (2023), [Erratum: Phys.Rev.D 110, 119904 (2024)], arXiv:2108.10835 [hep-lat]
-
[19]
W. Detmold, M. Illa, W. I. Jay, A. Parre˜ no, R. J. Perry, P. E. Shanahan, and M. L. Wagman (NPLQCD), Phys. Rev. D111, 114501 (2025), arXiv:2404.12039 [hep-lat]
-
[20]
Hadron-Hadron Interactions from Imaginary-time Nambu-Bethe-Salpeter Wave Function on the Lattice
N. Ishii, S. Aoki, T. Doi, T. Hatsuda, Y. Ikeda, T. Inoue, K. Murano, H. Nemura, and K. Sasaki (HAL QCD), Phys. Lett. B712, 437 (2012), arXiv:1203.3642 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[21]
Two-Baryon Potentials and H-Dibaryon from 3-flavor Lattice QCD Simulations
T. Inoue, S. Aoki, T. Doi, T. Hatsuda, Y. Ikeda, N. Ishii, K. Murano, H. Nemura, and K. Sasaki (HAL QCD), Nucl. Phys. A881, 28 (2012), arXiv:1112.5926 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[22]
T. Iritani, S. Aoki, T. Doi, T. Hatsuda, Y. Ikeda, T. In- oue, N. Ishii, H. Nemura, and K. Sasaki, Phys. Rev. D 96, 034521 (2017), arXiv:1703.07210 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[23]
S. Aoki and T. Doi, Front. in Phys.8, 307 (2020), arXiv:2003.10730 [hep-lat]
-
[24]
Bulavaet al.(BaSc), (2025), arXiv:2505.05547 [hep- lat]
J. Bulavaet al.(BaSc), (2025), arXiv:2505.05547 [hep- lat]
-
[25]
Resonance Scattering Phase Shifts on a Non-Rest Frame Lattice
K. Rummukainen and S. A. Gottlieb, Nucl. Phys. B450, 397 (1995), arXiv:hep-lat/9503028
work page internal anchor Pith review Pith/arXiv arXiv 1995
- [26]
- [27]
-
[28]
Two Particle States in an Asymmetric Box
X. Li and C. Liu, Phys. Lett. B587, 100 (2004), arXiv:hep-lat/0311035
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[29]
Electroweak Matrix Elements in the Two-Nucleon Sector from Lattice QCD
W. Detmold and M. J. Savage, Nucl. Phys. A743, 170 (2004), arXiv:hep-lat/0403005
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[30]
X. Feng, X. Li, and C. Liu, Phys. Rev. D70, 014505 (2004), arXiv:hep-lat/0404001
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[31]
N. H. Christ, C. Kim, and T. Yamazaki, Phys. Rev. D 72, 114506 (2005), arXiv:hep-lat/0507009
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[32]
C. h. Kim, C. T. Sachrajda, and S. R. Sharpe, Nucl. Phys. B727, 218 (2005), arXiv:hep-lat/0507006
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[33]
Resonance properties from the finite-volume energy spectrum
V. Bernard, M. Lage, U.-G. Meissner, and A. Rusetsky, JHEP08, 024 (2008), arXiv:0806.4495 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[34]
S. Bour, S. Koenig, D. Lee, H. W. Hammer, and U.-G. Meissner, Phys. Rev. D84, 091503 (2011), arXiv:1107.1272 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[35]
Improving the Volume Dependence of Two-Body Binding Energies Calculated with Lattice QCD
Z. Davoudi and M. J. Savage, Phys. Rev. D84, 114502 (2011), arXiv:1108.5371 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[36]
L. Leskovec and S. Prelovsek, Phys. Rev. D85, 114507 (2012), arXiv:1202.2145 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[37]
Scattering phases for meson and baryon resonances on general moving-frame lattices
M. Gockeler, R. Horsley, M. Lage, U. G. Meiss- ner, P. E. L. Rakow, A. Rusetsky, G. Schierholz, and J. M. Zanotti, Phys. Rev. D86, 094513 (2012), arXiv:1206.4141 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[38]
Derivation of L\"uscher's finite size formula for $N\pi$ and $NN$ system
N. Ishizuka, PoSLA T2009, 119 (2009), arXiv:0910.2772 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[39]
R. A. Briceno, Z. Davoudi, and T. C. Luu, Phys. Rev. D88, 034502 (2013), arXiv:1305.4903 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[40]
S. He, X. Feng, and C. Liu, JHEP07, 011 (2005), arXiv:hep-lat/0504019
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[41]
R. A. Briceno, Phys. Rev. D89, 074507 (2014), arXiv:1401.3312 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[42]
Z.-C. Huet al.(CLQCD), Phys. Rev. D109, 054507 (2024), arXiv:2310.00814 [hep-lat]
-
[43]
H.-Y. Duet al.(CLQCD), Phys. Rev. D111, 054504 (2025), arXiv:2408.03548 [hep-lat]
-
[44]
A novel quark-field creation operator construction for hadronic physics in lattice QCD
M. Peardon, J. Bulava, J. Foley, C. Morningstar, J. Dudek, R. G. Edwards, B. Joo, H.-W. Lin, D. G. Richards, and K. J. Juge (Hadron Spectrum), Phys. Rev. D80, 054506 (2009), arXiv:0905.2160 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2009
- [45]
- [46]
-
[47]
C. Morningstar, J. Bulava, B. Singha, R. Brett, J. Fallica, A. Hanlon, and B. H¨ orz, Nucl. Phys. B924, 477 (2017), arXiv:1707.05817 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[48]
M. Padmanath and S. Prelovsek, Phys. Rev. Lett.129, 032002 (2022), arXiv:2202.10110 [hep-lat]
-
[49]
S. Prelovsek, S. Collins, D. Mohler, M. Padmanath, and S. Piemonte, JHEP06, 035 (2021), arXiv:2011.02542 [hep-lat]
-
[50]
S. Piemonte, S. Collins, D. Mohler, M. Padmanath, and S. Prelovsek, Phys. Rev. D100, 074505 (2019), arXiv:1905.03506 [hep-lat]
- [51]
- [52]
- [53]
- [54]
-
[55]
Chiral effective field theory and nuclear forces
R. Machleidt and D. R. Entem, Phys. Rept.503, 1 (2011), arXiv:1105.2919 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2011
- [56]
-
[57]
J.-J. Wu, T. S. H. Lee, A. W. Thomas, and R. D. Young, Phys. Rev. C90, 055206 (2014), arXiv:1402.4868 [hep- lat]. 12 APPENDIX The total energyEin the lab frame for a specific two-particle interacting state with total momentumP= 2π L dis determined from the time dependence ofλ n(t) in Eq.7. The fitting process is detailed in Figs. 8, 9, and 10. In our anal...
work page internal anchor Pith review Pith/arXiv arXiv 2014
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.