Causal Viscous Fluids and Non-Singular Cosmological Bounces

Technical Summary: Causal Viscous Fluids and the Realization of Non-Singular Bounces

Problem Statement
The resolution of the initial cosmological singularity remains a fundamental challenge in theoretical physics. Within standard General Relativity (GR) coupled to classical matter, the Hawking–Penrose singularity theorems dictate that homogeneous and isotropic Friedmann–Lemaître–Robertson–Walker (FLRW) universes inevitably encounter a past singularity. Non-singular bouncing cosmologies offer an alternative, requiring the violation of the Null Energy Condition (NEC), expressed as $\rho + p_{\text{eff}} < 0$ at the bounce minimum.

While dissipative processes in the early Universe could theoretically modify effective pressure, the standard first-order Eckart formulation of relativistic hydrodynamics fails to support bounces in spatially flat universes. In the Eckart theory, the bulk viscous pressure is algebraically tied to the Hubble parameter ($\Pi = -3\zeta H$). Since the Hubble parameter must vanish ($H=0$) at the bounce point, the viscous pressure also vanishes ($\Pi_b = 0$). Consequently, for ordinary matter satisfying $\rho_b + p_b \geq 0$, the NEC cannot be violated, rendering the Eckart theory incapable of generating non-singular bounces in flat FLRW cosmologies.

Methodology
The authors investigate the realization of non-singular bounces by replacing the Eckart formulation with the causal, second-order Israel–Stewart (MIS) theory of relativistic thermodynamics. The methodology proceeds as follows:

1. Theoretical Framework: The study utilizes the MIS transport equation, which promotes the viscous pressure $\Pi$ to an independent dynamical variable with a finite relaxation time $\tau$. This restores hyperbolicity to the fluid equations, ensuring finite signal propagation speeds and thermodynamic stability.
2. Truncation and Validity: The authors analyze the truncated MIS equation ($\tau \dot{\Pi} + \Pi = -3\zeta H$) versus the full equation. They derive a dimensionless parameter $\epsilon(t)$ to quantify the error of the truncation near the bounce. They establish that for transport families where the relaxation time scales as $\tau(\rho) \propto \rho^p$ with $p > 1/2$, the truncation is valid near the bounce ($t \to 0$), as the omitted higher-order terms vanish.
3. Analytical Construction: Using a parametrized scale factor ansatz $a(t) = a_b(1 + t^2/t_0^2)^n$, the authors derive exact near-bounce solutions. They determine the necessary initial conditions for $\Pi_b$ to satisfy the Raychaudhuri equation while maintaining $H_b = 0$ and $\dot{H}_b > 0$.
4. Extension to Modified Gravity: The framework is extended to two alternative gravity theories:
   • $f(R)$ Gravity: Incorporating higher-curvature corrections where the effective scalar degree of freedom contributes to the bounce dynamics.
   • Loop Quantum Cosmology (LQC): Embedding causal viscous fluids within the effective dynamics of LQC, where quantum geometry effects naturally induce a bounce at a critical density $\rho_c$.
5. Numerical Verification: The authors perform numerical integrations of the dimensionless dynamical system using adaptive Runge–Kutta algorithms to verify analytical predictions and explore parameter spaces for viscosity coefficients ($\zeta$) and relaxation times ($\tau$).

Key Contributions and Results

• Resolution of the Eckart No-Go Theorem: The paper demonstrates that the Israel–Stewart formulation permits a controlled violation of the NEC. Unlike the Eckart theory, $\Pi$ is not algebraically constrained to vanish when $H=0$. Instead, $\Pi$ evolves dynamically, allowing for a finite, negative viscous pressure at the bounce ($\Pi_b < 0$) that drives $\dot{H}_b > 0$.
• Exact Bounce Solutions in GR: In spatially flat GR, the authors construct exact bounce solutions where the energy density vanishes at the bounce ($\rho_b = 0$) and the viscous pressure takes the specific value $\Pi_b = -\alpha/(4\pi G)$, where $\alpha = 2n/t_0^2$ represents the bounce slope. The existence of these solutions is mathematically guaranteed by the Picard–Lindelöf theorem under the conditions $\zeta \geq 0$, $\tau > 0$, and $p > 1/2$.
• Thermodynamic and Causality Constraints: The study establishes a rigorous "Admissibility Region" for model parameters. To ensure positive entropy production, subluminal signal propagation, and the validity of the truncated MIS approximation, the parameters must satisfy:
  • $\zeta \geq 0$ and $\tau > 0$.
  • The relaxation time exponent $p > 1/2$ (ensuring the truncation error $\epsilon \to 0$ as $t \to 0$).
  • The viscosity exponent $q < 1$ (ensuring regularity).
  • The effective sound speed $0 \leq c_s^2 \leq 1$.
• Role in Extended Gravity:
  • In $f(R)$ gravity, higher-curvature terms (specifically $\ddot{F}_b$) complement the viscous pressure. This allows for bounces with non-zero matter density ($\rho_b > 0$), expanding the parameter space for viable solutions compared to GR.
  • In LQC, the bounce is guaranteed by quantum geometry effects at $\rho = \rho_c$. Here, causal viscosity does not trigger the bounce but modifies the post-bounce dynamics and the slope $\dot{H}_b$. The authors show that for typical parameters, the viscous relaxation prevents $\Pi_b$ from becoming so negative that it would violate the LQC bounce condition ($\rho_c + p_b + \Pi_b > 0$).

Significance and Claims
The paper claims to establish a unified description in which bulk viscosity, modeled via the causal Israel–Stewart theory, provides a physically consistent mechanism for singularity resolution across different theories of gravity.

• Physical Consistency: The authors argue that causal viscosity is not an ad hoc modification but a natural consequence of dissipative processes in the hot, dense early Universe. It resolves the singularity problem without introducing exotic matter fields or ghost condensates that often suffer from instabilities.
• Robustness: The results highlight that viscosity-induced bounces are robust features that persist when embedded in both higher-curvature corrections ($f(R)$) and quantum geometry effects (LQC).
• Phenomenological Implications: The study suggests that the relaxation dynamics of the viscous pressure control the duration of the NEC-violating phase. This mechanism may leave distinctive imprints on the primordial power spectrum and relic abundances, offering a potential bridge between microscopic transport physics and cosmological observations.

The authors conclude that once causal transport processes are accounted for, non-singular bounces may be a generic outcome of early-universe physics, offering a new paradigm where transport processes are indispensable to resolving the initial singularity.