Speed-of-Sound Parametrization

• Speed-of-sound parametrization is a formal framework that expresses the variation of sound speed with macroscopic variables, linking microphysics and observable phenomena.
• It encompasses methods like constant speed-of-sound (CSS) and density-dependent models to efficiently simulate phase transitions and neutron star equation-of-state behavior.
• These parametrizations ensure analytic tractability and thermodynamic consistency, making them essential for computational studies in astrophysics, plasma physics, and condensed matter.

A speed-of-sound parametrization is any formal scheme, analytic or algorithmic, for expressing the variation of the sound speed $c_s$ (or, more commonly, its square $c_s^2\equiv\partial p/\partial\varepsilon$ at fixed entropy) as a function of macroscopic variables—such as density, energy, temperature, or composition—within a given physical system. Such parametrizations play a central role across astrophysics, nuclear matter, plasma physics, and condensed matter, as they form the core link between microphysics and observable bulk dynamics, compressibility, and wave propagation properties.

1. General Structure of Speed-of-Sound Parametrizations

Speed-of-sound parametrizations appear in two primary forms: (i) direct expressions $c_s^2(x)$, with $x$ a physical variable (density, energy density, pressure, temperature, etc.), and (ii) indirect relations where $c_s^2$ is derived from a parametrized equation of state (EoS) via thermodynamic or hydrodynamic derivatives. The aim is to capture essential physics (e.g., phase transitions, stiffening/softening at high density, crossovers) in a compact, predictive, and computationally tractable manner, often subject to physical requirements such as causality ($c_s^2\leq 1$ in units where $c=1$), monotonicity, and stability (positivity).

Several research domains have developed customized parametrizations; prominent among these are:

• The constant speed of sound (CSS) frameworks for high-density matter and hybrid stars (Pal et al., 3 Apr 2025, Christian et al., 2017).
• Density-dependent analytic schemes for baryonic matter (Pal et al., 3 Jul 2025).
• Piecewise and spectral fits for neutron-star EoSs ensuring $c_s^2$ continuity (Servignat et al., 2023).
• Scaling laws for simple liquids in terms of temperature, density, and freezing point (Khrapak, 22 Jul 2025).
• Polynomial/virial expansions for fluids and mixtures (Lozano-Martín et al., 2024).
• Hydrodynamic modeling in high-energy collisions via direct extraction from data (Gao et al., 2015).

In all cases, the parameterization acts as both a technical tool for efficient computation and as a vehicle for mapping micro-physical uncertainty (e.g., quark matter EoS, plasma correlations) onto testable macroscopic observables.

2. The Constant Speed-of-Sound (CSS) Framework in Compact Stars

The CSS parametrization is widely used for modeling quark matter or generic high-density phases in neutron stars, especially when interfaced with nucleonic EoSs via a first-order (Maxwell) transition.

In its canonical form (Pal et al., 3 Apr 2025, Christian et al., 2017): $$p(\varepsilon) = \begin{cases} p_{\rm had}(\varepsilon), & \varepsilon < \varepsilon_{\rm tr} \ p_{\rm had}(\varepsilon_{\rm tr}) + c_s^2 (\varepsilon - \varepsilon_{\rm tr}) - \Delta\varepsilon, & \varepsilon \geq \varepsilon_{\rm tr} \end{cases}$$ where

• $c_s^2$ is the constant squared speed of sound in the high-density (e.g., quark-matter) phase.
• $\varepsilon_{\rm tr}$ is the transition energy density (associated with transition baryon density $\rho_{\rm tr}$).
• $\Delta\varepsilon$ is the latent heat or energy density jump at the interface.

Physical and observational properties of compact objects, such as existence of stable hybrid branches ("twin stars"), radius–mass relations, and maximum mass, are tightly constrained by the choice of $(\rho_{\rm tr},\ \Delta\varepsilon,\ c_s^2)$. The Seidov limit specifies the minimum energy jump that produces instability at the interface: $$\Delta\varepsilon_{\rm cr} = \frac{1}{2}\varepsilon_{\rm tr} + \frac{3}{2}p_{\rm tr}$$

The parameter space can be rapidly explored to match observed neutron star properties and gravitational-wave constraints, as demonstrated in studies of low-mass hybrid candidates (e.g., HESS J1731-347) (Pal et al., 3 Apr 2025). The CSS approach thus forms the backbone of many modern neutron-star inference pipelines.

3. Density-Dependent and Flexible Parametrizations for Neutron-Star Matter

To interpolate between soft and stiff hadronic EoSs and enable efficient scans over broad physical assumptions, continuous density-dependent forms are favored. An example is the three-parameter model (Pal et al., 3 Jul 2025): $$c_s^2(\varepsilon) = a\left[1 - \exp\left(-b(\varepsilon/\varepsilon_0)^c\right)\right]$$ with the parameters $a$, $b$, $c$ controlling stiffness, saturation, and curvature, respectively, and $\varepsilon_0$ a reference (saturation) energy density. Thermodynamic integrals of $c_s^2$ yield $p(\varepsilon)$, which can then be used to solve the Tolman–Oppenheimer–Volkoff equations.

This approach allows the delineation of permissible ($a$, $b$, $c$) domains consistent with neutron-star mass–radius landscape, tidal deformability bounds from gravitational-wave events, and requirements of causality.

4. Analytic and Spectral EoS Fits with Continuous Sound Speed

For neutron-star evolution, oscillation, and binary-merger simulations, the analytic continuity of $c_s^2$ is numerically and physically crucial. Two families of parametrizations predominate (Servignat et al., 2023):

• Pseudo-polytrope fits: Use polynomials of $\ln n_B$ modulated by a leading power-law index to fit specific internal energy per baryon, guaranteeing analytic expressions for $p(n_B)$ and $c_s^2(n_B)$.
• Potekhin–Pearson spectral representations: Fit $\log_{10}p$ vs.\ $\log_{10}e$ using sigmoid transitions (soft-steps), allowing high-fidelity spectral accuracy.

Both types are coupled to multi-segment “crust-stitching” algorithms that enforce continuity of $p$, $e$, and $c_s^2$ at EoS domain boundaries. Validation against known EoSs demonstrates percent-level reproduction of macroscopic observables, and these schemes are readily generalized to multi-parameter (composition-dependent) forms.

5. Performance and Physical Interpretation in Diverse Physical Systems

Speed-of-sound parametrizations are not limited to neutron stars or quark matter. They appear broadly:

• Dense fluids, plasmas, and liquids: Empirical scaling laws provide $c_s$ in terms of temperature, thermal velocity, heat-capacity ratio, and freezing line—e.g.,

$$\frac{c_s}{v_T} \simeq \sqrt{\gamma} + \alpha \left(\frac{T_{\rm fr}}{T}\right)^\beta$$

for Lennard-Jones fluids and simple atomic liquids, with parameters supplied by fit or simulation (Khrapak, 22 Jul 2025).

• Yukawa one-component plasmas: The sound speed $s(\kappa, \Gamma)$ is mapped across five physical regimes through closed-form expressions with domain-dependent accuracy (~5%) (Silvestri et al., 2019).
• Mixtures and chemical fluids: Acoustic virial expansions and polynomial fits link $c_s$ to composition, pressure, and temperature (e.g., for hydrogen–methane mixtures) with precisely fitted virial coefficients and adiabatic exponents (Lozano-Martín et al., 2024).
• High-energy collisions: Phenomenological extractions from (pseudo)rapidity spectra provide an effective $c_s^2$ in multiparticle fireballs, linked directly to the width of the central Gaussian rapidity distribution and revealing energy-insensitive “perfect liquid” sound speeds at midrapidity (Gao et al., 2015).

6. Mathematical and Algorithmic Aspects

Most parametrizations are constructed with a view to analytic tractability of derivatives, invertibility (for microphysical inference), and ease of integration in simulations:

• Piecewise-linear, spectral, and polynomial forms allow for fast evaluation and guarantee smooth thermodynamics.
• In multi-phase or phase-transition scenarios, parametrizations are explicitly matched at the phase boundary by enforcing continuity of pressure and chemical potential (Maxwell construction), with parameters such as latent heat/jump ($\Delta\varepsilon$) playing a critical role.
• Fitting and inversion strategies (e.g., for Arctic ocean acoustics) exploit parametrization schemes with minimal degrees of freedom for robust and computationally efficient inference (Weng et al., 10 Aug 2025).

7. Physical and Observational Constraints

Speed-of-sound parametrizations are tightly constrained by both theoretical bounds and observational data:

• Causality: $c_s^2\leq 1$ (natural units), though some models permit saturation at the causal limit for stiffest possible EoS (Christian et al., 2017).
• Stability: Avoidance of mechanical and convective instabilities (e.g., positivity and monotonicity of $c_s^2$).
• Mass–Radius and Tidal Deformability: Consistency with observed neutron-star properties imposes sharp upper and lower bounds on parameter choices.
• Wave Propagation and Spectral Data: For fluids/plasmas, direct comparison with molecular dynamics or laboratory data ensures global accuracy and fidelity.

In heavy-ion physics and cosmology, effective $c_s^2$ inferred from data serves not only as a diagnostic of underlying microphysics (phase transitions, degrees of freedom) but also as a phenomenological tool for hydrodynamic evolution.

---

Selected Table: CSS Parameter Effects on Hybrid Stars (Pal et al., 3 Apr 2025)

Parameter	Physical Meaning	Effect on $M$–$R$ Relation
$\rho_{\rm tr}$	Transition baryon density	Low $\rho_{\rm tr}$: early deconfinement, larger quark core in low $M$ stars; high $\rho_{\rm tr}$: delayed transition, smaller/no core
$\Delta\varepsilon$	Energy-density jump (latent heat)	Higher values soften post-transition branch, reduce $R$ for given $M$
$c_s^2$	Quark-phase stiffness	$>1/3$ allows higher $M_{\rm max}$, stiffer EoS; lower $c_s^2$ yields more compact stars

---

Parametrizations of the speed of sound are thus indispensable tools for bridging microscopic physics, equations of state, and macroscopic observables across astrophysics, plasma physics, condensed matter, and experimental high-energy physics. Their development continues to be driven by new observational constraints and the need for computationally robust, physically transparent models.