Cosmic Hysteresis in Reconstructed $f(T)$ Bounce Models A Torsion-Based Thermodynamic Perspective

1. Problem Statement

The standard Big Bang cosmological model, while observationally successful, suffers from a fundamental singularity at $t=0$ where General Relativity (GR) breaks down. To resolve this, bouncing and cyclic cosmologies have been proposed, where the universe transitions smoothly from contraction to expansion, avoiding the singularity.

However, a critical thermodynamic question remains: Is the evolution of a cyclic universe reversible? In standard GR with perfect matter, a cyclic universe would theoretically retrace its path exactly, performing zero net thermodynamic work over a cycle. Recent studies in curvature-based modified gravity (e.g., $f(R)$ and Einstein-Gauss-Bonnet) have shown that cosmic hysteresis—a non-vanishing thermodynamic work integral over a cycle—can arise due to asymmetric scalar field dynamics.

The Gap: It was unknown whether this phenomenon of cosmic hysteresis persists in teleparallel gravity, specifically $f(T)$ gravity, where gravitation is attributed to torsion ($T$) rather than curvature. The authors investigate whether torsion-induced geometric modifications generate irreversible thermodynamic behavior and a cosmological arrow of time in cyclic models.

2. Methodology

A. Theoretical Framework: $f(T)$ Gravity

The authors utilize $f(T)$ gravity, an extension of the Teleparallel Equivalent of General Relativity (TEGR).

• Geometry: Instead of the Levi-Civita connection (curvature), they use the Weitzenbäck connection (torsion-free curvature, non-zero torsion).
• Action: The gravitational action is modified from the torsion scalar $T$ to an arbitrary function $f(T)$, minimally coupled to a canonical scalar field $\phi$:
  $$S = \int d^4x \, e \left[ \frac{1}{2\kappa^2}f(T) - \frac{1}{2}g^{\mu\nu}\partial_\mu\phi\partial_\nu\phi - V(\phi) \right]$$
• Field Equations: The modified Friedmann equations are derived for a spatially flat FLRW universe. Crucially, unlike $f(R)$ gravity, $f(T)$ gravity retains second-order equations of motion.

B. Reconstruction Procedure

Rather than guessing an $f(T)$ function and solving for the scale factor $a(t)$, the authors employ a reconstruction method:

1. Prescribe $a(t)$: They define two specific analytic bouncing scale factors:
   • Exponential Bounce: $a(t) = a_0 \cosh(\alpha t) + a_b$ (Time-symmetric).
   • Power-Law Bounce: $a(t) = a_b \left[ 1 + (t/t_b)^{2n} \right]^{1/2n}$ (Flexible transition sharpness).
2. Derive $T(t)$: Using the relation $T = -6H^2$ (where $H = \dot{a}/a$), they express the torsion scalar as a function of time.
3. Solve for $f(T)$: They algebraically invert the modified Friedmann equations to find the specific $f(T)$ function that generates the prescribed bouncing evolution.

C. Numerical Simulation

To quantify hysteresis, the authors performed numerical simulations using a fourth-order Runge-Kutta scheme:

• System: Coupled evolution of the reconstructed $f(T)$ geometry and the scalar field $\phi$.
• Detection: Used local extrema detection (identifying $\dot{a}=0$ and $\ddot{a} \gtrless 0$) to precisely locate bounces (minima) and turnarounds (maxima), avoiding the ambiguity of threshold-based methods.
• Integration: Calculated the thermodynamic work integral $W = \oint p_\phi dV$ over multiple complete cycles.

3. Key Contributions

1. First Demonstration of Hysteresis in $f(T)$ Gravity: The paper establishes that cosmic hysteresis is not limited to curvature-based theories but is a generic feature of geometric modifications, including torsion-based formulations.
2. Mechanism of Irreversibility: The authors identify the microscopic origin of hysteresis in $f(T)$ gravity. It arises from the asymmetric coupling between the torsion scalar and the scalar field:
   • During contraction ($H < 0$), the friction term $3H\dot{\phi}$ in the Klein-Gordon equation acts as anti-friction, amplifying the scalar field's kinetic energy.
   • During expansion ($H > 0$), it acts as standard friction, damping the field.
   • Torsion Corrections: The modified Friedmann equations introduce additional geometric terms ($f_T, f_{TT}$) that amplify this asymmetry, causing the pressure $p_\phi$ to differ at the same scale factor $a$ depending on the direction of evolution.
3. Thermodynamic Memory: The study demonstrates that the universe possesses a "memory" of its expansion/contraction history encoded in the geometric state (torsion), leading to a non-zero work integral even in the absence of explicit dissipative matter terms.

4. Key Results

A. Non-Vanishing Thermodynamic Work

The numerical integration of the work integral $W = \oint p_\phi dV$ over six complete cycles yielded:

• Total Work: $W_{tot} \simeq -1.7 \times 10^{13}$ (in natural units).
• Average Work per Cycle: $\langle W \rangle \simeq -2.8 \times 10^{12}$.
• Significance: The non-zero value confirms irreversibility. The negative sign indicates a net transfer of energy from the scalar field to the geometric sector (torsion) over each cycle.

B. Hysteresis Loops

In the $(w_\phi, a)$ plane (Equation of State vs. Scale Factor), the evolution traces closed loops rather than a single curve.

• At a given scale factor $a$, the pressure $p_\phi$ (and thus $w_\phi$) is higher during contraction than during expansion.
• The area enclosed by these loops corresponds to the dissipated work.

C. Secular Evolution of Cycles

The cycles are not identical; they exhibit systematic drifts indicating a far-from-equilibrium state:

• Bounce Minima: Decrease by $\Delta a_{min} \simeq -5.96$ per cycle (universe contracts tighter).
• Turnaround Maxima: Increase by $\Delta a_{max} \simeq 4.97$ per cycle (expansion reaches larger volumes).
• Period Variability: Cycle periods fluctuate significantly (Coefficient of Variation $\approx 36.6\%$), ranging from 3.7 to 8.5, unlike the constant periods of equilibrium oscillators.

5. Significance and Implications

• Cosmological Arrow of Time: The study provides a novel mechanism for the arrow of time in cyclic universes. Unlike standard explanations relying on low-entropy initial conditions (Past Hypothesis), this arrow emerges geometrically from the asymmetric response of torsion-modified gravity to expansion vs. contraction.
• Universality of Hysteresis: The findings suggest that any modification of gravity (curvature, torsion, or non-metricity) that alters the coupling between geometry and matter in time-dependent backgrounds will generically induce hysteresis and irreversibility.
• Fate of Cyclic Universes: The cumulative dissipation implies that cyclic universes in $f(T)$ gravity do not settle into a stable periodic attractor. Instead, they undergo secular evolution, potentially leading to amplitude decay, thermalization, or transitions to different evolutionary regimes over infinite time.
• Observational Prospects: The systematic variation in cycle parameters could leave distinctive imprints on primordial observables (e.g., gravitational waves or CMB non-Gaussianities), offering potential future tests for bouncing scenarios.

In conclusion, the paper successfully bridges teleparallel gravity and thermodynamics, proving that torsion acts as an effective geometric dissipative channel, making cosmic hysteresis a robust and universal feature of modified gravity cyclic cosmologies.