Recognition: 2 theorem links
· Lean TheoremUnified Extraction of In-Medium Heavy Quark Potentials from RHIC to LHC Energies via Deep Learning
Pith reviewed 2026-05-10 17:19 UTC · model grok-4.3
The pith
Deep learning applied to bottomonium data from RHIC and LHC extracts an in-medium heavy quark potential whose real part stays near the vacuum Cornell form while the imaginary part drives most suppression.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that a simultaneous Bayesian extraction from bottomonium R_AA measured in Pb-Pb collisions at 5.02 TeV and 2.76 TeV and in Au-Au collisions at 200 GeV yields an in-medium heavy quark potential whose real part remains close to the vacuum Cornell form, implying a relatively weak screened Debye mass, while the imaginary part is more tightly constrained by the data and supplies the dominant contribution to suppression over the full energy range.
What carries the argument
The in-medium heavy quark potential V(T,r) with separate real and imaginary parts, inserted into a time-dependent Schrödinger equation for b b-bar dipole evolution inside a hydrodynamically expanding medium; parameters are learned by convolutional neural networks trained on simulated R_AA and inverted via stochastic gradient Langevin dynamics to obtain posterior distributions.
If this is right
- The real part's closeness to vacuum implies that color screening remains modest at the temperatures reached in these collisions.
- Dissociation driven by the imaginary part, rather than reduction in binding energy, is the main mechanism for the observed suppression.
- A single potential parametrization works across the density range spanned by RHIC and LHC data.
- The method provides a unified, data-driven route to the potential instead of separate fits at each energy.
- The extracted imaginary part can be directly inserted into calculations of other quarkonium states or heavy-flavor transport.
Where Pith is reading between the lines
- Models of in-medium quarkonium should emphasize medium-induced width over adjustments to the real binding potential.
- The weak screening result could guide targeted lattice QCD studies of the real part at finite temperature.
- The same extraction pipeline can be applied to charmonium data or to future higher-precision measurements to test consistency.
- If the hydrodynamic background is varied, the stability of the extracted imaginary part would indicate how much the conclusion depends on the medium evolution model.
Load-bearing premise
The true in-medium potential can be captured by the chosen functional form and that the hydrodynamic background used for the evolution is sufficiently accurate; if either is wrong the extracted parameters will not reflect the actual medium effects.
What would settle it
A new lattice QCD calculation of the heavy quark potential at the relevant temperatures that shows substantially stronger Debye screening than the weak value extracted here, or a new set of bottomonium R_AA measurements at an intermediate energy that cannot be reproduced by the same potential parameters.
read the original abstract
We use deep learning under Bayesian perspective to quantitatively extract the in-medium heavy quark (HQ) potential from bottomonium nuclear modification factors ($R_{AA}$) measured across multiple heavy ion collision systems at the Large Hadron Collider (LHC) and the Relativistic Heavy-Ion Collider (RHIC). The in-medium HQ potential, comprising both a real and imaginary part, is parameterized and incorporated into a time-dependent Schr\"odinger equation to model the wave function evolution of $b\bar{b}$ dipoles within a hydrodynamically evolving hot QCD medium. We construct Convolutional Neural Networks (CNNs) to capture the non-linear correspondence between the heavy quark potential $V(T,r)$ and the bottomonium $R_{AA}$ for Pb-Pb collisions at 5.02 TeV and 2.76 TeV, and Au-Au collisions at 200 GeV. Training datasets are generated by sampling the potential parameters and are further augmented using Principal Component Analysis (PCA) and Gaussian Process Regression (GPR). After validating the stability and correctness of the CNNs, we employ Stochastic Gradient Langevin Dynamics (SGLD) to perform a simultaneous Bayesian inverse extraction of the optimal potential parameters and their posterior distributions using experimental data of bottomonium $R_{AA}$ in both LHC and RHIC energies. Our joint multi-energy extraction suggests that, within the present parametrization and hydrodynamic background, the real part of the in-medium potential remains close to the vacuum Cornell form, corresponding to a relatively weak screened Debye mass across RHIC to LHC energies. By contrast, the imaginary part is more strongly constrained by the data and provides the dominant contribution to bottomonium suppression from RHIC to LHC energies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a Bayesian deep learning framework to extract parameters of a parameterized in-medium heavy quark potential (real and imaginary parts) from bottomonium R_AA data across RHIC (Au-Au at 200 GeV) and LHC (Pb-Pb at 2.76 and 5.02 TeV) energies. Training data are generated by solving the time-dependent Schrödinger equation for b b-bar wave functions in a hydrodynamically evolving medium with sampled potential parameters; CNNs learn the forward map from V(T,r) to R_AA, augmented via PCA and GPR, and SGLD performs the inverse extraction to obtain posterior distributions. The central result is that, within the chosen parametrization and background, Re(V) remains close to the vacuum Cornell form with weak Debye screening while Im(V) dominates suppression.
Significance. If robust, the work supplies a unified multi-energy, data-driven extraction of the in-medium potential, clarifying the relative roles of real and imaginary components in bottomonium suppression. It combines established components (TDSE evolution, hydrodynamics) with CNN regression and SGLD sampling for Bayesian inversion, providing a reproducible template for similar inverse problems in heavy-ion physics. The joint analysis across energies is a clear strength.
major comments (2)
- [Training data generation (methods)] The hydrodynamic background (temperature and flow profiles) is held fixed when generating training data for the CNNs. No marginalization or propagation of uncertainties arising from initial conditions, equation of state, shear viscosity, or medium lifetime is performed. Because any mismatch between the assumed evolution and reality is absorbed into the inferred potential parameters, the posteriors on Re(V) and Im(V) do not incorporate this systematic; this directly underpins the headline claim of weak screening in the real part.
- [CNN construction and validation] The validation of the CNNs is stated to confirm stability and correctness, yet no quantitative metrics are supplied (e.g., mean-squared error on held-out synthetic data, recovery accuracy for known input potentials, or calibration of predicted R_AA uncertainties). Without these, the reliability of the learned forward map prior to SGLD inversion cannot be assessed at the level required for the central extraction.
minor comments (2)
- [Abstract] The abstract refers to augmentation with PCA and GPR but does not detail the procedure, the resulting training-set size, or any checks that the augmentation preserves the physical mapping.
- [Potential parametrization] The explicit functional form and all free parameters of V(T,r) (real and imaginary) should be written out with equation numbers so that the dimensionality of the SGLD sampling is immediately clear.
Simulated Author's Rebuttal
We thank the referee for the thorough review and constructive feedback on our manuscript. We address each major comment below and have prepared revisions to strengthen the presentation of our methods and results.
read point-by-point responses
-
Referee: The hydrodynamic background (temperature and flow profiles) is held fixed when generating training data for the CNNs. No marginalization or propagation of uncertainties arising from initial conditions, equation of state, shear viscosity, or medium lifetime is performed. Because any mismatch between the assumed evolution and reality is absorbed into the inferred potential parameters, the posteriors on Re(V) and Im(V) do not incorporate this systematic; this directly underpins the headline claim of weak screening in the real part.
Authors: We agree that fixing the hydrodynamic background means the extracted posteriors are conditional on this choice and do not marginalize over uncertainties in initial conditions, EOS, viscosity, or lifetime. This is a standard limitation in such inverse problems given the prohibitive cost of regenerating large training datasets for varied hydro evolutions. The headline claim of weak screening in Re(V) is therefore presented as holding within the employed hydrodynamic model. In the revised manuscript we will add an explicit discussion of this systematic uncertainty, including a note that future extensions could incorporate marginalization over hydro parameters. revision: partial
-
Referee: The validation of the CNNs is stated to confirm stability and correctness, yet no quantitative metrics are supplied (e.g., mean-squared error on held-out synthetic data, recovery accuracy for known input potentials, or calibration of predicted R_AA uncertainties). Without these, the reliability of the learned forward map prior to SGLD inversion cannot be assessed at the level required for the central extraction.
Authors: We acknowledge that the manuscript describes CNN validation only qualitatively. In the revised version we will add quantitative metrics in the CNN construction section, including mean-squared error on a held-out test set of synthetic R_AA values, recovery accuracy when inverting known input potentials, and calibration checks on the predicted uncertainties. These additions will allow readers to assess the forward map reliability before the SGLD inversion step. revision: yes
Circularity Check
No significant circularity: standard Bayesian extraction from external data
full rationale
The paper parameterizes V(T,r), solves the TDSE in a fixed hydrodynamic background to generate training R_AA, trains CNN surrogates, and uses SGLD to infer posterior parameters from measured bottomonium R_AA at RHIC and LHC. The headline conclusions about Re(V) remaining close to the vacuum Cornell form and Im(V) dominating suppression are direct outputs of this fit to independent experimental data. No derivation step reduces the result to the inputs by construction, no self-citation is load-bearing for the central claim, and the hydro background plus parametrization are stated as explicit modeling choices rather than smuggled assumptions. The analysis is self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- potential parameters
axioms (2)
- domain assumption The evolution of b bbar dipoles is accurately described by the time-dependent Schrödinger equation inside a hydrodynamically evolving medium.
- domain assumption The hydrodynamic background provides a reliable temperature and flow profile for the medium.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclearVI =−i T a0 (a1r+a2 ra3), VR =−α (μD + e−μD r / r) + 2σ/μD − σ e−μD r (2 + μD r)/μD with μD = a4 T … ; Tsw switching to vacuum Cornell
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction uncleartime-dependent Schrödinger equation … coupled with the hydrodynamic model … CNNs … SGLD … joint multi-energy extraction
Forward citations
Cited by 1 Pith paper
-
Effects of event-by-event hydrodynamic fluctuations on bottomonium dynamics in Pb--Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV
Event-by-event hydrodynamic fluctuations have marginal effects on bottomonium R_AA and v2 in 5.02 TeV Pb-Pb collisions.
Reference graph
Works this paper leans on
-
[1]
A. Bazavov et al.,The chiral and deconfinement aspects of the QCD transition,Phys. Rev. D 85(2012) 054503 [1111.1710]. – 31 – [2]ALICEcollaboration,Elliptic flow of charged particles in Pb-Pb collisions at 2.76 TeV, Phys. Rev. Lett.105(2010) 252302 [1011.3914]. [3]STARcollaboration,Bulk Properties of the Medium Produced in Relativistic Heavy-Ion Collision...
- [2]
-
[3]
C. Gale, S. Jeon and B. Schenke,Hydrodynamic Modeling of Heavy-Ion Collisions,Int. J. Mod. Phys. A28(2013) 1340011 [1301.5893]
work page Pith review arXiv 2013
- [4]
-
[5]
Matsui and H
T. Matsui and H. Satz,J/ψSuppression by Quark-Gluon Plasma Formation,Phys. Lett. B 178(1986) 416
1986
-
[6]
T. Song and J. Zhao,Influence of heavy quark potential on quarkonium production in quark-gluon plasma,Phys. Rev. C113(2026) 024916 [2512.06339]
-
[7]
S. Chen and M. He,Heavy quarkonium dissociation by thermal gluons at next-to-leading order in the Quark–Gluon Plasma,Phys. Lett. B786(2018) 260 [1805.06346]
- [8]
-
[9]
M. Strickland and D. Bazow,Thermal Bottomonium Suppression at RHIC and LHC,Nucl. Phys. A879(2012) 25 [1112.2761]
-
[10]
X. Zhao and R. Rapp,Medium Modifications and Production of Charmonia at LHC,Nucl. Phys. A859(2011) 114 [1102.2194]
-
[11]
J.-P. Blaizot and M.A. Escobedo,Approach to equilibrium of a quarkonium in a quark-gluon plasma,Phys. Rev. D98(2018) 074007 [1803.07996]
- [12]
-
[13]
Karsch, M.T
F. Karsch, M.T. Mehr and H. Satz,Color Screening and Deconfinement for Bound States of Heavy Quarks,Z. Phys. C37(1988) 617
1988
-
[14]
L. Grandchamp, R. Rapp and G.E. Brown,In medium effects on charmonium production in heavy ion collisions,Phys. Rev. Lett.92(2004) 212301 [hep-ph/0306077]
- [15]
-
[16]
B. Krouppa, A. Rothkopf and M. Strickland,Bottomonium suppression using a lattice QCD vetted potential,Phys. Rev. D97(2018) 016017 [1710.02319]
- [17]
-
[18]
Statistical hadronization of heavy quarks in ultra-relativistic nucleus-nucleus collisions,
A. Andronic, P. Braun-Munzinger, K. Redlich and J. Stachel,Statistical hadronization of heavy quarks in ultra-relativistic nucleus-nucleus collisions,Nucl. Phys. A789(2007) 334 [nucl-th/0611023]. – 32 –
- [19]
-
[20]
N. Brambilla, M. ´A. Escobedo, M. Strickland, A. Vairo, P. Vander Griend and J.H. Weber, Bottomonium suppression in an open quantum system using the quantum trajectories method,JHEP05(2021) 136 [2012.01240]
-
[21]
S. Delorme, R. Katz, T. Gousset, P.B. Gossiaux and J.-P. Blaizot,Quarkonium dynamics in the quantum Brownian regime with non-abelian quantum master equations,JHEP06(2024) 060 [2402.04488]
-
[22]
M. Margotta, K. McCarty, C. McGahan, M. Strickland and D. Yager-Elorriaga,Quarkonium states in a complex-valued potential,Phys. Rev. D83(2011) 105019 [1101.4651]. [25]HotQCDcollaboration,Static quark-antiquark interactions at nonzero temperature from lattice QCD,Phys. Rev. D105(2022) 054513 [2110.11659]
- [23]
- [24]
- [25]
-
[26]
M. Omana Kuttan, K. Zhou, J. Steinheimer, A. Redelbach and H. Stoecker,An equation-of-state-meter for CBM using PointNet,JHEP21(2020) 184 [2107.05590]
- [27]
- [28]
- [29]
- [30]
- [31]
- [32]
- [33]
-
[34]
Satz,Colour deconfinement and quarkonium binding,J
H. Satz,Colour deconfinement and quarkonium binding,J. Phys. G32(2006) R25 [hep-ph/0512217]
-
[35]
Crank and P
J. Crank and P. Nicolson,A practical method for numerical evaluation of solutions of partial – 33 – differential equations of the heat-conduction type,Mathematical Proceedings of the Cambridge Philosophical Society43(1947) 50
1947
- [36]
- [37]
- [38]
-
[39]
L. Wen and B. Chen,Bottomonium sequential suppression and strong heavy-quark potential in heavy-ion collisions,Physics Letters B839(2023) 137774. [43]CMScollaboration,Suppression ofΥ(1S),Υ(2S)andΥ(3S)production in PbPb collisions at√sNN = 2.76 TeV,Phys. Lett. B770(2017) 357 [1611.01510]. [44]CMScollaboration,Measurement of nuclear modification factors ofΥ...
-
[40]
C.M.S. Collaboration,Measurement of the production cross section ratioσ(χ b2(1p)) / σ(χb1(1p)) in pp collisions at √s= 8 TeV,Physics Letters B743(2015) 383 [1409.5761]
-
[41]
J.-P. Lansberg,New observables in inclusive production of quarkonia,Physics Reports889 (2020) 1 [1903.09185]. [48]PHENIXcollaboration,υ(1S+ 2S+ 3S)production ind+Au andp+pcollisions at√sN N = 200GeV and cold-nuclear matter effects,Phys. Rev. C87(2013) 044909 [1211.4017]. [49]STARcollaboration,Suppression ofΥproduction in d+Au and Au+Au collisions at√sN N=...
- [42]
-
[43]
Aaijet al.(LHCb), JHEP07, 069 (2023), arXiv:2212.12664 [hep-ex]
M.L. Miller, K. Reygers, S.J. Sanders and P. Steinberg,Glauber modeling in high-energy nuclear collisions,Annual Review of Nuclear and Particle Science57(2007) 205. [52]LHCbcollaboration,Measurement ofΥproduction in pp collisions at √s= 5 TeV,JHEP 07(2023) 069 [2212.12664]. [53]ATLAScollaboration,Measurement of quarkonium production in proton–lead and pro...
-
[44]
Akamatsu,Heavy quark master equations in the Lindblad form at high temperatures,Phys
Y. Akamatsu,Heavy quark master equations in the Lindblad form at high temperatures, Phys. Rev. D91(2015) 056002 [1403.5783]
-
[45]
D. Lafferty and A. Rothkopf,Quarkonium Phenomenology from a Generalised Gauss Law, Universe5(2019) 119 [1906.04014]
-
[46]
Y. Burnier and A. Rothkopf,Complex heavy-quark potential and Debye mass in a gluonic medium from lattice QCD,Phys. Rev. D95(2017) 054511 [1607.04049]. – 34 –
-
[47]
B. Schenke, S. Jeon and C. Gale,Elliptic and triangular flow in event-by-event (3+1)D viscous hydrodynamics,Phys. Rev. Lett.106(2011) 042301 [1009.3244]
-
[48]
B. Schenke, S. Jeon and C. Gale,(3+1)D hydrodynamic simulation of relativistic heavy-ion collisions,Phys. Rev. C82(2010) 014903 [1004.1408]
-
[49]
J.E. Bernhard, J.S. Moreland, S.A. Bass, J. Liu and U. Heinz,Applying Bayesian parameter estimation to relativistic heavy-ion collisions: simultaneous characterization of the initial state and quark-gluon plasma medium,Phys. Rev. C94(2016) 024907 [1605.03954]. [60]HotQCDcollaboration,Equation of state in ( 2+1 )-flavor QCD,Phys. Rev. D90(2014) 094503 [1407.6387]
-
[50]
Loizides,Glauber modeling of high-energy nuclear collisions at the subnucleon level,Phys
C. Loizides,Glauber modeling of high-energy nuclear collisions at the subnucleon level,Phys. Rev. C94(2016) 024914 [1603.07375]
-
[51]
H.,Analysis of a complex of statistical variables into principal components.,Journal of Educational Psychology24(6)(1933) 417–441
H. H.,Analysis of a complex of statistical variables into principal components.,Journal of Educational Psychology24(6)(1933) 417–441
1933
-
[52]
J.E. Bernhard, P.W. Marcy, C.E. Coleman-Smith, S. Huzurbazar, R.L. Wolpert and S.A. Bass,Quantifying properties of hot and dense QCD matter through systematic model-to-data comparison,Phys. Rev. C91(2015) 054910 [1502.00339]
-
[53]
C.E. Rasmussen and C.K.I. Williams,Gaussian Processes for Machine Learning, The MIT Press (2006). [65]HotQCDcollaboration,Unscreened forces in the quark-gluon plasma?,Phys. Rev. D109 (2024) 074504 [2308.16587]
-
[54]
Y. Burnier and A. Rothkopf,Bayesian Approach to Spectral Function Reconstruction for Euclidean Quantum Field Theories,Phys. Rev. Lett.111(2013) 182003 [1307.6106]
-
[55]
A. Bzdak and V. Koch,Rapidity dependence of proton cumulants and correlation functions, Phys. Rev. C96(2017) 054905 [1707.02640]
-
[56]
B. Schenke, C. Shen and D. Teaney,Transverse momentum fluctuations and their correlation with elliptic flow in nuclear collision,Phys. Rev. C102(2020) 034905 [2004.00690]
- [57]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.