pith. sign in

arxiv: 2605.16977 · v1 · pith:Y2GDXQ6Cnew · submitted 2026-05-16 · ✦ hep-lat · hep-ph· nucl-th

Two-nucleon systems at m_(π)approx292 MeV from lattice QCD

Pith reviewed 2026-05-19 18:54 UTC · model grok-4.3

classification ✦ hep-lat hep-phnucl-th
keywords lattice QCDnucleon-nucleon scatteringvirtual statedeuterondi-neutronfinite volumeLuscher methodpion mass
0
0 comments X

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.

This paper computes finite-volume energies of two-nucleon systems in lattice QCD at a pion mass of roughly 292 MeV using three different spatial volumes. It converts those energies into scattering amplitudes in the ^3S1 and ^1S0 channels with Lüscher's method and locates a virtual-state pole in each channel. The poles lie on the unphysical sheet with binding energies of 6^{+5}_{-3} MeV and 11^{+6}_{-5} MeV. A reader cares because the result supplies a concrete datum on how nucleon-nucleon binding changes with light-quark mass, which is required to understand the physical deuteron and to guide extrapolations toward the real world.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2605.16977 by Chuan Liu, Jia-Jun Wu, Kang Yu, Kuan Zhang, Liuming Liu, Peng Sun, Ruilin Zhu, Yiqi Geng.

Figure 1
Figure 1. Figure 1: FIG. 1. The energy levels in the center-of-momentum frame [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Ω( [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. log [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The results for the parameters [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The distributions of binding energy in the [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Finite volume energies from the FVH method com [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Fitting details on for the ensemble C24P29, in the [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Details of the fitting for the ensemble C32P29, in the [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Details of the fitting process for the ensemble C48P29, in the [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [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.
  2. [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)
  1. [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.
  2. [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

2 responses · 1 unresolved

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
  1. 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

  2. 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

standing simulated objections not resolved
  • 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

0 steps flagged

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

0 free parameters · 2 axioms · 0 invented entities

The analysis rests on the validity of Lüscher's quantization condition for two-nucleon systems, the assumption that the chosen interpolating operators have good overlap with the ground states, and the neglect of higher partial waves at the energies considered. No new free parameters are introduced beyond the standard lattice QCD action parameters already fixed by the ensembles.

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.
    Invoked in the scattering analysis section when converting finite-volume energies to phase shifts.
  • domain assumption The distillation smearing and symmetric correlators sufficiently suppress excited-state contamination and ensure reliable ground-state energies.
    Stated in the methods when describing the operator construction.

pith-pipeline@v0.9.0 · 5798 in / 1519 out tokens · 30476 ms · 2026-05-19T18:54:54.445685+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

57 extracted references · 57 canonical work pages · 32 internal anchors

  1. [1]

    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

    + (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...

  2. [2]

    DP210103706 (DBL) and DP230101791 (AWT)

    It was also supported by the Australian Research Council through Grant Nos. DP210103706 (DBL) and DP230101791 (AWT)

  3. [3]

    Ab initio calculation of the neutron-proton mass difference

    S. Borsanyiet al.(BMW), Science347, 1452 (2015), arXiv:1406.4088 [hep-lat]

  4. [4]

    C. C. Changet al., Nature558, 91 (2018), arXiv:1805.12130 [hep-lat]

  5. [5]

    Luscher, Commun

    M. Luscher, Commun. Math. Phys.105, 153 (1986)

  6. [6]

    Luscher, Nucl

    M. Luscher, Nucl. Phys. B354, 531 (1991)

  7. [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]

  8. [8]

    S. R. Beaneet al.(NPLQCD), Phys. Rev. C88, 024003 (2013), arXiv:1301.5790 [hep-lat]

  9. [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]

  10. [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]

  11. [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]

  12. [12]

    Illaet al.(NPLQCD), Phys

    M. Illaet al.(NPLQCD), Phys. Rev. D103, 054508 (2021), arXiv:2009.12357 [hep-lat]

  13. [13]

    Orginos, A

    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. [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]

  15. [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]

  16. [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]

  17. [17]

    H¨ orzet al., Phys

    B. H¨ orzet al., Phys. Rev. C103, 014003 (2021), arXiv:2009.11825 [hep-lat]

  18. [18]

    Amarasinghe, R

    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. [19]

    Detmold, M

    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. [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]

  21. [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]

  22. [22]

    Are two nucleons bound in lattice QCD for heavy quark masses? -- Consistency check with L\"uscher's finite volume formula --

    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]

  23. [23]

    Aoki and T

    S. Aoki and T. Doi, Front. in Phys.8, 307 (2020), arXiv:2003.10730 [hep-lat]

  24. [24]

    Bulavaet al.(BaSc), (2025), arXiv:2505.05547 [hep- lat]

    J. Bulavaet al.(BaSc), (2025), arXiv:2505.05547 [hep- lat]

  25. [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

  26. [26]

    S. R. Beane, P. F. Bedaque, A. Parreno, and M. J. Savage, Nucl. Phys. A747, 55 (2005), arXiv:nucl- 11 th/0311027

  27. [27]

    S. R. Beane, P. F. Bedaque, A. Parreno, and M. J. Savage, Phys. Lett. B585, 106 (2004), arXiv:hep- lat/0312004

  28. [28]

    Two Particle States in an Asymmetric Box

    X. Li and C. Liu, Phys. Lett. B587, 100 (2004), arXiv:hep-lat/0311035

  29. [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

  30. [30]

    X. Feng, X. Li, and C. Liu, Phys. Rev. D70, 014505 (2004), arXiv:hep-lat/0404001

  31. [31]

    N. H. Christ, C. Kim, and T. Yamazaki, Phys. Rev. D 72, 114506 (2005), arXiv:hep-lat/0507009

  32. [32]

    C. h. Kim, C. T. Sachrajda, and S. R. Sharpe, Nucl. Phys. B727, 218 (2005), arXiv:hep-lat/0507006

  33. [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]

  34. [34]

    S. Bour, S. Koenig, D. Lee, H. W. Hammer, and U.-G. Meissner, Phys. Rev. D84, 091503 (2011), arXiv:1107.1272 [nucl-th]

  35. [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]

  36. [36]

    Scattering phase shifts for two particles of different mass and non-zero total momentum in lattice QCD

    L. Leskovec and S. Prelovsek, Phys. Rev. D85, 114507 (2012), arXiv:1202.2145 [hep-lat]

  37. [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]

  38. [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]

  39. [39]

    R. A. Briceno, Z. Davoudi, and T. C. Luu, Phys. Rev. D88, 034502 (2013), arXiv:1305.4903 [hep-lat]

  40. [40]

    S. He, X. Feng, and C. Liu, JHEP07, 011 (2005), arXiv:hep-lat/0504019

  41. [41]

    R. A. Briceno, Phys. Rev. D89, 074507 (2014), arXiv:1401.3312 [hep-lat]

  42. [42]

    Huet al.(CLQCD), Phys

    Z.-C. Huet al.(CLQCD), Phys. Rev. D109, 054507 (2024), arXiv:2310.00814 [hep-lat]

  43. [43]

    Duet al.(CLQCD), Phys

    H.-Y. Duet al.(CLQCD), Phys. Rev. D111, 054504 (2025), arXiv:2408.03548 [hep-lat]

  44. [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]

  45. [45]

    H. Yan, C. Liu, L. Liu, and Y. Meng, (2025), arXiv:2507.16070 [hep-lat]

  46. [46]

    Luscher and U

    M. Luscher and U. Wolff, Nucl. Phys. B339, 222 (1990)

  47. [47]

    Estimating the two-particle $K$-matrix for multiple partial waves and decay channels from finite-volume energies

    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]

  48. [48]

    Padmanath and S

    M. Padmanath and S. Prelovsek, Phys. Rev. Lett.129, 032002 (2022), arXiv:2202.10110 [hep-lat]

  49. [49]

    Prelovsek, S

    S. Prelovsek, S. Collins, D. Mohler, M. Padmanath, and S. Piemonte, JHEP06, 035 (2021), arXiv:2011.02542 [hep-lat]

  50. [50]

    Piemonte, S

    S. Piemonte, S. Collins, D. Mohler, M. Padmanath, and S. Prelovsek, Phys. Rev. D100, 074505 (2019), arXiv:1905.03506 [hep-lat]

  51. [51]

    Yu, G.-J

    K. Yu, G.-J. Wang, J.-J. Wu, and Z. Yang, JHEP04, 108 (2025), arXiv:2502.05789 [hep-lat]

  52. [52]

    L. Meng, E. Ortiz-Pacheco, V. Baru, E. Epelbaum, M. Padmanath, and S. Prelovsek, Phys. Rev. D111, 034509 (2025), arXiv:2411.06266 [hep-lat]

  53. [53]

    L. Meng, V. Baru, E. Epelbaum, A. A. Filin, and A. M. Gasparyan, Phys. Rev. D109, L071506 (2024), arXiv:2312.01930 [hep-lat]

  54. [54]

    Li, J.-j

    Y. Li, J.-j. Wu, D. B. Leinweber, and A. W. Thomas, Phys. Rev. D103, 094518 (2021), arXiv:2103.12260 [hep- lat]

  55. [55]

    Chiral effective field theory and nuclear forces

    R. Machleidt and D. R. Entem, Phys. Rept.503, 1 (2011), arXiv:1105.2919 [nucl-th]

  56. [56]

    Machleidt, K

    R. Machleidt, K. Holinde, and C. Elster, Phys. Rept. 149, 1 (1987)

  57. [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...