Expansion schemes for finite density QCD equation of state: which one is the best?
Abstract
The finite density equation of state of QCD is a fundamental quantity for understanding the phase structure of the theory, and also for the realistic modeling of the QCD medium in simulations. While direct lattice simulations at finite density are extremely challenging, several expansion schemes have been devised to extrapolate into the QCD phase diagram. All of these schemes are based on susceptibilities of conserved charges, i.e. derivatives of the pressure with respect to the chemical potential(s). In this work we set out to determine which of these schemes possesses the best convergence properties. To this end, we employ extreme-statistics lattice simulations on a 163×8 lattice, which allow us to obtain susceptibilities up to order eight with remarkably small errors. We compare in particular two methods, based on contours of constant (normalized) baryon density and of constant entropy density. We find that the latter converges more quickly, effectively requiring susceptibilities only up to order four to cover the chemical potential range up to at least μB ≈ 500 MeV. We also extend the comparison to the three dimensional space of chemical potentials, and find that the entropy-based method consistently shows superior convergence.
Related
- Quantum HHL Algorithm for Solving 1D Linear Advection Equation Using Explicit and Implicit Finite Difference Schemes
- NJL-Chiral Soliton and the Nucleon Equation of State at supra-saturation density: Impact of Chiral Symmetry Restoration
- Equation of state and transport coefficients of warm dense aluminum from mixed deterministic-stochastic density functional theory
- Nuclear equation-of-state at high density and multi-messenger astronomy: contribution of heavy-ion collisions
- Short-range correlated pairs from nucleon density profiles as a new probe of the nuclear equation of state
BibTeX
@misc{borsanyi2026expansion,
title = {{Expansion schemes for finite density QCD equation of state: which one is the best?}},
author = {Szabolcs Borsányi and Zoltán Fodor and Jana N. Guenther and Piyush Kumar and Paolo Parotto and Attila Pásztor and Chik Him Wong},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.25494},
eprint = {2610.25494},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.25494},
}
RIS
TY - UNPB
TI - Expansion schemes for finite density QCD equation of state: which one is the best?
AU - Szabolcs Borsányi
AU - Zoltán Fodor
AU - Jana N. Guenther
AU - Piyush Kumar
AU - Paolo Parotto
AU - Attila Pásztor
AU - Chik Him Wong
PY - 2026
DA - 2026/10/18
PB - arXiv
DO - 10.48550/arxiv.2610.25494
UR - https://arxiv.org/abs/2610.25494
AB - The finite density equation of state of QCD is a fundamental quantity for understanding the phase structure of the theory, and also for the realistic modeling of the QCD medium in simulations. While direct lattice simulations at finite density are extremely challenging, several expansion schemes have been devised to extrapolate into the QCD phase diagram. All of these schemes are based on susceptibilities of conserved charges, i.e. derivatives of the pressure with respect to the chemical potential(s). In this work we set out to determine which of these schemes possesses the best convergence properties. To this end, we employ extreme-statistics lattice simulations on a 16^3×8 lattice, which allow us to obtain susceptibilities up to order eight with remarkably small errors. We compare in particular two methods, based on contours of constant (normalized) baryon density and of constant entropy density. We find that the latter converges more quickly, effectively requiring susceptibilities only up to order four to cover the chemical potential range up to at least μ_B ≈ 500 MeV. We also extend the comparison to the three dimensional space of chemical potentials, and find that the entropy-based method consistently shows superior convergence.
ER -