Cosmological Stealth fields and Non-Equilibrium thermodynamics

Abstract: We investigate the connection between cosmological stealth scalar fields and non-equilibrium thermodynamics in a spatially flat Friedmann-Lemaître-Robertson-Walker (FLRW) background. We consider a non-minimally coupled scalar field whose energy-momentum tensor vanishes identically, allowing the field to evolve on a dissipative cosmological background without producing gravitational backreaction. We show that the stealth condition leads to a generalized Riccati equation for the scalar-field kinematics, where the dissipative pressure acts as a thermodynamic driving term. In terms of the variable $y=\dotφ/(Hφ)$, the system admits two thermodynamic branches: a stable attractor selected by entropy production and an unstable repeller. We also construct the corresponding phase-space structure in the $(y,ω_{\rm eff})$ plane and identify a near-critical regime associated with $ζ=1/4$. Finally, we reconstruct the stealth potential and present a bulk viscous realization in which irreversible entropy production drives the universe toward an asymptotic de Sitter state while the stealth field tracks the dissipative background. Our results suggest that stealth fields can be interpreted as dynamically non-trivial thermodynamic trackers of non-equilibrium cosmological evolution.

• The paper derives a Riccati equation showing that a non-minimally coupled stealth scalar tracks cosmological entropy production through the dissipative pressure while remaining gravitationally invisible.
• The analysis identifies a thermodynamically selected stable attractor, constrains strong dissipation to coupling values ζ ≥ 0.238, and permits effective phantom expansion without ghost-like kinetic terms.
• The paper reconstructs closed-form scalar profiles for bulk-viscous cosmologies, including evolution toward de Sitter space, while highlighting the need for causal thermodynamics and perturbation-level tests.

Overview

This paper establishes a connection between stealth scalar field configurations—non-trivial scalar profiles whose energy-momentum tensor vanishes identically—and the irreversible thermodynamics of dissipative cosmological fluids. The analysis is performed on a spatially flat FLRW background in geometrized units, with a scalar field non-minimally coupled to gravity through $-\frac{1}{2}\zeta R\phi^2$. The central result is that enforcing the stealth condition on a background containing a dissipative pressure $\Pi$ reduces the field's kinematics to a generalized Riccati equation in which $\Pi$ acts as a thermodynamic driving term. The field thereby evolves dynamically while remaining gravitationally invisible, functioning as a tracker of entropy production rather than as an additional dark-energy component.

Thermodynamic setup

The cosmic fluid is described by an effective pressure $P = p + \Pi$, where $\Pi$ encodes departures from local thermodynamic equilibrium. The Friedmann equations are modified accordingly, and the Gibbs relation yields an entropy production rate proportional to the dissipative pressure:

$\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$

Positive entropy production during expansion therefore requires $\Pi < 0$. Bulk viscosity is identified as the only dissipative mechanism compatible with FLRW symmetries, although particle creation is noted as an alternative entropy-generating channel. All entropy production in the model is encoded solely in $\Pi$; no additional fluid with independent thermodynamic degrees of freedom is introduced.

The stealth Riccati equation

Imposing $T^{(S)}_{\mu\nu} = 0$ for the non-minimally coupled field and eliminating the potential $V(\phi)$ between the time and spatial components of the stealth condition produces

$\Pi$0

where $\Pi$1. This is the paper's foundational equation: a nonlinear Riccati equation driven by the thermodynamic state of the universe. Introducing the dimensionless variable $\Pi$2 and a barotropic effective equation of state $\Pi$3, the system becomes autonomous in the variable $\Pi$4:

$\Pi$5

with effective coupling $\Pi$6. Rewriting the evolution in terms of entropy as the time variable, $\Pi$7 is directly proportional to the same quadratic polynomial divided by $\Pi$8, making explicit that the field cannot respond instantaneously to changes in the dissipation rate—a "dynamical inertia" regulated by $\Pi$9. In the adiabatic limit $\Pi$0, entropy ceases to be a suitable clock for the stealth dynamics, although evolution in $\Pi$1 remains well defined.

The potential is not specified a priori; it is reconstructed from the stealth condition itself,

$\Pi$2

and thus acts as an adaptive thermodynamic sink guaranteeing the field's invisibility.

Attractors, thermodynamic bounds, and phantom behavior

The fixed points $\Pi$3 are the roots of the quadratic polynomial, and their existence requires a non-negative discriminant,

$\Pi$4

Because the eigenvalues of the linearized dynamics are $\Pi$5 and $\Pi$6, the branch $\Pi$7 is always a stable attractor and $\Pi$8 an unstable repeller. The irreversible arrow of time therefore selects the physically admissible branch—a model-independent result, since the sign of $\Pi$9 is fixed by the second law.

For dissipative pressureless matter, the discriminant bound translates into a constraint between $P = p + \Pi$0 and $P = p + \Pi$1. Two concrete consequences are worth emphasizing. First, a de Sitter-like expansion driven purely by dissipation of matter ($P = p + \Pi$2) demands $P = p + \Pi$3; weakly coupled fields cannot track such strongly dissipative backgrounds while remaining stealth, and the conformal value $P = p + \Pi$4 is excluded. Second, an effective phantom regime (e.g., $P = p + \Pi$5) is achievable through macroscopic dissipation alone, requiring only $P = p + \Pi$6, without ghost-like negative kinetic terms. This offers a stable mechanism for phantom crossing in which the stealth field's Lagrangian remains regular.

Phase-space structure and the critical coupling $P = p + \Pi$7

Promoting $P = p + \Pi$8 to an independent dynamical variable with generic relaxation law $P = p + \Pi$9 yields a model-independent two-dimensional system in the $\Pi$0 plane. The Jacobian analysis confirms that $\Pi$1 is always stable along the stealth direction and becomes a global attractor when $\Pi$2; the branch selection is thus independent of the microscopic transport physics.

A structural transition occurs at the critical coupling $\Pi$3, where $\Pi$4 and the Riccati structure degenerates into a linear equation. In a near-critical example with $\Pi$5 and a linear relaxation toward $\Pi$6, the unstable branch diverges to $\Pi$7, leaving a single thermodynamically selected attractor. The phase portrait naturally separates decelerating ($\Pi$8), accelerating, and effective phantom ($\Pi$9) regimes, the last induced entirely by dissipation. Superimposing the reconstructed potential landscape shows that viable trajectories predominantly occupy the $\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$0 sector near criticality, with the potential's stationary point located at $\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$1 and the profile flattening as $\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$2.

Bulk viscous realizations

Two explicit implementations validate the framework. For a constant effective equation of state $\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$3, the Hubble parameter integrates to $\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$4 with $\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$5, and the Riccati equation reduces to an Euler–Cauchy problem after the standard substitution $\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$6. The scalar profile admits a closed-form reconstruction,

$\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$7

valid for $\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$8. In the Eckart model with $\frac{T}{V}\frac{dS}{dt} = -3H\Pi.$9 and barotropic $\Pi < 0$0, the background obeys a logistic equation approaching the de Sitter attractor $\Pi < 0$1, so that $\Pi < 0$2 at late times. The stealth fixed points reduce to $\Pi < 0$3, and the field evolves exponentially, $\Pi < 0$4, tracking the scale factor and maintaining the stealth condition throughout. Accelerated and phantom expansion conditions reduce to simple inequalities on $\Pi < 0$5.

Limitations and open questions

The paper concedes several restrictions on its results. The dynamical analysis relies on first-order (Eckart-type) thermodynamics, which suffers from well-known acausality and instability issues; the authors state that extension to the causal Israel–Stewart formalism is necessary. The bulk viscous realizations assume a constant viscosity coefficient; whether a density- or Hubble-dependent $\Pi < 0$6 produces richer phase-space topologies is left open. The entropy-based formulation also breaks down in the adiabatic limit, where entropy production no longer parametrizes the stealth dynamics. Finally, the analysis is purely at the background level: no perturbative observables, growth-of-structure predictions, or observational constraints are derived, so the phenomenological viability of the framework against precision cosmological data remains untested. The critical case $\Pi < 0$7 requires separate treatment and is not solved in closed form for the general dissipative background.

Conclusion

The paper demonstrates that the stealth condition for a non-minimally coupled scalar field on a dissipative FLRW background is equivalent to a Riccati equation driven by the dissipative pressure, linking the field's kinematics directly to entropy production. The irreversible arrow of time selects a unique stable attractor branch, a model-independent result, and the coupling bound $\Pi < 0$8 for strongly dissipative matter (excluding the conformal value $\Pi < 0$9) constitutes a sharp, quantitative constraint. The framework permits effective phantom expansion without ghost degrees of freedom, and explicit bulk viscous realizations confirm that the stealth profile can be reconstructed in closed form while tracking the background toward de Sitter. The main open questions concern causal thermodynamic extensions and observational confrontation of the framework.