Phase space analysis of Bianchi III Universe with $f(R,T)$ gravity theory

This study performs a phase space analysis of the Bianchi type III universe within three different $f(R,T)$ gravity models, revealing that while all models reproduce the standard cosmological evolution sequence, only two are physically viable as the third model suffers from unbounded energy density issues.

The Gist

Imagine the Universe as a giant, expanding balloon. For a long time, scientists thought this balloon was perfectly round and smooth, expanding at the exact same speed in every direction. This is the standard "ΛCDM" model, the gold standard of modern cosmology.

However, this paper asks a bold question: What if the balloon isn't perfectly round? What if it's slightly squashed or stretched in one direction?

The authors, Pranjal Sarmah and Umananda Dev Goswami, investigate this "squashed" universe using a specific shape called the Bianchi III model. Think of this not as a perfect sphere, but as a slightly lopsided, egg-shaped balloon. They want to see how this shape changes the history of the Universe when we use a new, fancy set of rules for gravity called $f(R, T)$ gravity.

The New Rules of Gravity ($f(R, T)$)

In our everyday life, gravity is described by Einstein's General Relativity. But Einstein's rules might be incomplete. Scientists have proposed "Modified Gravity" theories to explain things like Dark Energy (the force pushing the universe apart) without needing to invent invisible stuff.

The $f(R, T)$ theory is like a new recipe for gravity.

• $R$ represents the curvature of space (how bent the fabric of the universe is).
• $T$ represents the matter and energy inside that space.
• In standard gravity, these two are separate. In this new theory, they are mixed together in a smoothie. The theory suggests that the presence of matter ($T$) can actually change how gravity ($R$) works.

The Experiment: A Phase Space Analysis

To study this, the authors didn't just run a simulation; they treated the Universe like a dynamical system.

The Analogy: The Roller Coaster Map
Imagine the history of the Universe as a roller coaster track.

• The Track: Represents the path the Universe takes from the Big Bang to today.
• The Cars: Represent different eras:
  1. Radiation Era: The very hot, fast-moving beginning (like a car zooming at the top).
  2. Matter Era: When stars and galaxies formed (the car cruising along the middle).
  3. Dark Energy Era: The current era where expansion is speeding up (the car shooting up a steep hill).
• Phase Space: This is a giant map showing all possible paths the roller coaster could take.
• Fixed Points: These are the "stations" on the map where the car stops or settles. The scientists looked for these stations to see if the Universe naturally settles into the Radiation, Matter, and Dark Energy phases we observe.

They tested three different recipes (models) for the $f(R, T)$ gravity theory to see which one keeps the roller coaster on a safe, logical track.

The Results: Three Recipes, Three Outcomes

1. Recipe A: $f(R, T) = \alpha R + \beta f(T)$

• The Verdict: Pass.
• The Story: This recipe works beautifully. When they plotted the roller coaster map, the Universe followed the exact same path as the standard model: it started hot (Radiation), cooled down to form stars (Matter), and then started speeding up again (Dark Energy).
• The Twist: Because the Universe is "squashed" (anisotropic), the roller coaster doesn't quite reach the perfect "100% full" mark on the map. There's a tiny gap. The authors explain this gap as "Shear Energy."
  • Analogy: Imagine stretching a rubber sheet. If you stretch it unevenly, the sheet stores extra tension. That tension is the "Shear Energy." In this model, the Universe's lopsided shape stores a little bit of extra energy that the standard round-universe model doesn't account for.

2. Recipe B: $f(R, T) = R + 2f(T)$

• The Verdict: Pass.
• The Story: This is a simpler version of the first recipe. It behaves almost exactly the same way. It successfully guides the Universe through all three eras (Radiation $\to$ Matter $\to$ Dark Energy).
• The Twist: Just like the first recipe, the "squashed" shape of the universe causes the roller coaster to settle slightly below the perfect mark, again due to that stored "Shear Energy."

3. Recipe C: $f(R, T) = (\zeta + \eta \tau T)R$

• The Verdict: Fail.
• The Story: This recipe is a bit more complex and non-linear. When they tried to run the simulation, the roller coaster went off the rails.
• The Problem:
  • Negative Energy: At one point, the math suggested the energy density was negative. In physics, this is like saying a car has "negative weight." It's physically impossible and breaks the laws of the universe.
  • Weird Pressure: During the "Matter Era" (when stars should be forming), the math suggested the pressure was negative, which implies the universe should be accelerating right then. But we know the universe was slowing down during that time to let stars form.
• Conclusion: This specific recipe is broken. It cannot describe our Universe, even if the Universe is lopsided.

The Big Takeaway

1. Shape Matters, but Theory Matters More: The "squashed" shape of the Universe (Bianchi III) does change the details slightly (adding that "Shear Energy" gap), but it doesn't break the story. The Universe still evolves from hot to cold to speeding up.
2. Not All Gravity Theories Are Created Equal: Just because a theory sounds fancy doesn't mean it works. Two of the recipes they tested fit the data perfectly, while the third one led to nonsense (negative energy).
3. The Future: The authors suggest that while our current models work well, we need better telescopes (like the upcoming Thirty Meter Telescope) to look back at the very early Universe. We need to see if that "squashed" shape was more pronounced back then, which would help us confirm if these new gravity theories are truly the right ones.

In a nutshell: The authors took a lopsided universe model and tested it against new gravity rules. Two of the rules worked great, proving that even a "squashed" universe follows the same general story as our standard round one. One rule failed miserably, showing that not every new idea in physics is a good idea.

Technical Summary

1. Problem Statement

Standard cosmology relies on the Cosmological Principle, assuming the Universe is homogeneous and isotropic (described by the FLRW metric and the $\Lambda$CDM model). However, observational data (e.g., CMB anomalies, Hubble tension, large-scale structure) suggest potential deviations, including cosmic anisotropies. While modified gravity theories like $f(R, T)$ (where the gravitational Lagrangian depends on the Ricci scalar $R$ and the trace of the energy-momentum tensor $T$) have been extensively studied in isotropic backgrounds, their application to anisotropic backgrounds remains under-explored.

The specific problem addressed is the lack of dynamical system analysis for the Bianchi Type III (BIII) metric within the $f(R, T)$ framework. The BIII metric is particularly relevant as it includes an exponential term and allows for hyper-spherically curved geometry, making it a suitable candidate for studying anisotropic evolution. The authors aim to determine if specific $f(R, T)$ models can physically describe the evolution of an anisotropic Universe and whether they remain consistent with standard cosmological phases (radiation, matter, dark energy).

2. Methodology

The study employs a dynamical system analysis approach to investigate the stability and evolutionary tracks of the Universe.

• Metric and Framework: The authors utilize the Bianchi Type III metric:
  $$ds^2 = -dt^2 + a_1^2(t) dx^2 + a_2^2(t) e^{-2mx} dy^2 + a_3^2(t) dz^2$$
  where $m$ is a constant. They assume a perfect fluid energy-momentum tensor and consider the condition $m \neq 0$ with directional Hubble parameters $H_1 = H_2 \neq H_3$ to maintain the BIII geometry.
• Gravity Theory: The analysis is conducted within $f(R, T)$ gravity. The field equations are derived from the modified Einstein-Hilbert action. Crucially, the authors account for the non-conservation of the energy-momentum tensor ($\nabla_\mu T^{\mu\nu} \neq 0$) inherent to $f(R, T)$ theory, leading to a modified continuity equation.
• Models Investigated: Three specific forms of $f(R, T)$ are tested:
  1. Model I: $f(R, T) = \alpha R + \beta f(T)$ (with $f(T) = \lambda T$).
  2. Model II: $f(R, T) = R + 2f(T)$ (a reduced version of Model I).
  3. Model III (Special Case): $f(R, T) = (\zeta + \eta \tau T)R$.
• Dynamical Variables: The system is converted into an autonomous set of differential equations using dimensionless density parameters:
  • $X = \Omega_m$ (Matter density)
  • $Y = \Omega_r$ (Radiation density)
  • $Z = \Omega_\Lambda$ (Dark energy density)
  • $S$: A variable related to the shear scalar and anisotropy.
• Analysis Techniques:
  • Fixed Point Analysis: Identifying critical points ($P_i$) where the derivatives of the dynamical variables vanish.
  • Stability Analysis: Calculating eigenvalues of the Jacobian matrix at critical points to determine stability (saddle points, stable nodes, etc.).
  • Phase Portraits: Plotting trajectories ($X$ vs. $Y$) to visualize the heteroclinic connections between different cosmological eras.
  • Parameter Constraints: Using constrained values of model parameters derived from observational data (Hubble, BAO, Pantheon+, CMB) via Bayesian inference (referenced from their previous work [74]).

3. Key Contributions

• Extension to Anisotropic Backgrounds: The paper successfully constructs a dynamical system for the BIII metric within $f(R, T)$ gravity, a niche area compared to standard FLRW analyses.
• Quantification of Anisotropy vs. Model Effects: The study distinguishes between deviations caused by the specific $f(R, T)$ model parameters and those caused by the intrinsic anisotropy of the BIII metric.
• Identification of Physical Viability: The authors provide a rigorous test for the physical viability of specific $f(R, T)$ models in anisotropic settings, highlighting that not all mathematically valid models are physically consistent.
• Shear Energy Interpretation: The paper interprets the deviation of the sum of density parameters from unity ($X+Y+Z \neq 1$) as "shear energy" stored in the spacetime deformation due to anisotropy.

4. Results

Model I ($f(R, T) = \alpha R + \beta f(T)$) and Model II ($f(R, T) = R + 2f(T)$)

• Evolutionary Path: Both models successfully reproduce the standard cosmological sequence: Radiation-dominated $\to$ Matter-dominated $\to$ Dark Energy-dominated. This is evidenced by heteroclinic trajectories connecting the critical points $P_3 \to P_2 \to P_1$.
• Critical Points:
  • The fixed points are shifted from the standard unit values (e.g., $\Omega_m \approx 0.88$ instead of $1.0$ in the matter era).
  • This shift is attributed to the anisotropic background and the specific model parameters.
  • The sum of density parameters $X+Y+Z < 1$, with the deficit interpreted as shear energy.
• Equation of State: The effective equation of state ($\omega_{eff}$) in the matter-dominated phase is slightly negative (e.g., $\approx -0.027$), indicating a minute influence of dark energy or model effects, but generally consistent with standard expectations.
• Anisotropy vs. Model Role: By comparing anisotropic ($\gamma \neq 1$) and isotropic ($\gamma = 1$) cases, the authors find that the model parameters have a more dominant role in the evolution than the anisotropy itself, although anisotropy causes a measurable shift in the density parameter values.

Model III ($f(R, T) = (\zeta + \eta \tau T)R$)

• Physical Inconsistency: This model fails the dynamical system test for the BIII metric.
  • Negative Energy Density: In the radiation-dominated phase ($P_{c3}$), the total contribution $Z_t$ (Dark energy + anisotropy) becomes negative, which is physically impossible.
  • Unphysical Equation of State: In the matter-dominated phase ($P_{c2}$), $\omega_{eff}$ is negative, implying accelerated expansion during the matter era, which contradicts standard cosmological history.
• Conclusion: This model is deemed unsuitable for studying the evolution of the Universe with an anisotropic background.

5. Significance and Conclusion

• Validation of Models: The study confirms that linear $f(R, T)$ models (Model I and II) are robust enough to describe an anisotropic Universe, aligning with the $\Lambda$CDM evolutionary track while accommodating anisotropic effects.
• Constraint on Anisotropy: The research demonstrates that while anisotropy affects the specific values of density parameters (shifting them below unity), it does not disrupt the fundamental heteroclinic evolution path of the Universe.
• Model Selection Criteria: The paper establishes that non-linear $f(R, T)$ models (like Model III) must be scrutinized carefully, as they may lead to unphysical results (negative energy densities) in anisotropic geometries, even if they appear viable in isotropic limits.
• Future Outlook: The authors suggest that upcoming observational projects (TMT, ELT, CTA) will be crucial in detecting early-universe anisotropies to further validate or constrain these modified gravity models.

In summary, the paper provides a rigorous phase-space analysis demonstrating that while $f(R, T)$ gravity can accommodate an anisotropic Bianchi III Universe, the physical viability of the theory is highly dependent on the specific functional form of $f(R, T)$. Linear models show promise, whereas certain non-linear forms introduce physical inconsistencies.