{"id":"20fc0afd-9aac-42ef-979d-838d07d3c2b2","arxiv_id":"2412.07717","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Chiral perturbation theory predicts an exponentially suppressed dbar minus ubar asymmetry in the proton's TMD PDFs at large transverse separation.","lead":"The authors use chiral perturbation theory to compute the proton's dbar minus ubar asymmetry in transverse momentum dependent parton distributions, finding an exponential suppression at large transverse distances. If confirmed by future measurements, the result would give a first-principles description of the long-distance regime of TMD PDFs, a region usually handled with ad hoc models.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (13) replaces the pion's valence TMD PDF with its collinear limit; at bT ~ 1/m_pi, Lambda_QCD bT ~ 1, so this replacement is uncontrolled and the chiral K0 suppression is not established as the dominant long-distance mechanism.","rationale":"The paper is honest and exploratory: it provides explicit analytic formulas, a clear collinear comparison with E866/E906, and an openly stated uncertainty caveat. The collinear part after the ymax prescription is a reasonable phenomenological study. For the TMD central claim, however, the weakest step is not only the ymax delta-function expansion but also the replacement of the pion valence TMD PDF by its collinear limit in Eq. (13). This replacement is justified only if the pion TMD's intrinsic transverse momentum is parametrically Lambda_chi; physically, confinement suggests it is Lambda_QCD. At the bT values plotted, bT Lambda_QCD ~ O(1), so the neglected bT dependence is not a small correction. The authors themselves show in Fig. 4 that adding a Lambda_QCD-scale Gaussian changes the prediction dramatically, which undercuts the phrase 'without ad hoc model input.' A conditional acceptance is appropriate: the qualitative statement that chiral physics produces exponential suppression at large bT may survive, but the quantitative claim that this suppression is the chiral HDF's K0 tail requires a controlled treatment of the pion TMD. The reader's CONDITIONAL verdict already captures this, so no change in verdict is needed.","tokens_in":12920,"tokens_out":7368,"duration_ms":75266,"concrete_test":"Compute the pion valence TMD PDF at physical pion mass from lattice QCD (or extract its intrinsic width from pion-induced Drell-Yan data), then evaluate the full convolution in Eq. (9) using the TMD HDF in Appendix A. If the resulting bar-d(x,bT) - bar-u(x,bT) asymmetry at bT = 5 GeV^-1 differs from the Eq. (13) curve in Fig. 3 by more than a factor of 2, the central claim that the chiral HDF controls the long-distance shape is falsified; if the curves agree within uncertainties, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (13) is obtained from the exact TMD matching (9) by treating the pion valence TMD PDF qv_H(0)(alpha, pT) as a high-energy object with pT ~ Lambda_chi and expanding in Lambda_QCD/Lambda_chi, keeping only the collinear pion PDF. The central large-bT prediction is then essentially the TMD HDF's K0(bT m_pi) falloff. However, qv_H(0) is a nonperturbative pion TMD; its intrinsic transverse-momentum width is set by confinement, of order Lambda_QCD, not Lambda_chi. At bT ~ 1/m_pi ~ 5 GeV^-1, the relevant expansion parameter is bT Lambda_QCD ~ 1, not Lambda_QCD/Lambda_chi ~ 0.2, so the truncation in Eq. (13) is uncontrolled precisely in the regime where the chiral suppression is advertised. The paper's own Fig. 4 demonstrates the sensitivity: replacing the delta-function approximation by a Gaussian with kappa = 200-400 MeV changes the asymmetry at large bT by roughly an order of magnitude. Thus the 'chiPT only' and 'chiPT x pQCD' curves in Fig. 3 conflate the chiral HDF with the unknown intrinsic bT dependence of the pion TMD, and the claim that the chiral TMD HDF determines the long-distance TMD shape without ad hoc model input is not established. The ymax/delta-function issue identified by the reader is also real, but for the TMD-specific central claim this pion-TMD truncation is at least as load-bearing.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript uses SU(2) chiral perturbation theory to compute the pion/baryon hadronic distribution functions (HDFs) of the proton, first for collinear PDFs and then for TMD PDFs. It introduces a momentum-fraction cutoff ymax above which the pion HDF is asserted to vanish, fits ymax to the E866/E906 dbar-ubar data, and uses the same prescription to predict the bT-space dbar-ubar asymmetry. The paper finds that the TMD HDF falls off like K0(bT m_pi) at large bT, which produces exponential suppression relative to a fixed-order pQCD matching computation, and it closes by showing that a Gaussian model for the pion valence TMD changes this prediction substantially.","tokens_in":13303,"tokens_out":8445,"duration_ms":79528,"significance":"If the central chain of reasoning were established, the paper would offer a QCD-based, regulator-independent computation of the long-distance TMD shape and a concrete experimental target for separating chiral and confinement effects. The analytic results in Appendix A are explicit, the K0(bT m_pi) asymptotics of the TMD HDF is a genuine mathematical feature of the calculation, and the use of dimensional regularization with MS for the collinear HDFs is a cleaner framework than the cutoffs used in earlier chiral studies. These are real strengths. However, as detailed below, the advertised prediction is built on an uncontrolled replacement of the pion TMD PDF by its collinear limit and on a ymax prescription that is tuned to the data it validates, so the significance of the result is currently more exploratory than definitive.","major_comments":[{"comment":"The replacement of the pion valence TMD PDF q_v,H^(0)(alpha,pT) by its collinear limit is not controlled in the bT region that carries the paper's central claim. The argument that pT ~ Lambda_chi is an assumption; the intrinsic transverse-momentum scale of a confined pion is of order Lambda_QCD. At bT ~ 1/m_pi ~ 5 GeV^{-1}, where the advertised K0 suppression sets in, Lambda_QCD bT is of order one, so the O(Lambda_QCD/Lambda_chi) truncation in Eq. (13) is not justified. The paper's own Fig. 4 makes the impact concrete: a Gaussian in Eq. (17) with kappa = 200-400 MeV changes the asymmetry at large bT by roughly an order of magnitude. The 'chiPT only' and 'chiPT x pQCD' curves in Fig. 3 therefore include unknown intrinsic bT dependence of the pion TMD, and the abstract's claim that the effective theory gives a natural exponential suppression of the TMD PDF is not established.","section":"TMD matching, Eqs. (12)-(13) and Fig. 4"},{"comment":"The assertion f_pi p(y)=0 for y>ymax rests on the replacement delta(k+ - yP+) -> delta(yP+), but this is not a valid distribution identity; for a soft test function f(k+) the exact expression gives f(yP+), not zero. The vanishing is an operator-level assumption about the support of the EFT matrix elements, and it is load-bearing for both Eq. (8) and the TMD formulas (13)-(15). In addition, ymax is a free parameter chosen so the collinear calculation reproduces E866/E906, and the same tuned interval (0.3-0.35 in Fig. 3) is then used for the TMD predictions. The calculation is regulator-independent in the loop sense, but it is not parameter-free; the uncertainty bands in Fig. 3 underestimate the model uncertainty because they only vary ymax over a narrow range.","section":"Collinear HDFs, Eqs. (5)-(8)"},{"comment":"The conclusion that chiPT produces a suppression at large bT that pQCD cannot replicate is stronger than what is computed. The green curve is built from fixed-order matching coefficients at mu^2 = zeta^2 = 54 GeV^2; it is not a complete pQCD TMD prediction, which would include resummed Sudakov evolution and would itself suppress the distribution at large bT. The comparison should either include the full perturbative TMD prediction at the same accuracy or be worded as a comparison with fixed-order matching only.","section":"Numerical results, Eq. (15) and Fig. 3"}],"minor_comments":[{"comment":"The notation '0.25 ~ 2m_pi/M and 0.4 ~ 3m_pi/M' is awkward; it should read ymax in [0.25, 0.4], corresponding to roughly 2-3 m_pi/M.","section":"Page 3, near Fig. 2"},{"comment":"The coupling gNDelta = 3 sqrt(2) gA/5 and the sign of f_pi- Delta++ relative to f_pi+ Delta0 should be stated more explicitly, since the cancellation in Eq. (A1) is important for the final asymmetry.","section":"Appendix A, Eqs. (A3) and (A6)"},{"comment":"The uncertainty 'from varying the result by O(m_pi/Lambda_chi)' is not defined; please specify how that variation is implemented.","section":"Fig. 2 caption and text"},{"comment":"There are several typos and notation issues, e.g., 'caluclations' in the paragraph containing Eq. (4), 'Fig. (4)' for Fig. 4, and 'bT -> 1/Lambda_QCD ~ 5 GeV [-1]' should be '5 GeV^{-1}'.","section":"Throughout"},{"comment":"The lower limit of the first integral is x, but for x > ymax this integral is empty; a brief comment clarifying the intended domain would help the reader.","section":"Eq. (7)"}],"recommendation":"major_revision","confidential_remarks":"Editor-only remarks: The overlap with the authors' own Ref. [14] is considerable; the TMD convolution framework and the pi+n TMD HDF are taken from that preprint, and the new elements here are the ymax prescription and the numerical survey. That is a legitimate development, but the incremental novelty should be stated explicitly in the introduction. I also note that the ymax choice is effectively data-informed, which should be disclosed wherever the word 'prediction' is used."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know before citing this as evidence that chiPT fixes the long-distance TMD shape: the qualitative exponential suppression the authors find is real, but the quantitative claim rests on a power-counting step that is uncontrolled in the regime where they advertise it. The paper is still worth reading and refereeing, just not as the last word.\n\nWhat is actually new: they give the first chiral-EFT calculation of the TMD dbar-ubar asymmetry, and they do it with dimensional regularization rather than the phenomenological cutoffs used in earlier collinear work. The ymax prescription—separating small-y from large-y pion momentum fractions—is a genuinely new power-counting scheme, and it makes the collinear comparison with E866/E906 work without a fitted cutoff. The analytic HDF formulas in Appendix A are clean, and the K0(bT m_pi) falloff is a real mathematical feature of the TMD HDF, not a fit.\n\nSoft spots, in proportion. The ymax value is tuned to the same collinear data used to validate the framework, which is moderate circularity. The delta-function expansion for y > ymax is heuristic, though the power-counting argument in Appendix B gives it some support. More concerning is the step leading to Eq. (13): they replace the pion's valence TMD PDF with its collinear limit, assuming pT ~ Lambda_chi. But at bT ~ 1/m_pi ~ 5 GeV^-1, the relevant expansion parameter is bT Lambda_QCD ~ 1, not Lambda_QCD/Lambda_chi ~ 0.2. The truncation is therefore uncontrolled precisely where the chiral suppression is supposed to dominate. Their own Fig. 4 makes this vivid: putting back a Gaussian with kappa = 200-400 MeV changes the large-bT asymmetry by roughly an order of magnitude. So the claim that the chiral TMD HDF determines the long-distance TMD shape without ad hoc model input is not established.\n\nThe authors are candid about underestimated uncertainties, and they frame the TMD part as exploratory. I read them as aware of the fragility. This is not a fatal flaw, but it means the central quantitative prediction should be viewed as a motivated hypothesis, not a derivation.\n\nFor whom: people working on TMD phenomenology and proton structure will get value from the new HDF expressions and the ymax scheme. The paper deserves a serious referee; a good referee should push hard on the expansion in Eq. (13) and ask for a more explicit justification or a reduced claim. I would engage with it and would cite it, with the caveat above.","headline":"A worthwhile exploratory paper: the K0 suppression at large bT is a real and interesting chiral effect, but the quantitative TMD claim rests on a power-counting step that is uncontrolled at the advertised bT values, so treat the central prediction as provisional.","tokens_in":13803,"tokens_out":1573,"would_cite":true,"duration_ms":16327,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Chiral perturbation theory predicts that the proton's $\\bar{d}-\\bar{u}$ TMD asymmetry falls off exponentially at large transverse separation, with the pion mass setting the decay scale, and does so without fitted shape parameters.","keywords":["transverse momentum dependent PDFs","chiral perturbation theory","dbar-ubar asymmetry","hadronic distribution functions","TMD factorization","pion cloud","exponential suppression","power counting"],"falsifier":"A measurement of the proton's $\\bar{d}-\\bar{u}$ TMD asymmetry at $x\\approx 0.1$–$0.15$ and $b_T$ from 2 to 10 GeV$^{-1}$ that does not follow the exponentially suppressed $\\chi$PT curve of Fig. 3 would falsify the claim that chiral dynamics fixes the long-distance tail. In particular, a falloff that matches the unsuppressed pQCD-only curve, or a Gaussian falloff with width 200–400 MeV (as in the paper's $\\kappa$ model), would indicate that confinement-scale effects beyond $\\chi$PT dominate. Alternatively, a lattice QCD computation of the TMD HDF $\\tilde{f}_{\\pi p}(y,b_T)$ at $b_T\\approx 1/m_\\pi$ that deviates from $K_0(b_T m_\\pi)$ scaling would directly contradict the mechanism.","tokens_in":12720,"feed_emoji":"⚛️","tokens_out":8746,"duration_ms":67850,"temperature":0.7,"pith_summary":"The paper establishes that the long-distance behavior of the proton's transverse momentum dependent parton distributions—specifically the $\\bar{d}-\\bar{u}$ asymmetry—can be computed from chiral perturbation theory without ad hoc model input. The key move is to power-count the pion's longitudinal momentum fraction as small, of order $m_\\pi/M$, so that pions carrying a large fraction of the proton's momentum are excluded from the convolution. With that restriction, the collinear $\\bar{d}-\\bar{u}$ asymmetry agrees with the E866 and E906 data, and the chiral TMD hadronic distribution functions produce a natural exponential suppression at large $b_T$, asymptotically $\\sim K_0(b_T m_\\pi)$. Because the result is regulator-independent, the long-distance tail of the TMD asymmetry is a parameter-free prediction of chiral symmetry, which distinguishes it from the unsuppressed pQCD-only matching.","feed_headline":"Pion cloud sets proton's TMD tail with no free parameters","feed_subtitle":"The d-bar minus u-bar asymmetry falls exponentially at large transverse separation, something perturbative QCD cannot produce.","key_machinery":"The central objects are the hadronic distribution functions (HDFs) $f_{Hp}(y)$ and their transverse-momentum dependent counterparts $f_{Hp}(y,k_T)$, which encode the proton's fluctuations into pion–baryon states. In position space the TMD HDF for the pion satisfies $\\tilde{f}_{\\pi p}(y,b_T)\\sim K_0(b_T m_\\pi)$ at large $b_T$, and this modified Bessel function carries the exponential suppression of the TMD asymmetry. The companion machinery is the $y_\\text{max}$ prescription: the lightcone delta function $\\delta(k^+-yP^+)$ is expanded as $\\delta(yP^+)=0$ for $y>y_\\text{max}\\sim m_\\pi/M$, which restricts the convolution to the region where chiral power counting holds.","core_discovery":"The paper's central claim is that matching TMD PDFs onto chiral perturbation theory yields a cutoff-independent prediction for the $\\bar{d}-\\bar{u}$ asymmetry at large $b_T$, where the chiral hadronic distribution functions fall off as $\\sim K_0(b_T m_\\pi)$. This exponential suppression is generated by the pion-cloud dynamics themselves and is absent in the pQCD-only matching, which shows no large-$b_T$ suppression. In the collinear sector, the paper shows that previously noted disagreements with data in dimensionally regularized $\\chi$PT were due to allowing pions to carry momentum fractions $y \\gg m_\\pi/M$; imposing the $y_\\text{max}$ prescription restores agreement with both E866 and E906 data. The paper further compares three matching schemes and shows that the chiral and pQCD options are distinguishable in the large-$b_T$ region, and that additional confinement-scale physics, modeled as a Gaussian in the pion TMD, would produce a visibly different falloff.","pith_inferences":["If the measured tail falls faster than $K_0(b_T m_\\pi)$—e.g., like a Gaussian with a 200–400 MeV scale—the clean separation between chiral long-distance physics and confinement effects would break down, and the full convolution with pion TMDs (Eq. (9)) would be needed to interpret the data.","The $y_\\text{max}$ prescription is a matching condition rather than a fitted cutoff; a first-principles derivation of its value (for instance from heavy-baryon $\\chi$PT power counting) would shrink the uncertainty band and make the approach fully ab initio.","The same Bessel-function suppression mechanism should appear in other pion-cloud-sensitive observables, such as the $d/u$ ratio at large $x$ or generalized parton distributions at small skewness, offering cross-checks of the chiral TMD framework."],"forward_implications":["The $\\bar{d}-\\bar{u}$ TMD asymmetry is predicted to fall exponentially at $b_T\\gtrsim 2$ GeV$^{-1}$, with the pion mass setting the decay scale; this is a concrete, testable signal for semi-inclusive deep inelastic scattering or Drell-Yan measurements.","Because the $\\chi$PT result has no cutoff or regulator parameter, the same $y_\\text{max}$ prescription that fits the collinear data fixes the TMD prediction, leaving no free shape parameters.","The three matching prescriptions ('$\\chi$PT only', '$\\chi$PT$\\times$pQCD', and 'pQCD only') produce distinct $b_T$ curves, so future data can determine which description is realized.","If Eq. (13) describes the data, the framework extends to other sea-quark TMD asymmetries, such as $s-\\bar{s}$, and to valence TMDs, giving a systematic chiral treatment of TMDs."],"supporting_citations":[{"why":"Establishes the TMD hadronic distribution function convolution formalism that this paper applies to the $\\bar{d}-\\bar{u}$ asymmetry.","marker":"[14]"},{"why":"Derives the collinear convolution formula between valence PDFs and hadronic distribution functions that underlies Eq. (1).","marker":"[15]"},{"why":"Provides the first $\\chi$PT calculation of the collinear $\\bar{d}-\\bar{u}$ asymmetry with a cutoff regulator; this paper redoes it with dimensional regularization and the $y_\\text{max}$ prescription.","marker":"[11]"},{"why":"NuSea E866 measurement of the collinear $\\bar{d}-\\bar{u}$ asymmetry, used as the primary benchmark for the collinear prediction.","marker":"[3]"},{"why":"SeaQuest E906 data, the second benchmark that the collinear prediction must simultaneously describe.","marker":"[34]"},{"why":"Source of the pQCD matching coefficients used for the 'pQCD only' and '$\\chi$PT$\\times$pQCD' TMD matching prescriptions.","marker":"[29]"}],"fun_headline_variants":["Pion cloud predicts exponential drop in proton's sea quark asymmetry","Chiral theory yields exponential TMD asymmetry with no free parameters","Exponential sea asymmetry from pion cloud: chiral prediction","Nonperturbative TMD: pion cloud predicts exponential dbar-ubar falloff","Chiral PT predicts exponential tail for proton's sea asymmetry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction rests on treating the pion's longitudinal momentum fraction $y$ as a small quantity of order $m_\\pi/M$, so that pions with $y$ above about 0.3 are assumed to contribute nothing because the lightcone delta function can be expanded to zero; if real pions with larger momentum fractions contribute measurably, the collinear and TMD predictions lose their foundation.","fun_headline_variants_meta":{"raw":{"variants":["Pion cloud predicts exponential drop in proton's sea quark asymmetry","Chiral theory yields exponential TMD asymmetry with no free parameters","Exponential sea asymmetry from pion cloud: chiral prediction","Nonperturbative TMD: pion cloud predicts exponential dbar-ubar falloff","Chiral PT predicts exponential tail for proton's sea asymmetry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000483,"raw_usage":{"total_tokens":2348,"prompt_tokens":872,"completion_tokens":1476,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":488,"completion_tokens_details":{"reasoning_tokens":1387}},"tokens_in":488,"tokens_out":1476,"duration_ms":10586,"temperature":1.0,"reasoning_tokens":1387,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:34:04.023572+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement of the proton's $\\bar{d}-\\bar{u}$ TMD asymmetry at $x\\approx 0.1$–$0.15$ and $b_T$ from 2 to 10 GeV$^{-1}$ that does not follow the exponentially suppressed $\\chi$PT curve of Fig. 3 would falsify the claim that chiral dynamics fixes the long-distance tail. In particular, a falloff that matches the unsuppressed pQCD-only curve, or a Gaussian falloff with width 200–400 MeV (as in the paper's $\\kappa$ model), would indicate that confinement-scale effects beyond $\\chi$PT dominate. Alternatively, a lattice QCD computation of the TMD HDF $\\tilde{f}_{\\pi p}(y,b_T)$ at $b_T\\approx 1/m_\\pi$ that deviates from $K_0(b_T m_\\pi)$ scaling would directly contradict the mechanism.","supporting_citations":[{"cited_title":"Measurement of flavor asymmetry of light-quark sea in the proton with Drell-Yan dimuon production in $p+p$ and $p+d$ collisions at 120 GeV","cited_arxiv_id":"2212.12160","evidence_quote":"SeaQuest E906 data, the second benchmark that the collinear prediction must simultaneously describe."}],"review_version":1}