Polynomial Equation of State

• Polynomial Equation of State is a thermodynamic framework that expresses pressure and energy density as polynomial functions of key physical parameters.
• It facilitates accurate modeling in systems ranging from finite-density QCD and classical fluids to solid-state physics, with validations through lattice simulations and experimental data.
• Its flexible, analytic form extends to cosmological dark energy and exotic matter in general relativity, offering transparent parameter control and systematic convergence.

A polynomial equation of state (EoS) is a thermodynamic model expressing pressure, energy density, or fluid variables as polynomials in relevant physical parameters, such as chemical potential, density, or composite functionals. Such polynomial representations span a range of contexts: from lattice QCD at finite baryon density, to classical fluids, condensed-matter systems, cosmological dark energy, and exotic matter in general relativity. Their utility lies in analytic tractability, systematic convergence, empirical flexibility, and the transparent incorporation of physical constraints and observational data.

1. General Theory and Formal Structure

The essence of a polynomial EoS is to represent the pressure or other thermodynamic potential as a finite polynomial in one or more variables characteristic of the system. The general form encountered in quantum chromodynamics (QCD), condensed matter, and cosmology is

• As a Taylor expansion (QCD/statistical physics):

$P(T, X) = \sum_{n=0}^N c_n(T)\,X^n$

where $X$ may be the baryon chemical potential $\mu_B/T$ (Huovinen et al., 2014, Huovinen et al., 2012), molar volume, or redshift filter functions (Rudra et al., 24 Oct 2025).

• As a polynomial in physical quantities (e.g., energy, pressure-volume):

$P = \sum_{n=1}^{N} \alpha_n (PV)^n$

(Gordon et al., 2016).

Polynomials are preferred when physical observables demonstrate analytic or quasi-analytic dependence on the expansion variable, or when systematic improvements are needed up to a desired order of accuracy. In classical fluids, cubic and quartic polynomial roots naturally result from physical approximations, as in the van der Waals equation (Prodanov, 2022) and virial expansions (Liu, 2020).

2. Polynomial EoS in QCD and Hot/Dense Matter

The polynomial EoS is a foundational tool for modeling finite-density QCD matter as explored via lattice gauge theory. In this context, the pressure at temperature $T$ and baryon chemical potential $\mu_B$ is expanded in powers of $\mu_B/T$:

$$\frac{P(T, \mu_B)}{T^4} = \sum_{n=0}^6 c_n(T)\left(\frac{\mu_B}{T}\right)^n$$

where the coefficients $c_n(T)$ are Taylor coefficients computed from lattice simulations. The approach is constrained by charge conjugation symmetry ($c_1 = 0$), and even powers dominate (Huovinen et al., 2014, Huovinen et al., 2012).

Key features:

• Second-order coefficients ($c_2$) are obtained via continuum-extrapolated lattice actions (stout, HISQ).
• Fourth- and sixth-order ($c_4, c_6$) coefficients require careful treatment of lattice artefacts—corrected by a temperature shift to match the hadron resonance gas (HRG) in the confined regime.
• Each $c_n(T)$ is parametrize by an “inverse-polynomial plus Stefan–Boltzmann” form:

$$c_n(T) = c_n^{\rm SB} + \sum_{k=1}^{m_n} \frac{a_{n,k}}{\left[\widehat T\right]^{p_{n,k}}}$$

with appropriately chosen shifts and scaling for HRG matching.

• The series converges for physical scenarios with entropy-to-baryon ratio $s/n_B \gtrsim 40$, i.e., essentially all heavy-ion collision scenarios above $\sqrt{s_{NN}}\gtrsim 17$ GeV (Huovinen et al., 2014).

Once the pressure is known, all other thermodynamic quantities are derived by differentiation:

$$n_B = \frac{\partial P}{\partial\mu_B}, \quad s = \frac{\partial P}{\partial T}, \quad \varepsilon = -P + Ts + \mu_B n_B$$

This EoS is matched to the HRG at low $T$ by enforcing the continuity of the $c_n$ and their derivatives at switching temperatures (e.g., $T_{\rm SW} = 155\text{–}160$ MeV) for thermodynamic consistency (Huovinen et al., 2012).

3. Classical Fluid and Solid-State Applications

In classical thermodynamics, polynomial equations of state often emerge as:

• Cubic polynomial form for real gases (van der Waals EoS):

$pV^3 - (RT + b p)V^2 + a V - a b = 0$

where $a, b, R$ are material constants, $p$ is pressure, $V$ is molar volume, and $T$ is temperature (Prodanov, 2022).

• For condensed matter (empirical polynomial EOS):

$$P = \alpha_1 (P V) + \alpha_2 (P V)^2 + \alpha_3 (P V)^3$$

with material-dependent fitting coefficients, achieving sub-percent accuracy across broad pressure ranges and multiple phase transitions (Gordon et al., 2016). This approach leverages the observation that, for solids, $P$ as a function of $P V$ yields a near-linear trend with minor curvature, justifying low-order truncation.

The fitting procedure involves ordinary least squares on experimental $P$–$V$ data. Notably, the polynomial $PV$ model regularly outperforms conventional forms (Birch–Murnaghan, Vinet) in phase-transitional and high-pressure regimes.

Table: Comparison of Polynomial EoS Accuracy in Solids

Model	System	Max % Error
Cubic PV-polynomial	MgO	0.6
Birch–Murnaghan (3rd order)	MgO	∼1.0
Vinet	MgO	∼0.8
Cubic PV-polynomial	Te	5.5

4. Polynomial Parametrizations in Cosmology

Modern cosmology utilizes polynomial EoS forms to encode the dark energy equation of state as a function of redshift, $z$, or suitable surrogates to avoid divergences at $z \to -1$. A recently developed class parametrizes $w(z)$, the ratio $p/\rho$, via smooth redshift filters:

$$w(z) = -1 + c_1\,g_1[x(z)] + c_2\,g_2[x(z)], \quad x(z) = 1+\frac{z}{1+z^2}$$

where $g_1, g_2$ are polynomial basis functions (e.g., monomials, Legendre, Chebyshev) and $c_1, c_2$ are fit to current cosmological data (Rudra et al., 24 Oct 2025). Features include:

• Regularity for all $z \in [-1, \infty)$; no future singularities.
• Performance benchmarking via MCMC likelihood fits (Hubble, BAO, DESI) and algorithmic regression (ANN, SVR), with the Legendre polynomial basis statistically optimal under reduced $\chi^2$ and RMSE.
• Example best-fit for Legendre basis: $w_0 \approx -0.97 \pm 0.06$, compatible with $\Lambda$CDM at $1\sigma$.

5. Polynomial-Based Virial Expansions and Carnahan–Starling-Type EoS

In the statistical mechanics of hard-particle fluids, the polynomial ansatz extends to the virial expansion, expressing the virial coefficients $B_n$ as polynomials in $n$ (order):

$$B_n = a_0 + a_1 n + \cdots + a_{D-1} n^{D-1}, \quad n \geq 4$$

which, summed analytically, yields Carnahan–Starling-type EoS:

$Z_D(\eta) = \frac{N_D(\eta)}{(1-\eta)^D}$

where $Z = p/(k_B T \rho)$ is the compressibility factor, $\eta$ is the packing fraction, and $N_D(\eta)$ is a system-specific polynomial (Liu, 2020). For 3D hard spheres:

$$Z_{D=3}(\eta) = \frac{1 + \eta + \eta^2 - 0.6155\eta^3 - 1.000\eta^4 + 0.500\eta^5}{(1-\eta)^3}$$

Such representations attain AAD (average-absolute deviation) in $Z$ below 0.004% over the equilibrium fluid range, paralleling Padé approximants and exceeding the original Carnahan–Starling form in accuracy.

A universal consequence of the polynomial $B_n$ ansatz is a pole at $\eta = 1$, irrespective of the close-packing fraction, reflecting the limits of polynomial representations at extremely high densities.

6. Polynomial EoS in General Relativity and Exotic Fluids

In relativistic wormhole and exotic matter modeling, polynomial EoS facilitate the exploration of fluids violating the standard energy conditions. The general mixed (polynomial) EoS:

$$p(\rho) = \omega \rho + \sum_{i=1}^N \omega_i \rho^{n_i}, \quad n_i > 1$$

expands around the usual barotropic term, with higher-order “polytropic” corrections (Parsaei et al., 2019). In such geometries, this structure allows:

• Localization of null energy condition (NEC) violation to arbitrary small zones near the wormhole throat.
• Parameter control over the extent of exotic matter, redshift profile, and asymptotic matching to dark energy ($w(\infty) \to \omega$, $-1 < \omega < -1/3$).

In particular, for quadratic correction ($N=1, n=2$), the radius of the exotic region scales analytically with EoS parameters and the geometry exponent $\alpha$.

7. Strengths, Limitations, and Extensions

Strengths:

• Analytical tractability: All thermodynamic derivatives follow directly.
• Empirical flexibility: Fit systematically to high-order data across phases and temperature/pressure/density regimes (Gordon et al., 2016).
• Systematic convergence: Direct control over truncation order and error estimates (Huovinen et al., 2014).

Limitations:

• Breakdown at critical points or near physical singularities; e.g., radius of convergence in $\mu_B$ or divergence at packing fraction $\eta \to 1$ (Liu, 2020).
• Physical content depends on the polynomial variable choice; care is needed when extending to non-analytic or highly nonlinear regimes.

Extensions:

• Inclusion of temperature or other thermodynamic variables as arguments or parametric coefficients; e.g., $\alpha_i(T)$ in the solid-state EOS (Gordon et al., 2016).
• Application to mixtures, anisotropic matter, composite fluids, or gravity-modified gravity/cosmology contexts (Rudra et al., 24 Oct 2025, Parsaei et al., 2019).

---

Summary Table: Representative Polynomial Equations of State

Physical System	Polynomial Variable	Canonical Form/Truncation	Reference
Hot QCD (finite μ_B)	$\mu_B/T$	$P/T^4 = \sum_{n=0}^6 c_n(T)(\mu_B/T)^n$	(Huovinen et al., 2014)
Real gas (van der Waals)	$V$	$pV^3 - (RT + b p)V^2 + a V - a b = 0$	(Prodanov, 2022)
Condensed Matter (solids)	$P V$	$$P = \alpha_1 (P V) + \alpha_2 (PV)^2 + \alpha_3 (PV)^3$$	(Gordon et al., 2016)
Virial expansions (hard core)	$n$ (virial index)	$B_n = \sum_{k=0}^{D-1} a_k n^k$	(Liu, 2020)
Dark energy EoS	Filtered z ($x(z)$)	$w(z) = -1 + c_1 g_1[x(z)] + c_2 g_2[x(z)]$	(Rudra et al., 24 Oct 2025)
Exotic matter (GR wormholes)	$\rho$	$p(\rho) = \omega \rho + \sum \omega_i \rho^{n_i}$	(Parsaei et al., 2019)

Polynomial equations of state thus provide a unifying framework for modeling thermodynamic, cosmological, and gravitational phenomena across a spectrum of physical regimes. Their efficacy is rooted in the balance between analytical transparency and empirical adequacy, as demonstrated in the most advanced current research across high-energy, condensed matter, and cosmological systems.