Polarized Deep-Inelastic Scattering with Spin Correlations in Herwig 7
Pith reviewed 2026-06-26 07:44 UTC · model grok-4.3
The pith
Born-level spin-density initialization visibly affects shower-sensitive observables in NLO-matched polarized DIS, while real-emission spin propagation does not for longitudinal cases.
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 next-to-leading-order treatment of polarized DIS in Herwig 7 using the POWHEG matching scheme, including neutral-current γ/Z exchange and charged-current scattering. It constructs a spin-correlation treatment for NLO-matched events by initializing the shower from the spin-density matrix of the polarized DIS Born process and then propagating the spin-density matrix of the accepted real-emission configuration without changing the POWHEG event weight or hardest-emission choice. Validation confirms agreement with fixed-order cross sections. Parton-level comparisons isolate the accepted POWHEG real-emission kinematics. Shower-level observables show that the Born-level spin-de
What carries the argument
spin-density matrix of the polarized DIS Born process, initialized into the shower and then propagated from the accepted real-emission configuration under fixed POWHEG weights
If this is right
- NLO-matched cross sections remain consistent with fixed-order calculations for both neutral-current and charged-current polarized DIS.
- Shower evolution starting from the Born spin-density matrix produces measurable shifts in observables sensitive to parton branching.
- The real-emission spin-density propagation leaves no visible imprint on the longitudinally polarized shower-level distributions studied.
- The treatment preserves the original POWHEG event weights and hardest-emission kinematics while adding spin information.
- Parton-level tests without showering isolate the kinematics of the accepted real-emission configuration.
Where Pith is reading between the lines
- The visible role of Born-level spin suggests that future generators should prioritize accurate Born spin-density initialization for polarized processes.
- If transversely polarized observables or different kinematic regions are examined, real-emission spin effects might become resolvable without changing the matching procedure.
- The same initialization-plus-propagation approach could be tested on other NLO-matched processes to check whether Born spin suffices more generally.
- Maintaining fixed POWHEG weights while adding spin correlations keeps the method modular and avoids reweighting overhead.
Load-bearing premise
The method assumes that initializing the shower from the Born spin-density matrix and propagating the real-emission spin-density matrix without altering the POWHEG weight or hardest-emission choice correctly captures the relevant spin correlations in NLO-matched events.
What would settle it
A statistically significant difference appearing in shower-sensitive observables when the real-emission spin-density matrix is propagated versus when it is omitted would falsify the claim that this additional information is not resolved.
read the original abstract
We investigate polarized deep-inelastic scattering (DIS) in the context of fully exclusive Monte Carlo simulations for the Electron-Ion Collider (EIC). We present a next-to-leading-order (NLO) treatment of polarized DIS in Herwig 7 using the POWHEG matching scheme, including neutral-current $\gamma/Z$ exchange and charged-current scattering. We also construct a spin-correlation treatment for NLO-matched events, first by initializing the shower from the spin-density matrix of the polarized DIS Born process and then by propagating the spin-density matrix of the accepted real-emission configuration without changing the POWHEG event weight or hardest-emission choice. We validate the integrated cross sections against fixed-order calculations, use parton-level comparisons without subsequent shower evolution to isolate the accepted POWHEG real-emission kinematics, and employ shower-level observables to test the numerical importance of the Born-level and real-emission spin information. We find that the Born-level spin-density initialization can have a visible impact on shower-sensitive observables, while the additional spin information from the accepted real-emission configuration is not resolved for the longitudinally polarized observables considered here.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript describes the implementation of an NLO POWHEG-matched treatment of polarized deep-inelastic scattering (neutral-current γ/Z and charged-current) in Herwig 7. It constructs a spin-correlation treatment for NLO-matched events by initializing the parton shower from the Born spin-density matrix and propagating the spin-density matrix of the accepted real-emission configuration without modifying the POWHEG event weight or hardest-emission choice. Integrated cross sections are validated against fixed-order calculations; parton-level comparisons isolate real-emission kinematics; and shower-level observables are used to assess the numerical impact of Born-level versus real-emission spin information. The central finding is that Born-level spin-density initialization visibly affects shower-sensitive observables, while additional real-emission spin information is not resolved for the longitudinally polarized observables considered.
Significance. If the implementation is robust, the work supplies a practical tool for generating fully exclusive polarized DIS events at the EIC that incorporate spin correlations in the shower. The explicit validation of integrated cross sections against independent fixed-order results and the use of parton-level isolation to separate kinematics from shower effects are concrete strengths that support the basic implementation claim.
major comments (1)
- [Description of the NLO spin-correlation treatment and associated numerical studies] The spin-correlation procedure for NLO events (initialization from the Born spin-density matrix followed by propagation of the accepted real-emission spin-density matrix) is performed without altering the POWHEG event weight or the choice of hardest emission. This construction assumes that spin correlations factorize from the NLO matching kinematics selection. If spin-dependent matrix elements would modify the real-emission probability or kinematics inside POWHEG, the reported conclusion that real-emission spin information is not resolved could be an artifact of the enforced decoupling rather than a physical result. The integrated cross-section validation against fixed order does not test this differential consistency.
Simulated Author's Rebuttal
We thank the referee for the careful reading of the manuscript and the positive assessment of the implementation and validation strategy. We address the major comment below.
read point-by-point responses
-
Referee: The spin-correlation procedure for NLO events (initialization from the Born spin-density matrix followed by propagation of the accepted real-emission spin-density matrix) is performed without altering the POWHEG event weight or the choice of hardest emission. This construction assumes that spin correlations factorize from the NLO matching kinematics selection. If spin-dependent matrix elements would modify the real-emission probability or kinematics inside POWHEG, the reported conclusion that real-emission spin information is not resolved could be an artifact of the enforced decoupling rather than a physical result. The integrated cross-section validation against fixed order does not test this differential consistency.
Authors: The POWHEG matching is performed using the polarized matrix elements, which determine both the event weights and the selection of the hardest emission kinematics. The spin-density matrix treatment is applied as an additional layer after the matching step: the shower is initialized from the Born spin-density matrix, and the real-emission spin-density matrix is propagated for the accepted configuration. This does not modify the already-determined POWHEG weights or emission choice. This separation follows the existing structure of spin correlations in the Herwig parton shower for NLO-matched events. The integrated cross sections are validated against fixed-order results using the same polarized matrix elements, and parton-level comparisons (without shower) isolate the real-emission kinematics. We acknowledge that a fully coupled implementation, in which spin information would influence the real-emission probability inside POWHEG, lies beyond the current scope and could in principle affect differential results. However, within the present construction the numerical studies indicate that the additional real-emission spin information has negligible impact on the longitudinally polarized observables considered. We will add an explicit discussion of this approximation and its limitations in the revised manuscript. revision: partial
Circularity Check
No significant circularity; implementation and numerical tests are self-contained
full rationale
The paper describes an NLO POWHEG-matched implementation of polarized DIS in Herwig 7, with spin-density initialization from the Born process and propagation of real-emission spin information without altering weights or kinematics. Validation is performed against independent fixed-order calculations for integrated cross sections, plus parton-level and shower-level comparisons. No derivation reduces a result to its own inputs by construction, no fitted parameters are relabeled as predictions, and no load-bearing uniqueness theorems or ansatzes are imported via self-citation. The reported findings (Born-level spin effects visible, real-emission spin not resolved for longitudinal observables) are numerical outcomes of the described procedure, not tautological.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Standard QCD factorization and parton-shower approximations remain valid when extended to polarized DIS at NLO.
Reference graph
Works this paper leans on
-
[1]
A. Accardi et al.,Electron Ion Collider: The Next QCD Frontier: Understanding the glue that binds us all,Eur. Phys. J. A52(2016) 268 [1212.1701]
Pith/arXiv arXiv 2016
-
[2]
R. Abdul Khalek et al.,Science Requirements and Detector Concepts for the Electron-Ion Collider: EIC Yellow Report,Nucl. Phys. A1026(2022) 122447 [2103.05419]
Pith/arXiv arXiv 2022
-
[3]
Altarelli and G
G. Altarelli and G. Parisi,Asymptotic Freedom in Parton Language,Nucl. Phys. B126 (1977) 298
1977
-
[4]
Bardeen, A.J
W.A. Bardeen, A.J. Buras, D.W. Duke and T. Muta,Deep Inelastic Scattering Beyond the Leading Order in Asymptotically Free Gauge Theories,Phys. Rev. D18(1978) 3998
1978
-
[5]
S. Moch, J.A.M. Vermaseren and A. Vogt,The three-loop splitting functions in QCD: The nonsinglet case,Nucl. Phys. B688(2004) 101 [hep-ph/0403192]
Pith/arXiv arXiv 2004
-
[6]
A. Vogt, S. Moch and J.A.M. Vermaseren,The three-loop splitting functions in QCD: The singlet case,Nucl. Phys. B691(2004) 129 [hep-ph/0404111]
Pith/arXiv arXiv 2004
-
[7]
Zijlstra and W.L
E.B. Zijlstra and W.L. van Neerven,Order-α 2 s corrections to the polarized structure function g1px, Q2q,Nucl. Phys. B417(1994) 61
1994
-
[8]
R. Mertig and W.L. van Neerven,The Calculation of the two loop spin splitting functions P(ij)(1)(x),Z. Phys. C70(1996) 637 [hep-ph/9506451]
Pith/arXiv arXiv 1996
-
[9]
Vogelsang,A Rederivation of the spin-dependent next-to-leading order splitting functions, Phys
W. Vogelsang,A Rederivation of the spin-dependent next-to-leading order splitting functions, Phys. Rev. D54(1996) 2023 [hep-ph/9603366]. – 49 –
Pith/arXiv arXiv 1996
-
[10]
S.-O. Moch, J.A.M. Vermaseren and A. Vogt,The three-loop splitting functions in QCD: The helicity-dependent case,Nucl. Phys. B889(2014) 351 [1409.5131]
Pith/arXiv arXiv 2014
- [11]
-
[12]
S. Catani and M.H. Seymour,A general algorithm for calculating jet cross-sections in NLO QCD,Nucl. Phys. B485(1997) 291 [hep-ph/9605323]
Pith/arXiv arXiv 1997
- [13]
- [14]
-
[15]
Mankiewicz, A
L. Mankiewicz, A. Sch¨ afer and M. Veltri,PEPSI: A Monte Carlo generator for polarized leptoproduction,Comput. Phys. Commun.71(1992) 305
1992
-
[16]
G. Ingelman, A. Edin and J. Rathsman,LEPTO 6.5: A Monte Carlo generator for deep inelastic lepton-nucleon scattering,Comput. Phys. Commun.101(1997) 108 [hep-ph/9605286]
Pith/arXiv arXiv 1997
-
[17]
Charchula, G.A
K. Charchula, G.A. Schuler and H. Spiesberger,Combined QED and QCD radiative effects in deep inelastic lepton-proton scattering: The Monte Carlo generator DJANGO6,Comput. Phys. Commun.81(1994) 381
1994
-
[18]
E.C. Aschenauer, T. Burton, T. Martini, H. Spiesberger and M. Stratmann,Prospects for charged current deep-inelastic scattering off polarized nucleons at a future electron-ion collider,Phys. Rev. D88(2013) 114025 [1309.5327]
Pith/arXiv arXiv 2013
-
[19]
S. Frixione and B.R. Webber,Matching NLO QCD Computations and Parton Shower Simulations,JHEP06(2002) 029 [hep-ph/0204244]
Pith/arXiv arXiv 2002
-
[20]
P. Nason,A New Method for Combining NLO QCD with Shower Monte Carlo Algorithms, JHEP11(2004) 040 [hep-ph/0409146]
Pith/arXiv arXiv 2004
-
[21]
S. Frixione, P. Nason and C. Oleari,Matching NLO QCD Computations with Parton Shower Simulations: the POWHEG method,JHEP11(2007) 070 [0709.2092]
Pith/arXiv arXiv 2007
-
[22]
T. Carli, T. Gehrmann and S. H¨ oche,Hadronic final states in deep-inelastic scattering with Sherpa,Eur. Phys. J. C67(2010) 73 [0912.3715]
Pith/arXiv arXiv 2010
-
[23]
S. H¨ oche, F. Krauss, M. Sch¨ onherr and F. Siegert,Automating the POWHEG method in Sherpa,JHEP04(2011) 024 [1008.5399]
Pith/arXiv arXiv 2011
-
[24]
M. D’Errico and P. Richardson,A Positive-Weight Next-to-Leading-Order Monte Carlo Simulation of Deep Inelastic Scattering and Higgs Boson Production via Vector Boson Fusion in Herwig++,Eur. Phys. J. C72(2012) 2042 [1106.2983]
Pith/arXiv arXiv 2012
-
[25]
S. H¨ oche, S. Kuttimalai and Y. Li,Hadronic final states in DIS at NNLO QCD with parton showers,Phys. Rev. D98(2018) 114013 [1809.04192]
Pith/arXiv arXiv 2018
- [26]
- [27]
-
[28]
I. Borsa and B. J¨ ager,Parton-shower effects in polarized deep inelastic scattering,JHEP (2024) 177 [2404.07702]. – 50 –
arXiv 2024
-
[29]
L. Buonocore, G. Limatola, P. Nason and F. Tramontano,An event generator for Lepton-Hadron deep inelastic scattering at NLO+PS with POWHEG including mass effects, JHEP08(2024) 083 [2406.05115]
arXiv 2024
-
[30]
S. Ferrario Ravasio, R. Gauld, B. J¨ ager, A. Karlberg and G. Zanderighi,An event generator for neutrino-induced deep inelastic scattering and applications to neutrino astronomy,Eur. Phys. J. C85(2025) 888 [2407.03894]
arXiv 2025
-
[31]
P. Meinzinger, D. Reichelt and F. Silvetti,Event generation at MEPS@NLO accuracy in neutral and charged current DIS at the EIC,Phys. Rev. D112(2025) 074039 [2506.08994]
arXiv 2025
-
[32]
I. Helenius, J.O. Laulainen and C.T. Preuss,Multi-jet production in deep inelastic scattering with Pythia,JHEP05(2025) 153 [2410.20950]
arXiv 2025
-
[33]
I. Helenius, J. Laulainen and C.T. Preuss,Multiplicative matching of neutral current deep-inelastic scattering processes at next-to-leading order in PYTHIA 8,2605.00502
-
[34]
Bahr et al.,Herwig++ Physics and Manual,Eur
M. Bahr et al.,Herwig++ Physics and Manual,Eur. Phys. J.C58(2008) 639 [0803.0883]
Pith/arXiv arXiv 2008
-
[35]
Gieseke et al.,Herwig++ 2.5 Release Note,1102.1672
S. Gieseke et al.,Herwig++ 2.5 Release Note,1102.1672
-
[36]
Arnold et al.,Herwig++ 2.6 Release Note,1205.4902
K. Arnold et al.,Herwig++ 2.6 Release Note,1205.4902
-
[37]
Bellm et al.,Herwig++ 2.7 Release Note,1310.6877
J. Bellm et al.,Herwig++ 2.7 Release Note,1310.6877
-
[38]
Bellm et al.,Herwig 7.0/Herwig++ 3.0 release note,Eur
J. Bellm et al.,Herwig 7.0/Herwig++ 3.0 release note,Eur. Phys. J.C76(2016) 196 [1512.01178]
Pith/arXiv arXiv 2016
-
[39]
Bellm et al.,Herwig 7.1 Release Note,1705.06919
J. Bellm et al.,Herwig 7.1 Release Note,1705.06919
-
[40]
Bellm et al.,Herwig 7.2 release note,Eur
J. Bellm et al.,Herwig 7.2 release note,Eur. Phys. J. C80(2020) 452 [1912.06509]
arXiv 2020
-
[41]
Bewick et al.,Herwig 7.3 Release Note,2312.05175
G. Bewick et al.,Herwig 7.3 Release Note,2312.05175
-
[42]
Bellm et al.,The Physics of Herwig 7,2512.16645
J. Bellm et al.,The Physics of Herwig 7,2512.16645
-
[43]
P. Richardson and S. Webster,Spin Correlations in Parton Shower Simulations,Eur. Phys. J. C80(2020) 83 [1807.01955]
Pith/arXiv arXiv 2020
-
[44]
V. Barone, A. Drago and P.G. Ratcliffe,Transverse polarisation of quarks in hadrons,Phys. Rept.359(2002) 1 [hep-ph/0104283]
Pith/arXiv arXiv 2002
-
[45]
Jaffe and A
R.L. Jaffe and A. Manohar,NUCLEAR GLUONOMETRY,Phys. Lett. B223(1989) 218
1989
-
[46]
S. Kumano and Q.-T. Song,Gluon transversity in polarized proton-deuteron Drell-Yan process,Phys. Rev. D101(2020) 054011 [1910.12523]
arXiv 2020
-
[47]
De Causmaecker, R
P. De Causmaecker, R. Gastmans, W. Troost and T.T. Wu,Multiple bremsstrahlung in gauge theories at high energies. 1. General formalism for quantum electrodynamics,Nucl. Phys. B 206(1982) 53
1982
-
[48]
Mertig, M
R. Mertig, M. Boehm and A. Denner,FEYN CALC: Computer algebraic calculation of Feynman amplitudes,Comput. Phys. Commun.64(1991) 345
1991
-
[49]
V. Shtabovenko, R. Mertig and F. Orellana,New Developments in FeynCalc 9.0, 1601.01167
-
[50]
V. Shtabovenko, R. Mertig and F. Orellana,FeynCalc 9.3: New features and improvements, Comput. Phys. Commun.256(2020) 107478 [2001.04407]
Pith/arXiv arXiv 2020
-
[51]
V. Shtabovenko, R. Mertig and F. Orellana,FeynCalc 10: Do multiloop integrals dream of computer codes?,Comput. Phys. Commun.306(2025) 109357 [2312.14089]. – 51 –
arXiv 2025
-
[52]
Seymour,Matrix element corrections to parton shower algorithms,Comput
M.H. Seymour,Matrix element corrections to parton shower algorithms,Comput. Phys. Commun.90(1995) 95 [hep-ph/9410414]
Pith/arXiv arXiv 1995
-
[53]
G. Corcella and M.H. Seymour,Initial state radiation in simulations of vector boson production at hadron colliders,Nucl. Phys. B565(2000) 227 [hep-ph/9908388]
Pith/arXiv arXiv 2000
-
[54]
’t Hooft and M.J.G
G. ’t Hooft and M.J.G. Veltman,Regularization and Renormalization of Gauge Fields,Nucl. Phys. B44(1972) 189
1972
-
[55]
Breitenlohner and D
P. Breitenlohner and D. Maison,Dimensional Renormalization and the Action Principle, Commun. Math. Phys.52(1977) 11
1977
-
[56]
Larin,The Renormalization of the axial anomaly in dimensional regularization,Phys
S.A. Larin,The Renormalization of the axial anomaly in dimensional regularization,Phys. Lett. B303(1993) 113 [hep-ph/9302240]
Pith/arXiv arXiv 1993
-
[57]
Kodaira,QCD higher-order effects in polarized electroproduction: Flavour singlet coefficient functions,Nucl
J. Kodaira,QCD higher-order effects in polarized electroproduction: Flavour singlet coefficient functions,Nucl. Phys. B165(1980) 129
1980
-
[58]
D. de Florian and R. Sassot,O (alpha-s) spin dependent weak structure functions,Phys. Rev. D51(1995) 6052 [hep-ph/9412255]
Pith/arXiv arXiv 1995
-
[59]
C. Bierlich, A. Buckley, J.M. Butterworth, C. Gutschow, L. Lonnblad, T. Procter et al., Robust independent validation of experiment and theory: Rivet version 4 release note, SciPost Phys. Codeb.36(2024) 1 [2404.15984]
arXiv 2024
-
[60]
M. Cacciari, G.P. Salam and G. Soyez,The anti-k t jet clustering algorithm,JHEP04(2008) 063 [0802.1189]
Pith/arXiv arXiv 2008
-
[61]
M. Cacciari, G.P. Salam and G. Soyez,FastJet User Manual,Eur. Phys. J. C72(2012) 1896 [1111.6097]. [62]NNPDFcollaboration,The path to proton structure at 1% accuracy,Eur. Phys. J. C82 (2022) 428 [2109.02653]
Pith/arXiv arXiv 2012
-
[62]
J. Cruz-Martinez, T. Hasenack, F. Hekhorn, G. Magni, E.R. Nocera, T.R. Rabemananjara et al.,NNPDFpol2.0: a global determination of polarised PDFs and their uncertainties at next-to-next-to-leading order,JHEP07(2025) 168 [2503.11814]
arXiv 2025
-
[63]
A. Buckley, J. Ferrando, S. Lloyd, K. Nordstr¨ om, B. Page, M. R¨ ufenacht et al.,LHAPDF6: parton density access in the LHC precision era,Eur. Phys. J. C75(2015) 132 [1412.7420]
Pith/arXiv arXiv 2015
-
[64]
Butterworth et al.,PDF4LHC recommendations for LHC Run II,J
J. Butterworth et al.,PDF4LHC recommendations for LHC Run II,J. Phys. G43(2016) 023001 [1510.03865]
Pith/arXiv arXiv 2016
- [65]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.