pith. sign in

arxiv: 2605.22199 · v1 · pith:QFXPBRAAnew · submitted 2026-05-21 · ✦ hep-ph · nucl-th

Equation of State at High Baryon Densities from a Thermodynamically Informed Neural Network

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

classification ✦ hep-ph nucl-th
keywords equation of stateneural networkthermodynamic consistencyheavy-ion collisionsbaryon densitylattice QCDhadron resonance gasphysics-informed machine learning
0
0 comments X

The pith

A neural network trained on thermodynamic constraints produces a four-dimensional equation of state that matches hadron resonance gas and lattice QCD results while extending into high baryon densities.

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

The paper constructs an equation of state for strongly interacting matter using a deep neural network that operates across temperature and three conserved charge densities. It is trained to reproduce hadron resonance gas thermodynamics at typical particlization conditions and to stay consistent with lattice QCD results at low baryon chemical potential. A physics-informed loss function enforces the required thermodynamic relations everywhere in the phase space so that the network can extrapolate into the high-density region relevant for beam-energy-scan experiments. The resulting equation of state is then inserted into an integrated hydrokinetic model at zero net strangeness and fixed electric-to-baryon ratio.

Core claim

A deep neural network supplied with a loss term that penalizes violations of thermodynamic identities yields an equation of state valid over the full range of temperature and chemical potentials; this equation of state reproduces known hadron resonance gas results at low to moderate densities, remains consistent with lattice QCD at small baryon chemical potential, and supplies thermodynamically stable extrapolations into the high-baryon-density domain that neither lattice calculations nor hadron resonance gas models can reach.

What carries the argument

A deep neural network whose training loss explicitly enforces thermodynamic consistency relations among pressure, entropy, and number densities across four dimensions.

If this is right

  • The equation of state can be inserted directly into hybrid hydrodynamic models for relativistic heavy-ion collisions.
  • Thermodynamic consistency is preserved across the entire temperature and chemical-potential plane.
  • The construction supplies controlled access to the density region targeted by RHIC BES, FAIR, HADES, and CBM experiments.
  • The same network architecture can be retrained at different values of net strangeness or electric charge.

Where Pith is reading between the lines

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

  • The method could be adapted to produce equations of state at the still higher densities and lower temperatures relevant for neutron-star cores.
  • Systematic variation of the network architecture or loss-function weights would quantify how much the high-density behavior depends on the choice of consistency constraints.
  • Embedding the network inside a larger Bayesian analysis could propagate its uncertainty into final hydrodynamic predictions.

Load-bearing premise

Training a neural network only on low-density data together with thermodynamic consistency constraints will produce physically reliable extrapolations at high baryon densities without introducing uncontrolled artifacts or violations of causality or stability.

What would settle it

A direct comparison of the neural-network predictions against new lattice QCD results at moderate-to-high baryon chemical potential or against measured particle yields and flow observables from heavy-ion runs at FAIR or CBM would show whether the high-density extrapolation is accurate.

Figures

Figures reproduced from arXiv: 2605.22199 by Musfer Adzhymambetov.

Figure 1
Figure 1. Figure 1: FIG. 1: Schematic of the hybrid EoS construction in [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Relative difference between the neural-network [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Contours of constant [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Adiabatic trajectories of constant [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
read the original abstract

We present a four-dimensional equation of state for strongly interacting matter at finite temperature and conserved charge densities, constructed using a deep neural network. It is designed for direct use in hybrid models of relativistic heavy-ion collisions: it reproduces hadron resonance gas thermodynamics at typical particlization scales, is consistent with lattice QCD at low baryon chemical potential, and extrapolates into the high-density region inaccessible to either approach, which is precisely the regime targeted by RHIC BES, FAIR, HADES, and CBM. Thermodynamic consistency throughout the full phase space is enforced via a physics-informed loss function. We demonstrate the developed equation of state by implementing it at zero net strangeness and fixed electric-to-baryon charge ratio within the integrated hydrokinetic model.

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

2 major / 2 minor

Summary. The manuscript presents a four-dimensional equation of state for strongly interacting matter constructed via a deep neural network. The network is trained on hadron resonance gas thermodynamics at low densities and lattice QCD results at small baryon chemical potential, with a physics-informed loss function enforcing thermodynamic consistency relations such as s = (∂p/∂T)_μ and n_B = (∂p/∂μ_B)_T. The resulting EoS is extrapolated to high baryon densities and demonstrated within the integrated hydrokinetic model at zero net strangeness and fixed electric-to-baryon ratio, targeting regimes relevant to RHIC BES, FAIR, HADES, and CBM.

Significance. If the extrapolation maintains physical properties including causality and stability, the approach offers a practical way to bridge low-density lattice and HRG results with the high-density domain inaccessible to direct computation, enabling consistent hydrodynamic modeling across the full phase space probed by current and future heavy-ion experiments. The physics-informed training is a clear methodological strength that could reduce discontinuities in hybrid models.

major comments (2)
  1. The physics-informed loss enforces first-derivative thermodynamic identities but provides no explicit constraint or post-training verification on second derivatives that control the speed of sound c_s^2 = (∂p/∂ε)_{s/n} or the sign of the compressibility. Without such checks or high-density anchor points, the extrapolation can satisfy the training constraints while producing regions with c_s^2 > 1 or indefinite pressure Hessian, undermining usability in hydrodynamics.
  2. Section describing the high-density extrapolation and model implementation: the central claim that the EoS is 'physically usable' in hybrid models requires demonstration that causality and mechanical stability hold throughout the extrapolated domain; the current results show implementation but do not report systematic scans or bounds on these quantities.
minor comments (2)
  1. Notation for the four conserved charges and the precise definition of the fixed electric-to-baryon ratio should be clarified in the methods section to aid reproducibility.
  2. Figure captions for the EoS slices could explicitly state the range of μ_B over which thermodynamic consistency was verified versus extrapolated.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The two major comments correctly identify a gap in our current presentation regarding verification of second-order thermodynamic properties in the extrapolated regime. We will revise the manuscript to address both points explicitly.

read point-by-point responses
  1. Referee: The physics-informed loss enforces first-derivative thermodynamic identities but provides no explicit constraint or post-training verification on second derivatives that control the speed of sound c_s^2 = (∂p/∂ε)_{s/n} or the sign of the compressibility. Without such checks or high-density anchor points, the extrapolation can satisfy the training constraints while producing regions with c_s^2 > 1 or indefinite pressure Hessian, undermining usability in hydrodynamics.

    Authors: We agree that the physics-informed loss function enforces only the first-order thermodynamic relations (entropy and baryon density from pressure derivatives) and does not directly constrain second derivatives. While the training data and architecture are chosen to favor smooth, physical behavior, we did not perform or report systematic post-training checks for causality (c_s^2 < 1) or positive compressibility throughout the high-density extrapolation. In the revised manuscript we will add a dedicated subsection with quantitative scans of c_s^2 and the pressure Hessian eigenvalues over the full (T, μ_B, μ_S, μ_Q) domain relevant to the hydrokinetic simulations, including explicit bounds and any regions where violations occur. revision: yes

  2. Referee: Section describing the high-density extrapolation and model implementation: the central claim that the EoS is 'physically usable' in hybrid models requires demonstration that causality and mechanical stability hold throughout the extrapolated domain; the current results show implementation but do not report systematic scans or bounds on these quantities.

    Authors: The manuscript currently demonstrates implementation of the EoS inside the integrated hydrokinetic model at zero net strangeness and fixed charge ratio, but indeed does not include systematic scans or quantitative bounds on causality and stability in the extrapolated region. We will revise the relevant section to include these checks (as described in the response to the first comment) and will update the abstract and conclusions to reflect the new verification results rather than asserting usability without supporting evidence. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper trains a four-dimensional neural network on low-density lattice QCD and hadron resonance gas data while using a physics-informed loss to enforce thermodynamic identities such as s = (∂p/∂T)_μ and n_B = (∂p/∂μ_B)_T. The high-baryon-density regime is produced by the trained network's extrapolation rather than by any algebraic reduction or redefinition that makes the output identical to the training inputs. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior work by the same authors are invoked to force the result. The construction is therefore a standard data-driven model with consistency constraints; the extrapolated values remain independent of the inputs by construction and can be tested against future high-density observables.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The model rests on the assumption that thermodynamic identities can be enforced through a loss term during training and that the resulting network will generalize outside its training domain. The only explicit free parameters are the neural-network weights themselves.

free parameters (1)
  • Neural-network weights and biases
    All internal parameters of the deep network are determined by minimizing the combined data-matching and thermodynamic-consistency loss.
axioms (1)
  • domain assumption Thermodynamic consistency (correct relations among pressure, energy density, entropy, and chemical potentials) must hold everywhere in the tabulated domain.
    This is invoked to justify the physics-informed loss function and is a standard requirement for any equation of state used in relativistic hydrodynamics.

pith-pipeline@v0.9.0 · 5654 in / 1437 out tokens · 54080 ms · 2026-05-22T05:10:08.664901+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

56 extracted references · 56 canonical work pages · 30 internal anchors

  1. [1]

    The equation of state in (2+1)-flavor QCD

    A. Bazavovet al.(HotQCD), Phys. Rev. D90, 094503 (2014), arXiv:1407.6387 [hep-lat]

  2. [2]

    The QCD equation of state with dynamical quarks

    S. Borsanyi, G. Endrodi, Z. Fodor, A. Jakovac, S. D. 8 Katz, S. Krieg, C. Ratti, and K. K. Szabo, JHEP11, 077, arXiv:1007.2580 [hep-lat]

  3. [3]

    QCD Equation of State and Hadron Resonance Gas

    P. Huovinen and P. Petreczky, Nucl. Phys. A837, 26 (2010), arXiv:0912.2541 [hep-ph]

  4. [4]

    Hadron Resonance Gas Equation of State from Lattice QCD

    V. Vovchenko, D. V. Anchishkin, and M. I. Gorenstein, Phys. Rev. C91, 024905 (2015), arXiv:1412.5478 [nucl- th]

  5. [5]

    Interacting hadron resonance gas meets lattice QCD

    A. Andronic, P. Braun-Munzinger, J. Stachel, and M. Winn, Phys. Lett. B718, 80 (2012), arXiv:1201.0693 [nucl-th]

  6. [6]

    Simulating QCD at finite density

    P. de Forcrand, PoSLA T2009, 010 (2009), arXiv:1005.0539 [hep-lat]

  7. [7]

    M. M. Aggarwalet al.(STAR), An Experimental Explo- ration of the QCD Phase Diagram: The Search for the Critical Point and the Onset of De-confinement (2010), arXiv:1007.2613 [nucl-ex]

  8. [8]

    The High-Acceptance Dielectron Spectrometer HADES

    G. Agakishievet al.(HADES), Eur. Phys. J. A41, 243 (2009), arXiv:0902.3478 [nucl-ex]

  9. [9]

    Challenges in QCD matter physics - The Compressed Baryonic Matter experiment at FAIR

    T. Ablyazimovet al.(CBM), Eur. Phys. J. A53, 60 (2017), arXiv:1607.01487 [nucl-ex]

  10. [10]

    Comparison of Chemical Freeze-Out Criteria in Heavy-Ion Collisions

    J. Cleymans, H. Oeschler, K. Redlich, and S. Wheaton, Phys. Rev. C73, 034905 (2006), arXiv:hep-ph/0511094

  11. [11]

    Decoding the phase structure of QCD via particle production at high energy

    A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel, Nature561, 321 (2018), arXiv:1710.09425 [nucl-th]

  12. [12]

    M. A. Stephanov, K. Rajagopal, and E. V. Shuryak, Phys. Rev. Lett.81, 4816 (1998), arXiv:hep-ph/9806219

  13. [13]

    Bzdak, S

    A. Bzdak, S. Esumi, V. Koch, J. Liao, M. Stephanov, and N. Xu, Phys. Rept.853, 1 (2020), arXiv:1906.00936 [nucl-th]

  14. [14]

    Chodos, R

    A. Chodos, R. L. Jaffe, K. Johnson, C. B. Thorn, and V. F. Weisskopf, Phys. Rev. D9, 3471 (1974)

  15. [15]

    Nuclei in a chiral SU(3) model

    P. Papazoglou, D. Zschiesche, S. Schramm, J. Schaffner- Bielich, H. Stoecker, and W. Greiner, Phys. Rev. C59, 411 (1999), arXiv:nucl-th/9806087

  16. [16]

    Motornenko, J

    A. Motornenko, J. Steinheimer, V. Vovchenko, S. Schramm, and H. Stoecker, Phys. Rev. C101, 034904 (2020), arXiv:1905.00866 [hep-ph]

  17. [17]

    The Equation of State of Deconfined Matter at Finite Chemical Potential in a Quasiparticle Description

    A. Peshier, B. Kampfer, and G. Soff, Phys. Rev. C61, 045203 (2000), arXiv:hep-ph/9911474

  18. [18]

    R. A. Schneider and W. Weise, Phys. Rev. C64, 055201 (2001), arXiv:hep-ph/0105242

  19. [19]

    Matching Excluded Volume Hadron Resonance Gas Models and Perturbative QCD to Lattice Calculations

    M. Albright, J. Kapusta, and C. Young, Phys. Rev. C 90, 024915 (2014), arXiv:1404.7540 [nucl-th]

  20. [20]

    Monnai, B

    A. Monnai, B. Schenke, and C. Shen, Phys. Rev. C100, 024907 (2019), arXiv:1902.05095 [nucl-th]

  21. [21]

    J. Gong, H. Roch, and C. Shen, Phys. Rev. C111, 044912 (2025), arXiv:2410.22160 [nucl-th]

  22. [22]

    Omana Kuttan, J

    M. Omana Kuttan, J. Steinheimer, K. Zhou, and H. Stoecker, Phys. Rev. Lett.131, 202303 (2023), arXiv:2211.11670 [hep-ph]

  23. [23]

    Liu, X.-Y

    F.-L. Liu, X.-Y. Wu, S. Cao, G.-Y. Qin, and X.-N. Wang, Phys. Lett. B848, 138355 (2024), arXiv:2304.08787 [hep- ph]

  24. [24]

    Ayriyan, D

    A. Ayriyan, D. Blaschke, J. P. Carlomagno, G. A. Con- trera, and A. G. Grunfeld, Universe11, 141 (2025), arXiv:2501.00115 [nucl-th]

  25. [25]

    S. Soma, L. Wang, S. Shi, H. St¨ ocker, and K. Zhou, JCAP 08, 071, arXiv:2201.01756 [hep-ph]

  26. [26]

    Morawski and M

    F. Morawski and M. Bejger, Astron. Astrophys.642, A78 (2020), arXiv:2006.07194 [astro-ph.HE]

  27. [27]

    Physics Informed Deep Learning (Part I): Data-driven Solutions of Nonlinear Partial Differential Equations

    M. Raissi, P. Perdikaris, and G. E. Karniadakis, Physics informed deep learning (part i): Data-driven solu- tions of nonlinear partial differential equations (2017), arXiv:1711.10561 [cs.AI]

  28. [28]

    Li, H.-L

    F.-P. Li, H.-L. L¨ u, L.-G. Pang, and G.-Y. Qin, Phys. Lett. B844, 138088 (2023), arXiv:2211.07994 [hep-ph]

  29. [29]

    Li, L.-G

    F.-P. Li, L.-G. Pang, and G.-Y. Qin, Phys. Lett. B868, 139692 (2025), arXiv:2501.10012 [nucl-th]

  30. [30]

    S. A. Bass and A. Dumitru, Phys. Rev. C61, 064909 (2000), arXiv:nucl-th/0001033

  31. [31]

    Fully integrated transport approach to heavy ion reactions with an intermediate hydrodynamic stage

    H. Petersen, J. Steinheimer, G. Burau, M. Bleicher, and H. St¨ ocker, Phys. Rev. C78, 044901 (2008), arXiv:0806.1695 [nucl-th]

  32. [32]

    Cooper and G

    F. Cooper and G. Frye, Phys. Rev. D10, 186 (1974)

  33. [33]

    Particlization in hybrid models

    P. Huovinen and H. Petersen, Eur. Phys. J. A48, 171 (2012), arXiv:1206.3371 [nucl-th]

  34. [34]

    Kapusta, T

    J. Kapusta, T. Welle, and C. Plumberg, Phys. Rev. C 106, 014909 (2022), arXiv:2112.07563 [nucl-th]

  35. [35]

    Parotto, M

    P. Parotto, M. Bluhm, D. Mroczek, M. Nahrgang, J. Noronha-Hostler, K. Rajagopal, C. Ratti, T. Sch¨ afer, and M. Stephanov, Phys. Rev. C101, 034901 (2020), arXiv:1805.05249 [hep-ph]

  36. [36]

    Noronha-Hostler, P

    J. Noronha-Hostler, P. Parotto, C. Ratti, and J. M. Stafford, Phys. Rev. C100, 064910 (2019), arXiv:1902.06723 [hep-ph]

  37. [37]

    Higher order fluctuations and correlations of conserved charges from lattice QCD

    S. Borsanyi, Z. Fodor, J. N. Guenther, S. K. Katz, K. K. Szabo, A. Pasztor, I. Portillo, and C. Ratti, JHEP10, 205, arXiv:1805.04445 [hep-lat]

  38. [38]

    Dashen, S.-K

    R. Dashen, S.-K. Ma, and H. J. Bernstein, Phys. Rev. 187, 345 (1969)

  39. [39]

    Venugopalan and M

    R. Venugopalan and M. Prakash, Nucl. Phys. A546, 718 (1992)

  40. [40]

    Karpenko, hadronSampler: a Cooper-Frye hadron sampler + resonance decayer + an interface to the UrQMD afterburner,https://github.com/yukarpenko/ hadronSampler, gitHub repository

    I. Karpenko, hadronSampler: a Cooper-Frye hadron sampler + resonance decayer + an interface to the UrQMD afterburner,https://github.com/yukarpenko/ hadronSampler, gitHub repository

  41. [41]

    Adzhymambetov, Neural Network Equation of State for Heavy-Ion Collisions,https://github.com/Musfer/ NNEoS(2026), gitHub repository

    M. Adzhymambetov, Neural Network Equation of State for Heavy-Ion Collisions,https://github.com/Musfer/ NNEoS(2026), gitHub repository

  42. [42]

    Self-consistent conversion of a viscous fluid to particles

    D. Molnar and Z. Wolff, Phys. Rev. C95, 024903 (2017), arXiv:1404.7850 [nucl-th]

  43. [43]

    Anchishkin, J

    D. Anchishkin, J. Phys. G49, 055109 (2022), arXiv:2303.06655 [nucl-th]

  44. [44]
  45. [45]

    Adzhymambetov and Y

    M. Adzhymambetov and Y. Sinyukov, Phys. Rev. C112, 034906 (2025), arXiv:2412.00458 [hep-ph]

  46. [46]

    G´ oes-Hirayama, J

    R. G´ oes-Hirayama, J. Egger, Z. Paul´ ınyov´ a, I. Karpenko, and H. Elfner, Phys. Rev. C113, 044906 (2026), arXiv:2507.19389 [hep-ph]

  47. [47]

    Systematic Investigation of Negative Cooper-Frye Contributions in Heavy Ion Collisions Using Coarse-grained Molecular Dynamics

    D. Oliinychenko, P. Huovinen, and H. Petersen, Phys. Rev. C91, 024906 (2015), arXiv:1411.3912 [nucl-th]

  48. [48]

    Y. M. Sinyukov, S. V. Akkelin, and Y. Hama, Phys. Rev. Lett.89, 052301 (2002), arXiv:nucl-th/0201015

  49. [49]

    PyTorch: An Imperative Style, High-Performance Deep Learning Library

    A. Paszke, S. Gross, F. Massa, A. Lerer, J. Brad- bury, G. Chanan, T. Killeen, Z. Lin, N. Gimelshein, L. Antiga, A. Desmaison, A. K¨ opf, E. Yang, Z. De- Vito, M. Raison, A. Tejani, S. Chilamkurthy, B. Steiner, L. Fang, J. Bai, and S. Chintala, Pytorch: An impera- tive style, high-performance deep learning library (2019), arXiv:1912.01703 [cs.LG]

  50. [50]

    K. He, X. Zhang, S. Ren, and J. Sun, Deep residual learning for image recognition (2015), arXiv:1512.03385 [cs.CV]

  51. [51]

    D. P. Kingma and J. Ba, Adam: A method for stochastic optimization (2017), arXiv:1412.6980 [cs.LG]

  52. [52]

    Sound velocity bound and neutron stars

    P. Bedaque and A. W. Steiner, Phys. Rev. Lett.114, 9 031103 (2015), arXiv:1408.5116 [nucl-th]

  53. [53]

    I. Tews, J. Carlson, S. Gandolfi, and S. Reddy, Astrophys. J.860, 149 (2018), arXiv:1801.01923 [nucl-th]

  54. [54]

    Annala, T

    E. Annala, T. Gorda, A. Kurkela, J. N¨ attil¨ a, and A. Vuorinen, Nature Phys.16, 907 (2020), arXiv:1903.09121 [astro-ph.HE]

  55. [55]

    Bleicheret al., J

    M. Bleicheret al., J. Phys. G25, 1859 (1999), arXiv:hep- ph/9909407

  56. [56]

    Rathod, Y

    N. Rathod, Y. Sinyukov, M. Adzhymambetov, and H. Zbroszczyk, Phys. Rev. C113, 024911 (2026), arXiv:2506.13944 [nucl-th]