Dynamical systems analysis of an Einstein-Cartan ekpyrotic nonsingular bounce cosmology

This paper presents a phenomenological Einstein-Cartan ekpyrotic model where a Weyssenhoff fluid's spin-torsion term, coupled to a scalar field with a steep exponential potential, dynamically damps shear and triggers a nonsingular bounce in a homogeneous universe, with numerical analysis confirming the stability of these trajectories against chaotic behavior.

The Gist

The Big Picture: Avoiding the "Big Crunch"

Imagine the history of our universe as a giant movie. The standard version (the Big Bang theory) starts with a scene where everything is squished into an infinitely small, infinitely hot point—a singularity. It's like a movie that starts with the camera lens shattering; we can't see what happened before that moment because the physics breaks down.

This paper asks a bold question: What if the universe didn't start with a shattering lens, but rather with a "bounce"?

Think of the universe like a giant rubber ball. In the standard story, the ball falls from the sky and smashes into the ground (the Big Bang). In this paper's story, the ball falls, hits the ground, squishes down, and then bounces back up without ever breaking. The author, Jackson Stingley, builds a mathematical model to show how this bounce could happen without the universe collapsing into a singularity.

The Cast of Characters

To make this bounce work, the author uses three main "actors" in his cosmic play:

1. The Scalar Field (The Director):
   Imagine a invisible force field that guides the universe. In this model, it's a "steep hill" that the universe rolls down.

   • The Ekpyrotic Phase: When the universe is shrinking (contracting), this field acts like a heavy weight on a wobbly table. It forces the universe to become perfectly smooth and flat, smoothing out any wrinkles or bumps (anisotropy) that might have formed. It's like a heavy hand pressing down on a crumpled piece of paper until it's perfectly flat again.

2. The Spin-Torsion Fluid (The Spring):
   This is the star of the show. In standard physics, matter just gets denser as you squeeze it. But in this model (Einstein-Cartan gravity), the "spin" of tiny particles acts like a super-spring.

   • As the universe shrinks, this spring gets tighter and tighter. Eventually, it gets so stiff that it pushes back harder than gravity pulls in. Instead of crushing the universe into a singularity, this springy force stops the collapse and pushes the universe back out.

3. The "Softening" Switch (The Brake):
   Here is the tricky part. If the "Director" (the scalar field) keeps rolling down the steep hill, it gets too heavy, and the "Spring" (torsion) can't push back hard enough to cause a bounce.

   • The author introduces a "softening" mechanism. Imagine the steep hill suddenly turns into a gentle slope or a flat plateau right before the bottom. This slows down the Director just enough so that the Spring can take over and trigger the bounce.

The Plot: How the Universe Bounces

The paper uses complex math (dynamical systems) to simulate the universe's life cycle. Here is the sequence of events:

1. Expansion: The universe starts big and expanding (like our current universe).
2. Contraction: Eventually, it stops expanding and starts shrinking.
3. The Smoothing: As it shrinks, the "Director" (scalar field) takes over. It acts like a cosmic iron, pressing out all the wrinkles and making the universe perfectly smooth and uniform. This is crucial because, in normal physics, a shrinking universe usually gets messy and chaotic (like a crumpled ball of paper).
4. The Handoff: As the universe gets very small and dense, the "Softening Switch" flips. The Director slows down.
5. The Bounce: Now, the "Spring" (spin-torsion) kicks in. It hits a density where it becomes repulsive. It pushes back against the shrinking universe.
6. The Rebound: The universe stops shrinking, hits a minimum size (but never zero!), and starts expanding again.

The "Dynamical Systems" Part (The Map)

The author didn't just guess this would work; he built a massive 6-dimensional map (a phase space) to track the universe's journey.

• Think of this map like a topographical map of a mountain range.
• The "valleys" on the map are stable paths the universe can follow.
• The "peaks" are unstable paths where the universe might crash.
• The author calculated the "Lyapunov exponents" (a fancy way of measuring chaos). He found that on this map, the paths are not chaotic. The universe doesn't get lost in a maze of random jumps; it follows a smooth, predictable track toward the bounce.

The "Basin of Viability" (The Goldilocks Zone)

One of the most interesting findings is that this bounce doesn't happen for every setting. It requires a specific "Goldilocks" tuning.

• If the "softening" happens too early, the spring never gets strong enough to bounce.
• If it happens too late, the universe collapses before the spring can push back.
• The author found a specific "basin" (a region on his map) where the settings are just right. It's not a single point, but a whole region where the bounce works. This suggests the model is robust, not just a fluke.

What This Means for Us

• No Singularity: The Big Bang wasn't a "beginning" from nothing, but a bounce from a previous shrinking phase.
• Smoothness: The model explains why our universe is so smooth and uniform today (the "Director" smoothed it out during the contraction).
• No Chaos: Even though the universe was shrinking, it didn't turn into a chaotic mess. It stayed orderly.

The Caveats (The "But...")

The author is very honest about what this model doesn't do yet:

• It's a Homogeneous Model: It assumes the universe is the same everywhere (like a smooth soup). It doesn't yet account for galaxies, stars, or clumps of matter.
• Entropy: It doesn't fully solve the "entropy problem" (why time moves forward and why the universe doesn't just get messier with every bounce).
• Microphysics: The "spring" is treated as a phenomenological fluid (a mathematical shortcut) rather than being derived from the deepest laws of particle physics.

Summary Analogy

Imagine a trampoline.

1. Standard Big Bang: You drop a bowling ball, and it smashes through the trampoline into a hole in the ground.
2. This Model: You drop the bowling ball. As it hits the trampoline, the fabric stretches. But right before it tears, the fabric turns into a super-rigid, springy material (Spin-Torsion). The ball stops, the fabric pushes back, and the ball bounces high into the air.
3. The "Director": Before the ball even hits the trampoline, a giant hand (the Scalar Field) smoothed out the wrinkles in the fabric so the ball lands perfectly flat.

This paper provides the mathematical proof that such a "trampoline bounce" is possible within the rules of Einstein-Cartan gravity, offering a potential alternative to the idea that the universe began with a singularity.

Technical Summary

1. Problem Statement

The paper addresses the cosmological singularity problem (the Big Bang) and the challenges associated with cyclic cosmological models. Specifically, it seeks to construct a nonsingular bounce within a homogeneous, nearly Friedmann–Lemaître–Robertson–Walker (FLRW) background using Einstein–Cartan (EC) gravity.

Key challenges in this domain include:

• The BKL Instability: In General Relativity (GR), contracting universes are generically unstable to anisotropic shear (the Belinski–Khalatnikov–Lifshitz instability), leading to a chaotic singularity.
• Entropy Accumulation: Cyclic models often struggle with the accumulation of entropy across cycles.
• The "Stiff" Problem: In standard EC models, the spin-torsion term (which scales as $a^{-6}$) acts as a repulsive force that can cause a bounce. However, if the contracting phase is dominated by a "stiff" fluid or kinetic energy (equation of state $w=1$), the shear and spin densities scale identically ($a^{-6}$), meaning the shear may dominate or remain comparable, potentially preventing a clean, isotropic bounce.
• Need for a Unified Framework: Previous works often focused on microscopic spinor condensates or specific perturbative solutions rather than the global phase-space structure of a combined scalar-fluid-spin system.

2. Methodology

The author constructs a phenomenological model called the Einstein–Cartan Ekpyrotic Model (ECEM) and analyzes it using dynamical systems theory.

• Theoretical Framework:

  • Gravity: Einstein–Cartan gravity, where intrinsic fermion spin sources spacetime torsion.
  • Matter Content: A canonical scalar field $\phi$ with a specific potential $V(\phi)$, a barotropic fluid (matter/radiation), and a phenomenological Weyssenhoff fluid representing the spin-torsion sector.
  • Spin-Torsion Scaling: The spin-torsion energy density is modeled as $\rho_s \propto a^{-6}$ (stiff component) but with a negative effective contribution to the Friedmann equation ($\rho_{eff} = \rho_{tot} - \rho_s$), providing repulsion at high densities.
  • Scalar Potential: A "softening" potential that interpolates between a steep negative exponential (driving ekpyrotic contraction, $w_\phi \gg 1$) and a positive plateau (allowing $w_\phi < 1$).

• Dynamical System Formulation:

  • The system is extended to a six-dimensional phase space with variables: $(x, y, z, \Sigma, \Omega_k, \Omega_s)$, representing normalized kinetic/potential energy, matter density, shear, curvature, and spin-torsion density.
  • The equations of motion are recast in terms of the deceleration parameter $q$ to create a compact, autonomous system.
  • The author derives the full $6 \times 6$ Jacobian matrix to perform linear stability analysis at fixed points.

• Numerical Analysis:

  • Fixed Points: Analyzed for kinetic, fluid, scalar-dominated, and scaling regimes.
  • Lyapunov Exponents: Computed using the Benettin algorithm to test for chaotic behavior (maximal Lyapunov exponent $\lambda_{max}$).
  • Trajectory Integration: Numerical integration performed in both e-fold time ($N = \ln a$) for phase-space structure and cosmic time ($t$) to resolve the bounce singularity where $H \to 0$.
  • Parameter Scans: A two-parameter scan over the softening scale $(\phi_b, \Delta)$ to identify the "basin of viability" for nonsingular bounces.

3. Key Contributions

1. Six-Dimensional Phase Space Extension: The paper extends the standard Copeland–Liddle–Wands (CLW) formalism to include shear, curvature, and spin-torsion, providing a global view of the Einstein–Cartan ekpyrotic dynamics.
2. Explicit Bounce Mechanism: It demonstrates how a single scalar field with a softening potential can drive the entire cycle: scalar-dominated expansion $\to$ fluid era $\to$ ekpyrotic contraction $\to$ torsion-induced bounce.
3. Shear Suppression Proof: It rigorously shows that the ekpyrotic phase ($w_\phi \gg 1$) acts as a strong attractor, exponentially damping homogeneous shear ($\Sigma \to 0$) before the bounce, thereby taming the anisotropy instability that plagues standard contracting GR models.
4. Lyapunov Stability: The computation of maximal Lyapunov exponents on the constrained phase space reveals negative values ($\lambda_{max} < 0$), indicating no chaotic mixing in the homogeneous truncation, even when curvature modes are included.
5. Basin of Viability: The paper identifies a finite region in the softening parameter space $(\phi_b, \Delta)$ where the spin-torsion term successfully overtakes the scalar field to trigger a bounce, rather than the universe collapsing into a singularity.

4. Key Results

• Ekpyrotic Shear Damping: During the contraction phase driven by the steep negative potential, the shear variable $\Sigma$ decays exponentially as $\Sigma \propto a^{q-2}$. Since $q \gg 2$ in the ekpyrotic regime, the universe becomes highly isotropic before reaching the high-density bounce.
• The Bounce Condition: A nonsingular bounce ($H=0, \dot{H}>0$) occurs when the spin-torsion density $\rho_s$ equals the total matter/scalar density $\rho_{tot}$, provided the total equation of state satisfies $w_{tot} < 1$. The softening of the potential is crucial here; if $w_\phi$ remained $>1$ (stiff) until the equality point, the bounce would not occur.
• Phase Space Structure:
  • Scalar-dominated point: A stable late-time attractor for $\lambda^2 < 2$ (accelerated expansion).
  • Scaling point: A saddle point unstable to curvature.
  • Ekpyrotic regime: A strong attractor for isotropy during contraction.
• Lyapunov Analysis: For the tested parameter ranges, the maximal Lyapunov exponent is negative, confirming that the homogeneous dynamics are non-chaotic. This suggests the model is robust against small perturbations in the homogeneous sector.
• Parameter Tuning: The bounce is not generic; it requires the softening transition to occur at a density scale comparable to where the spin-torsion term becomes dominant. The paper quantifies this "tuning" as a finite "bounce basin" in the $(\phi_b, \Delta)$ plane (see Fig. 5).
• Curvature and Cyclic Behavior:
  • For $k=0$ (flat), the model produces a single bounce.
  • For $k=+1$ (closed), the model supports repeated turnaround and bounce cycles (Fig. 1), though the paper notes this does not yet solve the entropy accumulation problem.

5. Significance and Limitations

Significance:

• The work provides a mathematically consistent, phenomenological framework for a nonsingular bounce without introducing new propagating degrees of freedom beyond GR and matter.
• It bridges the gap between microscopic spinor condensate models and macroscopic dynamical systems, offering a clear picture of how ekpyrotic contraction and torsion-driven bounces interact globally.
• It demonstrates that the anisotropy problem (a major hurdle for cyclic models) can be solved within EC gravity via the ekpyrotic mechanism, provided the initial shear is subdominant to the spin sector.

Limitations & Future Work:

• Homogeneous Truncation: The analysis is confined to homogeneous backgrounds. It does not address the full BKL (inhomogeneous) dynamics, which could reintroduce chaos.
• Entropy: The model does not resolve the thermodynamic entropy accumulation problem inherent in cyclic cosmologies.
• Phenomenological Tuning: The softening potential and the alignment of the spin-torsion scale with the scalar transition are treated as tuned parameters rather than derived from a specific ultraviolet (UV) completion.
• Perturbations: The paper does not evolve scalar or tensor perturbations through the bounce, so predictions for the Cosmic Microwave Background (CMB) power spectrum or gravitational waves are not yet available.

In conclusion, Stingley's paper establishes that a tuned Einstein–Cartan ekpyrotic model can successfully realize a nonsingular, isotropic bounce in a homogeneous universe, offering a viable alternative to the Big Bang singularity within the realm of effective field theory, pending further work on perturbations and inhomogeneities.