⚛️ general relativity

Dirac stars in modified $F(R)$ gravity

This paper investigates Dirac mini-stars in modified $R^2$ gravity sourced by linear and nonlinear spinor fields, demonstrating that increasing the nonlinearity parameter significantly alters the configurations' mass profiles, energy density, and metric functions compared to general relativity.

Original authors: Zhuldyz Alibekova, Vladimir Dzhunushaliev, Vladimir Folomeev, Sergei Odintsov

Published 2026-10-14

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Zhuldyz Alibekova, Vladimir Dzhunushaliev, Vladimir Folomeev, Sergei Odintsov

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the vast, quiet architecture of the universe, gravity is the invisible force that sculpts everything from falling apples to colliding galaxies. For nearly a century, our best description of this force has been Albert Einstein's theory of general relativity, which treats gravity not as a pull, but as a curvature in the fabric of space and time caused by matter and energy. However, this theory, while spectacularly successful, hits a wall when we look at the very beginning of the universe or the most extreme environments, like the centers of black holes. It also struggles to explain why the universe is expanding faster and faster without inventing a mysterious, unseen substance called dark energy. To solve these puzzles, physicists have proposed "modified gravity" theories, which tweak Einstein's equations to include extra geometric features. One such idea, known as $R^2$ gravity, suggests that the curvature of space itself has a squared component that becomes significant in high-energy environments, offering a way to explain cosmic acceleration without dark energy.

Parallel to this search for new gravity laws is the study of "Dirac stars." These are not stars made of burning gas like our Sun, but rather compact, self-gravitating clouds of fundamental particles called fermions, which make up all ordinary matter. In these theoretical objects, the particles are held together by their own gravity, balanced by the quantum pressure that prevents them from collapsing. While general relativity has long been used to model these objects, scientists have recently begun to ask: what happens if we build these stars using the rules of modified gravity instead? This question is the heart of a new study by researchers Zhuldyz Alibekova, Vladimir Dzhunushaliev, Vladimir Folomeev, and Sergei Odintsov, who set out to construct and analyze these exotic objects within the framework of $R^2$ gravity.

The researchers focused on a specific type of modified gravity where the mathematical description of space curvature includes a term proportional to the square of the curvature itself. They modeled a star composed of a spinor field—a mathematical description of particles like electrons or neutrons that includes both linear behavior and a "nonlinear" self-interaction. In simple terms, the particles in this model not only feel the pull of gravity but also interact with each other in a way that can be either attractive or repulsive, depending on a specific parameter in the equations. The team used powerful computer simulations to solve the complex equations governing these systems, comparing the results against the standard predictions of Einstein's general relativity. They varied the strength of the particle interactions to see how the star's shape, density, and total mass would change under different conditions.

The simulations revealed that the choice of gravity theory fundamentally alters the nature of these stars. When the researchers compared the modified gravity models to the standard ones, they found that for the same conditions, the total mass of a Dirac star in the modified theory was consistently lower than in general relativity. This difference became even more dramatic when they increased the strength of the nonlinear interaction between the particles. In the standard theory, increasing this interaction changed the star's profile only slightly, but in the modified gravity model, it caused significant shifts in the distribution of mass and the behavior of the particles themselves. The "density" of the star, which describes how much mass is packed into a given space, developed a pronounced peak at a certain distance from the center when the interaction was strong, a feature that was absent or much less distinct in the standard models.

Perhaps the most striking finding was how the modified gravity theory responded to the frequency at which the particles oscillate within the star. In general relativity, the relationship between the star's mass and this oscillation frequency follows a predictable curve. However, in the modified $R^2$ gravity, this relationship warped significantly as the nonlinearity parameter increased. The researchers observed that for very strong interactions, the entire structure of the solution branches changed, suggesting that the internal balance of forces holding the star together is governed by different rules. The study also confirmed that these modified stars remain stable and "asymptotically flat," meaning that far away from the star, space returns to a normal, flat state, just as it does in our everyday experience. This is a crucial check, as it ensures the models are physically realistic and do not produce impossible distortions in the universe at large distances.

The work also shed light on the role of the curvature of space itself. In standard general relativity, the curvature is directly tied to the amount of matter present. In the modified theory, the curvature is influenced by additional terms arising from the modification of the gravity equations. The researchers found that the curvature of space inside these stars reached its maximum value at the center and then smoothly decayed to zero as one moved outward, confirming the stability of their solutions. The study did not claim to have discovered a new type of star in the sky, but rather provided a rigorous theoretical map of how such objects would behave if the laws of gravity were slightly different from what Einstein proposed. By showing that modified gravity can significantly change the mass, density, and internal structure of these compact objects, the paper offers a new perspective on how fundamental particles might organize themselves under extreme conditions, providing a testing ground for theories that could eventually explain the deepest mysteries of our cosmos.

Technical Summary: Dirac Stars in Modified $R^2$ Gravity

Problem Statement
The paper investigates the properties of compact, self-gravitating configurations formed by spinor fields, known as Dirac stars (or Dirac mini-stars), within the framework of modified gravity. Specifically, the authors analyze these systems in $R^2$ gravity (a specific class of $F(R)$ gravity), where the gravitational action includes a quadratic curvature term $R^2$ in addition to the standard Einstein-Hilbert term. The study contrasts these modified gravity solutions with those derived in General Relativity (GR). The spinor fields are modeled as both linear and nonlinear, with the nonlinearity introduced via a power-law term in the Lagrangian. The primary objective is to elucidate how the additional geometric degrees of freedom in $R^2$ gravity, combined with spinor self-interaction, alter the structure, stability, and mass properties of these compact objects.

Methodology
The authors formulate the problem using a total action comprising the $F(R)$ gravitational term and the action for a classical spinor field $\psi$ . The gravitational sector is defined by $F(R) = R + \alpha R^2$ , where $\alpha$ is a coupling parameter. The spinor Lagrangian includes a mass term and a nonlinear interaction term $F(S) \propto -\lambda (\bar{\psi}\psi)^{k+1}$ , with the study focusing on the case $k=1$ (yielding a quartic interaction).

Key methodological steps include:

Mathematical Formulation: The authors derive the coupled system of equations consisting of the modified Einstein field equations and the curved-spacetime Dirac equation. To handle the mathematical consistency of classical spinors (which are Grassmann numbers), the spinor is decomposed into an ordinary complex-valued function and a Grassmann number, with integration over the Grassmann variables performed to yield a real-valued action.
Ansatz and Symmetry: A spherically symmetric metric is assumed, and a stationary ansatz for the spinor field is employed. This ansatz ensures a static spacetime with vanishing total spin and a spherically symmetric energy-momentum tensor.
Dimensionless Formulation: The system is reduced to a set of five coupled ordinary differential equations (ODEs) for the radial functions describing the spinor components ( $f, g$ ), the metric functions ( $m, \sigma$ ), and the scalar curvature ( $\Sigma$ ). These equations are rendered dimensionless using the Compton wavelength of the spinor mass.
Numerical Solution: The system is treated as a boundary-value problem. Regularity conditions are imposed at the origin ( $x=0$ ), and asymptotic flatness conditions are applied at infinity. The spinor frequency $\Omega$ is treated as an eigenvalue. The authors employ the shooting method to integrate the equations from a point near the origin to a finite boundary where the fields decay.
Parameter Space: Numerical simulations are conducted for various values of the nonlinearity parameter $\lambda$ (ranging from 0 to 20) and the modified gravity parameter $\alpha = -10^{13} \text{ cm}^2$ . The negative sign of $\alpha$ is chosen to ensure asymptotically decaying solutions for the metric and mass functions.

Key Contributions and Results
The study constructs two-dimensional (mass vs. frequency) and three-dimensional (mass vs. frequency vs. central spinor value) profiles for Dirac mini-stars in both GR and $R^2$ gravity. The main findings are:

Mass Reduction in Modified Gravity: For comparable frequencies and nonlinearity parameters ( $\lambda=0, 1$ ), the total mass of the configuration in $R^2$ gravity is consistently lower than in GR.
Impact of Nonlinearity: Increasing the nonlinearity parameter $\lambda$ significantly alters the system's properties. While differences between GR and $R^2$ gravity are modest for small $\lambda$ , the effects become pronounced for large values ( $\lambda=10, 20$ ).
Spinor and Metric Profiles:
- Spinor Functions: In $R^2$ gravity, increasing $\lambda$ leads to substantial changes in the shape and amplitude of the spinor functions $f(x)$ and $g(x)$ .
- Metric Functions: The mass function $m(x)$ and the metric function $\sigma(x)$ exhibit distinct profiles in $R^2$ gravity compared to GR. Within the modified gravity framework, higher $\lambda$ values induce noticeable deviations in these metric functions.
- Scalar Curvature: Unlike in GR, where the scalar curvature is algebraically related to the trace of the energy-momentum tensor, in $R^2$ gravity, the curvature is determined by differential equations involving additional terms from the modification. The scalar curvature $\Sigma$ peaks at the center and decays to zero at infinity, confirming asymptotic flatness.
Energy Density: The effective energy density $\rho_{\text{eff}}$ includes contributions from the spinor field and a geometric term arising from the $R^2$ modification. As $\lambda$ increases, the distribution of effective energy density changes significantly, developing a pronounced maximum at a finite distance from the center for large $\lambda$ .
Mass-Frequency Relations: The dependence of the total mass on the oscillation frequency ( $\bar{M}(\Omega)$ ) is altered by both the modification of gravity and the nonlinearity parameter. The 3D plots reveal that apparent overlaps in 2D projections are artifacts of the projection, and the solution branches are distinct in the full parameter space $(\Omega, g_c, \bar{M})$ .

Significance
The paper claims that the modification of gravity to $R^2$ form, combined with nonlinear spinor self-interactions, fundamentally alters the physical properties of Dirac stars. The study demonstrates that the additional geometric degrees of freedom in $F(R)$ gravity do not merely perturb GR solutions but can significantly change the profiles of spinor and metric functions, the effective energy density, and the mass-frequency relationship. The most dramatic effects are observed when the nonlinearity parameter is large. These results highlight the sensitivity of compact fermionic configurations to the underlying theory of gravity and the nature of spinor self-interactions, providing a theoretical basis for understanding how modified gravity might manifest in systems supported by fundamental fermion fields.

More like this