In-medium QCD splittings beyond the soft, large-N_c and harmonic-oscillator approximations all at once
Pith reviewed 2026-06-25 19:14 UTC · model grok-4.3
The pith
The BDMPS-Z equations for in-medium QCD radiation have been solved numerically in full generality for the first time.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The BDMPS-Z equations admit a stable, cutoff-free numerical solution that yields in-medium splitting functions including finite-energy effects, subleading-color contributions, and a realistic parton-medium interaction model; these functions deviate substantially from those obtained under the soft, large-Nc, and harmonic-oscillator approximations.
What carries the argument
The BDMPS-Z equations for the fully differential in-medium splitting probability, solved numerically across all phase space.
If this is right
- Jet quenching calculations that feed the new splitting functions into medium-modified parton showers will produce different predictions for observables such as jet shapes, fragmentation functions, and R_AA ratios.
- Constraints on the jet quenching parameter extracted from heavy-ion data will shift once the approximate splitting functions are replaced.
- The size of the reported deviations supplies a quantitative error budget for every earlier study that relied on the soft or large-Nc limits.
- The numerical framework can be extended to compute higher-order corrections or to treat more complex medium geometries without additional analytic approximations.
Where Pith is reading between the lines
- The same numerical engine could be applied to related evolution equations that appear in other transport problems, such as QED cascades or neutrino propagation in dense matter.
- Once the splitting functions are tabulated, they can be interfaced with existing event generators to produce updated Monte-Carlo predictions that can be confronted with upcoming LHC and RHIC datasets.
- The quantified deviations suggest that analytic resummations valid only in limited corners of phase space may need to be abandoned for precision work.
Load-bearing premise
The numerical method converges and remains stable for all relevant kinematics without hidden cutoffs or uncontrolled approximations.
What would settle it
A direct comparison between the new splitting functions and an independent analytic or Monte-Carlo implementation of the same BDMPS-Z equations in a controlled kinematic window where the approximations are known to fail.
Figures
read the original abstract
Nearly thirty years ago, Baier, Dokshitzer, Mueller, Peign\'e, Schiff, and Zakharov (BDMPS-Z) introduced a formalism to calculate the fully differential probability for a high-energy quark or gluon to radiate inside a finite-volume QCD plasma. We report on the first, complete numerical solution to the BDMPS-Z equations for in-medium QCD splittings. Our numerical routines are precise across phase-space, enabling a determination of the in-medium splitting functions that is significantly beyond the state-of-the-art, including finite-energy effects, subleading-color contributions, and a realistic model for parton-medium interactions. We quantify the uncertainties associated with standard approximations in the literature, revealing substantial deviations across phase-space. This work opens a path toward more precise calculations of jet observables and for powerful new constraints of medium parameters from high-energy heavy-ion collider data.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims to deliver the first complete numerical solution of the BDMPS-Z integro-differential equations for in-medium QCD splittings. It asserts that the routines achieve precision across phase space, incorporate finite-energy effects, subleading-color contributions, and a realistic parton-medium interaction model, and that they reveal substantial deviations from the soft, large-Nc, and harmonic-oscillator approximations used in the literature.
Significance. If the numerical implementation is shown to be stable and accurate, the work would constitute a genuine technical advance by removing three long-standing approximations simultaneously. This could tighten theoretical predictions for medium-induced radiation and thereby improve the extraction of medium parameters from jet observables at the LHC and RHIC.
major comments (1)
- [Abstract, §3] Abstract and §3 (numerical implementation): the central claim that the routines are 'precise across phase-space' and constitute a 'complete' solution is not accompanied by any convergence tests, grid-independence studies, residual-error metrics, or explicit comparisons against the known analytic limits (soft-gluon, large-Nc, harmonic-oscillator). Without these, the assertion that the results are 'significantly beyond the state-of-the-art' cannot be evaluated.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for highlighting the need for explicit validation of the numerical implementation. We address the major comment below.
read point-by-point responses
-
Referee: [Abstract, §3] Abstract and §3 (numerical implementation): the central claim that the routines are 'precise across phase-space' and constitute a 'complete' solution is not accompanied by any convergence tests, grid-independence studies, residual-error metrics, or explicit comparisons against the known analytic limits (soft-gluon, large-Nc, harmonic-oscillator). Without these, the assertion that the results are 'significantly beyond the state-of-the-art' cannot be evaluated.
Authors: We agree that the current manuscript does not present explicit convergence tests, grid-independence studies, residual-error metrics, or direct comparisons to the analytic limits. To substantiate the claims of precision and completeness, we will revise §3 to include: (i) results from multiple grid resolutions demonstrating numerical convergence, (ii) quantitative residual-error estimates across phase space, and (iii) side-by-side comparisons of the numerical solutions against the known soft-gluon, large-Nc, and harmonic-oscillator analytic limits. These additions will allow the reader to evaluate the accuracy of the implementation. revision: yes
Circularity Check
No circularity: numerical solution of established external equations
full rationale
The paper reports a numerical implementation of the BDMPS-Z equations, which originate from prior independent literature (Baier et al.). The claimed results follow directly from solving these integro-differential equations across phase space, with no evidence of self-definitional loops, fitted inputs renamed as predictions, or load-bearing self-citations that reduce the central claim to the paper's own inputs. The derivation chain is self-contained against the external BDMPS-Z benchmark and does not exhibit any of the enumerated circularity patterns.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption BDMPS-Z formalism correctly describes in-medium parton splittings
Reference graph
Works this paper leans on
-
[1]
L. Apolin´ ario, Y.-J. Lee, and M. Winn, Prog. Part. Nucl. Phys.127, 103990 (2022), arXiv:2203.16352 [hep-ph]
arXiv 2022
- [2]
-
[3]
Connors, C
M. Connors, C. Nattrass, R. Reed, and S. Salur, Rev. Mod. Phys.90, 025005 (2018)
2018
-
[4]
L. Cunqueiro and A. M. Sickles, Prog. Part. Nucl. Phys. 124, 103940 (2022), arXiv:2110.14490 [nucl-ex]
arXiv 2022
-
[5]
S. Acharyaet al.(ALICE), Eur. Phys. J. C84, 813 (2024), arXiv:2211.04384 [nucl-ex]
arXiv 2024
-
[6]
A. Hayrapetyanet al.(CMS), Phys. Rept.1115, 219 (2025), arXiv:2405.10785 [nucl-ex]
arXiv 2025
-
[7]
Z. Citronet al., CERN Yellow Rep. Monogr.7, 1159 (2019), arXiv:1812.06772 [hep-ph]
Pith/arXiv arXiv 2019
-
[8]
Mehtar-Tani, (2025), arXiv:2509.26394 [hep-ph]
Y. Mehtar-Tani, (2025), arXiv:2509.26394 [hep-ph]
arXiv 2025
-
[9]
R. Baier, Y. L. Dokshitzer, S. Peigne, and D. Schiff, Phys. Lett. B345, 277 (1995), arXiv:hep-ph/9411409
Pith/arXiv arXiv 1995
- [10]
- [11]
-
[12]
B. G. Zakharov, JETP Lett.63, 952 (1996), arXiv:hep- ph/9607440
arXiv 1996
-
[13]
B. G. Zakharov, JETP Lett.65, 615 (1997), arXiv:hep- ph/9704255
arXiv 1997
-
[14]
J.-P. Blaizot, F. Dominguez, E. Iancu, and Y. Mehtar- Tani, JHEP01, 143, arXiv:1209.4585 [hep-ph]
-
[15]
L. Apolinario, N. Armesto, and C. A. Salgado, Phys. Lett. B718, 160 (2012), arXiv:1204.2929 [hep-ph]
Pith/arXiv arXiv 2012
-
[16]
M. Gyulassy and X.-n. Wang, Nucl. Phys. B420, 583 (1994), arXiv:nucl-th/9306003
Pith/arXiv arXiv 1994
-
[17]
M. Gyulassy, P. Levai, and I. Vitev, Phys. Rev. D66, 014005 (2002), arXiv:nucl-th/0201078
Pith/arXiv arXiv 2002
-
[18]
P. B. Arnold, Phys. Rev. D80, 025004 (2009), arXiv:0903.1081 [nucl-th]
Pith/arXiv arXiv 2009
-
[19]
N. Armesto, C. A. Salgado, and U. A. Wiedemann, Phys. Rev. D69, 114003 (2004), arXiv:hep-ph/0312106
Pith/arXiv arXiv 2004
-
[20]
U. A. Wiedemann, Nucl. Phys. B588, 303 (2000), arXiv:hep-ph/0005129
Pith/arXiv arXiv 2000
-
[21]
M. Gyulassy, P. Levai, and I. Vitev, Phys. Rev. Lett.85, 5535 (2000), arXiv:nucl-th/0005032
Pith/arXiv arXiv 2000
-
[22]
M. Djordjevic and M. Gyulassy, Nucl. Phys. A733, 265 (2004), arXiv:nucl-th/0310076
Pith/arXiv arXiv 2004
-
[23]
G. Ovanesyan and I. Vitev, Phys. Lett. B706, 371 (2012), arXiv:1109.5619 [hep-ph]
Pith/arXiv arXiv 2012
-
[24]
P. B. Arnold, G. D. Moore, and L. G. Yaffe, JHEP06, 030, arXiv:hep-ph/0204343
-
[25]
S. Jeon and G. D. Moore, Phys. Rev. C71, 034901 (2005), arXiv:hep-ph/0309332
Pith/arXiv arXiv 2005
-
[26]
L. Apolin´ ario, N. Armesto, J. G. Milhano, and C. A. Salgado, JHEP02, 119, arXiv:1407.0599 [hep-ph]
-
[27]
F. Dom´ ınguez, J. G. Milhano, C. A. Salgado, K. Ty- woniuk, and V. Vila, Eur. Phys. J. C80, 11 (2020), arXiv:1907.03653 [hep-ph]
arXiv 2020
-
[28]
J. H. Isaksen and K. Tywoniuk, JHEP21, 125, arXiv:2107.02542 [hep-ph]
-
[29]
M. D. Sievert and I. Vitev, Phys. Rev. D98, 094010 (2018), arXiv:1807.03799 [hep-ph]
Pith/arXiv arXiv 2018
- [30]
-
[31]
J. Barata, M. Li, W. Qian, C. A. Salgado, and J. M. Silva, (2026), arXiv:2604.11616 [hep-ph]
Pith/arXiv arXiv 2026
-
[32]
X. Feal and R. Vazquez, Phys. Rev. D98, 074029 (2018), arXiv:1811.01591 [hep-ph]
Pith/arXiv arXiv 2018
-
[33]
X. Feal, C. A. Salgado, and R. A. Vazquez, Phys. Lett. B816, 136251 (2021), arXiv:1911.01309 [hep-ph]
arXiv 2021
-
[34]
Mehtar-Tani, JHEP07, 057, arXiv:1903.00506 [hep- ph]
Y. Mehtar-Tani, JHEP07, 057, arXiv:1903.00506 [hep- ph]
arXiv 1903
- [35]
- [36]
-
[37]
J. H. Isaksen, A. Takacs, and K. Tywoniuk, JHEP02, 156, arXiv:2206.02811 [hep-ph]
- [38]
-
[39]
M. V. Kuzmin, X. Mayo L´ opez, J. Reiten, and A. V. Sadofyev, Phys. Rev. D109, 014036 (2024), arXiv:2309.00683 [hep-ph]
arXiv 2024
-
[40]
M. V. Kuzmin and X. Mayo L´ opez, (2024), arXiv:2406.14628 [hep-ph]
arXiv 2024
-
[41]
S. Caron-Huot and C. Gale, Phys. Rev. C82, 064902 (2010), arXiv:1006.2379 [hep-ph]
Pith/arXiv arXiv 2010
-
[42]
A. Kovner and U. A. Wiedemann, , 192 (2003), arXiv:hep-ph/0304151
Pith/arXiv arXiv 2003
-
[43]
Wang and M
X.-N. Wang and M. Gyulassy, Phys. Rev. Lett.68, 1480 (1992)
1992
- [44]
-
[45]
J. H. Isaksen and K. Tywoniuk, JHEP09, 049, arXiv:2303.12119 [hep-ph]
-
[46]
A. G. Borisov and S. V. Shabanov, Journal of Computa- tional Physics216, 391 (2006)
2006
-
[47]
R. Diogo Soares and M. Schir` o, SciPost Physics17, 10.21468/scipostphys.17.5.128 (2024)
-
[48]
Andersson, G
B. Andersson, G. Gustafson, L. Lonnblad, and U. Pet- tersson, Z. Phys. C43, 625 (1989)
1989
-
[49]
Y. Mehtar-Tani and K. Tywoniuk, JHEP04, 125, arXiv:1610.08930 [hep-ph]
-
[50]
Y.-T. Chien and I. Vitev, Phys. Rev. Lett.119, 112301 (2017), arXiv:1608.07283 [hep-ph]
Pith/arXiv arXiv 2017
-
[51]
G. Milhano, U. A. Wiedemann, and K. C. Zapp, Phys. Lett. B779, 409 (2018), arXiv:1707.04142 [hep-ph]
Pith/arXiv arXiv 2018
- [52]
-
[53]
J. Casalderrey-Solana, G. Milhano, D. Pablos, and K. Rajagopal, JHEP01, 044, arXiv:1907.11248 [hep-ph]
arXiv 1907
- [54]
-
[55]
L. Wang, J.-W. Kang, Q. Zhang, S. Shen, W. Dai, B.- W. Zhang, and E. Wang, Chin. Phys. Lett.40, 032101 (2023), arXiv:2211.13674 [nucl-th]
arXiv 2023
-
[56]
Tachibanaet al.(JETSCAPE), Phys
Y. Tachibanaet al.(JETSCAPE), Phys. Rev. C110, 044907 (2024), arXiv:2301.02485 [hep-ph]
arXiv 2024
-
[57]
A. S. Kudinoor, D. Pablos, and K. Rajagopal, (2025), arXiv:2509.08881 [hep-ph]
arXiv 2025
-
[58]
L. Apolin´ ario, D. Costa, and A. Soto-Ontoso, (2026), arXiv:2601.13310 [hep-ph]
arXiv 2026
- [59]
- [60]
- [61]
-
[62]
Z. Yang, Y. He, I. Moult, and X.-N. Wang, Phys. Rev. Lett.132, 011901 (2024), arXiv:2310.01500 [hep-ph]
arXiv 2024
-
[63]
Y. Fu, B. M¨ uller, and C. Sirimanna, Phys. Rev. Lett. 135, 112302 (2025), arXiv:2411.04866 [nucl-th]
arXiv 2025
-
[64]
W.-J. Xing, S. Cao, G.-Y. Qin, and X.-N. Wang, Phys. Rev. Lett.134, 052301 (2025), arXiv:2409.12843 [hep- ph]
arXiv 2025
- [65]
- [66]
- [67]
- [68]
- [69]
- [70]
-
[71]
L. Apolin´ ario, R. Kunnawalkam Elayavalli, N. O. Madureira, J.-X. Sheng, X.-N. Wang, and Z. Yang, Phys. Rev. D112, 054018 (2025), arXiv:2502.11406 [hep-ph]
arXiv 2025
-
[72]
J. Barata, M. V. Kuzmin, I. Moult, A. V. Sadofyev, and J. M. Silva, (2026), arXiv:2604.21971 [hep-ph]
Pith/arXiv arXiv 2026
-
[73]
K. C. Zapp, F. Krauss, and U. A. Wiedemann, JHEP03, 080, arXiv:1212.1599 [hep-ph]
-
[74]
J. Casalderrey-Solana, D. C. Gulhan, J. G. Milhano, D. Pablos, and K. Rajagopal, JHEP10, 019, [Erratum: JHEP 09, 175 (2015)], arXiv:1405.3864 [hep-ph]
Pith/arXiv arXiv 2015
-
[75]
S. Caoet al.(JETSCAPE), Phys. Rev. C96, 024909 (2017), arXiv:1705.00050 [nucl-th]
Pith/arXiv arXiv 2017
-
[76]
P. Caucal, E. Iancu, A. H. Mueller, and G. Soyez, Phys. Rev. Lett.120, 232001 (2018), arXiv:1801.09703 [hep- ph]
Pith/arXiv arXiv 2018
-
[77]
R. Belmontet al., Nucl. Phys. A1043, 122821 (2024), arXiv:2305.15491 [nucl-ex]
arXiv 2024
-
[78]
S. Schlichting and I. Soudi, Phys. Rev. D105, 076002 (2022), arXiv:2111.13731 [hep-ph]
arXiv 2022
-
[79]
P. Aurenche, F. Gelis, and H. Zaraket, JHEP05, 043, arXiv:hep-ph/0204146
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.