Recognition: unknown
TMDs in the Lens of Generative AI: A Pixel-Based Approach to Partonic Imaging
Pith reviewed 2026-05-08 08:09 UTC · model grok-4.3
The pith
A nonparametric pixel-based framework with generative AI solves the TMD inverse problem for unbiased parton imaging.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper presents a novel nonparametric pixel-based framework for the Bayesian inference and imaging of TMD parton distributions. Built on a fully differentiable integration of TMD evolution equations, it enables simultaneous extraction of partonic distributions and the nonperturbative evolution kernel. Efficient sampling of the high-dimensional posterior is achieved via a hybrid normalizing flow-driven Metropolis-Hastings approach. Validation through multi-scale closure tests and SVD analysis reveals null TMDs in the kernel's null space, demonstrating how the synergy of machine learning and multi-scale data removes degeneracies for unbiased 3D partonic imaging.
What carries the argument
The pixel-based discretization of TMDs combined with the hybrid normalizing flow-driven Metropolis-Hastings sampler in a Bayesian context, which allows nonparametric reconstruction and identification of unconstrained null components via singular value decomposition.
If this is right
- The framework permits joint determination of TMDs and the evolution kernel from multi-scale observables.
- Closure tests confirm recovery of input distributions even in convoluted structure function scenarios.
- SVD decomposition exposes functional components of TMDs that produce no observable effects and remain undetermined.
- Degeneracies inherent in the TMD inverse problem are eliminated through this machine learning integration.
Where Pith is reading between the lines
- This method could generalize to other inverse problems in QCD, such as extracting generalized parton distributions from exclusive processes.
- Experiments could be optimized to target the constrained components of TMDs identified by this SVD analysis.
- The approach opens the possibility of incorporating more data sources into the Bayesian pixel framework to further constrain the distributions.
Load-bearing premise
The hybrid normalizing flow-driven Metropolis-Hastings sampler performs efficient and exact sampling of the high-dimensional posterior without introducing biases into the TMD reconstructions or null component identification.
What would settle it
If closure tests with known input TMDs show that the extracted distributions deviate significantly from the inputs or fail to correctly identify the null space components, the framework's ability to solve the inverse problem without bias would be falsified.
read the original abstract
This work introduces a novel, nonparametric pixel-based framework for the Bayesian inference and imaging of transverse momentum dependent (TMD) parton distributions. The methodology is built upon a fully differentiable framework that integrates TMD evolution with the Collins-Soper-Sterman formalism, enabling the simultaneous extraction of partonic distributions and the nonperturbative evolution kernel. To achieve efficient and exact sampling of the high-dimensional posterior, we leverage generative AI through a hybrid normalizing flow-driven Metropolis-Hastings approach. The framework is validated through multi-scale closure tests of increasing complexity, ranging from basic functional models to convoluted structure functions. Using singular value decomposition (SVD), we rigorously characterize the uncertainty of the reconstructed distributions and reveal the existence of null TMDs, which are functional components in the null space of the kernel that remain unconstrained by observables. The new framework provides the first integration of pixel-based discretization, generative AI, and SVD within a Bayesian context to solve the TMD inverse problem. This synergy between machine learning and multi-scale data removes inherent degeneracies and enables unbiased 3D partonic imaging.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces a pixel-based, nonparametric Bayesian framework for inferring transverse-momentum-dependent parton distributions (TMDs). It integrates a fully differentiable TMD evolution model based on the Collins-Soper-Sterman formalism to extract both the distributions and the nonperturbative kernel simultaneously. Posterior sampling of the high-dimensional parameter space is performed via a hybrid normalizing flow-driven Metropolis-Hastings algorithm. The approach is validated with multi-scale closure tests and uses singular value decomposition (SVD) to quantify uncertainties and identify null TMD components lying in the kernel's null space. The central claim is that this combination of pixel discretization, generative AI, and SVD removes inherent degeneracies and enables unbiased 3D partonic imaging.
Significance. If the hybrid sampler can be shown to deliver unbiased sampling and the pixel-based discretization remains truly nonparametric, the framework would be a notable advance for TMD phenomenology. It provides a systematic route to high-dimensional inverse problems in QCD, rigorous uncertainty quantification that includes null-space analysis, and a potential means to break degeneracies with multi-scale data. The explicit use of SVD to expose unconstrained functional components is a strength that could influence future Bayesian extractions.
major comments (3)
- [Methodology section on the hybrid sampler] The claim of unbiased imaging rests on the hybrid normalizing flow-driven Metropolis-Hastings sampler achieving efficient and exact sampling of the high-dimensional posterior (abstract and methodology section describing the generative AI component). Without reported diagnostics (e.g., effective sample size, Gelman-Rubin statistics, or direct bias quantification on reconstructed TMDs), residual approximation errors from the flow can propagate into the SVD null-space components and undermine the degeneracy-removal assertion.
- [SVD analysis and null TMDs section] The SVD analysis that reveals null TMDs is load-bearing for the nonparametric claim. The manuscript must explicitly construct the kernel matrix from the pixel discretization, compute its null space, and demonstrate that the identified null components produce zero contribution to the observables within data precision; otherwise the statement that the method 'removes inherent degeneracies' cannot be verified.
- [Closure tests section] Closure tests of increasing complexity are presented as validation, yet the text lacks details on how pixel resolution was chosen, how error propagation through the differentiable evolution is handled, and what safeguards prevent post-hoc data selection. These choices directly affect whether the extracted TMDs and null components are free of discretization bias.
minor comments (2)
- [Abstract] The abstract asserts a 'fully differentiable framework' but does not display the explicit form of the TMD evolution equations or the parameterization of the nonperturbative kernel; adding one or two key equations would improve immediate readability.
- [Results figures] Figures displaying reconstructed TMDs or structure functions should include full posterior uncertainty bands derived from the sampled ensemble rather than central values alone, to make the SVD-derived uncertainties visually apparent.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. We have revised the manuscript to incorporate additional diagnostics, explicit constructions, and methodological clarifications as requested. Below we respond point by point to the major comments.
read point-by-point responses
-
Referee: [Methodology section on the hybrid sampler] The claim of unbiased imaging rests on the hybrid normalizing flow-driven Metropolis-Hastings sampler achieving efficient and exact sampling of the high-dimensional posterior (abstract and methodology section describing the generative AI component). Without reported diagnostics (e.g., effective sample size, Gelman-Rubin statistics, or direct bias quantification on reconstructed TMDs), residual approximation errors from the flow can propagate into the SVD null-space components and undermine the degeneracy-removal assertion.
Authors: We agree that sampling diagnostics are essential to substantiate the unbiased nature of the posterior. In the revised manuscript we now report effective sample sizes, Gelman-Rubin statistics computed across independent chains, and direct bias quantification obtained by comparing the reconstructed TMDs against the known inputs in the closure tests. These metrics confirm that the hybrid sampler mixes efficiently and that any residual normalizing-flow approximation errors remain well below the level of the data uncertainties, thereby preserving the validity of the degeneracy-removal claim. revision: yes
-
Referee: [SVD analysis and null TMDs section] The SVD analysis that reveals null TMDs is load-bearing for the nonparametric claim. The manuscript must explicitly construct the kernel matrix from the pixel discretization, compute its null space, and demonstrate that the identified null components produce zero contribution to the observables within data precision; otherwise the statement that the method 'removes inherent degeneracies' cannot be verified.
Authors: We have expanded the SVD section to include the explicit construction of the kernel matrix from the pixel-discretized TMD representation, the full singular-value decomposition, and the numerical identification of the null-space basis vectors. Direct evaluation of the observables generated by these null components shows that their contributions lie within numerical zero (consistent with machine precision and data uncertainties), thereby verifying that the method isolates and exposes the unconstrained functional components. revision: yes
-
Referee: [Closure tests section] Closure tests of increasing complexity are presented as validation, yet the text lacks details on how pixel resolution was chosen, how error propagation through the differentiable evolution is handled, and what safeguards prevent post-hoc data selection. These choices directly affect whether the extracted TMDs and null components are free of discretization bias.
Authors: The revised closure-tests section now specifies that pixel resolution was determined through systematic convergence studies in which the number of pixels was increased until the posterior means and variances stabilized to within a chosen tolerance. Error propagation is performed via automatic differentiation through the fully differentiable Collins-Soper-Sterman evolution kernel. To guard against post-hoc selection we employed a pre-defined validation protocol using independent mock data sets generated before any analysis was performed, ensuring that the reported TMDs and null-space components are free of discretization bias. revision: yes
Circularity Check
No significant circularity detected in derivation chain
full rationale
The paper introduces a pixel-based Bayesian framework that integrates TMD evolution via the CSS formalism with a hybrid normalizing flow-driven Metropolis-Hastings sampler for posterior sampling, followed by SVD-based uncertainty characterization. The central claims of degeneracy removal and unbiased 3D imaging rest on the differentiability of the framework, multi-scale closure tests, and explicit identification of null-space components rather than on any self-definitional reduction, fitted input renamed as prediction, or load-bearing self-citation. Simultaneous extraction of distributions and nonperturbative kernel is presented as an enabled capability of the differentiable setup, not as a tautology. No equations or steps in the provided description reduce the output to the inputs by construction.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Collins and D.E
J.C. Collins and D.E. Soper,Back-To-Back Jets in QCD,Nucl. Phys. B193(1981) 381
1981
-
[2]
Collins and D.E
J.C. Collins and D.E. Soper,Parton Distribution and Decay Functions,Nucl. Phys. B194 (1982) 445
1982
-
[3]
Collins, D.E
J.C. Collins, D.E. Soper and G.F. Sterman,Transverse Momentum Distribution in Drell-Yan Pair and W and Z Boson Production,Nucl. Phys. B250(1985) 199
1985
-
[4]
A. Idilbi, X.-d. Ji, J.-P. Ma and F. Yuan,Collins-Soper equation for the energy evolution of transverse-momentum and spin dependent parton distributions,Phys. Rev. D70(2004) 074021 [hep-ph/0406302]
-
[5]
J.C. Collins and A. Metz,Universality of soft and collinear factors in hard-scattering factorization,Phys. Rev. Lett.93(2004) 252001 [hep-ph/0408249]
- [6]
-
[7]
J. Collins,Foundations of Perturbative QCD, vol. 32, Cambridge University Press (2011), 10.1017/9781009401845
-
[8]
M.G. Echevarria, A. Idilbi and I. Scimemi,Factorization Theorem For Drell-Yan At LowqT And Transverse Momentum Distributions On-The-Light-Cone,JHEP07(2012) 002 [1111.4996]
-
[9]
M.G. Echevarria, I. Scimemi and A. Vladimirov,Transverse momentum dependent fragmentation function at next-to–next-to–leading order,Phys. Rev. D93(2016) 011502 [1509.06392]
-
[10]
Piessens,The Hankel Transform, inThe Transforms and Applications Handbook, A.D
R. Piessens,The Hankel Transform, inThe Transforms and Applications Handbook, A.D. Poularikas, ed., (Boca Raton), CRC Press (2000)
2000
-
[11]
Tikhonov and V.Y
A.N. Tikhonov and V.Y. Arsenin,Solutions of Ill-Posed Problems, Winston & Sons, Washington, D.C. (1977)
1977
-
[12]
Groetsch,The Theory of Tikhonov Regularization for Fredholm Equations of the First Kind, Pitman, Boston (1984)
C.W. Groetsch,The Theory of Tikhonov Regularization for Fredholm Equations of the First Kind, Pitman, Boston (1984). – 46 –
1984
-
[13]
I. Scimemi and A. Vladimirov,Non-perturbative structure of semi-inclusive deep-inelastic and Drell-Yan scattering at small transverse momentum,JHEP06(2020) 137 [1912.06532]
- [14]
-
[15]
A. Bacchetta, F. Delcarro, C. Pisano, M. Radici and A. Signori,Extraction of partonic transverse momentum distributions from semi-inclusive deep-inelastic scattering, Drell-Yan and Z-boson production,JHEP06(2017) 081 [1703.10157]
-
[16]
A. Bacchetta, V. Bertone, C. Bissolotti, G. Bozzi, M. Cerutti, F. Piacenza et al.,Unpolarized transverse momentum distributions from a global fit of Drell-Yan and semi-inclusive deep-inelastic scattering data,JHEP10(2022) 127 [2206.07598]
- [17]
- [18]
-
[19]
S. Camarda, G. Ferrera and L. Rossi,Drell-Yan lepton pair production at low invariant masses: transverse-momentum resummation and non-perturbative effects in QCD,JHEP01 (2026) 150 [2508.06201]
-
[20]
A. Bacchetta, V. Bertone, C. Bissolotti, G. Bozzi, M. Cerutti, F. Delcarro et al.,Flavor dependence of unpolarized quark Transverse Momentum Distributions from a global fit,JHEP 08(2024) 232 [2405.13833]
-
[21]
A. Bacchetta, A. Bongallino, M. Cerutti, M. Radici and L. Rossi,Exploring the three-dimensional momentum distribution of longitudinally polarized quarks in the proton, Phys. Rev. Lett.134(2025) 121901 [2409.18078]
- [22]
-
[23]
M. Anselmino, M. Boglione, U. D’Alesio, A. Kotzinian, F. Murgia, A. Prokudin et al., Transversity and Collins functions from SIDIS and e+ e- data,Phys. Rev. D75(2007) 054032 [hep-ph/0701006]
-
[24]
M. Anselmino, M. Boglione, U. D’Alesio, S. Kotzinian, F. Murgia, A. Prokudin et al., Transversity and Collins functions from SIDIS ande+e− data: an update,Phys. Rev. D87 (2013) 094019 [1303.3822]
-
[25]
M. Radici and A. Bacchetta,First Extraction of Transversity from a Global Analysis of Electron-Proton and Proton-Proton Data,Phys. Rev. Lett.120(2018) 192001 [1802.05212]
-
[26]
J. Cammarota, L. Gamberg, Z.B. Kang, J.A. Miller, D. Pitonyak, A. Prokudin et al.,Origin of single transverse-spin asymmetries in high-energy collisions,Phys. Rev. D102(2020) 054002 [2002.08384]
-
[27]
U. D’Alesio, C. Flore and A. Prokudin,Role of the Soffer bound in the extraction of transversity and Collins functions,Phys. Lett. B803(2020) 135347 [2001.01573]
-
[28]
L. Gamberg et al.,Transversity distributions and tensor charges of the nucleon from a global analysis,Phys. Rev. D106(2022) 034014 [2205.00999]. – 47 –
-
[29]
C. Cocuzza, A. Metz, D. Pitonyak, A. Prokudin, N. Sato and R. Seidl,First simultaneous global QCD analysis of dihadron fragmentation functions and transversity parton distribution functions,Phys. Rev. D109(2024) 034024 [2308.14857]
-
[30]
C. Cocuzza, A. Metz, D. Pitonyak, A. Prokudin, N. Sato and R. Seidl,Transversity Distributions and Tensor Charges of the Nucleon: Extraction from Dihadron Production and Their Universal Nature,Phys. Rev. Lett.132(2024) 091901 [2306.12998]
-
[31]
M. Boglione, U. D’Alesio, C. Flore, J.O. Gonzalez-Hernandez, F. Murgia and A. Prokudin, Simultaneous reweighting of Transverse Momentum Dependent distributions,Phys. Lett. B 854(2024) 138712 [2402.12322]
- [32]
- [33]
-
[34]
M. Anselmino, M. Boglione, U. D’Alesio, F. Murgia and A. Prokudin,Study of the sign change of the Sivers function from STAR CollaborationW/Zproduction data,Eur. Phys. J. A53(2017) 45 [1612.06413]
-
[35]
A. Bacchetta, F. Delcarro, C. Pisano and M. Radici,The 3-dimensional distribution of quarks in momentum space,Phys. Lett. B827(2022) 136961 [2004.14278]
-
[36]
M.G. Echevarria, Z.B. Kang and J. Terry,Global analysis of the Sivers functions at NLO+NNLL in QCD,JHEP01(2021) 126 [2009.10710]
-
[37]
M. Boglione, U. D’Alesio, C. Flore and J.O. Gonzalez-Hernandez,Assessing signals of TMD physics in SIDIS azimuthal asymmetries and in the extraction of the Sivers function,JHEP 07(2018) 148 [1806.10645]
-
[38]
M. Cerutti, L. Rossi, S. Venturini, A. Bacchetta, V. Bertone, C. Bissolotti et al.,Extraction of Pion Transverse Momentum Distributions from Drell-Yan data,Phys. Rev. D107(2023) 014014 [2210.01733]
- [39]
-
[40]
M. Boglione and A. Simonelli,Factorization of e+e- to hX cross section, differential in zh, PT and thrust, in the 2-jet limit,JHEP02(2021) 076 [2011.07366]
-
[41]
M. Boglione, J.O. Gonzalez-Hernandez and A. Simonelli,Transverse momentum dependent fragmentation functions from recent BELLE data,Phys. Rev. D106(2022) 074024 [2206.01639]
-
[42]
M. Boglione and A. Simonelli,Full treatment of the thrust distribution in single inclusive e+e− processes,JHEP09(2023) 006 [2306.02937]
-
[43]
M. Anselmino, M. Boglione, U. D’Alesio, J.O. Gonzalez Hernandez, S. Melis, F. Murgia et al.,Collins functions for pions from SIDIS and newe+e− data: A first glance at their transverse momentum dependence,Phys. Rev. D92(2015) 114023 [1510.05389]
- [44]
- [45]
-
[46]
U. D’Alesio, F. Murgia and M. Zaccheddu,First extraction of the Lambda polarizing fragmentation function from Belle e+e- data,Phys. Rev. D102(2020) 054001 [2003.01128]
-
[47]
D. Callos, Z.-B. Kang and J. Terry,Extracting the transverse momentum dependent polarizing fragmentation functions,Phys. Rev. D102(2020) 096007 [2003.04828]
-
[48]
L. Gamberg, Z.-B. Kang, D.Y. Shao, J. Terry and F. Zhao,Transverse Lambda polarization in e+e- collisions,Phys. Lett. B818(2021) 136371 [2102.05553]
-
[49]
U. D’Alesio, L. Gamberg, F. Murgia and M. Zaccheddu,Transverse Lambda polarization in e+e- processes within a TMD factorization approach and the polarizing fragmentation function,JHEP12(2022) 074 [2209.11670]
-
[50]
K.-B. Chen, Z.-T. Liang, Y.-L. Pan, Y.-K. Song and S.-Y. Wei,Isospin Symmetry of Fragmentation Functions,Phys. Lett. B816(2021) 136217 [2102.00658]
-
[51]
I.P. Fernando and D. Keller,Extraction of the Sivers function with deep neural networks, Phys. Rev. D108(2023) 054007 [2304.14328]
-
[52]
A. Bacchetta, V. Bertone, C. Bissolotti, M. Cerutti, M. Radici, S. Rodini et al.,A Neural-Network Extraction of Unpolarised Transverse-Momentum-Dependent Distributions, Phys. Rev. Lett.135(2025) 021904 [2502.04166]
-
[53]
Deep Neural Network extraction of Unpolarized Transverse Momentum Distributions
I.P. Fernando and D. Keller,Deep Neural Network extraction of Unpolarized Transverse Momentum Distributions,2510.17243
work page internal anchor Pith review Pith/arXiv arXiv
-
[54]
Z.-B. Kang, L. Sellers, C. Zhang and C. Zhou,AI-assisted modeling and Bayesian inference of unpolarized quark transverse momentum distributions from Drell-Yan data,2604.14133
work page internal anchor Pith review Pith/arXiv arXiv
-
[55]
Deconvolution problem of deeply virtual Compton scattering,
V. Bertone, H. Dutrieux, C. Mezrag, H. Moutarde and P. Sznajder,Deconvolution problem of deeply virtual Compton scattering,Phys. Rev. D103(2021) 114019 [2104.03836]
- [56]
- [57]
-
[58]
A Higher Order Perturbative Parton Evolution Toolkit (HOPPET)
G.P. Salam and J. Rojo,A Higher Order Perturbative Parton Evolution Toolkit (HOPPET), Comput. Phys. Commun.180(2009) 120 [0804.3755]
work page Pith review arXiv 2009
-
[59]
APFEL: A PDF Evolution Library with QED corrections
V. Bertone, S. Carrazza and J. Rojo,APFEL: A PDF Evolution Library with QED corrections,Comput. Phys. Commun.185(2014) 1647 [1310.1394]
work page Pith review arXiv 2014
-
[60]
V. Bertone, R.F.d. Castillo, M.G. Echevarria, Ó. del Río and S. Rodini,One-loop evolution of twist-2 generalized parton distributions,Phys. Rev. D109(2024) 034023 [2311.13941]
-
[61]
L. Dinh, J. Sohl-Dickstein and S. Bengio,Density estimation using Real NVP,ICLR 2017 (2017) [1605.08803]
work page internal anchor Pith review arXiv 2017
-
[62]
R.C. Aster, B. Borchers and C.H. Thurber,Parameter Estimation and Inverse Problems, Elsevier / Academic Press, Amsterdam, 3rd ed. (2018), 10.1016/C2015-0-04300-3
-
[63]
P.C. Barry et al.,First simultaneous analysis of transverse momentum dependent and collinear parton distributions in the proton,2510.13771. – 49 –
- [64]
-
[65]
S.M. Aybat and T.C. Rogers,TMD Parton Distributions and TMD Factorization at Higher Order,Phys. Rev. D83(2011) 114042 [1101.5057]
-
[66]
Semi-inclusive deep inelastic scattering at small transverse momentum
A. Bacchetta, M. Diehl, K. Goeke, A. Metz, P.J. Mulders and M. Schlegel,Semi-inclusive deep inelastic scattering at small transverse momentum,JHEP02(2007) 093 [hep-ph/0611265]. – 50 –
work page Pith review arXiv 2007
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.