{"id":"4b12c04c-4625-490f-bad3-a6e237174e2a","arxiv_id":"2501.00283","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Using OPE matched to Pionless EFT, the authors reproduce the deuteron's high-momentum nucleon distribution from scattering-state inputs, validating the method against the AV18 potential.","lead":"This paper combines the operator product expansion with a low-energy nuclear theory to compute the high-momentum tail of the nucleon momentum distribution in the deuteron. It demonstrates the matching procedure on the AV18 potential and offers a route from short-range nuclear forces to scattering data or lattice QCD.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The OPE basis in Eqs. (71) and (85)-(88) omits the O2(SD)/tensor operator while the AV18 A00 input includes coupled-channel D-wave T-matrix contributions; the few-percent agreement with rho18 is therefore not yet established as an OPE prediction.","rationale":"Good-faith reading: the paper is a proof-of-principle that OPE plus Pionless EFT can compute high-momentum tails from scattering-state matching. The toy model is exact and the algebraic construction is coherent. The weakest point is the AV18 comparison: the calculation uses a truncated operator basis. My concern is not that the paper is internally inconsistent, but that the claimed few-percent accuracy is not supported without controlling the O2(SD)/tensor operator. The reader's weakest_assumption identified the same D-wave issue, so I agree with the conditional verdict. I would not reject the paper: the missing check is well-defined and the authors already have the machinery to perform it. The concrete test above would settle whether the omission matters. If inclusion of O2(SD) leaves the curves unchanged at the 1% level, the central claim stands as stated; if not, the claim needs qualification.","tokens_in":17314,"tokens_out":10301,"duration_ms":112188,"concrete_test":"Re-run the AV18 matching with O2(SD) included: add fW_SD(k) M_SD(p) to Eq. (87) (and to Eq. (71)), with M_SD(p) the Pionless 3S1-3D1 matrix element generated at leading nontrivial order by C2(SD), and add fW_SD(k)<d|O2(SD)|d> to the deuteron OPE. Keep p1=4 MeV, p2=5 MeV, Lambda=800 MeV. If the new curve shifts from Fig. 6 by more than ~1% for 250<k<500 MeV, the S-wave-only basis controls the headline claim; if the shift is below ~1%, the omission is benign.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (19) defines A_{l'l}(k,p) as a sum over l'' of T_{l'l''}T_{l''l}; for the AV18 3S1-3D1 channel, A00 receives D-wave contributions through T_{02} and T_{20}. The OPE side, Eq. (71), contains only Pionless 3S1 operators O0 and O2 (with the NLO mixing term fW0^{O2}); the O2(SD) operator of Eq. (27) is never included in the matching basis or in the deuteron expectation values (69)-(70), which are computed with a purely S-wave Pionless deuteron. Matching at p1=4 and p2=5 MeV (Eqs. (85)-(88)) can therefore absorb D-wave/tensor dynamics into fW0 and fW2, but those same coefficients are then multiplied by S-wave-only deuteron matrix elements, while the AV18 benchmark rho18(k) includes the D-state contribution Gamma_D^2/(Bd+k^2/m)^2. No power-counting estimate is given for the size of the O2(SD) contribution to rho(k) for k>250 MeV, so the quoted few-percent agreement could be partly accidental. The problem is compounded by p1,p2 << 1/a (~36 MeV), where M2^LO/M0^LO ~ p^2/a^{-2} ~ 10^-2, making the extraction of fW2 from small differences numerically fragile.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an operator-product-expansion (OPE) framework, combined with Pionless effective field theory, for the large-momentum part of the single-nucleon momentum distribution in the deuteron. Wilson coefficients for the local Pionless operators O0 and O2 are fixed by matching exact matrix elements of the nonlocal operator Omega(k) between NN scattering states computed in an underlying nuclear potential and in the EFT. The method is checked analytically in a separable toy model and then applied numerically to the AV18 potential. The authors report that the OPE/EFT reproduces the AV18 deuteron momentum distribution to a few percent for k > 250 MeV once NLO EFT corrections and the O2 operator are included.","tokens_in":17684,"tokens_out":5106,"duration_ms":58278,"significance":"If the central claim holds, the paper offers a systematic, state-independent framework for short-range correlations that connects short-distance nuclear dynamics to Wilson coefficients, and it provides a concrete alternative to generalized contact formalisms. The paper has several genuine strengths: the toy-model check is analytic and exact; the renormalized operator basis in Eqs. (58)-(62) is constructed explicitly; and the matching is performed with scattering states rather than by fitting the deuteron momentum distribution itself. The AV18 comparison is a meaningful benchmark. However, the few-percent claim for AV18 currently lacks a quantitative treatment of the D-wave/tensor content of the deuteron and of the sensitivity to the matching procedure, so the soundness of the central numerical result is not yet fully established.","major_comments":[{"comment":"The matching is performed at p1 = 4 MeV and p2 = 5 MeV, both much smaller than the 3S1 inverse scattering length a^{-1} ~ 36 MeV. Since M2^LO/M0^LO scales as p^2 a^2, the O2 Wilson coefficient is extracted from a very small difference between A00(k,p1) and A00(k,p2). The numerical fragility of this extraction is not discussed, and no stability test with respect to the choice of (p1,p2) is shown. Please demonstrate that fW2 and the resulting rho(k) are stable when the matching momenta are varied within the Pionless EFT range, or quantify the induced uncertainty.","section":"Sec. V B, Eqs. (85)-(88)"},{"comment":"The central claim of few-percent agreement is presented without uncertainty quantification. The cutoff Lambda = 800 MeV and the matching points p1 = 4 MeV, p2 = 5 MeV are single choices, and the residual cutoff dependence noted around Eq. (40) is not scanned. In addition, the numerical values of the Wilson coefficients fW0 and fW2 as functions of k are not reported, so the relative size of the O0 and O2 terms cannot be checked from the paper. Please provide numerical tables or curves for the Wilson coefficients and error bands obtained by varying Lambda and the matching momenta.","section":"Sec. V B and Figs. 6-7"}],"minor_comments":[{"comment":"The abstract contains the phrase “short-rang structure”; this should be “short-range structure.”","section":"Abstract"},{"comment":"The caption says “The rest are the results from the OPE-Pionless EFT” without identifying which curve corresponds to which OPE order; please clarify the legend or caption.","section":"Fig. 6 caption"},{"comment":"Because the same AV18 potential supplies both the scattering-state input and the deuteron benchmark, the test is not circular, but it validates state independence only within a single underlying model. Testing the matching procedure with a second potential, for example a chiral EFT potential, would strengthen the claim that the extracted Wilson coefficients are state independent.","section":"Sec. V B"}],"recommendation":"major_revision","confidential_remarks":"The paper is a serious contribution within the scope of the journal. My main reservation is the missing treatment, or at least a quantitative estimate, of the D-wave/tensor operator O2(SD) in the AV18 comparison. This is load-bearing for the few-percent claim and should be addressed before publication. The absence of numerical Wilson coefficients and of sensitivity tests to Lambda and (p1,p2) compounds the issue. I do not see grounds for rejection, but the required additions are substantial rather than purely editorial."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper develops a matching prescription for OPE Wilson coefficients using NN scattering states at low momenta, then applies it to the deuteron momentum distribution. What is new is the multiple-kinematic-point matching scheme (p1=4, p2=5 MeV) with NLO Pionless-EFT corrections and O2-to-O0 operator mixing; that is a genuine technical step beyond unitary-Fermi-gas OPE and SRG-evolved operators. The toy model in Sec. V A is the cleanest part: it is analytic, exact, and shows how the double expansion in 1/(ak) and 1/(aMhi) works term by term. The paper is clearly written and honest about the framework's scope.\n\nThe soft spot is real. The stress-test concern lands. Eq. (19) for A00 includes D-wave contributions through T02 and T20, but Eq. (71) contains only S-wave operators O0 and O2. Matching at two soft momenta can absorb D-wave dynamics into fW0 and fW2, yet the deuteron matrix elements in Eqs. (69)-(70) are computed with an S-wave-only Pionless deuteron. Since the AV18 rho(k) includes the D-state vertex contribution, the quoted few-percent agreement is not yet established as an OPE prediction unless the O2(SD) contribution is demonstrably negligible. The paper gives no power-counting or numerical estimate for that contribution, and the missing O2(SD) operator in the matching basis is never justified. Second, there is no sensitivity study on p1, p2, or Lambda. The extreme softness of p1,p2 (4-5 MeV, compared to a^{-1} ~ 36 MeV) makes the extraction of fW2 from M2/M0 ~ 10^{-2} numerically fragile. Third, the numerical Wilson coefficients and the actual curves behind Fig. 6 are not reported, so the few-percent claim is not reproducible from the text. The regulator check in Sec. IV is a good start but covers only one systematic.\n\nOverall, the framework is coherent, the toy model is exact, and the AV18 result is promising. But the central numerical claim needs a direct D-wave accounting (or an explicit demonstration that D-wave is irrelevant for k>250 MeV) plus a sensitivity analysis on the matching points and cutoff. This paper deserves a serious referee; I would send it out, expecting major revision before publication.","headline":"Solid OPE+Pionless-EFT matching scheme with a clean toy model, but the AV18 few-percent claim needs a D-wave accounting and sensitivity analysis before it is convincing.","tokens_in":18170,"tokens_out":6029,"would_cite":true,"duration_ms":64951,"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 paper shows that combining the operator product expansion with Pionless effective field theory reproduces the deuteron's high-momentum distribution to a few percent once Wilson coefficients are fixed by matching nucleon-nucleon…","keywords":["operator product expansion","nucleon momentum distribution","Pionless effective field theory","short-range correlations","deuteron","Wilson coefficients","AV18 potential","nuclear EFT"],"falsifier":"Repeat the matching while including the $^3D_1$ component of the deuteron or the $O_2(SD)$ mixing operator; if the OPE/EFT prediction for $\\rho(k)$ shifts by more than the claimed few percent over $k=250$-$600$ MeV, the S-wave-only operator basis is insufficient.","tokens_in":17120,"feed_emoji":"⚛️","tokens_out":5777,"duration_ms":55350,"temperature":0.7,"pith_summary":"The paper argues that the operator product expansion (OPE), combined with Pionless effective field theory, gives a systematic way to compute the large-momentum part of the single-nucleon momentum distribution in nuclei, the part governed by short-range nucleon-nucleon correlations. Its central move is to fix the OPE's Wilson coefficients not in the bound state but in low-energy nucleon-nucleon scattering states, where the matrix elements can be calculated both exactly in an underlying potential model and in the effective field theory. Matching the two at small momenta determines coefficients that are independent of the external state, so the same coefficients can later be used for other nuclei. Tested on the deuteron, the OPE/EFT expression reproduces the AV18 momentum distribution to within a few percent for $k>250$ MeV, and a separable toy model shows the expansion organizes into powers of $1/(ak)$ and $1/(aM_{\\mathrm{hi}})$. A reader should care because this offers a parameter-free route from soft nuclear interactions to hard-probe observables such as short-range correlations.","feed_headline":"OPE matches deuteron's high-momentum tail to a few percent","feed_subtitle":"Wilson coefficients fixed by nucleon-nucleon scattering states let Pionless EFT reach momenta above 250 MeV.","key_machinery":"The load-bearing object is the OPE of the nonlocal nucleon-pair operator $\\Omega(\\vec{k})$ whose expectation value defines the single-nucleon momentum distribution. For short separation, the paper expands this operator into local Pionless EFT operators, the leading ones being the one-body density $N^\\dagger N$ and the two-body contact operators $O_0(^3S_1)$ and $O_2(^3S_1)$, with Fourier-transformed Wilson coefficients $\\widetilde W_n(k)$. The coefficients are determined by matching the exact potential-model matrix elements $A(k,p)$ to the EFT matrix elements $M_n(p)$ through Eqs. (85)-(88), using the renormalized matrix elements $M_{R0}$ and $M_{R2}$ and the deuteron expectation values of the renormalized operators. The expansion is a double series in $1/(ak)$, controlled by higher-dimension operators, and $1/(aM_{\\mathrm{hi}})$, controlled by higher-order EFT corrections to the Wilson coefficients.","core_discovery":"The central claim is that the large-momentum single-nucleon momentum distribution $\\rho(k)$ can be written as an OPE of the nonlocal operator $\\Omega(\\vec{k})=\\int d^3r\\,e^{-i\\vec{k}\\cdot\\vec{r}}N^\\dagger(-\\vec{r}/2)N(\\vec{r}/2)$ into local Pionless EFT operators, with state-independent Wilson coefficients $\\widetilde W_n(k)$. These coefficients are fixed by matching the exact matrix elements of $\\Omega(\\vec{k})$ between $^3S_1$-$^3D_1$ nucleon-nucleon scattering states computed from the underlying potential against the same matrix elements computed in Pionless EFT, at the soft momenta $p_1=4$ MeV and $p_2=5$ MeV. The resulting OPE/EFT expression reproduces the AV18 deuteron momentum distribution at the few-percent level for $k>250$ MeV, and the numerical results show that the next-to-leading-order EFT correction to the leading operator matters more than adding the next operator $O_2$ in the kinematic region studied. The same matching logic is demonstrated analytically in a separable toy model, where the OPE/EFT reproduces the exact momentum distribution term by term as a double expansion in $1/(ak)$ and $1/(aM_{\\mathrm{hi}})$.","pith_inferences":["The same matching technology should extend to other short-range-dominated observables, such as the two-nucleon momentum distribution or generalized contact parameters, because the separation into state-independent Wilson coefficients and state-dependent EFT matrix elements is generic.","The S-wave-only assumption could be tested by adding the $^3D_1$ channel or the $O_2(SD)$ mixing operator; the claimed few-percent agreement would be reinforced if these extra operators barely change $\\rho(k)$, and undermined if they shift it significantly.","A direct experimental check is possible through $y$-scaling electron-scattering data, where the extracted deuteron momentum distribution should follow the OPE prediction and deviations would locate where the one-body Wilson coefficient stops being negligible."],"forward_implications":["Once matched in two-body scattering, the same Wilson coefficients apply to any nucleus: computing $\\rho_A(k)$ requires only the Pionless EFT matrix element $\\langle A|O_n|A\\rangle$.","At very large $k$, the one-term factorization holds and ratios such as $\\rho_A(k)/\\rho_d(k)$ become approximately constant; at intermediate $k$ the $O_2$ term introduces $k$-dependent corrections to those ratios.","For $k$ around the pion mass, including the NLO EFT interaction matters more than adding the $O_2$ operator, so the practical ordering is EFT order first, then operator dimension.","The matching procedure is ready to accept lattice-QCD inputs for the two-body matrix elements, replacing the phenomenological potential as the underlying theory."],"supporting_citations":[{"why":"Supplies the AV18 nucleon-nucleon potential used as the underlying theory for the central numerical comparison.","marker":"[38]"},{"why":"Provides the variational Monte Carlo benchmark momentum distributions for light nuclei with AV18 that motivate and anchor the comparison.","marker":"[36]"},{"why":"Introduces the separable toy potential whose exact scattering-state matrix elements give analytic Wilson coefficients.","marker":"[56, 57]"},{"why":"Supplies the Pionless EFT Lagrangian operators and the operator basis into which the OPE is expanded.","marker":"[51]"},{"why":"Provides the Pionless EFT framework and power counting used for the EFT side of the matching.","marker":"[50]"},{"why":"Gives the related treatment of short-range operators under unitary transformations, which the paper contrasts with its OPE matching and cites for the vanishing one-body coefficient.","marker":"[30, 31]"},{"why":"Establishes the unitary Fermi gas large-momentum behavior that the OPE coefficients reproduce in the small-$k$ limit.","marker":"[39, 40]"},{"why":"Provides the perturbative NLO EFT potentials and renormalization conditions used in the matching equations.","marker":"[53, 58]"}],"fun_headline_variants":["OPE matches deuteron's high-momentum tail to few percent","Wilson coefficients from NN scattering states pin down deuteron tail","Few-percent accuracy for deuteron momentum distribution via OPE","Operator expansion reproduces deuteron's short-range structure","Systematic OPE for large-momentum nuclear distributions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole comparison assumes that matching $^3S_1$ scattering states at $p_1=4$ MeV and $p_2=5$ MeV, using only the $^3S_1$ operators $O_0$ and $O_2$ with no $^3D_1$ or tensor-operator contributions, captures all the short-distance dynamics that determines the deuteron momentum distribution for $k>250$ MeV.","fun_headline_variants_meta":{"raw":{"variants":["OPE matches deuteron's high-momentum tail to few percent","Wilson coefficients from NN scattering states pin down deuteron tail","Few-percent accuracy for deuteron momentum distribution via OPE","Operator expansion reproduces deuteron's short-range structure","Systematic OPE for large-momentum nuclear distributions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000666,"raw_usage":{"total_tokens":3028,"prompt_tokens":924,"completion_tokens":2104,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":540,"completion_tokens_details":{"reasoning_tokens":2021}},"tokens_in":540,"tokens_out":2104,"duration_ms":39789,"temperature":1.0,"reasoning_tokens":2021,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:54:22.918321+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the matching while including the $^3D_1$ component of the deuteron or the $O_2(SD)$ mixing operator; if the OPE/EFT prediction for $\\rho(k)$ shifts by more than the claimed few percent over $k=250$-$600$ MeV, the S-wave-only operator basis is insufficient.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the variational Monte Carlo benchmark momentum distributions for light nuclei with AV18 that motivate and anchor the comparison."}],"review_version":1}