pith. sign in

arxiv: 2503.15719 · v2 · submitted 2025-03-19 · ❄️ cond-mat.stat-mech · hep-th· nlin.SI· nucl-th

Phase transitions and finite-size effects in integrable virial statistical models

Xin An , Francesco Giglio , Giulio Landolfi This is my paper

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

classification ❄️ cond-mat.stat-mech hep-thnlin.SInucl-th
keywords virial expansionphase transitionsfinite-size effectshydrodynamic PDEsshock wavesQCD phase diagramcritical pointsintegrable systems
0
0 comments X

The pith

Virial expansion models for fluids with finite particle number N are exactly solvable at any order via solutions to nonlinear C-integrable hydrodynamic PDEs, and phase transitions appear as shock waves when N goes to infinity.

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

The paper establishes that thermodynamic models built from a virial expansion of internal energy in terms of volume density admit exact solutions for any finite number of particles. Expectation values of observables such as volume density are obtained directly from solutions of the associated nonlinear PDEs of hydrodynamic type. In the thermodynamic limit the same equations produce phase transitions through the formation of classical shock waves. The framework is used to build a QCD phase diagram that locates critical points for the nuclear liquid-gas transition and the hadron gas to quark-gluon plasma transition, while showing how finite particle number smooths out the critical signatures.

Core claim

The models formulated for finite-size systems with N particles are exactly solvable to any expansion order, as expectation values of physical observables are determined from solutions to nonlinear C-integrable partial differential equations of hydrodynamic type. In the limit N to infinity, phase transitions emerge as classical shock waves in the space of thermodynamic variables. Near critical points the volume density exhibits a scaling behavior consistent with the Universality Conjecture in viscous transport PDEs. The same construction yields a global QCD phase diagram containing critical points for the nuclear liquid-gas transition and the hadron gas-quark-gluon plasma transition.

What carries the argument

Nonlinear C-integrable partial differential equations of hydrodynamic type that encode the virial expansion and supply the expectation values of thermodynamic observables for any finite N.

If this is right

  • Expectation values of observables remain computable exactly at arbitrary order in the virial expansion for any fixed particle number.
  • Phase transitions in the infinite-particle limit take the concrete form of shock discontinuities in thermodynamic variables.
  • Volume density obeys a definite scaling law near critical points that matches the Universality Conjecture for viscous transport equations.
  • A single global phase diagram can be drawn for nuclear and quark matter that places both the liquid-gas and hadron-quark-gluon critical points.
  • Finite particle number systematically rounds the critical signatures, altering their visibility in finite systems.

Where Pith is reading between the lines

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

  • Numerical integration of the finite-N PDEs could serve as a practical route to approximate the infinite-volume thermodynamics without summing the full virial series.
  • Heavy-ion collision experiments may need to vary system size explicitly to map how the apparent location of the QCD critical point shifts with N.
  • The same hydrodynamic integrability might extend to other statistical ensembles whose equations of state admit a virial-like density expansion.

Load-bearing premise

The system of hydrodynamic equations obtained from the virial expansion of the internal energy is C-integrable at every order in the density expansion.

What would settle it

An explicit low-order truncation of the virial series whose resulting hydrodynamic PDE fails to be C-integrable, or an experimental measurement in heavy-ion collisions showing that finite-size smearing of critical signatures is absent.

Figures

Figures reproduced from arXiv: 2503.15719 by Francesco Giglio, Giulio Landolfi, Xin An.

Figure 1
Figure 1. Figure 1: FIG. 1. Universal behavior and finite-size corrections in [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Isothermal curves predicted by the EOS ( [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
read the original abstract

We analyze thermodynamic models for fluid systems in equilibrium based on a virial expansion of the internal energy in terms of the volume density. We prove that the models, formulated for finite-size systems with $N$ particles, are exactly solvable to any expansion order, as expectation values of physical observables (e.g., volume density) are determined from solutions to nonlinear C-integrable partial differential equations (PDEs) of hydrodynamic type. In the limit $N\to \infty$, phase transitions emerge as classical shock waves in the space of thermodynamic variables. Near critical points, we argue that the volume density exhibits a scaling behavior consistent with the Universality Conjecture in viscous transport PDEs. As an application, we employ our framework to nuclear and quark matter, constructing a global quantum chromodynamics (QCD) phase diagram that reveals critical points for the nuclear liquid-gas transition and the hadron gas-quark-gluon plasma transition. We demonstrate how finite-size effects smear critical signatures, implying their potential impact on the search for the QCD critical point.

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 / 1 minor

Summary. The manuscript analyzes thermodynamic models based on a virial expansion of internal energy in volume density for finite-N particle systems. It claims to prove that these models are exactly solvable to any expansion order because expectation values of observables are obtained from solutions of nonlinear C-integrable PDEs of hydrodynamic type; in the N→∞ limit phase transitions appear as classical shock waves, with an application to constructing a global QCD phase diagram that includes critical points for nuclear liquid-gas and hadron-QGP transitions while noting finite-size smearing of signatures.

Significance. If the asserted C-integrability at arbitrary virial order can be established with a general structure, the framework would supply a new class of exactly solvable finite-size statistical models that connect virial expansions to hydrodynamic integrability and shock-wave descriptions of phase transitions, offering a concrete route to quantitative predictions for systems such as nuclear matter.

major comments (2)
  1. [derivation of the hydrodynamic PDEs from the virial expansion] The exact-solvability claim to arbitrary virial order rests on the hydrodynamic equations remaining C-integrable for any truncation of the virial series. The manuscript invokes this property directly from the form of the equations but does not exhibit an order-independent structure (Lax pair, recursion operator, or infinite family of conservation laws) that survives arbitrary coefficients; without such a demonstration the reduction of observables to solutions of the PDE system does not follow at every order.
  2. [statement of exact solvability for finite N] The finite-N exact solvability is stated to hold because expectation values are determined from the PDE solutions, yet no explicit error estimates or verification that the truncation error remains controlled under the C-integrability assumption are supplied; this step is load-bearing for the central assertion that the models are exactly solvable at any expansion order.
minor comments (1)
  1. [discussion near critical points] The Universality Conjecture for viscous transport PDEs is invoked near critical points without a specific reference or precise statement of the conjecture being used.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We appreciate the recognition of the potential significance of connecting virial expansions to C-integrable hydrodynamics. We address each major comment below and indicate where revisions will be made for clarification.

read point-by-point responses

  1. Referee: [derivation of the hydrodynamic PDEs from the virial expansion] The exact-solvability claim to arbitrary virial order rests on the hydrodynamic equations remaining C-integrable for any truncation of the virial series. The manuscript invokes this property directly from the form of the equations but does not exhibit an order-independent structure (Lax pair, recursion operator, or infinite family of conservation laws) that survives arbitrary coefficients; without such a demonstration the reduction of observables to solutions of the PDE system does not follow at every order.

    Authors: The PDEs are obtained directly by computing the time derivatives of the relevant moments from the finite-N virial expansion of the internal energy, yielding a closed quasilinear system of hydrodynamic type whose coefficients depend on the chosen truncation but whose structural form (conservation-law structure with velocity-dependent fluxes) remains identical at every order. Systems of this hydrodynamic type are C-integrable via the hodograph method or Riemann invariants whenever they are hyperbolic, a property independent of the specific numerical values of the coefficients. We will add a brief paragraph and reference to the general theory of C-integrable hydrodynamic systems to make this order-independent structural argument explicit. revision: partial

  2. Referee: [statement of exact solvability for finite N] The finite-N exact solvability is stated to hold because expectation values are determined from the PDE solutions, yet no explicit error estimates or verification that the truncation error remains controlled under the C-integrability assumption are supplied; this step is load-bearing for the central assertion that the models are exactly solvable at any expansion order.

    Authors: For any fixed truncation order the virial model is defined by that finite set of coefficients; the resulting PDE system is then closed and exact for the finite-N dynamics, so the expectation values are given precisely by the PDE solutions. The truncation error relative to a hypothetical infinite virial series is an approximation inherent to the model definition itself rather than an additional error inside the solvability claim. We will insert a clarifying sentence in the introduction and methods section to distinguish the exact solvability of the truncated model from the separate question of virial-series convergence. revision: partial

Circularity Check

0 steps flagged

No circularity: solvability follows from asserted C-integrability of virial-derived hydrodynamic PDEs without reduction to fitted inputs or self-citation chains

full rationale

The paper's central claim is that finite-N models based on virial expansion of internal energy yield expectation values determined by solutions of C-integrable hydrodynamic PDEs, with the N→∞ limit producing shock waves. No step in the provided abstract or described derivation reduces a result to its own inputs by construction (e.g., no fitted parameter renamed as prediction, no self-definitional loop, and no load-bearing self-citation of an unverified uniqueness theorem). The integrability property is invoked as following from the structure of the virial equations, but the abstract presents this as a proof rather than a tautology; absent any quoted equation that equates the output directly to the input by definition, the derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the assumption that the virial expansion produces C-integrable PDEs at every order and that the hydrodynamic limit correctly captures thermodynamic phase transitions; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption The virial expansion of internal energy in volume density yields a closed system of hydrodynamic PDEs that remain C-integrable to arbitrary order.
    Invoked to establish exact solvability for finite N.
  • domain assumption Classical shock waves in the hydrodynamic equations correspond to thermodynamic phase transitions in the N to infinity limit.
    Required for the phase-transition interpretation.

pith-pipeline@v0.9.0 · 5718 in / 1498 out tokens · 30695 ms · 2026-05-22T23:38:03.377107+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages · 14 internal anchors

  1. [1]

    Parisi, Nobel lecture: Multiple equilibria , Rev

    G. Parisi, Nobel lecture: Multiple equilibria , Rev. Mod. Phys. 95 (2023) 030501

    work page 2023

  2. [2]

    B.J. ˙I. Beycan ˙Ibrahimo˘ glu,Critical States at Phase Transitions of Pure Substances, Springer (2022)

    work page 2022

  3. [3]

    Cheng, Phase Transitions in Polymers: The Role of Metastable States , Elsevier (2008)

    S.Z. Cheng, Phase Transitions in Polymers: The Role of Metastable States , Elsevier (2008)

    work page 2008

  4. [4]

    Pelissetto and E

    A. Pelissetto and E. Vicari, Critical phenomena and renormalization-group theory, Phys. Rep. 368 (2002) 549

    work page 2002

  5. [5]

    Combescot and C.B

    M. Combescot and C.B. ` a la Guillaume,Two critical points for an electron-hole system , Phys. Rev. Lett. 44 (1980) 182

    work page 1980

  6. [6]

    Mukherjee, Emergence of a bicritical end point in the random-crystal-field blume-capel model, Phys

    Sumedha and S. Mukherjee, Emergence of a bicritical end point in the random-crystal-field blume-capel model, Phys. Rev. E 101 (2020) 042125 [1907.06454]

    work page arXiv 2020

  7. [7]

    Maritan, M

    A. Maritan, M. Cieplak, M.R. Swift, F. Toigo and J.R. Banavar, Random anisotropy Blume-Emery-Griffiths model, Phys. Rev. Lett. 69 (1992) 221

    work page 1992

  8. [8]

    F.S. P.G. Debenedetti and G. Zerze, Second critical point in two realistic models of water , Science 369 (2020) 289

    work page 2020

  9. [9]

    Gartner, P.M

    T.E. Gartner, P.M. Piaggi, R. Car, A.Z. Panagiotopoulos and P.G. Debenedetti, Liquid-liquid transition in water from first principles, Phys. Rev. Lett. 129 (2022) 255702 [2208.13633]

    work page arXiv 2022

  10. [10]

    Multiple Critical Points in Effective Quark Models

    L. Ferroni, V. Koch and M.B. Pinto, Multiple critical points in effective quark models , Phys. Rev. C 82 (2010) 055205 [ 1007.4721]

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  11. [11]

    Multiple critical end points in magnetized three flavor quark matter

    M. Ferreira, P. Costa and C. Providˆ encia,Multiple critical endpoints in magnetized three flavor quark matter, Phys. Rev. D 97 (2018) 014014 [1712.08378]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  12. [12]

    Shock dynamics of phase diagrams

    A. Moro, Shock dynamics of phase diagrams , Ann. Phys. 343 (2014) 49 [ 1307.7512]

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  13. [13]

    Exact solution of the van der Waals model in the critical region

    A. Barra and A. Moro, Exact solution of the van der Waals model in the critical region , Ann. Phys. 359 (2015) 290 [ 1412.1951]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  14. [14]

    Yang and T.D

    C.N. Yang and T.D. Lee, Statistical theory of equations of state and phase transitions. i. theory of condensation, Phys. Rev. 87 (1952) 404

    work page 1952

  15. [15]

    Ruelle, Statistical Mechanics, World Scientific, Singapore (1999)

    D. Ruelle, Statistical Mechanics, World Scientific, Singapore (1999)

    work page 1999

  16. [16]

    Classification of phase transitions in small systems

    P. Borrmann, O. M¨ ulken and J. Harting, Classification of phase transitions in small systems , Phys. Rev. Lett. 84 (2000) 3511 [cond-mat/9909184]

    work page internal anchor Pith review Pith/arXiv arXiv 2000

  17. [17]

    Bernhardt, C.S

    J. Bernhardt, C.S. Fischer, P. Isserstedt and B.-J. Schaefer, Critical endpoint of QCD in a finite volume, Phys. Rev. D 104 (2021) 074035 [2107.05504]

    work page arXiv 2021

  18. [18]

    Kov´ acs, P

    G. Kov´ acs, P. Kov´ acs, P.M. Lo, K. Redlich and G. Wolf, Sensitivity of finite size effects to the boundary conditions and the vacuum term , Phys. Rev. D 108 (2023) 076010 [ 2307.10301]

    work page arXiv 2023

  19. [19]

    Lorenzoni and A

    P. Lorenzoni and A. Moro, Exact analysis of phase transitions in mean-field potts models , Phys. Rev. E 100 (2019) 022103

    work page 2019

  20. [20]

    Callen, Thermodynamics and an Introduction to Thermostatistics, Wiley, Singapore (1991)

    H. Callen, Thermodynamics and an Introduction to Thermostatistics, Wiley, Singapore (1991)

    work page 1991

  21. [21]

    Charged AdS Black Holes and Catastrophic Holography

    A. Chamblin, R. Emparan, C. Johnson and R. Myers, Charged AdS black holes and catastrophic holography, Phys. Rev. D 60 (1999) 064018 [hep-th/9902170]

    work page internal anchor Pith review Pith/arXiv arXiv 1999

  22. [22]

    Oertel, M

    M. Oertel, M. Hempel, T. Kl¨ ahn and S. Typel, Equations of state for supernovae and compact stars, Rev. Mod. Phys. 89 (2017) 015007

    work page 2017

  23. [23]

    C. Chin, R. Grimm, P. Julienne and E. Tiesinga, Feshbach resonances in ultracold gases, Rev. Mod. Phys. 82 (2010) 1225–1286 [ 0812.1496]

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  24. [24]

    Kapusta and C

    J.I. Kapusta and C. Gale, Finite-temperature field theory: Principles and applications , Cambridge Monographs on Mathematical Physics, Cambridge University Press (2011)

    work page 2011

  25. [25]

    Le Bellac, F

    M. Le Bellac, F. Mortessagne and G. Batrouni, Equilibrium and Non-Equilibrium Statistical Thermodynamics, Cambridge University Press (2006)

    work page 2006

  26. [26]

    Exact equations of state for nematics

    G. De Matteis, F. Giglio and A. Moro, Exact equations of state for nematics , Ann. Phys. 396 (2018) 386 [ 1802.07246]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  27. [27]

    Calogero, Why are certain nonlinear PDEs both widely applicable and integrable? , in What Is Integrability?, (Berlin), pp

    F. Calogero, Why are certain nonlinear PDEs both widely applicable and integrable? , in What Is Integrability?, (Berlin), pp. 1–62, Springer (1991)

    work page 1991

  28. [28]

    G.B.Whitham, Linear and Nonlinear Waves , John Wiley & Sons, Inc. (1973)

    work page 1973

  29. [29]

    Thom, Structural Stability And Morphogenesis , 7 CRC Press (1989)

    R. Thom, Structural Stability And Morphogenesis , 7 CRC Press (1989)

    work page 1989

  30. [30]

    Integrable extended van der Waals model

    F. Giglio, G. Landolfi and A. Moro, Integrable extended van der Waals model , Physica D: Nonlinear Phenomena 333 (2016) 293 [1602.03975]

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  31. [31]

    Bressan, Hyperbolic systems of conservation laws: the one-dimensional Cauchy problem , Oxford University Press, USA (2000)

    A. Bressan, Hyperbolic systems of conservation laws: the one-dimensional Cauchy problem , Oxford University Press, USA (2000)

    work page 2000

  32. [32]

    On Hamiltonian perturbations of hyperbolic systems of conservation laws, II: universality of critical behaviour

    B. Dubrovin, On hamiltonian perturbations of hyperbolic systems of conservation laws, ii: Universality of critical behaviour , Comm. Math. Phys. 267 (2006) 117 [ math-ph/0510032]

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  33. [33]

    On critical behavior in nonlinear evolutionary PDEs with small viscosity

    B. Dubrovin and M. Elaeva, On the critical behavior in nonlinear evolutionary PDEs with small viscosity, Russ. J. Math. Phys. 19 (2012) 449 [1301.7216]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  34. [34]

    QCD phase diagram: an overview

    M. Stephanov, QCD phase diagram: An Overview , PoS LA T2006(2006) 024 [ hep-lat/0701002]

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  35. [35]

    Sorensen et al., Dense nuclear matter equation of state from heavy-ion collisions, Prog

    A. Sorensen et al., Dense nuclear matter equation of state from heavy-ion collisions , Prog. Part. Nucl. Phys. 134 (2024) 104080 [ 2301.13253]

    work page arXiv 2024

  36. [36]

    Pochodzalla, T

    J. Pochodzalla, T. M¨ ohlenkamp, T. Rubehn, A. Sch¨ uttauf, A. W¨ orner, E. Zude et al.,Probing the nuclear liquid-gas phase transition , Phys. Rev. Lett. 75 (1995) 1040

    work page 1995

  37. [37]

    Wang, Y.-G

    R. Wang, Y.-G. Ma, R. Wada, L.-W. Chen, W.-B. He, H.-L. Liu et al., Nuclear liquid-gas phase transition with machine learning , Phys. Rev. Res. 2 (2020) 043202 [ 2010.15043]

    work page arXiv 2020

  38. [38]

    Sorensen and V

    A. Sorensen and V. Koch, Phase transitions and critical behavior in hadronic transport with a relativistic density functional equation of state , Phys. Rev. C 104 (2021) 034904 [ 2011.06635]

    work page arXiv 2021

  39. [39]

    Omana Kuttan, J

    M. Omana Kuttan, J. Steinheimer, K. Zhou and H. Stoecker, QCD equation of state of dense nuclear matter from a bayesian analysis of heavy-ion collision data , Phys. Rev. Lett. 131 (2023) 202303 [ 2211.11670]

    work page arXiv 2023

  40. [40]

    The phase diagram of dense QCD

    K. Fukushima and T. Hatsuda, The phase diagram of dense QCD , Rep. Prog. Phys. 74 (2010) 014001 [1005.4814]

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  41. [41]

    L. Xu, P. Kumar, S.V. Buldyrev, S.-H. Chen, P.H. Poole, F. Sciortino et al., Relation between the widom line and the dynamic crossover in systems with a liquid–liquid phase transition , Proceedings of the National Academy of Sciences 102 (2005) 16558

    work page 2005

  42. [42]

    Guenther, Overview of the QCD phase diagram: Recent progress from the lattice , Eur

    J.N. Guenther, Overview of the QCD phase diagram: Recent progress from the lattice , Eur. Phys. J. A 57 (2021) 136 [ 2010.15503]

    work page arXiv 2021

  43. [43]

    STAR collaboration, Nonmonotonic Energy Dependence of Net-Proton Number Fluctuations , Phys. Rev. Lett. 126 (2021) 092301 [ 2001.02852]

    work page arXiv 2021

  44. [44]

    NA61/SHINE collaboration, Search for the QCD critical point by NA61/SHINE at the CERN SPS , in 10th International Conference on New Frontiers in Physics, 8, 2023 [ 2308.04254]

    work page arXiv 2023

  45. [45]

    X. An, D. Mesterh´ azy and M.A. Stephanov,On spinodal points and Lee-Yang edge singularities , J. Stat. Mech. 1803 (2018) 033207 [ 1707.06447]

    work page internal anchor Pith review Pith/arXiv arXiv 2018