arxiv: 2605.05314 · v1 · submitted 2026-05-06 · ✦ hep-ph · hep-th· nucl-th

Recognition: unknown

First next-to-next-to-leading-order extraction of fragmentation functions for Lambda hyperons

Valerio Bertone , Alessia Bongallino , Amedeo Chiefa , Miguel G. Echevarria , Gunar Schnell

Authors on Pith no claims yet

Pith reviewed 2026-05-08 16:27 UTC · model grok-4.3

classification ✦ hep-ph hep-thnucl-th

keywords fragmentation functionsLambda hyperonsnext-to-next-to-leading orderperturbative QCDsemi-inclusive deep inelastic scatteringsingle-inclusive annihilationneural networks

0 comments

The pith

The first global next-to-next-to-leading-order analysis extracts collinear unpolarised fragmentation functions for Lambda hyperons.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper carries out the first next-to-next-to-leading-order global fit of fragmentation functions that describe how Lambda hyperons form from partons. It combines single-inclusive electron-positron annihilation data with both neutral-current and charged-current semi-inclusive deep-inelastic scattering measurements. Neural-network parametrization and Monte Carlo sampling are used to determine seven independent parton-to-Lambda functions, separating valence-quark contributions for the first time. The resulting set supplies new information on the hadronisation of strange baryons and supplies a reference for later studies.

Core claim

We present MAPFF1.0_Lambda, the first global analysis at next-to-next-to-leading order in perturbative QCD of the collinear unpolarised fragmentation functions of Lambda hyperons. The fit is based on data from single-inclusive electron-positron annihilation, and from both neutral-current and -- for the first time -- charged-current semi-inclusive deep-inelastic scattering. We have adopted a statistical framework based on Monte Carlo sampling and parametrised fragmentation functions in terms of a neural network. The fragmentation function set comprises a total of seven independent parton flavours, allowing for the first independent determination of valence-quark distributions.

What carries the argument

Neural-network parametrization of seven independent parton-to-Lambda fragmentation functions fitted at NNLO in a Monte Carlo statistical framework to SIA and SIDIS data.

If this is right

Offers new insights into the hadronisation mechanism of strange baryons.
Establishes a baseline for future phenomenological and experimental investigations.
Enables the first independent determination of valence-quark distributions to Lambda hyperons.
Incorporates charged-current semi-inclusive deep-inelastic scattering data for the first time.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The set can be used to compute Lambda production cross sections in proton-proton collisions with reduced perturbative uncertainty.
It provides a starting point for global fits that combine hyperon data with other identified-hadron measurements.
Measurements at an electron-ion collider could tighten the valence distributions using the same NNLO framework.

Load-bearing premise

The neural-network parametrization together with the chosen SIA and SIDIS data sets is sufficient to determine seven independent parton-to-Lambda fragmentation functions at NNLO without large biases from missing higher-order terms, data inconsistencies, or parametrization flexibility.

What would settle it

A high-precision measurement of Lambda production rates in charged-current semi-inclusive deep-inelastic scattering that lies outside the uncertainty band of the fit would show that the extracted functions require revision.

read the original abstract

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper gives the first NNLO global extraction of Lambda hyperon fragmentation functions and adds charged-current SIDIS data for the first time.

read the letter

The main takeaway is that they have produced MAPFF1.0_Lambda, the first set of unpolarised Lambda fragmentation functions extracted at next-to-next-to-leading order. They fit a neural-network parametrisation to single-inclusive annihilation data plus both neutral-current and charged-current semi-inclusive DIS, with the charged-current channel appearing here for the first time. The fit treats seven independent parton flavours and uses Monte Carlo sampling to propagate uncertainties, which lets them separate valence-quark distributions more cleanly than earlier NLO work could manage. They include the usual fit-quality diagnostics and cross-checks, so the central values and error bands rest on explicit comparisons rather than just the abstract claim. That combination of higher order, new data, and independent valence extraction is the concrete advance. The methodology follows the standard perturbative framework for fragmentation functions, with NNLO coefficient functions applied consistently and no obvious internal contradictions in the setup. The paper does a reasonable job showing how the new data pulls on the strange-quark sector and how the neural net handles the extra freedom at NNLO. The soft spots are the familiar ones for any global fit of this type. The neural-network flexibility means the output is still shaped by data selection and regularisation choices, and the available Lambda data remain limited, so some flavours stay less constrained than others. There is the usual circularity that the functions are tuned to reproduce the input measurements, but the authors flag this and provide validation plots rather than over-claiming independence. Overall the evidence supports the claim that this is a usable new baseline rather than a revolutionary shift. This work is aimed at people doing phenomenology with strange baryons, whether at colliders or in hadronisation studies. A reader who needs updated Lambda FFs at NNLO accuracy or wants to test the impact of charged-current data will find it directly useful. It deserves a serious referee because the technical extension is clear, the diagnostics are present, and the gap it fills is real.

Referee Report

0 major / 2 minor

Summary. The manuscript presents MAPFF1.0_Lambda, the first global analysis at next-to-next-to-leading order (NNLO) in perturbative QCD of the collinear unpolarised fragmentation functions of Lambda hyperons. The extraction employs a neural-network parametrization with Monte Carlo sampling and is based on single-inclusive electron-positron annihilation data together with both neutral-current and, for the first time, charged-current semi-inclusive deep-inelastic scattering data. This yields seven independent parton-to-Lambda fragmentation functions, permitting the first separate determination of valence-quark distributions.

Significance. If the reported fit quality and uncertainty estimates hold, the result supplies the first NNLO-accurate Lambda fragmentation functions and incorporates new data channels that were previously unavailable. It thereby improves the theoretical baseline for predictions of strange-baryon production in high-energy collisions and offers a concrete reference for future phenomenological studies and experimental analyses.

minor comments (2)

The abstract would benefit from a concise statement of the achieved chi-squared per degree of freedom or the number of data points included, providing immediate context for the robustness of the NNLO extraction.
In the discussion of the neural-network architecture and hyper-parameter choices, a short table summarising the final network configuration and the Monte Carlo replica count would improve reproducibility for readers.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation for minor revision. The report correctly identifies MAPFF1.0_Lambda as the first NNLO global extraction of unpolarised Lambda hyperon fragmentation functions, employing neural-network parametrisation and Monte Carlo sampling on SIA together with both neutral- and charged-current SIDIS data. No specific major comments appear in the report, so we have no point-by-point rebuttals to provide. We remain prepared to incorporate any minor suggestions during revision.

Circularity Check

0 steps flagged

No significant circularity: standard data-driven extraction

full rationale

The manuscript describes a global fit of seven independent Lambda fragmentation functions at NNLO using SIA and SIDIS data, neural-network parametrization, and Monte Carlo sampling. No load-bearing step reduces a claimed prediction or first-principles result to its own inputs by construction. The central output is the fitted set MAPFF1.0_Lambda itself, obtained via explicit minimization against external experimental measurements; coefficient functions, data selection, and uncertainty propagation follow standard perturbative QCD without self-referential closure. Self-citations, if present, are not invoked to justify uniqueness or forbid alternatives. The analysis is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard perturbative QCD assumptions and a data-driven neural-network fit; no new physical entities are introduced.

free parameters (1)

Neural network parameters
The fragmentation functions are represented by a neural network whose weights and biases are fitted to data, introducing a large number of free parameters.

axioms (2)

domain assumption Collinear factorization applies at NNLO for the SIA and SIDIS processes considered
This is the standard assumption that allows fragmentation functions to be extracted from the cited data.
domain assumption The Monte Carlo sampling framework reliably quantifies uncertainties in the global fit
Invoked to handle the statistical treatment of the fit.

pith-pipeline@v0.9.0 · 5444 in / 1454 out tokens · 32558 ms · 2026-05-08T16:27:22.801730+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

43 extracted references · 39 canonical work pages

[1]

Parton distribution and decay functions,

J. C. Collins and D. E. Soper, “Parton distribution and decay functions,”Nucl. Phys. B194 (1982) 445–492

1982
[2]

Collins,Foundations of perturbative QCD, vol

J. Collins,Foundations of perturbative QCD, vol. 32. Cambridge University Press, 11, 2013

2013
[3]

Parton fragmentation functions,

A. Metz and A. Vossen, “Parton fragmentation functions,”Prog. Part. Nucl. Phys.91(2016) 136–202,arXiv:1607.02521 [hep-ex]

work page arXiv 2016
[4]

Factorization of Hard Processes in QCD

J. C. Collins, D. E. Soper, and G. F. Sterman, “Factorization of Hard Processes in QCD,” Adv. Ser. Direct. High Energy Phys.5(1989) 1–91,arXiv:hep-ph/0409313

work page Pith review arXiv 1989
[5]

O(α 2 s) contributions to the longitudinal fragmentation function ine +e− annihilation,

P. J. Rijken and W. L. van Neerven, “O(α 2 s) contributions to the longitudinal fragmentation function ine +e− annihilation,”Phys. Lett. B386(1996) 422–428,arXiv:hep-ph/9604436

work page arXiv 1996
[6]

Higher order QCD corrections to the transverse and longitudinal fragmentation functions in electron–positron annihilation,

P. J. Rijken and W. L. van Neerven, “Higher order QCD corrections to the transverse and longitudinal fragmentation functions in electron–positron annihilation,”Nucl. Phys. B487 (1997) 233–282,arXiv:hep-ph/9609377

work page arXiv 1997
[7]

O(α 2 s) contributions to the asymmetric fragmentation function ine +e− annihilation,

P. J. Rijken and W. L. van Neerven, “O(α 2 s) contributions to the asymmetric fragmentation function ine +e− annihilation,”Phys. Lett. B392(1997) 207–215,arXiv:hep-ph/9609379

work page arXiv 1997
[8]

O(α2 s) timelike Wilson coefficients for parton-fragmentation functions in Mellin space,

J. Bl¨ umlein and V. Ravindran, “O(α2 s) timelike Wilson coefficients for parton-fragmentation functions in Mellin space,”Nucl. Phys. B749(2006) 1–24,arXiv:hep-ph/0604019

work page arXiv 2006
[9]

QCD corrections to semi-inclusive hadron production in electron–positron annihilation at two loops,

A. Mitov and S.-O. Moch, “QCD corrections to semi-inclusive hadron production in electron–positron annihilation at two loops,”Nucl. Phys. B751(2006) 18–52, arXiv:hep-ph/0604160

work page arXiv 2006
[10]

Semi-Inclusive Deep-Inelastic Scattering at Next-to-Next-to-Leading Order in QCD,

L. Bonino, T. Gehrmann, and G. Stagnitto, “Semi-Inclusive Deep-Inelastic Scattering at Next-to-Next-to-Leading Order in QCD,”Phys. Rev. Lett.132no. 25, (2024) 251901, arXiv:2401.16281 [hep-ph]

work page arXiv 2024
[11]

Next-to-Next-to-Leading Order QCD Corrections to Semi-Inclusive Deep-Inelastic Scattering,

S. Goyal, S.-O. Moch, V. Pathak, N. Rana, and V. Ravindran, “Next-to-Next-to-Leading Order QCD Corrections to Semi-Inclusive Deep-Inelastic Scattering,”Phys. Rev. Lett.132 no. 25, (2024) 251902,arXiv:2312.17711 [hep-ph]

work page arXiv 2024
[12]

Neutral and charged current semi-inclusive deep-inelastic scattering at NNLO QCD,

L. Bonino, T. Gehrmann, M. L¨ ochner, K. Sch¨ onwald, and G. Stagnitto, “Neutral and charged current semi-inclusive deep-inelastic scattering at NNLO QCD,”JHEP10(2025) 016,arXiv:2506.19926 [hep-ph]

work page arXiv 2025
[13]

NNLO QCD corrections to unpolarized and polarized SIDIS,

S. Goyal, R. N. Lee, S.-O. Moch, V. Pathak, N. Rana, and V. Ravindran, “NNLO QCD corrections to unpolarized and polarized SIDIS,”Phys. Rev. D111no. 9, (2025) 094007, arXiv:2412.19309 [hep-ph]

work page arXiv 2025
[14]

Czakon, T

M. Czakon, T. Generet, A. Mitov, and R. Poncelet, “Identified Hadron Production at Hadron Colliders in Next-to-Next-to-Leading-Order QCD,”Phys. Rev. Lett.135no. 17, (2025) 17,arXiv:2503.11489 [hep-ph]

work page arXiv 2025
[15]

Fragmentation functions for pions, kaons, and protons at next-to-leading order

B. A. Kniehl, G. Kramer, and B. Potter, “Fragmentation functions for pions, kaons, and protons at next-to-leading order,”Nucl. Phys. B582(2000) 514–536, arXiv:hep-ph/0010289

work page arXiv 2000
[16]

Fragmentation functions from flavor inclusive and flavor taggede +e− annihilations,

S. Kretzer, “Fragmentation functions from flavor inclusive and flavor taggede +e− annihilations,”Phys. Rev. D62(2000) 054001,arXiv:hep-ph/0003177. – 25 –

work page arXiv 2000
[17]

Next-to-leading order determination of fragmentation functions,

L. Bourhis, M. Fontannaz, J. P. Guillet, and M. Werlen, “Next-to-leading order determination of fragmentation functions,”Eur. Phys. J. C19(2001) 89–98, arXiv:hep-ph/0009101

work page arXiv 2001
[18]

Fragmentation functions from semi-inclusive DIS pion production and implications for the polarized parton densities,

S. Kretzer, E. Leader, and E. Christova, “Fragmentation functions from semi-inclusive DIS pion production and implications for the polarized parton densities,”Eur. Phys. J. C22 (2001) 269–276,arXiv:hep-ph/0108055

work page arXiv 2001
[19]

Fragmentation functions for light charged hadrons with complete quark flavor separation,

S. Albino, B. A. Kniehl, and G. Kramer, “Fragmentation functions for light charged hadrons with complete quark flavor separation,”Nucl. Phys. B725(2005) 181–206, arXiv:hep-ph/0502188

work page arXiv 2005
[20]

Resummation of soft gluon logarithms in the DGLAP evolution of fragmentation functions,

S. Albino, B. A. Kniehl, G. Kramer, and W. Ochs, “Resummation of soft gluon logarithms in the DGLAP evolution of fragmentation functions,”Phys. Rev. D73(2006) 054020, arXiv:hep-ph/0510319

work page arXiv 2006
[21]

de Florian, R

D. de Florian, R. Sassot, and M. Stratmann, “Global analysis of fragmentation functions for pions and kaons and their uncertainties,”Phys. Rev. D75(2007) 114010, arXiv:hep-ph/0703242

work page arXiv 2007
[22]

Determination of fragmentation functions and their uncertainties,

M. Hirai, S. Kumano, T. H. Nagai, and K. Sudoh, “Determination of fragmentation functions and their uncertainties,”Phys. Rev. D75(2007) 094009,arXiv:hep-ph/0702250

work page arXiv 2007
[23]

Fragmentation functions at next-to-next-to-leading order accuracy,

D. P. Anderle, F. Ringer, and M. Stratmann, “Fragmentation functions at next-to-next-to-leading order accuracy,”Phys. Rev. D92no. 11, (2015) 114017, arXiv:1510.05845 [hep-ph]

work page arXiv 2015
[24]

Parton-to-pion fragmentation reloaded,

D. de Florian, R. Sassot, M. Epele, R. J. Hern´ andez-Pinto, and M. Stratmann, “Parton-to-pion fragmentation reloaded,”Phys. Rev. D91no. 1, (2015) 014035, arXiv:1410.6027 [hep-ph]

work page arXiv 2015
[25]

First Monte Carlo analysis of fragmentation functions from single-inclusivee +e− annihilation,

N. Sato, J. J. Ethier, W. Melnitchouk, M. Hirai, S. Kumano, and A. Accardi, “First Monte Carlo analysis of fragmentation functions from single-inclusivee +e− annihilation,”Phys. Rev. D94no. 11, (2016) 114004,arXiv:1609.00899 [hep-ph]. [26]NNPDFCollaboration, V. Bertone, S. Carrazza, N. P. Hartland, E. R. Nocera, and J. Rojo, “A determination of the fragme...

work page arXiv 2016
[26]

Parton-to-kaon fragmentation revisited,

D. de Florian, M. Epele, R. J. Hern´ andez-Pinto, R. Sassot, and M. Stratmann, “Parton-to-kaon fragmentation revisited,”Phys. Rev. D95no. 9, (2017) 094019, arXiv:1702.06353 [hep-ph]

work page arXiv 2017
[27]

First QCD analysis of charged hadron fragmentation functions and their uncertainties at next-to-next-to-leading order,

M. Soleymaninia, M. Goharipour, and H. Khanpour, “First QCD analysis of charged hadron fragmentation functions and their uncertainties at next-to-next-to-leading order,”Phys. Rev. D98no. 7, (2018) 074002,arXiv:1805.04847 [hep-ph]. [29]xFITTERDevelopers’ TeamCollaboration, H. Abdolmaleki, M. Soleymaninia, H. Khanpour, S. Amoroso, F. Giuli, A. Glazov, A. Lu...

work page arXiv 2018
[28]

Towards a Global QCD Analysis of Fragmentation Functions at Next-to-Next-to-Leading Order Accuracy,

I. Borsa, R. Sassot, D. de Florian, M. Stratmann, and W. Vogelsang, “Towards a Global QCD Analysis of Fragmentation Functions at Next-to-Next-to-Leading Order Accuracy,” Phys. Rev. Lett.129no. 1, (2022) 012002,arXiv:2202.05060 [hep-ph]

work page arXiv 2022
[29]

Global analysis of fragmentation functions to charged hadrons with high-precision data from the LHC,

J. Gao, C. Liu, X. Shen, H. Xing, and Y. Zhao, “Global analysis of fragmentation functions to charged hadrons with high-precision data from the LHC,”Phys. Rev. D110no. 11, (2024) 114019,arXiv:2407.04422 [hep-ph]

work page arXiv 2024
[30]

Analysis ofπ 0,K 0 S, andηfragmentation functions using single-inclusive annihilation data including BESIII measurements,

M. Li, D. P. Anderle, H. Xing, and Y. Zhao, “Analysis ofπ 0,K 0 S, andηfragmentation functions using single-inclusive annihilation data including BESIII measurements,”Phys. Rev. D111no. 3, (2025) 034030,arXiv:2404.11527 [hep-ph]

work page arXiv 2025
[31]

standard

D. de Florian, M. Stratmann, and W. Vogelsang, “QCD analysis of unpolarized and polarized Λ-baryon production in leading and next-to-leading order,”Phys. Rev. D57(1998) 5811–5824,arXiv:hep-ph/9711387

work page arXiv 1998
[32]

Fragmentation functions forK 0 S and Λ with complete quark flavor separation,

S. Albino, B. A. Kniehl, and G. Kramer, “Fragmentation functions forK 0 S and Λ with complete quark flavor separation,”Nucl. Phys. B734(2006) 50–61,arXiv:hep-ph/0510173

work page arXiv 2006
[33]

AKK update: Improvements from new theoretical input and experimental data,

S. Albino, B. A. Kniehl, and G. Kramer, “AKK update: Improvements from new theoretical input and experimental data,”Nucl. Phys. B803(2008) 42–104,arXiv:0803.2768 [hep-ph]

work page arXiv 2008
[34]

Global analysis of fragmentation functions to light neutral hadrons,

J. Gao, C. Liu, M. Li, X. Shen, H. Xing, Y. Zhao, and Y. Zhou, “Global analysis of fragmentation functions to light neutral hadrons,”Phys. Rev. D112no. 5, (2025) 054045, arXiv:2503.21311 [hep-ph]

work page arXiv 2025
[35]

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–1668,arXiv:1310.1394 [hep-ph]

work page Pith review arXiv 2014
[36]

Bertone, PoSDIS2017, 201 (2018), 1708.00911

V. Bertone, “APFEL++: A new PDF evolution library in C++,”PoSDIS2017(2018) 201, arXiv:1708.00911 [hep-ph]. [42]H1Collaboration, F. D. Aaronet al., “Strangeness production at lowQ 2 in deep-inelasticep scattering at HERA,”Eur. Phys. J. C61(2009) 185–205,arXiv:0810.4036 [hep-ex]. [43]DELPHICollaboration, P. Abreuet al., “Charged and identified particles in ...

work page arXiv 2018
[37]

Λ Production ine +e− Annihilation at 29 GeV,

C. de la Vaissiereet al., “Λ Production ine +e− Annihilation at 29 GeV,”Phys. Rev. Lett. 54(1985) 2071–2074. [Erratum: Phys.Rev.Lett. 55, 263 (1985)]

1985
[38]

Λ Production ine +e− Annihilations at 29 GeV: A Comparison With Lund Model Predictions,

P. S. Baringeret al., “Λ Production ine +e− Annihilations at 29 GeV: A Comparison With Lund Model Predictions,”Phys. Rev. Lett.56(1986) 1346. [55]SLDCollaboration, K. Abeet al., “Production ofπ +,K +,K 0,K ∗0,ϕ,pand Λ 0 in hadronicZ 0 decays,”Phys. Rev. D59(1999) 052001,arXiv:hep-ex/9805029. [56]H1Collaboration, S. Aidet al., “Strangeness production in de...

work page arXiv 1986
[39]

Muoproduction of Neutral Strange Hadrons at 225 GeV,

R. G. Hickset al., “Muoproduction of Neutral Strange Hadrons at 225 GeV,”Phys. Rev. Lett.45(1980) 765–768. [59]E665Collaboration, M. R. Adamset al., “Production of neutral strange particles in muon–nucleon scattering at 490 GeV,”Z. Phys. C61(1994) 539–550. [60]European MuonCollaboration, M. Arneodoet al., “Studies of quark and diquark fragmentation into i...

work page arXiv 1980
[40]

The path to proton structure at 1% accuracy,

G. D’Agostini,Bayesian Reasoning in Data Analysis. World Scientific, Singapore, 2003. [65]NNPDFCollaboration, R. D. Ballet al., “The path to proton structure at 1% accuracy,” Eur. Phys. J. C82no. 5, (2022) 428,arXiv:2109.02653 [hep-ph]

work page arXiv 2003
[41]

Ceres Solver,

S. Agarwal, K. Mierle, and T. C. S. Team, “Ceres Solver,” 10, 2023. https://github.com/ceres-solver/ceres-solver. – 28 –

2023
[42]

On the derivatives of feed-forward neural networks,

R. Abdul Khalek and V. Bertone, “On the derivatives of feed-forward neural networks,” arXiv:2005.07039 [physics.comp-ph]. [68]NNPDFCollaboration, R. D. Ballet al., “Parton distributions from high-precision collider data,”Eur. Phys. J. C77no. 10, (2017) 663,arXiv:1706.00428 [hep-ph]

work page arXiv 2005
[43]

A Variable number flavor scheme for charged current heavy flavor structure functions,

R. S. Thorne and R. G. Roberts, “A Variable number flavor scheme for charged current heavy flavor structure functions,”Eur. Phys. J. C19(2001) 339–349,arXiv:hep-ph/0010344. – 29 –

work page arXiv 2001