Relativistic stellar modeling with perfect fluid core and anisotropic envelope fluid

This study investigates the stability of relativistic core-envelope stellar models with anisotropic fluids and demonstrates that the distortion energy accumulated through anisotropy-induced density perturbations can reach magnitudes comparable to gamma-ray bursts, thereby suggesting a potential physical connection between starquakes in self-bound compact stars and gamma-ray bursts.

The Gist

Imagine a super-dense star, such as a neutron star, not as a solid, homogeneous rock ball, but as a multi-layered dessert. This article treats the star as having a core (the dense, homogeneous center) surrounded by a crust (a slightly different, more complex outer layer).

Here is the story of what the authors found, explained simply:

1. The star is like a pressurized balloon with a twist

Normally, scientists imagine the pressure inside a star pushing outward equally in all directions, like air in a perfectly round balloon. However, this article suggests that in the outer layer (the crust) of these super-dense stars, the pressure is anisotropic.

Think of it like a rubber band wrapped around a ball. When you squeeze the ball, the rubber band pushes back harder in the direction it is wrapped (tangential) than in the direction you are pressing (radial). The authors propose that the outer shell of these stars acts like this rubber band, with the "sideways" pressure being slightly stronger than the "up-and-down" pressure.

2. The concept of "cracking"

The authors use a concept called "cracking" to investigate whether the star is stable. Imagine a patch of dry mud. When the mud dries unevenly, cracks form because different parts shrink or expand at different rates.

In the star, this means that if the "sideways" pressure and the "up-and-down" pressure behave differently when the star wobbles or changes density, a situation arises where the material wants to move in opposite directions.

• The analogy: Imagine two people holding a heavy rope. If one pulls slightly harder than the other, the rope snaps or slips. In the star, this means that if the "sound waves" (which transmit pressure) travel in different directions at different speeds, a "crack" or fault forms in the star's crust.

3. The "starquake" and energy release

The article suggests that these stars are like overstretched springs.

• Because the outer layer has this additional "sideways" pressure, energy builds up in the crust, just as tension builds in a stretched rubber band or a fault in the Earth's crust before an earthquake.
• The authors calculated that if this tension is suddenly released (a starquake), an enormous amount of energy could be released.
• The scale: They found that even a tiny pressure difference (so small it is almost invisible) could release an amount of energy equivalent to $10^{50}$ ergs. To put this in perspective, the article notes that this is roughly the amount of energy released in a Gamma-Ray Burst (GRB) or a massive outburst from a magnetar. It is as if the Sun were to release all the energy it will produce over its entire 10-billion-year life in just a few seconds.

4. How they did it

The researchers used a mathematical model (the TRV model) to simulate a star with a perfect fluid core and a "rubber-band-like" anisotropic crust.

• They checked the "speed of sound" inside the star. If sound travels sideways faster than up-and-down, the star is potentially unstable and prone to cracking.
• They found that their model classifies the star as potentially stable (it will not collapse immediately), but it does build up tension.
• They calculated that if the star "cracks" (an earthquake), the released energy matches the massive energy outbursts we observe coming from the depths of space.

5. The conclusion

The article proposes a new way to understand why some stars suddenly flash with intense gamma rays.

• The cause: A tiny imbalance in pressure between the "sideways" and "up-and-down" directions in the outer shell of the star.
• The effect: This imbalance stores deformation energy. When the star finally "cracks" or reorganizes (a starquake), this stored energy is released in a massive outburst.
• The connection: This mechanism could explain the origin of some of the most energetic events in the universe, such as Gamma-Ray Bursts, linking the tiny physics of pressure inside a star to the colossal explosions we see throughout the galaxy.

In short: The authors propose that these super-dense stars are like tension-filled balloons that, when they finally burst or crack due to internal pressure differences, release enough energy to illuminate the entire universe for a brief moment.

Technical Summary

Technical Summary: Relativistic Star Modeling with a Perfect Fluid Core and an Anisotropic Envelope Fluid

Problem Description
The article addresses the stability of compact star structures and specifically investigates the effects of density perturbations and local anisotropy on self-bound stars within the framework of General Relativity. A primary focus lies on the potential for "cracking" (overturning) in anisotropic matter configurations, a phenomenon proposed as a mechanism for starquakes. These starquakes are hypothesized as precursors to high-energy astrophysical events such as Soft Gamma Repeaters (SGRs) and Gamma-Ray Bursts (GRBs), which release energies comparable to the entire 10-billion-year lifespan of the Sun. The study aims to determine whether anisotropic pressure in the envelope of a compact star can accumulate sufficient strain energy to trigger such events.

Methodology
The authors employ a core-envelope model for superdense stars, specifically the Thomas-Ratanpal-Vinodkumar (TRV) model. The methodology comprises the following steps:

1. Geometric Framework: The interior of the star is modeled using a spherically symmetric line element in Schwarzschild coordinates. The core ($0 \le r \le R_C$) is treated as a perfect fluid (isotropic pressure), while the envelope ($R_C < r \le R_E$) is modeled with an anisotropic fluid ($p_r \neq p_\perp$).
2. Field Equations: Einstein's field equations are solved using an ansatz with pseudo-spheroidal geometry. The authors derive generalized Tolman-Oppenheimer-Volkoff (TOV) equations that incorporate an anisotropic stress tensor.
3. Equation of State (EoS): A geometrically derived equation of state in the quadratic form $p = \gamma + \alpha\rho + \beta\rho^2$ is used and fitted to density variation parameters ($\lambda$).
4. Stability Analysis (Cracking): The study utilizes the concept of "cracking" introduced by Herrera et al. Stability is assessed by analyzing the difference between the squares of the tangential ($v_{s\perp}$) and radial ($v_{sr}$) sound speeds. The condition $-1 \le v_{s\perp}^2 - v_{sr}^2 \le 1$ is used to identify potentially stable ($v_{s\perp}^2 \ge v_{sr}^2$) and unstable ($v_{s\perp}^2 < v_{sr}^2$) regions.
5. Numerical Estimation: The authors numerically integrate the TOV equations to calculate stellar parameters (mass, radius, gravitational energy, moment of inertia) for varying anisotropy magnitudes ($\Delta$) in the range of $10^{-8}$ to $10^{-6}$. They compute the differences in these parameters between isotropic and anisotropic cases to estimate the energy release during a hypothetical starquake.

Main Contributions and Results

• Stability via Sound Speeds: The study confirms that for the TRV model with the chosen parameters, the condition $-1 \le v_{s\perp}^2 - v_{sr}^2 \le 0$ is satisfied throughout the entire anisotropic envelope region. This suggests that the configuration is potentially stable against cracking under the specific perturbations considered, although the presence of anisotropy creates a region where strain energy accumulates.
• Estimation of Energy Release: The authors calculate the difference in gravitational energy ($\Delta E_g$) and the difference in the moment of inertia ($\Delta I$) resulting from changes in the anisotropy parameter. They note that for a self-bound star with a mass of approximately $1.43 M_\odot$ and a tangential pressure slightly exceeding the radial pressure (anisotropy parameter $\Delta \approx 10^{-6}$), the energy released during a quake-like event can reach $\sim 10^{50}$ erg.
• Correlation with Observations: The calculated energy scale ($10^{50}$ erg) proves to be of the same order of magnitude as the energy associated with giant gamma-ray bursts and SGR superflares (referenced in Table 1 of the article, e.g., GRB 070714b at $1.2 \times 10^{51}$ erg).
• Structural Properties: The model yields mass-radius relations consistent with observational constraints from binary radio pulsars, the GW170817 event, and NICER measurements of PSR J0030+0451. The study notes that at lower masses, the radius decreases monotonically, a characteristic feature of strange/quark stars.

Significance and Claims
The article claims that the proposed core-envelope model with an anisotropic envelope offers a viable theoretical framework for understanding the origin of GRBs and SGRs. The authors posit that:

1. Strain Accumulation: Anisotropy in the envelope leads to the accumulation of strain energy. When the tangential pressure is somewhat more significant than the radial pressure, this stored energy can be released during a starquake.
2. Precursor Mechanism: Self-bound compact stars with anisotropic envelopes are potential precursors for starquakes. The release of accumulated stress energy could explain the high-energy bursts observed in magnetars and GRBs.
3. Glitches and Spin Dynamics: The study suggests that changes in the anisotropy parameter could lead to changes in the star's moment of inertia, potentially causing pulsar glitches ($\Delta \Omega/\Omega$) and influencing spin dynamics.

The authors conclude that while the model suggests a mechanism for energy release comparable to observed astrophysical phenomena, the stability of the configuration depends on the specific relationship between radial and tangential sound speeds, which in their model indicates a potentially stable structure capable of storing the necessary energy for such events.