{"id":"cfe3ade2-59ef-4b7f-abdc-1a91fc46f0fc","arxiv_id":"2502.07438","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"The first lattice QCD calculation of single-channel DD scattering finds a weakly repulsive S-wave isovector interaction and a slightly attractive P-wave isoscalar interaction.","lead":"This paper reports the first lattice QCD calculation of DD meson scattering, finding a weak repulsion in the S-wave isovector channel and a slight attraction in the P-wave isoscalar channel. It matters because model-independent DD scattering input is key to studies of the Tcc tetraquark and the predicted DDK three-body state.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Physical r0 from two-point linear m_pi^2 extrapolation lacks any model/systematic error; alternative forms shift it by ~1 fm, underlining that Eq. (12) should be treated as preliminary.","rationale":"The paper's primary contribution is a first single-channel lattice QCD determination of DD scattering. The energy-level data are plausible: the A1+ levels are systematically above the free DD threshold, giving a negative S-wave scattering length, and the T1- levels give a positive P-wave phase shift. These qualitative findings are consistent across both pion masses and are not called into question by my analysis. The quantitative headline, however, rests on the chiral extrapolation in Eq. (11). With two data points the linear ansatz cannot be tested, and the paper itself notes this limitation. My concern is sharper for r0 than for a0: a0 is basically flat (-0.24 to -0.25) so the physical intercept is stable, whereas r0 moves by a factor of about 2.4 between the two masses and the extrapolated value is far outside the data. The effective range is a derivative-like quantity and is notoriously more sensitive to the choice of extrapolation. My own quick recomputation with m_pi and m_pi^3 gives r0_phys = -6.1 and -5.1 fm, respectively, compared to -5.5 fm with m_pi^2; this spread of about 1 fm is close to the reported 2.1 fm statistical error, and a chiral-log form could shift it further. Thus the error budget in Eq. (12) is incomplete. The abstract/full-text inconsistency (a0 = -0.25(8)(12) vs -0.26(5), r0 = -5.7(45)(17) vs -5.5(21)) is real but is a reporting issue; it does not change the physics. For these reasons I would keep the reader's CONDITIONAL verdict: the qualitative claim is credible, but the physical-pion-mass values should be labeled as model-dependent until a third pion mass or second lattice spacing is available.","tokens_in":13902,"tokens_out":16066,"duration_ms":151680,"concrete_test":"Using the Z=1 values in Table IV (m_pi=305 MeV: a0=-0.24(3), r0=-1.81(21); m_pi=207 MeV: a0=-0.25(3), r0=-4.3(14)), recompute the physical-pion-mass intercepts from Eq. (11) under at least three alternative two-parameter ansatze: linear in m_pi, linear in m_pi^2, and linear in m_pi^3, plus an NLO ChPT-inspired form with a pion-loop contribution. If the spread of r0_phys across these forms exceeds the quoted +/-2.1 fm (or a0_phys exceeds +/-0.05 fm), the central values in Eq. (12) are model-dominated and must carry an explicit systematic error; if the spread is much smaller, the reported precision is acceptable for a first calculation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative result, Eq. (12), is obtained by a two-point linear extrapolation in m_pi^2 (Eq. 11) using the Z=1 rows of Table IV. Since there are exactly two pion masses, the fit has zero degrees of freedom and the quoted uncertainties are purely statistical; no curvature, chiral-log, continuum, or finite-volume error is included. This is acknowledged in Section VI, but the numbers are still quoted as central results. The S-wave scattering length a0 is stable (about -0.24 to -0.26 fm across masses), so the qualitative 'weak repulsive' claim is robust. However, r0 changes from -1.81(21) fm at 305 MeV to -4.3(14) fm at 207 MeV, and the extrapolated r0_phys=-5.5(21) fm lies beyond the data range and depends entirely on the assumed linearity in m_pi^2. Repeating the same two-point fit with m_pi or m_pi^3 as the variable shifts r0_phys by roughly 1 fm, i.e., a model uncertainty comparable to the quoted statistical error. With one lattice spacing, no continuum estimate is possible, and the paper does not quote a systematic error.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a lattice QCD calculation of single-channel DD scattering in the isovector S-wave I(J^P)=1(0+) and isoscalar P-wave 0(1-) channels, using 2+1 flavor Wilson-Clover ensembles at a≈0.077 fm with m_pi≈305 and 207 MeV. The authors extract finite-volume energy levels with distillation and GEVP, apply a Lüscher quantization condition modified for the non-continuum D-meson dispersion relation (Z≠1), fit S- and P-wave effective-range expansion parameters, and extrapolate the S-wave scattering length and effective range to the physical pion mass via a linear form in m_pi^2, obtaining a0=-0.26(5) fm and r0=-5.5(21) fm. The paper concludes that the S-wave interaction is weakly repulsive and the P-wave slightly attractive.","tokens_in":14183,"tokens_out":6541,"duration_ms":55565,"significance":"If the results hold, this is the first lattice QCD determination of single-channel DD scattering in these channels and provides useful model-independent input for the Tcc system and DDK three-body calculations. The treatment of the D-meson dispersion relation in the Lüscher formula is a methodological plus, and the paper is transparent about the limitations of a two-point chiral extrapolation with one lattice spacing. The central qualitative conclusion—weak repulsive S-wave isovector interaction—is stable across pion masses and across the Z=1/Z≠1 treatments. The quantitative physical-pion-mass values, however, currently carry only statistical uncertainties and should be regarded as preliminary.","major_comments":[{"comment":"The abstract quoted in the submission metadata reports a0=-0.25±0.08±0.12 fm and r0=-5.7±4.5±1.7 fm, while the body abstract and Eq. (12) report a0=-0.26±0.05 fm and r0=-5.5±2.1 fm. These are two different versions of the central quantitative result, with different central values and different error decompositions; they must be reconciled before publication.","section":"Abstract and Sec. VI, Eq. (12)"},{"comment":"The physical-pion-mass values in Eq. (12) are obtained from a two-point linear fit in m_pi^2 with zero degrees of freedom, and the paper explicitly notes that the uncertainties are purely statistical. The extrapolated r0=-5.5(21) fm lies outside the range of the input points and depends entirely on the assumed linearity. Please add a model-variation estimate (e.g., comparing fits in m_pi, m_pi^2, and m_pi^3, or including a curvature term), or state clearly that this value is a preliminary extrapolation rather than a final result.","section":"Sec. VI, Eqs. (11)-(12)"},{"comment":"The evidence for a slightly attractive P-wave interaction is marginal: the P-wave scattering length is a1=1.7(21) fm^3 at 305 MeV and 0.9(17) fm^3 at 207 MeV, each within about 1.5 standard deviations of zero. Please quantify the significance of the attraction or temper the abstract and conclusion wording accordingly.","section":"Sec. V, Table IV"}],"minor_comments":[{"comment":"The first sentence contains a typo: 'Qauntum Chromodynamics' should be 'Quantum Chromodynamics'.","section":"Sec. I"},{"comment":"In the paragraph after Table IV, 'isoscaclar' should be 'isoscalar', and 'differences of of the ERE parameters' has a duplicated 'of'.","section":"Sec. V"},{"comment":"The phrase 'the Z values in Table II bare uncertainties' should read 'bear uncertainties' or 'carry uncertainties'.","section":"Sec. V"},{"comment":"The conclusions state 'at two different pion masses 207 and 350 MeV', but Table I lists 305 MeV; also the channel '0(1+)' should be '0(1-)'.","section":"Sec. VII"},{"comment":"The figure labels 'P303' and 'P210' appear inconsistent with the ensemble names P30 and P21; please clarify or rename the labels.","section":"Figs. 4 and 5"},{"comment":"The horizontal axes in Fig. 7 are labeled 'm [GeV]' without a subscript; use 'm_pi [GeV]' for clarity.","section":"Fig. 7"},{"comment":"The zeta function Zjs is used in Eqs. (8) and (A10) but defined only in Eq. (A11); define it at first use.","section":"Appendix A, Eq. (A10)"}],"recommendation":"major_revision","confidential_remarks":"The paper's qualitative conclusions are defensible, but the headline numbers in Eq. (12) are not yet controlled predictions because the two-point extrapolation has no systematic or model uncertainty. The discrepancy between the metadata abstract and the body abstract is a consistency problem that must be fixed. Given the first-of-its-kind nature of the calculation, major revision is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me cut to it: this is the first LQCD calculation of single-channel DD scattering in I=1(0+) and I=0(1-), and that alone makes it worth reading. The qualitative result—weak repulsion in the S-wave, weak attraction in the P-wave—is stable across pion masses and across the Z=1 vs Z≠1 treatments. The methodology is sound and honestly executed: distillation, GEVP, modified Lüscher for the noncontinuum dispersion relation, and two volumes at each mass to get more kinematic points. The paper also clearly states its own limitations at the end.\n\nThe soft spots are exactly where the reader's report puts them. The headline numbers at the physical pion mass, Eq. (12), come from a two-point linear fit in m_pi^2. With only two pion masses and one lattice spacing, the system is fully constrained; there is no way to test curvature, chiral logs, or continuum effects. The paper acknowledges this, but still presents a0=-0.26(5) and r0=-5.5(21) as central results with statistical errors only. The stress-test point is fair: changing the variable to m_pi or m_pi^3 shifts r0 by roughly 1 fm. So those numbers are preliminary, not determinations.\n\nThere are also small internal inconsistencies that a referee would want cleaned up: the abstract's quoted central values with two uncertainties don't match the full-text one-sigma values; the conclusions mention m_pi=350 MeV and I=0(1+) which contradict 305 MeV and 0(1-) elsewhere; and 'isoscaclar' is a typo. More substantively, F32P21's energy levels are excluded from the ERE fits, so the m_pi=207 MeV column in Table IV rests on one volume, which further thins the chiral extrapolation.\n\nWho is this for? People working on Tcc widths and DDK three-body systems. It gives them a first-principles input where previously they had models. The citation pattern looks clean; the prior lattice work on DD* and coupled channels is properly credited. My summary: the qualitative physics is credible, the methodology is right, but the quantitative physical-pion values should be reported with a model-dependence caveat, and ideally with more ensembles before being used as central inputs.\n\nI'd send this to a serious referee; it's a legitimate first result with truthful limitations, and the remaining issues are fixable in revision. If I were citing it, I'd cite it for the qualitative S-wave repulsion and as a first calculation, not for the precise r0.","headline":"First lattice QCD calculation of single-channel DD scattering, with a robust qualitative message but quantitative physical-pion values that rest on a two-point extrapolation.","tokens_in":14770,"tokens_out":1821,"would_cite":true,"duration_ms":16123,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The first lattice QCD calculation of single-channel $DD$ scattering finds a weakly repulsive isovector S-wave, with physical-pion scattering length $a_0=-0.26\\pm0.05$ fm, and a slightly attractive isoscalar P-wave.","keywords":["lattice QCD","DD scattering","Lüscher finite-volume method","scattering length","effective range expansion","charmed mesons","Tcc(3875)+","Wilson-Clover fermions"],"falsifier":"Compute the same $A_1^+$ spectrum on a third ensemble with $m_\\pi\\approx 250$ MeV, or at a second lattice spacing, using the identical analysis; if the resulting $a_0$ does not fall on the straight line through the 207 and 305 MeV values within statistical errors, the linear extrapolation used for Eq. (12) is contradicted.","tokens_in":13638,"feed_emoji":"⚛️","tokens_out":10497,"duration_ms":85892,"temperature":0.7,"pith_summary":"This paper reports the first lattice QCD calculation of single-channel $DD$ scattering for the quantum numbers $I(J^P)=1(0^+)$ and $0(1^-)$, using the Lüscher finite-volume method on $2+1$ flavor Wilson-Clover ensembles with a lattice spacing $a\\simeq0.077$ fm and pion masses near 207 and 305 MeV. The central result is that the isovector S-wave interaction is weakly repulsive, with a physical-pion scattering length $a_0=-0.26\\pm0.05$ fm and effective range $r_0=-5.5\\pm2.1$ fm, while the isoscalar P-wave channel is slightly attractive. If correct, these are the first model-independent numbers for $DD$ final-state interactions below the $D^*D^*$ threshold, and they provide input for calculations of the $T_{cc}^+$ width and for predictions of a possible $DDK$ three-body bound state. The authors state that a more careful chiral and continuum extrapolation will require additional pion masses and lattice spacings.","feed_headline":"First lattice-QCD D-D scattering result: weak repulsion","feed_subtitle":"Physical-pion S-wave scattering length lands at -0.26 fm, tightening Tcc and DDK predictions.","key_machinery":"The load-bearing object is Lüscher's finite-volume quantization condition, which relates two-particle energy levels in a box to continuum phase shifts; the paper modifies the standard condition by replacing the continuum dispersion relation with $E^2=m_D^2+Zp^2$ for each particle, introducing the factor $(\\omega_1+\\omega_2)/(Z_1\\omega_2+Z_2\\omega_1)$ in the zeta-function matrix. The energy levels themselves are obtained with the generalized eigenvalue problem applied to interpolating operators projected onto the $A_1^+$, $E^+$, $T_2^+$ (isovector) and $T_1^-$, $T_2^-$, $A_2^-$ (isoscalar) irreps, and the phase shifts are parameterized by the effective-range expansion $\\cot\\delta_l=p^{-2l-1}(1/a_l+\\tfrac12 r_l p^2)$, which is fitted directly to the spectra. This combination turns a small number of lattice energy levels into scattering lengths and effective ranges.","core_discovery":"Using interacting energy levels extracted from correlation functions in the $A_1^+$ and $T_1^-$ irreps and fitted with the effective-range expansion, the paper determines the S-wave isovector and P-wave isoscalar $DD$ scattering parameters at two pion masses. It finds negative $a_0$ and $r_0$, indicating weak repulsion, and positive $a_1$ and $r_1$, indicating slight attraction; the near-threshold energy levels are close to the free two-$D$-meson levels, so higher partial waves are consistent with negligible interaction. Assuming the continuum dispersion relation, a linear extrapolation in $m_\\pi^2$ gives $a_0^{\\mathrm{phy}}=-0.26(5)$ fm and $r_0^{\\mathrm{phy}}=-5.5(21)$ fm at the physical pion mass, with the same $a_0$ within errors at the chiral limit. The paper does not extrapolate the P-wave parameters because their statistical uncertainties are large. It also notes that a single lattice spacing prevents a continuum extrapolation.","pith_inferences":["Because the extrapolation is a straight line through two points, the agreement between the chiral-limit value $a_0^\\chi=-0.26(6)$ fm and the physical value is essentially built into the fit; a third pion mass would test whether the true curvature is small enough for this conclusion to survive.","The pattern $|a_0|\\ll|r_0|$ matches what is seen in $I=2$ $\\pi\\pi$ and $I=1$ $KK$ scattering, suggesting a common near-threshold behavior for repulsive hadron-hadron channels that could be probed by collecting more repulsive S-wave systems.","One testable extension is to repeat the extraction with the P-wave parameters included in a coupled-channel $DD$-$D^*D^*$ analysis; the present single-channel values would then serve as a constrained starting point rather than final answers.","A natural next computation is the same analysis on a finer lattice at a different lattice spacing; if the extracted $Z$ factor changes the energy levels significantly, the role of the modified Lüscher formula will be directly visible."],"forward_implications":["If the central values hold, the $I=1$ S-wave $DD$ interaction is weakly repulsive, so no $DD$ bound state appears in this channel near threshold.","The physical-pion $a_0$ and $r_0$ provide direct input for quantifying the effect of $DD$ final-state interactions on the predicted $T_{cc}^+$ width.","The same parameters serve as two-body input in continuum and finite-volume studies of a possible $DDK$ three-body bound state.","The modified finite-volume quantization condition, which accounts for a lattice dispersion relation with $Z\\neq 1$, can be carried over to other heavy-meson scattering problems where such artifacts matter.","Because only two pion masses and one lattice spacing are used, additional ensembles are needed before the physical-pion result can be considered reliable, as the paper itself notes."],"supporting_citations":[{"why":"Introduces the finite-volume quantization condition that connects lattice spectra to scattering phase shifts.","marker":"[1]"},{"why":"Extends the quantization condition used by this paper to two-particle systems.","marker":"[2]"},{"why":"Provides the partial-wave reduction and zeta-function form of the quantization condition employed in the fits.","marker":"[3]"},{"why":"Supplies the $2+1$ flavor Wilson-Clover gauge ensembles used for the calculation.","marker":"[46]"},{"why":"Describes the distillation method used to build precise interpolating operators and correlation functions.","marker":"[54]"},{"why":"Sets up the generalized eigenvalue problem used to extract the finite-volume energy levels.","marker":"[55]"},{"why":"Reports the analogous $|a_0|\\ll|r_0|$ pattern in $I=2$ $\\pi\\pi$ scattering, used for comparison.","marker":"[10]"},{"why":"Gives the same $|a_0|\\ll|r_0|$ pattern in $I=1$ $KK$ scattering, used for comparison.","marker":"[16]"}],"fun_headline_variants":["First lattice QCD DD scattering: weak repulsion","Lattice QCD finds weak repulsion in DD scattering","New lattice QCD calculation of DD scattering shows repulsion","DD scattering in lattice QCD: first result is repulsive"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quoted physical-pion values assume that $a_0$ and $r_0$ are exactly linear functions of $m_\\pi^2$ between 207 and 305 MeV; with only two pion masses and one lattice spacing, any curvature, chiral logarithm, or lattice artifact would shift the results beyond the reported uncertainties.","fun_headline_variants_meta":{"raw":{"variants":["First lattice QCD DD scattering: weak repulsion","Lattice QCD finds weak repulsion in DD scattering","New lattice QCD calculation of DD scattering shows repulsion","DD scattering in lattice QCD: first result is repulsive"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000533,"raw_usage":{"total_tokens":2559,"prompt_tokens":938,"completion_tokens":1621,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":554,"completion_tokens_details":{"reasoning_tokens":1554}},"tokens_in":554,"tokens_out":1621,"duration_ms":10831,"temperature":1.0,"reasoning_tokens":1554,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T12:45:30.437326+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same $A_1^+$ spectrum on a third ensemble with $m_\\pi\\approx 250$ MeV, or at a second lattice spacing, using the identical analysis; if the resulting $a_0$ does not fall on the straight line through the 207 and 305 MeV values within statistical errors, the linear extrapolation used for Eq. (12) is contradicted.","supporting_citations":[{"cited_title":"L¨ uscher, Commun","cited_arxiv_id":null,"evidence_quote":"Introduces the finite-volume quantization condition that connects lattice spectra to scattering phase shifts."},{"cited_title":"L¨ uscher and U","cited_arxiv_id":null,"evidence_quote":"Sets up the generalized eigenvalue problem used to extract the finite-volume energy levels."},{"cited_title":"The I=2 pipi S-wave Scattering Phase Shift from Lattice QCD","cited_arxiv_id":"1107.5023","evidence_quote":"Reports the analogous $|a_0|\\ll|r_0|$ pattern in $I=2$ $\\pi\\pi$ scattering, used for comparison."},{"cited_title":"Interactions of two and three mesons including higher partial waves from lattice QCD","cited_arxiv_id":"2106.05590","evidence_quote":"Gives the same $|a_0|\\ll|r_0|$ pattern in $I=1$ $KK$ scattering, used for comparison."}],"review_version":1}