Boson Stars in $D \ge 4$ Dimensions: Stability, Oscillation Frequencies, and Dynamical Evolutions

Basic Information

• Paper ID: 2510.13988
• Title: Boson Stars in $D \ge 4$ Dimensions: Stability, Oscillation Frequencies, and Dynamical Evolutions
• Authors: Gareth Arturo Marks, Abdullah Al Zaif (Cambridge University)
• Classification: gr-qc (General Relativity and Quantum Cosmology)
• Publication Date: October 15, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2510.13988

Abstract

This paper constructs spherically symmetric boson star solutions in $D \in \{4,5,6\}$ dimensional spacetimes, accounting for quartic self-interaction terms and solitonic potentials. Through perturbative analysis, the pulsation equations are generalized to arbitrary dimensions and potential functions, establishing the existence of radially stable higher-dimensional boson star solutions. The study complements linear results with spherically symmetric nonlinear dynamical evolution, employing dimensional reduction techniques that preserve the full gauge freedom of standard numerical relativity methods while evolving spacetimes of arbitrary background dimensionality. The evolution results demonstrate that solutions identified as perturbatively stable are indeed generically stable under nonlinear spherically symmetric dynamics.

Research Background and Motivation

Core Problem

The core problem addressed by this research is the dynamical stability of higher-dimensional boson stars. Conventionally, it has been believed that in $D > 4$ dimensional spacetimes, due to the characteristic decay of the gravitational field as $1/r^{D-3}$, equilibrium between gravitational attraction and scalar field kinetic energy may be unattainable, leading to inevitable fragmentation instability of boson stars.

Significance

1. Theoretical Importance: Boson stars as horizonless compact object models provide essential tools for studying gravitational wave sources, dark matter candidates, and black hole alternatives
2. Numerical Applications: In four dimensions, boson stars have been employed to study critical phenomena and test the hoop conjecture; the existence of higher-dimensional stable solutions will open new avenues for investigating higher-dimensional gravitational theories
3. Fundamental Physics: Understanding gravitationally bound states in higher-dimensional spacetimes is crucial for comprehending the nature of gravity itself

Existing Limitations

1. Previous studies indicated that higher-dimensional mini boson stars (without self-interaction) are perturbatively unstable for $D > 4$
2. Lack of compelling evidence for dynamical stability of higher-dimensional asymptotically flat boson stars
3. While linear stability analysis is important, it is insufficient for definitive conclusions; nonlinear dynamical evolution verification is necessary

Core Contributions

1. Demonstrated the existence of stable higher-dimensional boson stars: Through introduction of sufficiently strong self-interaction (quartic terms) or solitonic potentials, radially stable boson star solutions are constructed in $D = 5, 6$ dimensions
2. Generalized the pulsation equations: Extended the calculation method for radial oscillation frequencies to arbitrary dimensions and arbitrary potential functions, representing the most general form in the literature
3. Developed a higher-dimensional numerical relativity framework: Based on modified Cartoon methods, implemented BSSN and CCZ4 formulations capable of evolving spacetimes of arbitrary dimensionality
4. Verified consistency between linear and nonlinear theory: Extensive dynamical evolutions confirm that linear stability analysis predictions remain valid in the nonlinear regime

Detailed Methodology

Theoretical Framework

The study employs $D$-dimensional general relativity with a minimally coupled complex scalar field $\phi$, with action: $S = \int d^Dx \sqrt{-g}\left\{\frac{R}{16\pi} - \frac{1}{2}\left[g^{\mu\nu}\nabla_\mu\bar{\phi}\nabla_\nu\phi + V(|\phi|^2)\right]\right\}$

The spherically symmetric metric is adopted as: $ds^2 = -\alpha^2(r)dt^2 + X^2(r)dr^2 + r^2d\Omega^2_{D-2}$

Potential Function Selection

The study considers two classes of potential functions:

1. Quartic Self-Interaction Potential: $V_{\text{massive}}(|\phi|^2) = \mu^2|\phi|^2 + \frac{\lambda}{2}|\phi|^4$
2. Solitonic Potential: $V_{\text{solitonic}}(|\phi|^2) = \mu^2|\phi|^2\left(1 - 2\frac{|\phi|^2}{\sigma_0^2}\right)^2$

Perturbative Analysis Method

By decomposing the complex scalar field into two real fields and introducing perturbations: $\psi_1 = A_b(1 + \delta\psi_1), \quad \psi_2 = A_b\delta\psi_2$$\alpha = \alpha_b\left(1 + \frac{1}{2}\delta\alpha\right), \quad X = X_b\left(1 + \frac{1}{2}\delta X\right)$

The pulsation equations applicable to arbitrary dimensions and potential functions are derived:

$f'' = \left[\left(\frac{2\omega^2 - \chi^2}{\alpha_b^2} + 2A_b^2\frac{d^2V}{d(|\phi|^2)^2} + 2\frac{dV}{d|\phi|^2}\right) + \text{other terms}\right]X_b^2f + \cdots$

where $\chi$ denotes the radial oscillation frequency, and $f$ and $g$ are perturbation functions.

Numerical Evolution Scheme

The dimensional reduction technique based on modified Cartoon methods projects physical quantities onto specific coordinate axes. The evolution system is based on BSSN and CCZ4 formulations:

• Conformal Factor: $\chi = \gamma^{1/(D-1)}$
• Conformal Metric: $\tilde{\gamma}_{ij} = \chi\gamma_{ij}$
• Mean Curvature: $K = \gamma^{ij}K_{ij}$
• Traceless Extrinsic Curvature: $\tilde{A}_{ij} = \chi(K_{ij} - \frac{1}{D-1}K\gamma_{ij})$

Experimental Setup

Parameter Space Exploration

• Dimensions: $D \in \{4,5,6\}$
• Quartic Coupling: $\hat{\lambda} := \lambda/\mu^2 \in [0, 4000]$
• Solitonic Parameter: $\sigma_0 \in [0.06, 0.5]$
• Central Amplitude: $A_0$ range adjusted according to different models

Numerical Methods

• Fourth-order spacetime discretization
• Runge-Kutta time integration
• Courant factor: 0.4 for $D=4,5$, 0.25 for $D=6$
• Grid resolution: $\Delta x \in \{1/8, 1/12, 1/18\}$

Stability Criteria

• Linear Stability: Fundamental radial mode $\chi_0^2 > 0$
• Dynamical Stability: Central amplitude variation less than 10% over evolution time $t \sim 10^4 M^{1/(D-3)}$
• Instability Timescale: Time for central amplitude to deviate 10% from initial value

Experimental Results

Main Findings

1. Critical Parameter Values

2. Oscillation Frequency Verification

Power spectrum analysis reveals that oscillation frequencies in numerical evolution are highly consistent with perturbative calculation results, validating the accuracy of theoretical predictions.

3. Dynamical Evolution Results

• Stable Models: Successfully identified families of dynamically stable boson stars in $D=5,6$ dimensions
• Unstable Models: Display different behaviors according to binding energy sign:
  • $E_B < 0$: Migration to stable branch or collapse to black hole
  • $E_B > 0$: Dispersion or collapse

4. Six-Dimensional Peculiarities

In $D=6$ dimensions, stable quartic interaction models do not include the limit solution at $A_0 = 0$, with a minimum achievable compactness.

Constraint Verification

The CCZ4 formulation demonstrates superior constraint control compared to BSSN in long-term evolution, with both Hamiltonian and momentum constraint violations converging to finite values.

Related Work

Historical Development

• Kaup (1968) and Ruffini & Bonazzola (1969): First introduced the boson star concept
• Lee & Pang (1989): Established the theoretical foundation for four-dimensional stability
• Franzin (2024): Proved perturbative instability of higher-dimensional mini boson stars

Relationship to This Work

This paper provides the first compelling evidence for dynamical stability of higher-dimensional asymptotically flat boson stars, filling an important gap in the field.

Conclusions and Discussion

Main Conclusions

1. Self-interaction can stabilize higher-dimensional boson stars: In $D=5,6$ dimensions, sufficiently strong quartic self-interaction or appropriate solitonic potentials can produce radially stable solutions
2. Linear and nonlinear theory are consistent: Predictions from perturbative stability analysis are verified in nonlinear spherically symmetric dynamics
3. Binding energy is not a necessary and sufficient condition for stability: Models with positive binding energy but dynamical stability are discovered

Limitations

1. Spherical Symmetry Only: Non-radial instabilities are not explored
2. Dimensional Restriction: Investigation limited to $D=6$; behavior at higher dimensions remains unknown
3. Numerical Challenges: Higher-dimensional evolution requires more stringent numerical conditions

Future Directions

1. Non-Spherically Symmetric Analysis: Study stability under $SO(D-2)$ or $SO(D-3)$ symmetry
2. Other Compact Objects: Extension to other models such as Proca stars
3. Higher-Curvature Gravity: Consider modified theories such as Gauss-Bonnet gravity
4. Large $D$ Limit: Explore universal behavior at arbitrarily high dimensions

In-Depth Evaluation

Strengths

1. Theoretical Completeness: Complete research chain from linear analysis to nonlinear verification
2. Numerical Innovation: The developed dimensional reduction method has broad applicability
3. Result Reliability: Convergence and stability verified through multiple numerical tests
4. Physical Insight: Reveals new physics of higher-dimensional gravitationally bound states

Weaknesses

1. Computational Cost: Full-dimensional non-spherically symmetric evolution remains computationally infeasible
2. Parameter Space: Pathological behavior exists in certain parameter regions (e.g., $D=6$ solitonic models)
3. Physical Interpretation: Deeper understanding of stabilization mechanisms is needed

Impact

1. Pioneering: First demonstration of stable higher-dimensional boson star existence
2. Methodology: Provided numerical framework applicable to other higher-dimensional gravitational studies
3. Application Prospects: Provides tools for higher-dimensional critical phenomena and hoop conjecture research

Applicable Scenarios

• Higher-dimensional gravitational theory research
• Compact object astrophysical modeling
• Numerical relativity method development
• Nonlinear field theory stability analysis

References

This paper cites 77 important references spanning multiple fields including boson star theory, numerical relativity, and higher-dimensional gravity, providing a solid theoretical foundation for the research.