pith. machine review for the scientific record. sign in

arxiv: 2604.09071 · v1 · submitted 2026-04-10 · ✦ hep-ph · nucl-th

Recognition: no theorem link

A numerical implementation of the NLO DIS structure functions in the dipole picture

Henri H\"anninen , Heikki M\"antysaari , Jani Penttala
Authors on Pith no claims yet

Pith reviewed 2026-05-10 17:44 UTC · model grok-4.3

classification ✦ hep-ph nucl-th
keywords NLODISdipole picturestructure functionsnumerical implementationmassive quarksimpact factors
0
0 comments X

The pith

Numerical program computes NLO DIS structure functions stably in the dipole picture with massive quarks

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

This paper introduces a numerical implementation for calculating deep inelastic scattering structure functions at next-to-leading order accuracy within the dipole picture. The key advance is rewriting the NLO impact factors involving massive quarks into a form that permits stable numerical evaluation of the cross sections. A reader would care because this provides a practical tool for high-precision predictions in QCD processes at small x, where the dipole model is commonly applied. The work focuses on making theoretical expressions computationally reliable without altering their underlying physics content.

Core claim

We present a numerical program that evaluates deep inelastic scattering (DIS) structure functions at next-to-leading order (NLO) accuracy in the dipole picture. In this numerical implementation the NLO DIS impact factors with massive quarks are written in a form that ensures a stable numerical evaluation of the DIS cross sections.

What carries the argument

The reformulated NLO DIS impact factors with massive quarks, which carry the argument by enabling stable numerical integration while preserving the original physical expressions.

If this is right

  • The code allows for the numerical computation of DIS cross sections at NLO including quark mass effects.
  • It facilitates comparisons between dipole model predictions and experimental data at higher perturbative orders.
  • The method ensures numerical stability for the integrals involved in the structure functions.
  • This implementation can be extended to other observables in the dipole framework at NLO.

Where Pith is reading between the lines

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

  • Such a stable code could be integrated into global fits of parton distributions or dipole amplitudes to extract more accurate parameters from data.
  • It opens the door to NLO calculations of related processes like exclusive vector meson production or diffractive DIS in the same framework.
  • Future work might test the approach against full NLO QCD calculations in other schemes to verify equivalence.

Load-bearing premise

The rewritten form of the NLO impact factors with massive quarks is mathematically equivalent to the original expressions and achieves numerical stability without introducing new approximations.

What would settle it

If running the program on a known test case at NLO yields results that disagree with independent analytic calculations or established Monte Carlo simulations beyond numerical error margins, the implementation's correctness would be questioned.

Figures

Figures reproduced from arXiv: 2604.09071 by Heikki M\"antysaari, Henri H\"anninen, Jani Penttala.

Figure 1
Figure 1. Figure 1: Example diagrams contributing to the elastic [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
read the original abstract

We present a numerical program that evaluates deep inelastic scattering (DIS) structure functions at next-to-leading order (NLO) accuracy in the dipole picture. In this numerical implementation the NLO DIS impact factors with massive quarks are written in a form that ensures a stable numerical evaluation of the DIS cross sections.

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.

Referee Report

1 major / 1 minor

Summary. The manuscript presents a numerical program for evaluating deep inelastic scattering (DIS) structure functions at next-to-leading order (NLO) accuracy in the dipole picture. The central contribution is a reformulation of the NLO DIS impact factors involving massive quarks into an alternative expression that enables stable numerical evaluation of the cross sections, avoiding instabilities from cancellations between real and virtual contributions.

Significance. If the reformulation is algebraically equivalent to the standard expressions and the implementation is verified, the work supplies a practical computational tool for NLO dipole-model calculations of DIS structure functions with heavy quarks. This addresses a known numerical challenge in the field and could support more reliable phenomenological analyses of small-x data, including heavy-flavor contributions relevant to HERA measurements and future EIC studies. The provision of a dedicated numerical code is a concrete asset for reproducibility.

major comments (1)
  1. The central claim rests on the assertion that the rewritten NLO massive-quark impact factors are exactly equivalent to the original expressions while being numerically stable. The manuscript should therefore include explicit numerical cross-checks (e.g., direct comparison of the new and standard forms) in kinematic regions where the original expressions are tractable, such as moderate Q² and x values away from the endpoints. Without such benchmarks, the equivalence cannot be confirmed and remains an unverified assumption.
minor comments (1)
  1. The abstract states that a numerical program is presented but provides no information on code availability, programming language, or installation instructions; these details should be added for usability.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and for the constructive comment. We agree that explicit numerical verification of the equivalence is important and have revised the manuscript accordingly to include the requested benchmarks.

read point-by-point responses

  1. Referee: The central claim rests on the assertion that the rewritten NLO massive-quark impact factors are exactly equivalent to the original expressions while being numerically stable. The manuscript should therefore include explicit numerical cross-checks (e.g., direct comparison of the new and standard forms) in kinematic regions where the original expressions are tractable, such as moderate Q² and x values away from the endpoints. Without such benchmarks, the equivalence cannot be confirmed and remains an unverified assumption.

    Authors: We agree that direct numerical cross-checks are required to substantiate the claim of algebraic equivalence. In the revised manuscript we have added a new subsection (Section 4.1) that presents explicit comparisons between the original and reformulated NLO massive-quark impact factors. These benchmarks are performed at moderate kinematics (Q² = 10–50 GeV², x = 0.01–0.1) where the original expressions remain numerically stable. The two forms agree to better than 0.1 % relative difference, consistent with the Monte-Carlo integration precision. We also include a supplementary figure showing the cancellation pattern in the original expressions and the absence of such cancellations in the new form. These additions directly address the referee’s concern and confirm both equivalence and improved stability. revision: yes

Circularity Check

0 steps flagged

Numerical implementation of existing NLO formulas shows no circularity in derivation chain

full rationale

The paper presents a numerical program implementing pre-existing NLO DIS impact factors (with massive quarks) in the dipole picture, reformulated only for numerical stability. No first-principles derivation, parameter fitting, or predictions are claimed; the work is purely an implementation and algebraic rearrangement of known expressions. No self-citations, self-definitional steps, or reductions of results to inputs by construction appear in the abstract or context. The central claim reduces to providing stable code for established formulas, which is self-contained against external benchmarks and carries no circularity burden.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No information available from the abstract to identify specific free parameters, axioms, or invented entities; the contribution is framed as a numerical implementation of established dipole-picture concepts.

pith-pipeline@v0.9.0 · 5342 in / 1098 out tokens · 38230 ms · 2026-05-10T17:44:49.787412+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Confronting Color Glass Condensate at next-to-leading order with HERA data

    hep-ph 2026-04 unverdicted novelty 6.0

    A Bayesian global fit at full NLO+NLL accuracy extracts the posterior distribution for the non-perturbative initial condition of the NLO Balitsky-Kovchegov equation from HERA inclusive and charm data.

Reference graph

Works this paper leans on

44 extracted references · 41 canonical work pages · cited by 1 Pith paper

  1. [1]

    Abramowicz, et al., Eur

    H. Abramowicz et al., Combination of measurements of inclusive deep inelastic e^ p scattering cross sections and QCD analysis of HERA data , Eur. Phys. J. C 75(12), 580 (2015), doi:10.1140/epjc/s10052-015-3710-4, 1506.06042

    work page doi:10.1140/epjc/s10052-015-3710-4 2015

  2. [2]

    The Color glass condensate and high-energy scattering in QCD

    E. Iancu and R. Venugopalan, The Color glass condensate and high-energy scattering in QCD , pp. 249--3363, World Scientific , doi:10.1142/9789812795533_0005, In Quark-gluon plasma 4 https://doi.org/10.1142/7588 (eds. R.C. Hwa and X.-N. Wang) (2003), hep-ph/0303204

    work page doi:10.1142/9789812795533_0005 2003

  3. [3]

    Garcia-Montero and S

    O. Garcia-Montero and S. Schlichting, Effective theories for nuclei at high energies , Eur. Phys. J. A 61(3), 54 (2025), doi:10.1140/epja/s10050-025-01523-7, 2502.09721

    work page doi:10.1140/epja/s10050-025-01523-7 2025

  4. [4]

    Morreale and F

    A. Morreale and F. Salazar, Mining for Gluon Saturation at Colliders , Universe 7(8), 312 (2021), doi:10.3390/universe7080312, 2108.08254

    work page doi:10.3390/universe7080312 2021

  5. [5]

    Abdul Khalek, et al., Nucl

    R. Abdul Khalek et al., Science Requirements and Detector Concepts for the Electron-Ion Collider : EIC Yellow Report , Nucl. Phys. A 1026, 122447 (2022), doi:10.1016/j.nuclphysa.2022.122447, 2103.05419

    work page doi:10.1016/j.nuclphysa.2022.122447 2022

  6. [6]

    Lappi and H

    T. Lappi and H. M\"antysaari, Next-to-leading order Balitsky-Kovchegov equation with resummation , Phys. Rev. D 93(9), 094004 (2016), doi:10.1103/PhysRevD.93.094004, 1601.06598

    work page doi:10.1103/physrevd.93.094004 2016

  7. [7]

    Balitsky and G

    I. Balitsky and G. A. Chirilli, Next-to-leading order evolution of color dipoles , Phys. Rev. D 77, 014019 (2008), doi:10.1103/PhysRevD.77.014019, 0710.4330

    work page doi:10.1103/physrevd.77.014019 2008

  8. [8]

    Iancu, J

    E. Iancu, J. D. Madrigal, A. H. Mueller, G. Soyez and D. N. Triantafyllopoulos, Resumming double logarithms in the QCD evolution of color dipoles , Phys. Lett. B 744, 293 (2015), doi:10.1016/j.physletb.2015.03.068, 1502.05642

    work page doi:10.1016/j.physletb.2015.03.068 2015

  9. [9]

    Iancu, J

    E. Iancu, J. D. Madrigal, A. H. Mueller, G. Soyez and D. N. Triantafyllopoulos, Collinearly-improved BK evolution meets the HERA data , Phys. Lett. B 750, 643 (2015), doi:10.1016/j.physletb.2015.09.071, 1507.03651

    work page doi:10.1016/j.physletb.2015.09.071 2015

  10. [10]

    A. H. Mueller and B. Patel, Single and double BFKL pomeron exchange and a dipole picture of high-energy hard processes , Nucl. Phys. B 425, 471 (1994), doi:10.1016/0550-3213(94)90284-4, hep-ph/9403256

    work page doi:10.1016/0550-3213(94)90284-4 1994

  11. [11]

    G. Beuf, T. Lappi and R. Paatelainen, Massive quarks in NLO dipole factorization for DIS: Transverse photon , Phys. Rev. D 106(3), 034013 (2022), doi:10.1103/PhysRevD.106.034013, 2204.02486

    work page doi:10.1103/physrevd.106.034013 2022

  12. [12]

    G. Beuf, T. Lappi and R. Paatelainen, Massive quarks in NLO dipole factorization for DIS: Longitudinal photon , Phys. Rev. D 104(5), 056032 (2021), doi:10.1103/PhysRevD.104.056032, 2103.14549

    work page doi:10.1103/physrevd.104.056032 2021

  13. [13]

    G. Beuf, T. Lappi and R. Paatelainen, Massive Quarks at One Loop in the Dipole Picture of Deep Inelastic Scattering , Phys. Rev. Lett. 129(7), 072001 (2022), doi:10.1103/PhysRevLett.129.072001, 2112.03158

    work page doi:10.1103/physrevlett.129.072001 2022

  14. [14]

    Balitsky and G

    I. Balitsky and G. A. Chirilli, Photon impact factor in the next-to-leading order , Phys. Rev. D 83, 031502 (2011), doi:10.1103/PhysRevD.83.031502, 1009.4729

    work page doi:10.1103/physrevd.83.031502 2011

  15. [15]

    Beuf, NLO corrections for the dipole factorization of DIS structure functions at low x , Phys

    G. Beuf, NLO corrections for the dipole factorization of DIS structure functions at low x , Phys. Rev. D 85, 034039 (2012), doi:10.1103/PhysRevD.85.034039, 1112.4501

    work page doi:10.1103/physrevd.85.034039 2012

  16. [16]

    Balitsky and G

    I. Balitsky and G. A. Chirilli, Photon impact factor and k_T -factorization for DIS in the next-to-leading order , Phys. Rev. D 87(1), 014013 (2013), doi:10.1103/PhysRevD.87.014013, 1207.3844

    work page doi:10.1103/physrevd.87.014013 2013

  17. [17]

    Beuf, Dipole factorization for DIS at NLO: Loop correction to the ^*_ T,L q q light-front wave functions , Phys

    G. Beuf, Dipole factorization for DIS at NLO: Loop correction to the ^*_ T,L q q light-front wave functions , Phys. Rev. D 94(5), 054016 (2016), doi:10.1103/PhysRevD.94.054016, 1606.00777

    work page doi:10.1103/physrevd.94.054016 2016

  18. [18]

    H \"a nninen, T

    H. H \"a nninen, T. Lappi and R. Paatelainen, One-loop corrections to light cone wave functions: the dipole picture DIS cross section , Annals Phys. 393, 358 (2018), doi:10.1016/j.aop.2018.04.015, 1711.08207

    work page doi:10.1016/j.aop.2018.04.015 2018

  19. [19]

    Beuf, Dipole factorization for DIS at NLO: Combining the q q and q q g contributions , Phys

    G. Beuf, Dipole factorization for DIS at NLO: Combining the q q and q q g contributions , Phys. Rev. D 96(7), 074033 (2017), doi:10.1103/PhysRevD.96.074033, 1708.06557

    work page doi:10.1103/physrevd.96.074033 2017

  20. [20]

    a nninen, H. M \

    H. H \"a nninen, H. M \"a ntysaari, R. Paatelainen and J. Penttala, Proton Structure Functions at Next-to-Leading Order in the Dipole Picture with Massive Quarks , Phys. Rev. Lett. 130(19), 192301 (2023), doi:10.1103/PhysRevLett.130.192301, 2211.03504

    work page doi:10.1103/physrevlett.130.192301 2023

  21. [21]

    Bertilsson, T

    M. Bertilsson, T. Lappi, H. M \"a ntysaari and X.-B. Tong, The DIS dipole picture cross section in exact kinematics (2026), 2601.07302

    work page arXiv 2026

  22. [22]

    Duclou \'e , H

    B. Duclou \'e , H. H \"a nninen, T. Lappi and Y. Zhu, Deep inelastic scattering in the dipole picture at next-to-leading order , Phys. Rev. D 96(9), 094017 (2017), doi:10.1103/PhysRevD.96.094017, 1708.07328

    work page doi:10.1103/physrevd.96.094017 2017

  23. [23]

    Y. V. Kovchegov, Small x F(2) structure function of a nucleus including multiple pomeron exchanges , Phys. Rev. D 60, 034008 (1999), doi:10.1103/PhysRevD.60.034008, hep-ph/9901281

    work page doi:10.1103/physrevd.60.034008 1999

  24. [24]

    Operator expansion for high-energy scattering

    I. Balitsky, Operator expansion for high-energy scattering , Nucl. Phys. B 463, 99 (1996), doi:10.1016/0550-3213(95)00638-9, hep-ph/9509348

    work page doi:10.1016/0550-3213(95)00638-9 1996

  25. [25]

    Jalilian-Marian, A

    J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, The BFKL equation from the Wilson renormalization group , Nucl. Phys. B 504, 415 (1997), doi:10.1016/S0550-3213(97)00440-9, hep-ph/9701284

    work page doi:10.1016/s0550-3213(97)00440-9 1997

  26. [26]

    Iancu, A

    E. Iancu, A. Leonidov and L. D. McLerran, Nonlinear gluon evolution in the color glass condensate. 1. , Nucl. Phys. A 692, 583 (2001), doi:10.1016/S0375-9474(01)00642-X, hep-ph/0011241

    work page doi:10.1016/s0375-9474(01)00642-x 2001

  27. [27]

    A. H. Mueller, A Simple derivation of the JIMWLK equation , Phys. Lett. B 523, 243 (2001), doi:10.1016/S0370-2693(01)01343-0, hep-ph/0110169

    work page doi:10.1016/s0370-2693(01)01343-0 2001

  28. [28]

    Lappi and H

    T. Lappi and H. M \"a ntysaari, On the running coupling in the JIMWLK equation , Eur. Phys. J. C 73(2), 2307 (2013), doi:10.1140/epjc/s10052-013-2307-z, 1212.4825

    work page doi:10.1140/epjc/s10052-013-2307-z 2013

  29. [29]

    S. Cali, K. Cichy, P. Korcyl, P. Kotko, K. Kutak and C. Marquet, On systematic effects in the numerical solutions of the JIMWLK equation , Eur. Phys. J. C 81(7), 663 (2021), doi:10.1140/epjc/s10052-021-09380-6, 2104.14254

    work page doi:10.1140/epjc/s10052-021-09380-6 2021

  30. [30]

    Beuf, Improving the kinematics for low- x QCD evolution equations in coordinate space , Phys

    G. Beuf, Improving the kinematics for low- x QCD evolution equations in coordinate space , Phys. Rev. D 89(7), 074039 (2014), doi:10.1103/PhysRevD.89.074039, 1401.0313

    work page doi:10.1103/physrevd.89.074039 2014

  31. [31]

    a nninen, T. Lappi and H. M \

    G. Beuf, H. H \"a nninen, T. Lappi and H. M \"a ntysaari, Color Glass Condensate at next-to-leading order meets HERA data , Phys. Rev. D 102, 074028 (2020), doi:10.1103/PhysRevD.102.074028, 2007.01645

    work page doi:10.1103/physrevd.102.074028 2020

  32. [32]

    a nninen and H. M \

    C. Casuga, H. H \"a nninen and H. M \"a ntysaari, Initial condition for the Balitsky-Kovchegov equation at next-to-leading order , Phys. Rev. D 112(3), 034003 (2025), doi:10.1103/54zd-hyvg, 2506.00487

    work page doi:10.1103/54zd-hyvg 2025

  33. [33]

    J. L. Albacete, N. Armesto, J. G. Milhano, P. Quiroga-Arias and C. A. Salgado, AAMQS: A non-linear QCD analysis of new HERA data at small-x including heavy quarks , Eur. Phys. J. C 71, 1705 (2011), doi:10.1140/epjc/s10052-011-1705-3, 1012.4408

    work page doi:10.1140/epjc/s10052-011-1705-3 2011

  34. [34]

    Hänninen, H

    H. Hänninen, H. Mäntysaari and J. Penttala, nlodis , https://github.com/hejajama/nlodis (2026)

    2026

  35. [35]

    Hänninen, H

    H. Hänninen, H. M\"antysaari and J. Penttala, NLODIS: A numerical implementation of the NLO DIS structure functions in the dipole picture , doi:10.5281/zenodo.19367635

    work page doi:10.5281/zenodo.19367635

  36. [36]

    anninen and H. M\

    C. Casuga, H. H\"anninen and H. M\"antysaari, Posterior Samples from Bayesian fit of Initial condition to the BK Equation, repo: camcasuga/bayesian-nlodisfit-dipole: v1.0.0 doi:10.5281/zenodo.15552940

    work page doi:10.5281/zenodo.15552940

  37. [37]

    Galassi, J

    M. Galassi, J. Davies, J. Theiler, B. Gough, G. Jungman, P. Alken, M. Booth, F. Rossi and R. Ulerich, GNU Scientific Library, Free Software Foundation, 2.8, Reference manual (2024)

    2024

  38. [38]

    Hahn, CUBA: A Library for multidimensional numerical integration , Comput

    T. Hahn, CUBA: A Library for multidimensional numerical integration , Comput. Phys. Commun. 168, 78 (2005), doi:10.1016/j.cpc.2005.01.010, hep-ph/0404043

    work page doi:10.1016/j.cpc.2005.01.010 2005

  39. [39]

    Berger and A

    J. Berger and A. Stasto, Numerical solution of the nonlinear evolution equation at small x with impact parameter and beyond the LL approximation , Phys. Rev. D 83, 034015 (2011), doi:10.1103/PhysRevD.83.034015, 1010.0671

    work page doi:10.1103/physrevd.83.034015 2011

  40. [40]

    M \"a ntysaari, J

    H. M \"a ntysaari, J. Penttala, F. Salazar and B. Schenke, Finite-size effects on small-x evolution and saturation in proton and nuclear targets , Phys. Rev. D 111(5), 054033 (2025), doi:10.1103/PhysRevD.111.054033, 2411.13533

    work page doi:10.1103/physrevd.111.054033 2025

  41. [41]

    Lappi and H

    T. Lappi and H. M \"a ntysaari, Single inclusive particle production at high energy from HERA data to proton-nucleus collisions , Phys. Rev. D 88, 114020 (2013), doi:10.1103/PhysRevD.88.114020, 1309.6963

    work page doi:10.1103/physrevd.88.114020 2013

  42. [42]

    Casuga, M

    C. Casuga, M. Karhunen and H. M \"a ntysaari, Inferring the initial condition for the Balitsky-Kovchegov equation , Phys. Rev. D 109(5), 054018 (2024), doi:10.1103/PhysRevD.109.054018, 2311.10491

    work page doi:10.1103/physrevd.109.054018 2024

  43. [43]

    H. A. Bethe, Zur Theorie der Metalle. i. Eigenwerte und Eigenfunktionen der linearen Atomkette , Zeit. f \"u r Phys. 71 , 205 (1931), 10.1007\

    1931

  44. [44]

    Ginsparg, It was twenty years ago today

    P. Ginsparg, It was twenty years ago today... , http://arxiv.org/abs/1108.2700

    work page arXiv