REVIEW 3 major objections 5 minor 51 references
Lattice QCD study of $\Lambda_c \Lambda_c$ scattering
T0 review · 3 major / 5 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read This paper reports the first lattice QCD calculation of $\Lambda_c\Lambda_c$ scattering in the $I(J^P)=0(0^+)$ channel and finds the interaction is repulsive, with scattering length $a_0 = -0.21(4)(8)$ fm and no bound state in the studied…
desk verdict First lattice QCD study of Lambda_c Lambda_c scattering: a competent calculation with an honest error budget, but the repulsive conclusion leans on a qualitative no-mixing check rather than a quantitative coupled-channel bound. 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 method is L\"uscher's finite-volume formalism: the relation $k\cot\delta(k) = \frac{2\sqrt{\pi}}{L} Z_{00}(1;q^2)$ converts each finite-volume energy into the infinite-volume s-wave phase shift $\delta(k)$, with $Z_{00}$ the generalized zeta function encoding the box geometry. The phase shift is then parameterized by the effective range expansion $k\cot\delta = \frac{1}{a_0} + \frac{1}{2} r_0 k^2$. Energy shifts are extracted from ratios of the two-baryon correlation function to single-baryon correlators, and the free energies are evaluated with the continuum dispersion relation, a choice that reduces the impact of the observed roughly 4% deviation of the lattice speed-of-light parameter $c$ from unity.
What would settle it
A higher-statistics GEVP with additional $\Xi_{cc}N$ operators that reveals level repulsion between $\Lambda_c\Lambda_c$ and $\Xi_{cc}N$ states would refute the no-mixing assumption; alternatively, a second lattice spacing at the same pion mass that moves $a_0$ to a positive value would indicate that the repulsion is a discretization artifact.
Extended reading notes
Core claim
The central claim is that the $\Lambda_c\Lambda_c$ interaction in the $I(J^P)=0(0^+)$ channel is repulsive near threshold. From seven finite-volume energy levels across two lattice volumes, the s-wave scattering length and effective range are determined as $a_0 = -0.21(4)(8)$ fm and $r_0 = -0.05(13)(25)$ fm, where the first error is statistical and the second systematic; no pole appears in the amplitude over the fitted energy range. The $\Xi_{cc}N$ channel is omitted because a GEVP including both $\Lambda_c\Lambda_c$ and $\Xi_{cc}N$ operators shows no mixing between the levels, and the $\Sigma_c\Sigma_c$ channel is omitted because the analyzed energies lie below its threshold. This is presented as the first lattice QCD constraint on double-charm baryon scattering.
Load-bearing premise
The extraction assumes the measured finite-volume levels are purely $\Lambda_c\Lambda_c$ states and are not secretly mixed with the $\Xi_{cc}N$ channel; if a weak coupling is hiding below the statistical noise, the quoted scattering length would be biased and the repulsive conclusion could change.
Editorial extensions
If this is right
- At $m_\pi \sim 303$ MeV, the single-channel $\Lambda_c\Lambda_c$ system in this quantum number is unbound, so any double-charm dibaryon in this channel would have to arise from coupled-channel effects rather than from the $\Lambda_c\Lambda_c$ attraction alone.
- The quoted $a_0$ and $r_0$ give the first lattice-QCD numbers for double-charm baryon scattering, giving model calculations a concrete target to reproduce.
- The observed absence of mixing between $\Lambda_c\Lambda_c$ and $\Xi_{cc}N$ operators supports the single-channel analysis only within current statistical precision; a stronger coupled-channel study would be needed near the $\Xi_{cc}N$ threshold.
- Because the calculation uses one lattice spacing and one pion mass, these parameters are not yet the physical-point values; the paper quotes systematic errors from the effective-range expansion and fit-range choices but not a continuum-extrapolation error.
Reading between the lines
- If the repulsion persists toward the physical pion mass, the most plausible route to a double-charm dibaryon in this channel runs through the $\Sigma_c\Sigma_c$ channel, and a future calculation with levels above that threshold would test that scenario.
- The roughly 4% deviation of the lattice dispersion coefficient $c$ from 1 suggests charm-quark discretization effects; a second lattice spacing could turn the quoted one-point result into a continuum-extrapolated value.
- The same correlation functions could be re-analyzed with a coupled-channel L\"uscher treatment or an alternative potential method; agreement between independent extractions would strengthen confidence in the repulsion, while disagreement would expose a systematic not covered by the quoted errors.
- The no-mixing observation is consistent with a genuinely weak coupling or with operators that do not overlap the mixed state; adding higher-momentum $\Xi_{cc}N$ operators and more statistics would discriminate between these readings.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents a lattice QCD study of ΛcΛc scattering in the I(J^P)=0(0+) channel using two N_f=2+1 Wilson-Clover ensembles with pion mass mπ≈303 MeV and lattice spacing a=0.07746 fm (L/a=32 and 48). The authors compute the finite-volume spectrum from GEVP and diagonal correlation functions, extract energy shifts from ratio fits, and apply Lüscher's finite-volume method with an effective-range expansion to obtain a0=-0.21(4)(8) fm and r0=-0.05(13)(25) fm. They interpret the negative scattering length as evidence for a repulsive interaction and thus the absence of a ΛcΛc bound state in this channel. The coupled channels ΞccN and ΣcΣc are argued to be negligible based on the observed absence of operator mixing and the energy window below the ΣcΣc threshold, respectively.
Significance. If the result holds, this is the first lattice QCD constraint on double-charm baryon scattering and provides a useful benchmark for phenomenological models, which currently disagree on whether ΛcΛc is bound. The paper benefits from a transparent analysis: the GEVP fits, the ratio method for energy shifts, the covariance-matrix fits, and the systematic-error estimates from fit-range variations and the ERE are all described in detail, and the absence of input-output circularity is clear. The main sources of concern are the unquantified effect of the open ΞccN channel and the mixing of continuum and lattice dispersion conventions in Eq. (12); both bear directly on the central quantitative claim. The single-lattice-spacing and single-pion-mass setup is acknowledged and limits the conclusions to the simulated quark masses.
major comments (3)
- [Sec. III.B and Figs. 2, 5] The argument for neglecting the open ΞccN channel is not quantitative. The observation that the ΛcΛc and ΞccN operators have small overlap with the eigenstates does not imply that the ΛcΛc → ΞccN scattering coupling is negligible, because the energy shift of a finite-volume level depends on the coupled-channel T-matrix and the phase space, not directly on the operator overlaps. Since mΞcc + mN ≈ 4.817 GeV is about 9 MeV below 2mΛc ≈ 4.826 GeV (Table II), ΞccN is open throughout the fitted energy range, and a weak coupling could shift the single-channel levels by more than the statistical errors. Please provide a quantitative bound, for example by performing a two-channel Lüscher fit or by inserting a test coupling and showing that the resulting shifts are within the quoted errors.
- [Sec. IV, Eq. (12) and Table II] The free energies used in Eq. (12) are computed with the continuum dispersion relation, while the measured Λc dispersion relation has c ≈ 0.99 (Table II). Because the ΔEn values in Eq. (11) are extracted from the ratio to the single-baryon correlators at the actual lattice momenta, the combination in Eq. (12) mixes the lattice and continuum dispersion conventions. This mismatch enters the Lüscher analysis and can bias a0 and r0. The expectation that this effect is smaller than the statistical error should be demonstrated numerically, e.g., by repeating the fits using the measured dispersion relation in the quantization condition or by assigning a systematic error from the difference.
- [Sec. IV, Eq. (16) and Table III] The evidence for a repulsive interaction, which is the central quantitative claim, has limited statistical strength when the quoted systematic error is included: a0 = -0.21(4)(8) fm is only about 2.3σ from zero, and the ground-state energy shift on F32P30 is 0.00085(47) in lattice units, i.e., about 1.8σ. Please state the significance of the repulsive sign explicitly, for example by giving the probability that a0 is negative, and temper the abstract and conclusions accordingly if that probability is not high.
minor comments (5)
- [Abstract and Sec. V] The abstract and summary should explicitly state that the repulsive interaction and the scattering parameters are obtained at mπ ≈ 303 MeV, since no chiral extrapolation is performed.
- [Eqs. (7)-(8)] Equations (7) and (8) contain an extra closing parenthesis in the definitions of OΞccN and OΣcΣc.
- [Sec. III.B] The notation 'p2' in the operator labels is ambiguous; using 'n' or '|p|^2' would be clearer.
- [Fig. 6] The axis label of Fig. 6 appears truncated ('0/') and should be completed, e.g., 'δ0 (rad)' or 'δ0/π'.
- [Sec. IV] The statement that the effect of the c≈0.99 dispersion should be 'much smaller than the statistical error' is an estimate; please either provide a numerical estimate or move this caveat to the systematic-error discussion.
Circularity Check
No circularity: the scattering parameters are extracted from independent lattice finite-volume spectra via the external Lüscher relation; no fitted input is fed back into the spectrum.
full rationale
The paper's central result, a0 = -0.21(4)(8) fm and r0 = -0.05(13)(25) fm, is obtained by (i) computing finite-volume two-baryon energies from lattice correlation functions (Sec. III), and (ii) inverting the standard Lüscher quantization condition (Eq. 13) with an effective-range parameterization (Eq. 14) via a chi-square fit to those energies (Eq. 15). Nothing in the input encodes the sign or magnitude of the scattering length: the input is the set of measured energy shifts in Table III, and the output is a fit to those independent lattice data. The use of the continuum dispersion relation in Eq. (12) is an explicit consistency choice, not a fitted parameter recycled as a prediction. The omission of the Xi_ccN and Sigma_cSigma_c channels is justified by observed level shifts and thresholds, and the paper explicitly acknowledges in Sec. V that a quantitative coupled-channel study would require more operators and three-body treatments; that is a systematic-uncertainty concern, not circularity. The CLQCD ensemble citation is data provenance, and the Lüscher and distillation references are external standard methods. No self-citation is load-bearing, and no derivation step reduces by construction to its inputs.
Assumptions & free parameters
free parameters (3)
- a0 (s-wave scattering length) =
-0.21(4)(8) fm
- r0 (effective range) =
-0.05(13)(25) fm
- Energy window cut (aE <= 1.98) =
1.98 (lattice units)
assumptions (6)
- standard math Lüscher's finite volume formula (Eq. 13) with the A1+ irrep and neglect of partial waves l >= 4 connects the finite-volume spectrum to the infinite-volume phase shift.
- domain assumption The CLQCD ensembles provide Nf=2+1 dynamical QCD with the stated bare quark masses, lattice spacing a = 0.07746(18) fm, and pion mass 303 MeV.
- domain assumption The valence charm quark mass is tuned to reproduce the spin-averaged eta_c and J/psi mass, as stated in Sec. II.
- domain assumption The effective range expansion truncated at O(k^2) describes the phase shift in the fitted energy window.
- ad hoc to paper The continuum dispersion relation E^2 = m^2 + p^2 is used to build finite-volume energies in Eq. (12), even though the measured Lambda_c dispersion has c ~ 0.99.
- domain assumption The off-diagonal elements of the Lambda_c Lambda_c correlation matrix are negligible, so the F48P30 spectrum can be extracted from diagonal correlators alone.
Cite this review
Pith. "Pith review of Lattice QCD study of $\Lambda_c \Lambda_c$ scattering." pith.science (2026). https://pith.science/paper/JQK7HFM4
@misc{pith2026250205546,
author = {Pith},
title = {Pith review of: Lattice QCD study of $\Lambda_c \Lambda_c$ scattering},
year = {2026},
howpublished = {\url{https://pith.science/paper/JQK7HFM4}},
note = {Machine review of arXiv:2502.05546}
}
abstract
We present the first lattice result of the near threshold $\Lambda_c\Lambda_c$ scattering with $I(J^P) = 0(0^+)$. The calculation is performed on two $N_f = 2+1$ Wilson-Clover ensembles with pion mass $m_\pi \sim 303$\,MeV and lattice spacing $a = 0.07746$\,fm. The L\"uscher's finite volume method is utilized to extract the scattering parameters from the finite-volume spectrum. The coupled channel $\Xi_{cc}N$ is ignored in the scattering analysis based on the observation that the energy levels computed from the $\Lambda_c\Lambda_c$ and $\Xi_{cc}N$ operators do not mix. The $\Sigma_c\Sigma_c$ channel is not included either since the energy range explored in this study is well below its threshold. Our results indicate that the interaction in the $\Lambda_c\Lambda_c$ single channel is repulsive, and the scattering length is determined to be $a_0 = -0.21(4)(8)$\,fm, where the first error is the statistical error and the second is the systematic error.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[1]
R. L. Jaffe, Phys. Rev. Lett. 38, 195 (1977), [Erratum: Phys.Rev.Lett. 38, 617 (1977)]
1977
-
[2]
S. R. Beane, E. Chang, S. D. Cohen, W. Detmold, H. W. Lin, T. C. Luu, K. Orginos, A. Parreno, M. J. Savage, and A. Walker-Loud (NPLQCD), Phys. Rev. D 87, 034506 (2013), 1206.5219
arXiv 2013
-
[3]
S. R. Beane et al. (NPLQCD), Phys. Rev. C 88, 024003 (2013), 1301.5790
arXiv 2013
-
[4]
E. Berkowitz, T. Kurth, A. Nicholson, B. Joo, E. Rinaldi, M. Strother, P. M. Vranas, and A. Walker-Loud, Phys. Lett. B 765, 285 (2017), 1508.00886
arXiv 2017
-
[5]
K. Orginos, A. Parreno, M. J. Savage, S. R. Beane, E. Chang, and W. Detmold, Phys. Rev. D 92, 114512 (2015), [Erratum: Phys.Rev.D 102, 039903 (2020)], 1508.07583
arXiv 2015
-
[6]
M. L. Wagman, F. Winter, E. Chang, Z. Davoudi, W. Detmold, K. Orginos, M. J. Savage, and P. E. Shanahan, Phys. Rev. D 96, 114510 (2017), 1706.06550
arXiv 2017
- [7]
-
[8]
S. Amarasinghe, R. Baghdadi, Z. Davoudi, W. Detmold, M. Illa, A. Parreno, A. V. Pochinsky, P. E. Shanahan, and M. L. Wagman, Phys. Rev. D 107, 094508 (2023), 2108.10835
arXiv 2023
Show all 51 references
-
[9]
J. R. Green, A. D. Hanlon, P. M. Junnarkar, and H. Wittig, Phys. Rev. Lett. 127, 242003 (2021), 2103.01054
2021 arXiv
-
[10]
Francis, J
A. Francis, J. R. Green, P. M. Junnarkar, C. Miao, T. D. Rae, and H. Wittig, Phys. Rev. D 99, 074505 (2019), 1805.03966
2019 arXiv
- [11]
- [12]
-
[13]
S. R. Beane et al., Mod. Phys. Lett. A 26, 2587 (2011), 1103.2821
2011 arXiv
-
[14]
P. E. Shanahan, A. W. Thomas, and R. D. Young, Phys. Rev. Lett. 107, 092004 (2011), 1106.2851
2011 arXiv
-
[15]
Dong, F.-K
X.-K. Dong, F.-K. Guo, and B.-S. Zou, Commun. Theor. Phys. 73, 125201 (2021), 2108.02673
2021 arXiv
-
[16]
R. Chen, A. Hosaka, and X. Liu, Phys. Rev. D 96, 116012 (2017), 1707.08306
2017 arXiv
-
[17]
Lee, Z.-G
N. Lee, Z.-G. Luo, X.-L. Chen, and S.-L. Zhu, Phys. Rev. D 84, 014031 (2011), 1104.4257
2011 arXiv
- [18]
- [19]
-
[20]
Oka, Nucl
M. Oka, Nucl. Phys. A 914, 447 (2013)
2013
- [21]
- [22]
-
[23]
S. M. Gerasyuta and E. E. Matskevich, Int. J. Mod. Phys. E 21, 1250058 (2012), 1109.2338
2012 arXiv
-
[24]
Lu, L.-S
J.-X. Lu, L.-S. Geng, and M. P. Valderrama, Phys. Rev. D 99, 074026 (2019), 1706.02588
2019 arXiv
-
[25]
Y. Lyu, H. Tong, T. Sugiura, S. Aoki, T. Doi, T. Hatsuda, J. Meng, and T. Miyamoto, Phys. Rev. Lett. 127, 072003 (2021), 2102.00181
2021 arXiv
-
[26]
Mathur, M
N. Mathur, M. Padmanath, and D. Chakraborty, Phys. Rev. Lett. 130, 111901 (2023), 2205.02862
2023 arXiv
-
[27]
P. M. Junnarkar and N. Mathur, Phys. Rev. D 106, 054511 (2022), 2206.02942
2022 arXiv
-
[28]
P. M. Junnarkar and N. Mathur, Phys. Rev. D 111, 014512 (2025), 2410.08519
2025 arXiv
-
[29]
Junnarkar and N
P. Junnarkar and N. Mathur, Phys. Rev. Lett. 123, 162003 (2019), 1906.06054
2019 arXiv
-
[30]
T. F. Carames and A. Valcarce, Phys. Rev. D 92, 034015 (2015), 1507.08278
2015 arXiv
-
[31]
L¨ uscher, Commun
M. L¨ uscher, Commun. Math. Phys.104, 177 (1986)
1986
-
[32]
L¨ uscher, Commun
M. L¨ uscher, Commun. Math. Phys.105, 153 (1986). 10
1986
-
[33]
L¨ uscher, Nucl
M. L¨ uscher, Nucl. Phys. B354, 531 (1991)
1991
- [34]
-
[35]
H. Yan, C. Liu, L. Liu, Y. Meng, and H. Xing (2024), 2404.13479
2024 arXiv
- [36]
-
[37]
H. Liu, L. Liu, P. Sun, W. Sun, J.-X. Tan, W. Wang, Y.-B. Yang, and Q.-A. Zhang, Phys. Lett. B 841, 137941 (2023), 2303.17865
2023 arXiv
-
[38]
H. Xing, J. Liang, L. Liu, P. Sun, and Y.-B. Yang (2022), 2210.08555
2022 arXiv
-
[39]
H. Liu, J. He, L. Liu, P. Sun, W. Wang, Y.-B. Yang, and Q.-A. Zhang, Sci. China Phys. Mech. Astron. 67, 211011 (2024), 2207.00183
2024 arXiv
- [40]
-
[41]
Meng, J.-L
Y. Meng, J.-L. Dang, C. Liu, Z. Liu, T. Shen, H. Yan, and K.-L. Zhang, Phys. Rev. D 109, 074511 (2024), 2401.13475
2024 arXiv
-
[42]
H. Yan, M. Garofalo, M. Mai, U.-G. Meißner, and C. Urbach (2024), 2407.16659
2024 arXiv
-
[43]
Peardon, J
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. D 80, 054506 (2009), 0905.2160
2009 arXiv
-
[44]
J. J. Dudek, R. G. Edwards, M. J. Peardon, D. G. Richards, and C. E. Thomas, Phys. Rev. D82, 034508 (2010), 1004.4930
2010 arXiv
-
[45]
Padmanath and S
M. Padmanath and S. Prelovsek, Phys. Rev. Lett. 129, 032002 (2022), 2202.10110
2022 arXiv
-
[46]
Prelovsek, S
S. Prelovsek, S. Collins, D. Mohler, M. Padmanath, and S. Piemonte, JHEP 06, 035 (2021), 2011.02542
2021 arXiv
-
[47]
Piemonte, S
S. Piemonte, S. Collins, D. Mohler, M. Padmanath, and S. Prelovsek, Phys. Rev. D 100, 074505 (2019), 1905.03506
2019 arXiv
-
[48]
R. G. Edwards and B. Joo (SciDAC, LHPC, UKQCD), Nucl. Phys. B Proc. Suppl. 140, 832 (2005), hep-lat/0409003
2005 arXiv
-
[49]
M. A. Clark, R. Babich, K. Barros, R. C. Brower, and C. Rebbi (QUDA), Comput. Phys. Commun. 181, 1517 (2010), 0911.3191
2010 arXiv
-
[50]
Babich, M
R. Babich, M. A. Clark, B. Joo, G. Shi, R. C. Brower, and S. Gottlieb (QUDA), in International Conference for High Performance Computing, Networking, Storage and Analysis (2011), 1109.2935
2011 arXiv
-
[51]
M. A. Clark, B. Jo´ o, A. Strelchenko, M. Cheng, A. Gambhir, and R. C. Brower (QUDA), in International Conference for High Performance Computing, Networking, Storage and Analysis (2016), 1612.07873
2016 arXiv
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.