⚛️ general relativity

Assessing the stability of ultracompact spinning boson stars with nonlinear evolutions

Using fully nonlinear numerical relativity simulations with two different Einstein equation formulations, the study finds no evidence of instability in ultracompact spinning boson stars with stable light rings over timescales of $10^4$ scalar mass units, even when subjected to various perturbations that decay slowly without triggering immediate collapse.

Original authors: Tamara Evstafyeva, Nils Siemonsen, William E. East

Published 2026-02-18

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Tamara Evstafyeva, Nils Siemonsen, William E. East

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

Imagine the universe is filled with mysterious, ultra-dense objects. For a long time, we thought the only things this dense were Black Holes—cosmic vacuum cleaners with such strong gravity that not even light can escape. But recently, scientists have started wondering: "What if there are other objects just as dense, but without the 'event horizon' (the point of no return)?"

These hypothetical objects are called Ultracompact Objects. One specific type is the Boson Star. Think of a Boson Star not as a ball of gas or rock, but as a giant, swirling cloud of invisible, heavy particles (like a cosmic cloud of dust) held together by its own gravity.

The Big Question: Are They Stable?

In 2023, a study suggested that these spinning Boson Stars might be unstable. The theory was that because these stars are so dense, they trap light in a "traffic circle" around them (called a Light Ring). The researchers thought that this trapped light would slowly build up energy, like water filling a bathtub with a slow leak, until the whole star eventually exploded or collapsed into a black hole.

This paper says: "Hold on, let's check that again."

The Experiment: A Cosmic Stress Test

The authors of this paper, Tamara, Nils, and William, decided to run the ultimate stress test. They used supercomputers to simulate these stars in a virtual universe.

The Setup: They created digital models of these spinning Boson Stars.
The Nudge: To see if the stars would fall apart, they gave them a little push.
- Sometimes, the "push" was just the tiny, unavoidable errors that happen when you do math on a computer (like rounding errors).
- Other times, they manually added a "kick" to the star, shaking it up or adding a bump to its shape.
The Watch: They watched these stars for a very long time (in computer time, this is about 10,000 "heartbeat" cycles of the star).

The Results: The Stars Are Tougher Than We Thought

Here is what they found:

No Collapse: Even after giving the stars a good shake, they didn't collapse. They didn't explode. They didn't turn into black holes.
The "Wobble" that Fades: When they shook the stars, the stars did wobble and vibrate. But instead of the vibrations getting louder and louder (which would mean instability), they slowly died down. It's like plucking a guitar string; it vibrates for a bit, but eventually, the sound fades away, and the string goes back to being still.
The "Fake" Instability: The researchers also discovered a tricky trap. In their computer simulations, sometimes the math itself started acting weird. It looked like the star was exploding, but it was actually just a glitch in the coordinate system (like a map that gets distorted as you zoom in). They had to be very careful to tell the difference between a real physical explosion and a computer glitch.

The Analogy: The Spinning Top

Imagine a spinning top.

The Old Theory: Someone claimed that if you spin a top fast enough, the air trapped around it would eventually push it over, making it fall.
This Paper's Finding: The authors spun the top (the Boson Star) thousands of times. They poked it, they shook it, and they let the computer's tiny errors nudge it. The top wobbled, but it kept spinning. It didn't fall over.

Why Does This Matter?

Black Hole Mimickers: If these stars are stable, it means they could actually exist in our universe. If they do, they might look exactly like black holes to our telescopes, but they would be different inside. This changes how we interpret the signals from gravitational waves (the "ripples" in space-time).
The "Light Ring" Mystery: The idea that trapped light causes an explosion was a hot topic. This paper suggests that while trapped light does exist, it doesn't necessarily lead to a disaster. The universe might be more stable than we feared.
Better Math: The paper also taught other scientists how to avoid "fake" computer errors when simulating these extreme objects, which helps everyone do better science in the future.

The Bottom Line

The authors didn't find the "smoking gun" that proves these stars are unstable. Instead, they found that Ultracompact Boson Stars are surprisingly resilient. They can take a beating, wobble for a while, and then settle back down, staying exactly as they were.

So, while the universe is full of mysteries, these specific cosmic clouds seem to be holding their ground, refusing to collapse under their own weight—at least for now.

1. Problem Statement

The paper addresses the dynamical stability of ultracompact spinning boson stars (BSs) that possess stable light rings (LRs).

Context: Ultracompact objects (UCOs) are proposed alternatives to black holes. Unlike black holes, which have only an unstable light ring, UCOs possess an inner stable light ring where null geodesics are trapped.
The Hypothesis: It has been conjectured that stable trapping of massless modes in the inner LR could lead to a nonlinear "light ring" instability. Theoretically, weak nonlinear effects could accumulate over time, causing the object to either collapse into a black hole or migrate to a less compact state.
Previous Work: A previous study (Ref. [14]) reported such an instability in ultracompact scalar boson stars on a timescale of $t\mu \sim 5000$ . However, other recent works found no evidence of such instability in spherically symmetric or axisymmetric cases.
Goal: To rigorously test the stability of these objects using fully nonlinear 3+1 and 2+1 numerical relativity simulations over extended timescales ( $t\mu \sim 10^4$ ), distinguishing between physical instabilities and numerical artifacts.

2. Methodology

The authors employed a robust, multi-faceted approach using two independent numerical relativity codes and different formulations of the Einstein equations.

A. Numerical Codes and Formulations

ExoZvezda:
- Based on the CCZ4 (Conformal and Covariant Z4) formulation.
- Uses the Moving Puncture gauge.
- Employs Adaptive Mesh Refinement (AMR) via the Chombo library.
- Utilizes 4th-order finite differencing and 6th-order Kreiss-Oliger dissipation.
GHC (Generalized Harmonic Code):
- Based on the Generalized Harmonic formulation.
- Uses a Stationary gauge (fixing source functions $H_a$ ).
- Includes constraint damping terms.
- Utilizes a compactified coordinate system to include spatial infinity.

B. Models and Initial Data

Primary Model (S005): An ultracompact spinning boson star ( $m=1$ , $\sigma_0=0.05$ ) previously reported as unstable in Ref. [14]. It has a compactness $C=0.33$ and two light rings.
Secondary Model (S02): A model with $m=3$ , $\sigma_0=0.2$ , and $C=0.38$ .
Perturbation Strategies:
1. Discretization Error: Evolving the star with no explicit perturbation, allowing numerical truncation error to act as the seed.
2. Azimuthal Perturbation: Modifying the scalar field as $\phi \to \phi e^{-2\epsilon \cos(n\phi)}$ to excite specific non-axisymmetric modes ( $n=1, 2$ ).
3. Radial/LR Perturbation: Adding a Gaussian perturbation centered on the stable light ring location: $|\phi| \to |\phi| + \epsilon e^{-(r-r_0)^2/\Delta^2}$ .

C. Diagnostics

To assess stability, the authors monitored:

Stationarity Measures: $\Phi = \max|\partial_t |\phi|^2|$ and $G = \max|\partial_t g_{tt}|$ .
Gauge-Invariant Quantity: $\Phi_{\text{max}} = \max(|\phi|) / \max(|\phi(t=0)|)$ .
Mode Decomposition: Azimuthal decomposition of $\Phi$ to track specific modes ( $\tilde{m}$ ).
Constraint Violations: Monitoring Hamiltonian and momentum constraint violations to distinguish physical instabilities from numerical gauge artifacts.

3. Key Contributions and Findings

A. Main Stability Results

No Evidence of Instability: For the S005 model, the authors found no evidence of the conjectured light ring instability on timescales up to $t\mu \sim 10^4$ .
Perturbation Behavior:
- When perturbed by discretization error, the system remained stationary.
- When subjected to explicit perturbations (azimuthal or radial), the perturbations exhibited slow decay or remained stable, even for amplitudes just below the threshold for immediate gravitational collapse.
- Black hole formation only occurred for "extortionately large" perturbations ( $\epsilon \ge 0.1$ for azimuthal, $\epsilon \ge 0.01$ for radial), which induced prompt collapse rather than a slow nonlinear instability.
Conclusion: The slow decay of linear modes does not lead to an eventual nonlinear instability within the simulated timescales.

B. Identification of Numerical Artifacts

A significant portion of the paper is dedicated to identifying and mitigating unphysical instabilities that could be mistaken for physical phenomena:

Gauge Instabilities:
- Using the Generalized Harmonic formulation with a stationary gauge, the authors observed an exponential growth in gauge-dependent quantities ( $\Phi$ ) for the S02 model.
- Crucial Distinction: While $\Phi$ grew exponentially, the gauge-invariant quantity $\Phi_{\text{max}}$ remained constant. Furthermore, constraint violations converged to zero with resolution. This confirmed the growth was a gauge artifact, not a physical instability.
- This instability was found to scale with the compactness of the star and is specific to the stationary gauge in the generalized harmonic formulation.
Constraint Damping:
- Simulations using the BSSNOK formulation (without constraint damping) showed a gradual, secular accumulation of constraint violations, leading to numerical drift.
- The CCZ4 formulation (with constraint damping) successfully suppressed these violations, proving that constraint damping is essential for long-term stability of ultracompact objects.

4. Significance

Challenging the Paradigm: The results directly contradict the findings of Ref. [14], suggesting that the previously reported instability may have been a numerical artifact or specific to a different setup.
Robustness of UCOs: The study implies that ultracompact spinning boson stars with stable light rings may be dynamically stable against the specific nonlinear "light ring" instability mechanism, at least on timescales of $10^4$ scalar mass units.
Methodological Guidance: The paper provides critical "recipes" for numerical relativity simulations of extreme spacetimes:
- The necessity of using gauge-invariant diagnostics to confirm physical instabilities.
- The superiority of CCZ4 with constraint damping over BSSNOK for these specific models.
- The danger of stationary gauges in the Generalized Harmonic formulation for highly compact objects.
Future Directions: The authors note that while no instability was found on $t\mu \sim 10^4$ , the possibility of a very weak instability on much longer timescales (due to slow logarithmic decay) remains open. They suggest future work should explicitly introduce test fields trapped in the light ring to study their gravitational backreaction.

Summary

This paper utilizes high-resolution, fully nonlinear numerical relativity to re-evaluate the stability of ultracompact spinning boson stars. By employing multiple codes, formulations, and rigorous diagnostic checks, the authors conclude that no nonlinear light ring instability exists on the timescales tested. They attribute previous claims of instability to numerical artifacts (gauge instabilities and lack of constraint damping) and provide a framework for correctly simulating these extreme objects in the future.