{"id":"6ed28c07-58c8-41ca-8835-e6b343748771","arxiv_id":"2412.18324","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Including new COMPASS 2022 deuteron SIDIS data in a TMD global fit markedly reduces the uncertainty on d and dbar quark Sivers and transversity distributions and on the tensor charge.","lead":"Physicists updated a global fit of quark transverse-momentum distributions, adding new COMPASS deuteron and STAR Drell-Yan data. The new data shrink the uncertainties on d and dbar quark Sivers functions, transversity, and the nucleon tensor charge.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Replica generation in Appendix A treats the common scale uncertainty as independent per data point, so the reported precision gain from COMPASS 2024 may be an artifact of the error-propagation scheme.","rationale":"The reader's weakest_assumption focuses on parametrization rigidity and the zeta-prescription evolution. That is a reasonable model-dependence concern, but it is not the most load-bearing issue. The paper's precision claim is quantified entirely by replica uncertainties, and Appendix A implements the correlated scale uncertainty incorrectly: it draws an independent Gaussian for every data point instead of one common shift per experiment. This directly contradicts the covariance matrix in Eq. (2), which the paper itself uses for chi^2 and which contains off-diagonal terms. Because the new COMPASS data are a single high-statistics set with a common target-polarization uncertainty, treating that uncertainty as uncorrelated makes the 38 points look like many independent constraints, artificially inflating the apparent reduction in d/bar-d uncertainties and the tensor-charge error. This is an internal inconsistency, not a matter of external consensus, and it cuts exactly at the central claim. The concrete test of re-running the replica generation with a common scale factor would settle whether the precision improvement survives. I keep the reader's CONDITIONAL verdict because the issue is fixable in principle, but the condition should be sharpened: the replica procedure must be corrected and the precision claim re-evaluated before acceptance. I disagree with the reader's choice of weakest assumption because the more immediate threat to the claim is the error-propagation bug, not the choice of TMD parametrization.","tokens_in":15838,"tokens_out":7005,"duration_ms":67553,"concrete_test":"Regenerate the 1000 replicas using one common Gaussian scale factor per experiment per replica: m_rep^i = (1+lambda_exp)(m_i + delta_i), with lambda_exp ~ N(0, sigma_cor) held constant over all points of that data set and delta_i ~ N(0, sigma_uncor^i), then refit. Compare the 68% bands for h_1^d, h_1^bar-d, f_{1T}^{perpendicular d}, f_{1T}^{perpendicular bar-d}, and g_T with and without COMPASS [44]. If the reduction between blue and red bands largely disappears, the central claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that including the 2024 COMPASS deuteron data 'notably improves' the precision of d/bar-d Sivers and transversity distributions and the tensor charge. The evidence is the shrinkage of the 68% replica bands in Figs. 2-6 and the change g_T = 0.80^{+1.69}_{-0.31} versus 1.73^{+6.38}_{-1.23}. However, the replica generation in Appendix A is internally inconsistent with the covariance matrix in Eq. (2). Equation (2) contains off-diagonal terms (sigma_cor)^2 m_i m_j, representing an overall relative scale uncertainty that should shift all points of a data set together. Yet Eq. (A.4) shifts each point by an independent draw c_rep^i = Random(0, sigma_cor), so the replicas have no off-diagonal correlations. For the 38 new COMPASS points, which share one target-polarization/dilution scale uncertainty, independent shifts largely cancel in the aggregate and artificially narrow the red bands. The apparent reduction in d/bar-d uncertainties and in the tensor-charge error may therefore be a property of the incorrect replica scheme rather than genuine constraining power of the new data.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a global Monte Carlo fit of TMD Sivers functions, transversity distributions, and Collins fragmentation functions using the ζ-prescription TMD evolution. It combines SIDIS data from HERMES, COMPASS (including the 2022 deuteron run), and JLab, Drell-Yan and W/Z data from COMPASS and STAR, and e+e− Collins data from BELLE, BABAR, and BESIII. The central claim is that the new COMPASS deuteron data markedly reduce the uncertainty on the d and anti-d quark Sivers and transversity distributions and on the tensor charge, with g_T changing from 1.73(+6.38/-1.23) without these data to 0.80(+1.69/-0.31) with them. The paper reports χ²/N values near one for most data sets, and the fit and uncertainty estimates follow the parametrization introduced in the authors' previous work [43].","tokens_in":16099,"tokens_out":8735,"duration_ms":82123,"significance":"If the error estimates are reliable, the result is significant: the 2022 COMPASS deuteron data are high-statistics and provide a strong direct handle on sea-quark transverse spin, so a demonstrated improvement in d/anti-d and tensor-charge precision would be a useful step for the TMD program. The global data set is broad, the fitting procedure is standard, and the reported χ²/N values suggest reasonable overall consistency. The paper's quantitative contribution, however, is essentially contained in the error bars: the claimed improvement is the shrinkage of the replica bands in Figs. 2-6 and the g_T interval. Consequently, the reliability of the replica-generation scheme is not a technical detail but the central load-bearing ingredient of the paper.","major_comments":[{"comment":"The replica generation is inconsistent with the stated covariance matrix. Equation (2) contains an off-diagonal term (σ_cor)^2 m_i m_j, which represents a common scale shift applied to all points of a data set. In Eq. (A.4), however, each point receives its own independent draw c_rep^i = Random(0,σ_cor), so the replicas have no common shift and the off-diagonal correlation is lost. For the 38 new COMPASS points, this makes the effective scale uncertainty cancel in the aggregate roughly as σ_cor/sqrt(38), artificially narrowing the red bands in Figs. 2-6 and the reported tensor-charge interval. Since the paper's central claim is precisely this narrowing, the uncertainty estimates as presented are not established. Please regenerate the replicas with a single common c_rep per correlated data set, or sample directly from the covariance matrix in Eq. (2), and re-evaluate all reported error bands and g_T.","section":"Appendix A, Eq. (A.4); §3, Eq. (2)"},{"comment":"The numerical values of the correlated scale uncertainties σ_cor are not reported, and the text does not state whether one common value is used for all SIDIS data sets or whether each experiment has its own value. Because the different data sets have different normalization and dilution-factor uncertainties, and because the central improvement claim depends on how the 2022 COMPASS points are correlated, these inputs must be listed explicitly. Without them, the covariance matrix in Eq. (2) and the replica recipe cannot be reproduced or checked.","section":"§3, Tables 1-3 and Appendix A"},{"comment":"The text says that the functional forms for f_1T^perp, h_1, and H_1^perp are provided in [43] but does not restate them. This matters because the claim of improved precision is conditional on the flexibility of that specific parametrization; a reader cannot tell whether the uncertainty reduction comes from the new data or from model rigidity. Please include the explicit functional forms, the allowed parameter ranges, and the number of free parameters in an appendix, or state more precisely which parameters were re-fitted in the present analysis.","section":"§3"}],"minor_comments":[{"comment":"In the COMPASS 2022 row, clarify whether N=38 refers to the h+ and h- samples combined or to each charge separately; the two reaction lines currently share one N entry.","section":"Table 1"},{"comment":"The caption reads 'compared with with Zeng et al.' and should be corrected to 'compared with Zeng et al.'.","section":"Fig. 10 caption"},{"comment":"The phrase 'the studies are slighted different' should read 'the studies are slightly different'.","section":"§3"},{"comment":"The Sivers χ²/N of 2.26 for the COMPASS proton data [39] is noticeably higher than the other data sets; a brief comment on whether this indicates tension or a model limitation would help the reader assess the fit quality.","section":"Table 1"},{"comment":"The notation is inconsistent: Eq. (2) uses σ_cor. and σ_uncor., while Appendix A uses σ_cor. and σ_uncor. without the same periods; please unify.","section":"Eq. (2) and Appendix A"},{"comment":"The wording 'as recently reported by COMPASS and STAR' is imprecise because the STAR W± data [46] are from 2016, while the STAR Z measurement [47] is from 2024; please specify the measurement years.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"To the editor: The central new result is the claimed reduction in uncertainty. The inconsistency between Eq. (2) and Appendix A means the current error estimates are not trustworthy as presented. This is fixable by rerunning the fit with a common scale shift per data set or by sampling from the stated covariance matrix, and I would not reject solely on this basis. The revised version must contain the corrected uncertainty analysis before the precision claim can be evaluated. I see no citation or scope problem; the paper is a standard continuation of the authors' earlier work."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the inclusion of the 2024 COMPASS high-statistics deuteron SIDIS data in a global TMD fit for Sivers, transversity, and Collins functions. The deuteron target should constrain d and dbar at moderate x, and the comparison of bands with and without the new data shows a clear narrowing. The paper is also honest that the STAR W/Z points add little. The χ²/N values are reasonable and the framework follows the authors' previous work, so the physics direction is not the issue.\n\nThe problem is the error propagation. Equation (2) defines a covariance matrix with off-diagonal correlated terms for a common relative scale uncertainty. Appendix A, however, generates replicas by drawing an independent Gaussian shift per data point for that common scale (c_rep^i). Those two procedures are not equivalent. For the 38 COMPASS deuteron points, which share one target-polarization/dilution scale uncertainty, independent shifts will largely cancel across the set, artificially narrowing the red uncertainty bands in Figs. 2–6. The central quantitative claim—that precision on d and dbar Sivers/transversity and the tensor charge is 'notably improved'—rests on those bands. The improvement may be real, but the quoted numbers (e.g., gT = 0.80^{+1.69}_{-0.31} vs 1.73^{+6.38}_{-1.23}) are likely too optimistic as they stand.\n\nOther soft spots are minor by comparison: the functional forms are deferred to [43], making the fit less self-contained, and no code or tabulated results are provided, so independent verification is hard. The shift in the tensor charge central value from 1.73 to 0.80 is within the old uncertainty but deserves a sentence of discussion.\n\nWho benefits: hadron physicists working on TMD extraction and nucleon spin structure. The qualitative conclusion—new deuteron data pin down sea-quark transverse spin—is probably right, but the error analysis needs to be corrected or defended.\n\nRecommendation: do not desk reject. Send to peer review with a request that the authors fix the replica scheme to match Eq. (2), or show that the independent-draw approximation is numerically irrelevant for these data. The uncertainty bands and tensor-charge errors should then be re-evaluated.","headline":"A useful update of a TMD global fit with new COMPASS deuteron data, but the replica generation understates correlated scale uncertainties, so the headline precision gain may be overstated.","tokens_in":16643,"tokens_out":4646,"would_cite":true,"duration_ms":44748,"reading_group":"yes","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A global TMD fit including the 2022 COMPASS deuteron data narrows the d and dbar quark Sivers and transversity distributions and tightens the nucleon tensor charge.","keywords":["Sivers function","transversity","Collins fragmentation function","TMD factorization","transverse single-spin asymmetry","tensor charge","semi-inclusive deep inelastic scattering","Drell-Yan"],"falsifier":"A refit of the same datasets with a more flexible x-dependence for the d and dbar Sivers and transversity distributions (for example, a neural-network or multi-parameter form): if the 68% uncertainty bands return to their pre-COMPASS width, the reported precision gain is an artifact of the chosen parametrization. A direct experimental check would be a future independent high-statistics measurement of the deuteron Collins and Sivers asymmetries at comparable Q2; if its data fall systematically outside the red error bands of this fit, the extraction is inconsistent.","tokens_in":1925,"feed_emoji":"⚛️","tokens_out":2476,"duration_ms":64459,"temperature":0.7,"pith_summary":"This paper tries to establish that the newest COMPASS measurements of transverse single-spin asymmetries on a polarized deuteron target, combined with existing world data, sharply improve what we know about the transverse spin structure of sea quarks. The analysis works inside transverse-momentum-dependent (TMD) factorization, where the Sivers function, the transversity distribution, and the Collins fragmentation function encode how a quark's transverse motion and spin are correlated with the nucleon spin. Fitting these functions to semi-inclusive deep-inelastic scattering, Drell-Yan, and W/Z production data, the authors report that adding the 2022 COMPASS deuteron data notably reduces the uncertainty on the d and dbar Sivers and transversity distributions. The same data tighten the tensor charge, the integral of the isovector combination of transversity, which matters for searches for new physics beyond the Standard Model. A sympathetic reader would care because this is a concrete step toward a three-dimensional picture of the nucleon's spin structure and a sharper handle on sea-quark spin.","feed_headline":"New COMPASS data sharply cut sea-quark spin-map errors","feed_subtitle":"A global TMD fit narrows the d/dbar Sivers and transversity bands and tightens the tensor charge.","key_machinery":"The machinery is the ζ-prescription for TMD evolution, combined with a fixed $b$-space parametrization of the three nonperturbative functions: the Sivers function $f_{1T}^{\\perp}(x,b)$, the transversity $h_1(x,b)$, and the Collins fragmentation function $H_1^{\\perp}(z,b)$. This parametrization, carried over from the authors' previous study [43], is inserted into the TMD factorization formulas for SIDIS, Drell-Yan, and W/Z production. A $\\chi^2$ with correlated scale uncertainties is minimized, and 1000 replicas of the world data are generated to propagate uncertainties into the extracted functions. The first transverse moments are integrated up to a cutoff $k_T^{\\rm cut} = Q \\times \\delta_{\\rm cut}$ with $\\delta_{\\rm cut} = 1$, and the tensor charge is obtained from the $x$-integrals of the valence combinations of transversity.","core_discovery":"The central claim is that a single global TMD fit, including the high-statistics COMPASS 2022 data [44], yields markedly better-determined Sivers, transversity, and Collins functions for up, down, and sea-quark flavors, with the largest gain in the down and anti-down sectors. The paper reports the isovector tensor charge as $g_T = 0.80^{+1.69}_{-0.31}$ with the new COMPASS data, compared with $1.73^{+6.38}_{-1.23}$ without them. It also states that the fit supports a nonzero sea-quark transversity distribution, a conclusion reached in the authors' earlier work [43]. The Drell-Yan and W/Z data, while included, give only a marginal improvement over what SIDIS already provides.","pith_inferences":["Because the improvement is driven by a single high-statistics deuteron dataset, a re-analysis with a less rigid functional form for the x-dependence would tell whether the quoted compression of the error bands is physical or an artifact of the parametrization.","If the sharper isovector tensor charge survives, it narrows the allowed parameter space for tensor-type interactions in precision low-energy probes of the nucleon.","The observed marginal role of W/Z and Drell-Yan data suggests that testing the predicted sign change of the Sivers function between SIDIS and Drell-Yan will need substantially more DY statistics than exist today.","One could extend the same framework to test flavor dependence by separately fitting the new COMPASS deuteron data and the proton-target data, checking that the combined improvement in d and dbar comes from the deuteron kinematics rather than from a global rescaling."],"forward_implications":["The d and dbar quark Sivers and transversity distributions can now be quoted with substantially smaller uncertainties, which directly improves predictions for future SIDIS and Drell-Yan experiments.","The tensor charge $g_T = 0.80^{+1.69}_{-0.31}$ becomes a sharper benchmark for lattice QCD calculations and for beyond-the-Standard-Model searches that depend on the tensor current.","Because SIDIS data dominate the fit, the marginal impact of current Drell-Yan and W/Z data means that upcoming high-statistics Drell-Yan measurements are the next lever for testing the Sivers sign-change prediction.","The nonzero sea-quark transversity implied by the fit, if confirmed, would show that the proton's transverse spin is carried not only by valence quarks but also by the sea."],"supporting_citations":[{"why":"Supplies the high-statistics 2022 COMPASS deuteron SIDIS data whose inclusion produces the reported uncertainty reduction.","marker":"[44]"},{"why":"Provides the b-space parametrization of Sivers, transversity, and Collins functions that the fit adopts.","marker":"[43]"},{"why":"Supplies the ζ-prescription TMD evolution framework used to evolve the functions to the data scales.","marker":"[52]"},{"why":"HERMES re-binned proton-target SIDIS data provide the high-x constraint in the fit.","marker":"[37]"},{"why":"COMPASS deuteron-target SIDIS data anchor the low-x sea-quark extraction before the 2024 update.","marker":"[38]"},{"why":"COMPASS proton-target data constrain the up-quark sector and help separate flavor dependence.","marker":"[39]"},{"why":"Drell-Yan data from COMPASS test the Sivers function in a process with a time-reversed gauge link.","marker":"[45]"},{"why":"STAR W± data constrain the large-x sea-quark Sivers function in hadron-hadron collisions.","marker":"[46]"},{"why":"STAR Z0 data add a single high-scale point for the Sivers function at much higher Q.","marker":"[47]"}],"fun_headline_variants":["Global TMD fit tightens sea-quark spin distributions","COMPASS data refine d-quark Sivers and transversity","Tensor charge precision boosted by TMD fit with COMPASS","Sea-quark transversity evidence strengthened in global fit","New COMPASS input cuts sea-quark spin-map errors"],"cache_read_input_tokens":18688,"weakest_assumption_plain":"The load-bearing premise is that the functional form chosen for the Sivers, transversity, and Collins functions, together with the ζ-prescription evolution, is flexible enough that the improved precision on d and dbar comes from the new COMPASS data rather than from the model's fixed shape.","fun_headline_variants_meta":{"raw":{"variants":["Global TMD fit tightens sea-quark spin distributions","COMPASS data refine d-quark Sivers and transversity","Tensor charge precision boosted by TMD fit with COMPASS","Sea-quark transversity evidence strengthened in global fit","New COMPASS input cuts sea-quark spin-map errors"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000583,"raw_usage":{"total_tokens":2664,"prompt_tokens":788,"completion_tokens":1876,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":404,"completion_tokens_details":{"reasoning_tokens":1795}},"tokens_in":404,"tokens_out":1876,"duration_ms":15172,"temperature":1.0,"reasoning_tokens":1795,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:48:08.933502+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A refit of the same datasets with a more flexible x-dependence for the d and dbar Sivers and transversity distributions (for example, a neural-network or multi-parameter form): if the 68% uncertainty bands return to their pre-COMPASS width, the reported precision gain is an artifact of the chosen parametrization. A direct experimental check would be a future independent high-statistics measurement of the deuteron Collins and Sivers asymmetries at comparable Q2; if its data fall systematically outside the red error bands of this fit, the extraction is inconsistent.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the high-statistics 2022 COMPASS deuteron SIDIS data whose inclusion produces the reported uncertainty reduction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the b-space parametrization of Sivers, transversity, and Collins functions that the fit adopts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the ζ-prescription TMD evolution framework used to evolve the functions to the data scales."},{"cited_title":"Alekseev, et al., Collins and Sivers asymmetries for pions and kaons in muon-deuteron DIS, Phys.Lett.B673(2009)127--135","cited_arxiv_id":null,"evidence_quote":"COMPASS deuteron-target SIDIS data anchor the low-x sea-quark extraction before the 2024 update."},{"cited_title":"Adolph, et al., Collins and Sivers asymmetries in muonproduc- tion of pions and kaons o ff transversely polarised protons, Phys.Lett","cited_arxiv_id":null,"evidence_quote":"COMPASS proton-target data constrain the up-quark sector and help separate flavor dependence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Drell-Yan data from COMPASS test the Sivers function in a process with a time-reversed gauge link."},{"cited_title":"Adamczyk, et al., Measurement of the transverse single-spin asym- metry in p↑ + p→ W±/Z0 at RHIC, Phys.Rev.Lett.116(13)(2016) 132301","cited_arxiv_id":null,"evidence_quote":"STAR W± data constrain the large-x sea-quark Sivers function in hadron-hadron collisions."},{"cited_title":"Collaboration, Measurements of the Z0/γ∗ cross section and trans- verse single spin asymmetry in 510 GeV p + p collisions, Phys.Lett","cited_arxiv_id":null,"evidence_quote":"STAR Z0 data add a single high-scale point for the Sivers function at much higher Q."}],"review_version":1}