---
title: Polynomial Equation of State
url: https://www.emergentmind.com/topics/polynomial-equation-of-state
type: topic
---

# Polynomial Equation of State

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 \) [1407.8532, 1202.3104], molar volume, or redshift filter functions [2510.21892].
- As a polynomial in physical quantities (e.g., energy, pressure-volume):
  $$
  P = \sum_{n=1}^{N} \alpha_n (PV)^n
  $$
  [1612.05559].

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 [2201.04009] and virial expansions [2010.14357].

## 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 [1407.8532, 1202.3104].

**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 [1407.8532].

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 [1202.3104].

## 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 [2201.04009].

- 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 [1612.05559]. 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 [2510.21892].  
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 [2010.14357].  
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 [1909.09899]. 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 [1612.05559].
- Systematic convergence: Direct control over truncation order and error estimates [1407.8532].

**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 \) [2010.14357].
- 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 [1612.05559].
- Application to mixtures, anisotropic matter, composite fluids, or gravity-modified gravity/cosmology contexts [2510.21892, 1909.09899].

---

**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 \)           | [1407.8532]       |
| Real gas (van der Waals)      | \( V \)                     | \( pV^3 - (RT + b p)V^2 + a V - a b = 0 \)             | [2201.04009]      |
| Condensed Matter (solids)     | \( P V \)                   | \( P = \alpha_1 (P V) + \alpha_2 (PV)^2 + \alpha_3 (PV)^3 \) | [1612.05559]      |
| Virial expansions (hard core) | \( n \) (virial index)      | \( B_n = \sum_{k=0}^{D-1} a_k n^k \)                   | [2010.14357]      |
| Dark energy EoS               | Filtered z (\( x(z) \))     | \( w(z) = -1 + c_1 g_1[x(z)] + c_2 g_2[x(z)] \)        | [2510.21892]      |
| Exotic matter (GR wormholes)  | \( \rho \)                  | \( p(\rho) = \omega \rho + \sum \omega_i \rho^{n_i} \) | [1909.09899]      |

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.

Source: https://www.emergentmind.com/topics/polynomial-equation-of-state