Dense Matter and Equations of State

An equation of state connects microscopic matter to relativistic stress, wave propagation, and stellar structure. Its values, derivatives, phase construction, composition, and uncertainty must all be consistent. This chapter begins with thermodynamic potentials, develops the local expansion around nuclear saturation, treats possible phase changes and hot composition, imposes causal high-density bounds, and finally maps laboratory and astronomical evidence to the density domains it can actually inform.

Helpful background. Hilbert Stress–Energy from Metric Variation supplies pressure and energy density as stress components; The Einstein Equation maps the equation of state into stellar gravity; and Relativistic Matter and Compact Objects gives the volume map.

Thermodynamics before interpolation

Thermodynamic Consistency of Equations of State starts from a Helmholtz free-energy density $\mathcal F(T,n_i)$. Entropy, chemical potentials, pressure, and energy follow by differentiation and Legendre transformation. The first law and Maxwell identities then become nontrivial checks, not independent definitions that may drift apart.

This matters most for tables. Interpolating pressure, energy, and entropy independently can preserve smooth values while violating derivatives between nodes. The page specifies analytic classical, degenerate, and radiation benchmarks; nested axis refinement; first-law and Maxwell residuals; and a versioned table environment. The CompOSE data model and interpolation conventions are described in CompOSE Core Team 2022, Reference Manual v3.01, §§3.5–3.6, §4.1.2, and Appendix A.1 (PDF).

Stability is multidimensional. Positive heat capacity, compressibility, and composition susceptibilities constrain different directions of the free-energy Hessian. A beta-equilibrated one-dimensional slice can hide an unstable composition mode, and a smooth spline across a phase boundary can manufacture a derivative that no equilibrium state possesses.

The local expansion around saturation

Nuclear Matter Near Saturation organizes cold matter by density displacement $x=(n-n_0)/(3n_0)$ and isospin asymmetry $\delta=1-2Y_p$. Symmetric matter has a minimum at $n_0$, so its pressure vanishes there and its leading compression term is quadratic. The symmetry energy begins with $J$ and slope $L$; near $n_0$, neutron-rich pressure is strongly correlated with $L$ under the quadratic-asymmetry approximation.

These coefficients are local derivatives and correlated inference products. Nuclear masses, collective modes, weak form factors, neutron skins, and microscopic neutron matter constrain different combinations. Chiral effective field theory supplies controlled nucleonic ensembles near saturation with widening uncertainty toward higher density Drischler, Holt, and Wellenhofer 2021, §§ 2–4. It does not justify extrapolating a low-order saturation polynomial through several times $n_0$ or across new degrees of freedom.

Beta equilibrium adds chemical-potential and charge-neutrality constraints. With trapped neutrinos, the equilibrium path differs from cold transparent matter. Effective masses and composition can change thermal and transport behavior even when two models have similar cold pressure.

Phase structure and hot composition

Phase Transitions and Hybrid Matter distinguishes a continuous crossover from first-order coexistence. Maxwell construction assumes locally neutral bulk phases separated sharply; Gibbs construction permits globally neutral charged mixtures. Surface tension, Coulomb screening, domain geometry, and conversion kinetics decide which bulk idealization is closer to a star. A large energy-density jump can destabilize the onset of a dense core, while sufficiently stiff high-density matter can restore stability and produce a third family. Oertel and collaborators review exotic constituents and phase constructions in Oertel et al. 2017, §§ III.D–IV.

The frontier synthesis is dated through 2026-08-18. Low-density theory and asymptotic QCD leave a wide middle-density gap; laboratory collisions probe different temperature and composition; present stellar observations restrict pressure but do not uniquely identify hyperons, quarks, crossover, or a strong first-order transition. Claims of twin stars therefore require stability, prior, and multimessenger checks rather than one mass–radius contour.

Finite Temperature and Composition replaces a cold barotrope with a state space such as $(n_B,T,Y_e)$. Nuclear statistical equilibrium, dissolution of nuclei, leptons, photons, and neutrino trapping must be included or excluded explicitly. The page checks $T\to0$, dilute ideal gas, radiation, and degenerate-Fermi limits, then tests table derivatives and composition sums. A table that already contains electrons or photons cannot be combined with a module that adds them again.

Causal completion and extremal bounds

Causality, Stability, and Maximum Stiffness uses

$$0\leq\frac{dp}{d\varepsilon}\leq1$$

for a cold perfect-fluid barotrope. The lower bound prevents local mechanical instability; the upper bound keeps the sound cone causal. It is distinct from the dominant energy condition and from full multidimensional thermodynamic convexity.

Above a matching state, a constant-sound-speed continuation

$$p=p_m+c_*^2(\varepsilon-\varepsilon_m)$$

with $c_*^2=1$ gives the maximally stiff causal envelope in this model class. Baryon density must also be integrated from the first law if chemical potential or baryon mass is needed. Lattimer and Prakash show how causal stellar bounds depend on the trusted low-density segment and matching point in Lattimer and Prakash 2016, §2.2, Eqs. (6)–(12), and §4.2. Success of the extremal curve proves possibility under assumptions, not that real matter reaches the causal limit.

Evidence by domain, version, and likelihood

Dense-Matter Uncertainty and Inference Domains treats inference as a forward problem. Nuclear observables pass through a functional; radio masses pass through timing and binary models; NICER pulse profiles pass through instrument response, background, atmosphere, and hotspot geometry; gravitational-wave tides pass through detector calibration, waveform models, and spin and mass priors. Each has a versioned dataset and a density-response range rather than one exact density.

The inventory is current through 2026-08-18 and names the expanded 2018–2022 NICER J0740 dataset, the GWOSC GW170817-v1 release, refined radio timing, CREX, and theory ensembles. It also states selection effects and shared information. Multiplying the J0740 radio likelihood after it already served as the NICER mass prior double counts the same data; multiplying posterior samples without dividing out their multidimensional priors is likewise invalid.

Kumar and collaborators compare present laboratory, pulsar, NICER, and gravitational-wave constraints and their modeling boundaries in Kumar et al. 2024, §§7–8 and Tables 9–11. Robust joint inference should repeat the fit across equation-of-state function priors and remove each evidence family in turn. A narrow posterior inherited from a restrictive prior is not the same as a likelihood-driven constraint.

Choosing and testing an equation of state

Begin with a declared potential, particle inventory, units, and rest-mass convention. Check first-law and Maxwell identities before using stellar data. Keep saturation expansions inside their calibrated density range. Treat coexistence with an explicit phase construction. Validate hot tables in cold, classical, radiation, and degenerate limits. Enforce stability and causality along the actual dynamical constraint.

Then version every theory ensemble and dataset, reproduce its original likelihood inputs, and preserve correlations when combining messengers. Report the density-response map and leave-one-family-out changes. Finally, separate conditional bounds from composition claims: a mass–radius curve constrains stress, but inferring which particles carry that stress requires additional model-dependent evidence.

Review the chapter

A cold table matches pressure and energy values at every node, but its independently interpolated pressure gives $dp/d\varepsilon\gt1$ between two nodes. Is the table acceptable because the nodes are causal?

Solution

No. Evolution samples interpolated states, and the local acoustic derivative controls its characteristic speed. The superluminal interval may be an interpolation artifact, but it is still part of the implemented equation of state. Interpolate a generating potential or use a derivative-controlled scheme, refine the relevant axis, verify thermodynamic identities, and recheck the causal derivative before use.

References

• CompOSE Core Team. CompOSE Reference Manual, Version 3.01 (2022; the official download retains v3.00 in its filename). Official PDF.
• Drischler, Christian, Jeremy W. Holt, and Christian Wellenhofer. “Chiral Effective Field Theory and the High-Density Nuclear Equation of State.” Annual Review of Nuclear and Particle Science 71, 403–432 (2021). Open preprint.
• Kumar, Rajesh, et al., and the MUSES Collaboration. “Theoretical and Experimental Constraints for the Equation of State of Dense and Hot Matter.” Living Reviews in Relativity 27, 3 (2024). DOI.
• Lattimer, James M., and Madappa Prakash. “The Equation of State of Hot, Dense Matter and Neutron Stars.” Physics Reports 621, 127–164 (2016). Open preprint.
• Oertel, Micaela, Matthias Hempel, Thomas Klähn, and Stefan Typel. “Equations of State for Supernovae and Compact Stars.” Reviews of Modern Physics 89, 015007 (2017). DOI.