Causal Bulk Viscous Cosmologies

• Causal bulk viscous cosmologies are models that incorporate dynamic bulk viscous pressure through extended irreversible thermodynamics to preserve causality and prevent singularities.
• They utilize hyperbolic transport equations, such as the Müller–Israel–Stewart framework, to drive accelerated expansion and modulate structure formation.
• These models offer a versatile framework for unifying dark sector dynamics and provide observational signatures that distinguish them from standard ΛCDM cosmology.

A causal bulk viscous cosmology generalizes the standard Friedmann framework by incorporating dissipative effects arising from the bulk viscosity of cosmic fluids, modeled within extended irreversible thermodynamics. Unlike the acausal and unstable Eckart theory, modern causal bulk viscous models employ the full Müller–Israel–Stewart (MIS) or related hyperbolic transport formulations, where the bulk viscous stress is treated as a dynamical variable with finite relaxation time. This approach preserves causality, ensures thermodynamic consistency, and enables a diverse phenomenology, including singularity avoidance, dynamical dark energy, modified structure growth, and non-singular bounces across different gravitational settings. The following sections provide a thorough technical overview of the modeling frameworks, analytic and numerical solution structures, physical consequences, and observational constraints of causal bulk viscous cosmologies.

1. Fundamental Principles and Theoretical Frameworks

The starting point for causal bulk viscous cosmologies is the spatially flat Friedmann–Robertson–Walker (FRW) or Friedmann–Lemaître–Robertson–Walker (FLRW) metric,

$ds^2 = -dt^2 + a^2(t) d\vec{x}^2,$

with scale factor $a(t)$ and Hubble parameter $H = \dot{a}/a$.

An imperfect (viscous) fluid is modeled via the energy–momentum tensor,

$$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$

where $\rho$ is the energy density, $p$ the equilibrium pressure, and $\Pi \leq 0$ the bulk viscous pressure (dissipative correction).

The key innovation in causal theory is the relaxation-type evolution law for $\Pi$. The full Israel–Stewart (IS) equation in homogeneous cosmologies takes the form,

$$\tau \, \dot{\Pi} + \Pi = -3 \zeta H - \frac{1}{2} \tau \Pi [3H + \frac{\dot{\tau}}{\tau} - \frac{\dot{\zeta}}{\zeta} - \frac{\dot{T}}{T}],$$

with $\zeta(\rho)$ the bulk viscosity coefficient, $a(t)$0 the relaxation time, and $a(t)$1 the temperature. The truncated MIS (Maxwell–Cattaneo) form omits the bracketed second-order term.

Physical closure necessitates equations of state of the form $a(t)$2, along with power-law ansätze for transport coefficients: $a(t)$3, $a(t)$4, $a(t)$5 (Tawfik et al., 2010).

Causality and stability require that the total characteristic speed—combining sound and viscous signal propagation—remain subluminal. This is quantified by

$a(t)$6

constrained by $a(t)$7 (Gavassino et al., 2023).

2. Dynamical Equations and Analytic Solution Structures

The Einstein equations in the causal viscous fluid context are

$a(t)$8

Substituting the viscous transport law and closure relations yields a second-order nonlinear ODE for $a(t)$9 (see (Tawfik et al., 2010, Cornejo-Pérez et al., 2012, D et al., 2017, Chimento et al., 2012)):

$H = \dot{a}/a$0

where the coefficients $H = \dot{a}/a$1 are determined by the equation-of-state and viscous parameters.

Special cases admit closed-form analytic or parametric solutions:

• QGP-motivated early universe: Under ultra-relativistic equation of state ($H = \dot{a}/a$2) and transport powers $H = \dot{a}/a$3, the solution for the Hubble parameter avoids the Big Bang singularity, yielding finite $H = \dot{a}/a$4, $H = \dot{a}/a$5, $H = \dot{a}/a$6, $H = \dot{a}/a$7, and $H = \dot{a}/a$8 (Tawfik et al., 2010).
• Factorization and Lie symmetries: For specific $H = \dot{a}/a$9 pairs, exact parametric and scaling solutions exist, including special inflationary or non-singular branches. At $$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$0, the system admits both power-law (scaling) and exponential expansion with explicit analytic forms for all state variables (Cornejo-Pérez et al., 2012, Chimento et al., 2012).
• Unified dark matter: Nonlinear extensions of full MIS theory allow the coexistence of de Sitter, quintessence, and phantom attractors depending on the exponent $$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$1, with late-time acceleration arising for $$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$2 and phantom behavior for $$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$3 (Palma et al., 13 Oct 2025).

The general solution space is richer than non-causal or algebraic models and supports aperiodic, singularity-free, or attractor-driven cosmological evolutions.

3. Physical Implications: Singularity Avoidance, Acceleration, and Cosmological Bounces

Causal bulk viscous cosmologies yield several phenomena absent in standard perfect-fluid treatments:

• Singularity avoidance: The presence of relaxation terms and dynamically generated negative pressure can prevent the initial curvature singularity, with all thermodynamic quantities remaining finite at $$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$4 in QGP-era models (Tawfik et al., 2010, Eshaghi et al., 2015). Non-singular bounces are accessible both in general relativity and in extensions such as $$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$5 gravity and loop quantum cosmology, subject to explicit conditions on the dissipative parameters ensuring positive entropy production and causal stability (Yildiz et al., 26 Aug 2025).
• Transient (self-inflation): The bulk viscosity can drive a brief era of accelerated expansion, with deceleration parameter $$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$6 and effective equation-of-state $$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$7, eventually relaxing to radiation- or matter-dominated deceleration. Such transitions occur naturally in models with QGP or stiff-fluid initial conditions (Tawfik et al., 2010, D et al., 2017).
• Late-time acceleration and cosmic speed-up: In unified dark sector scenarios, the MIS theory generates a robust mechanism for late-time cosmic acceleration without requiring a cosmological constant. For $$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$8 and suitably chosen $$T_{\mu\nu} = [\rho + p + \Pi] u_\mu u_\nu + [p + \Pi] g_{\mu\nu},$$9, the late-time attractor is either de Sitter (quintessence-like) or phantom, with basin-structure separating the two classes of solutions (Palma et al., 13 Oct 2025). The de Sitter attractor appears generically for $\rho$0, and cosmic acceleration arises even for small bulk viscosity, circumventing fine-tuning issues (Palma et al., 13 Oct 2025).

4. Structure Formation, Linear Perturbations, and Observational Signatures

Perturbation theory within causal bulk viscous cosmologies exhibits features not present in non-causal or perfect-fluid analogs. The Israel–Stewart (IS) relaxation equation leads to a third-order ODE for the density contrast $\rho$1 in the Newtonian limit:

$\rho$2

where the structure of $\rho$3 depends on the background cosmology, viscous exponent $\rho$4, sound speed $\rho$5, and other transport coefficients (Acquaviva et al., 2016, Acquaviva et al., 2018).

Distinctive consequences include:

• Suppression or enhancement of clustering: For standard parameter choices, bulk viscosity damps the growth of large-scale structure, leading to reduced matter power at small scales and potentially reconciling observational tension in $\rho$6 between CMB and low-$\rho$7 surveys. However, for $\rho$8 and large $\rho$9, causal corrections can enhance clustering relative to $p$0CDM, an effect tied to the sign and magnitude of causal terms in perturbation coefficients (Acquaviva et al., 2016, Acquaviva et al., 2018). This challenges the naive expectation that viscosity always dampens structure.
• Viable observational windows: Causal theory allows close tracking of $p$1CDM background and perturbations provided the viscous sound speed $p$2 lies within $p$3 (truncated MIS), producing compatible gravitational potentials and matter growth up to current measurement precision (Piattella et al., 2011).
• Redshift-space distortions and ISW: The redshift evolution of $p$4 and integrated Sachs–Wolfe effect are sensitive to the detailed viscous dynamics, offering avenues for empirical constraints through large-scale structure and CMB cross-correlation data (Acquaviva et al., 2018).

Physically realistic regimes must satisfy both causality and positive entropy production, bounding the parameter space and ensuring phenomenological viability.

5. Generalizations: Particle Creation, Anisotropic Cosmologies, and Nonlinear Rheology

Causal bulk viscosity models accommodate significant generalizations:

• Particle creation: Open-system cosmologies treating particle creation as an irreversible process introduce an additional "creation pressure" to the stress tensor. Explicit solutions demonstrate that bulk viscosity plus even modest particle creation rates soften the initial singularity and permit smooth transitions from inflation to radiation and late-time acceleration without requiring a cosmological constant (Eshaghi et al., 2015, Zeyauddin et al., 2013).
• Anisotropic frameworks: In Bianchi type V universes, combined effects of anisotropy, causal viscosity, and particle creation allow for both singular and non-singular exact cosmological models, with direct calculation of all physical observables and explicit algebraic distinction from Eckart and truncated theories (Zeyauddin et al., 2013).
• Nonlinear rheology: Recent advances employ full nonlinear causal transport equations capable of capturing pseudoplastic (shear-thinning) and dilatant (shear-thickening) bulk response, enabling the modeling of far-from-equilibrium cosmic epochs such as phase transitions, bounces, or viscoelastic dark sector physics (Gavassino et al., 2023).

These generalizations further underscore the flexibility and dynamical richness of the causal viscous approach.

6. Thermodynamics, Energy Conditions, and Bayesian Model Selection

The causal MIS framework guarantees that the entropy production rate,

$p$5

is strictly positive for all physical solutions, provided $p$6 and relaxation times are positive (Gavassino et al., 2023, D et al., 2017, Disconzi et al., 2014). This constraint excludes acausal models (e.g., Eckart), which can violate the local second law, and explains the preference for causal extensions in both early- and late-universe phenomenology.

Energy condition analysis reveals that strong transient or permanent violations of the strong or null energy conditions are possible, directly enabling non-singular bounces or late acceleration (Chimento et al., 2012, Yildiz et al., 26 Aug 2025). The stability of attractors and the approach to asymptotic states are determined by Lyapunov analyses in the relevant parameter regimes.

Comparative Bayesian inference using supernovae (e.g., Pantheon sample) indicates that while the causal viscous matter model can match $p$7CDM in $p$8 evolution and fit the data with similar likelihood, it is statistically disfavored, with a Bayes factor $p$9 (i.e., odds against the viscous model) (Shareef et al., 2021), and predicts a lower universe age ($\Pi \leq 0$0 Gyr) than $\Pi \leq 0$1CDM ($\Pi \leq 0$2 Gyr). Nevertheless, models remain viable for specific applications or as unified dark sector scenarios.

---

The causal bulk viscous cosmology framework thus represents a mathematically consistent and physically flexible extension of the standard cosmological model. It enables singularity-softening, dynamical acceleration, and scale-dependent structure growth, with predictive power across both background and perturbation evolution, provided the model parameters are chosen to ensure causality, stability, and agreement with empirical constraints (Tawfik et al., 2010, Palma et al., 13 Oct 2025, Eshaghi et al., 2015, D et al., 2017, Gavassino et al., 2023, Yildiz et al., 26 Aug 2025, Acquaviva et al., 2018, Piattella et al., 2011).